Conformational search using a molecular dynamics–minimization procedure: Applications to clusters of coulombic charges,Lennard–Jones particles,and waters

A new conformational search method, molecular dynamics–minimization (MDM), is proposed, which combines a molecular dynamics sampling strategy with energy minimizations in the search for low-energy molecular structures. This new method is applied to the search for low energy configurations of clusters of coulombic charges on a unit sphere, Lennard–Jones clusters, and water clusters. The MDM method is shown to be efficient in finding the lowest energy configurations of these clusters. A closer comparison of MDM with alternative conformational search methods on Lennard–Jones clusters shows that, although MDM is not as efficient as the Monte Carlo–minimization method in locating the global energy minima, it is more efficient than the diffusion equation method and the method of minimization from randomly generated structures. Given the versatility of the molecular dynamics sampling strategy in comparison to Monte Carlo in treating molecular complexes or molecules in explicit solution, one anticipates that the MDM method could be profitably applied to conformational search problems where the number of degrees of freedom is much greater. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 60–70, 1998     

含乙烯氧基柔链的苯并菲盘状液晶清亮点的平均场理论研究

我们前期的合成和实验发现, 对于一系列含端烯氧基醚链的苯并菲分子, 它们没有分子间氢键, 其生成的液晶的清亮点的高低与端烯氧基醚链的数目有关. 本工作采用Gay-Berne势模拟液晶分子的双体势, 用量子化学和分子力学相结合的ONIOM(MP2/6-31G*(0.25):UFF)方法优化含端烯氧基醚链的苯并菲分子的单体和双体的构型和能量, 从而获得Gay-Berne势的参数. 在此基础上, 根据平均场理论, 解一个关于序参量和温度的积分方程, 求出液晶的清亮点. 计算结果与实验趋势相同, 可以解释清亮点和分子端烯氧基醚链数目的关系.     

Simple computing of the viscosity of water–dioxane mixtures,according to a fluctuating SPC/E‐Ih interstitial model

In this contribution, we propose a new simple and fast computing approach to model liquid water and water‐1,4 dioxane mixture, which is intermediate between molecular dynamics calculation and lattice fluid H‐bonding theory in which water molecules are arranged according to a fluctuating hexagonal I_h network. It generalizes the interstitial mixture models of Samoilov and Pauling by increasing the number of spatial configurations into the space of phases due to the addition of rotations and vibrations of the molecules around their equilibrium positions. The potential energy landscape is determined by the generalization of the Madelung's approach to the molecular water network using simple point charge model (SPC/E partial charges). This approach allows useful predictions as the possibility of incorporation of a molecule of dioxane in the water I_h‐network or how some cooperative dynamic jumps (site ‐ hole) of water molecules, contribute to lower Eyring's activation energies and thus to facilitate the viscous flow. © 2017 Wiley Periodicals, Inc.     

Efficient free energy calculation of water across lipid membranes

An efficient free energy (FE) calculation of a water molecule to go across lipid membranes is presented. Both overlapping distribution and cavity insertion Widom methods are complementarily used. The former is useful for a dense region where water molecules are found, i.e., from the interfacial to bulk water region, while the latter works well in the low density region, i.e., the hydrocarbon region. Since both methods evaluate the excess chemical potential of water, the obtained FE profile is free from the fitting problem usually arisen when a FE difference method is used. A diphytanyl phosphatidylcholine bilayer is used for our test calculations. An excellent and fast convergence of the chemical potential is obtained when each method is applied for the appropriate region. The estimated FE barrier using the Ewald method for the electrostatic interaction is approximately 7.2 kcal/mol, which is higher than that using the interaction cutoff of 20 A by about 0.9 kcal/mol.     

The molecular-level relationship between the properties of liquid water molecules and orientational order

The relationship between the orientational (tetrahedral) order (q) of an individual liquid water molecule and its various properties such as Voronoi volume, potential energy, kinetic energy, and nearest neighbors was thoroughly examined using molecular dynamics simulations of TIP5P model at 278, 298, and 318 K. By constructing Voronoi polyhedra (VP), we found that the average volume of water molecules classified according to q decreased monotonically as q increased, while the surface of VP increased in the range of high q. Kinetic energy was almost invariant but potential energy decreased monotonously as q increased. The volumes of molecules having a very large q increased as temperature decreased, implying a possible density maximum phenomena. Using time correlation functions, it was shown that total energy rather than potential energy was a more significant factor in the determination of the orientational order. With varying temperature, the relation between the properties of central molecules and those of nearest neighbor molecules were investigated. It required a very systematic cooperative motion to obtain LDL (low-density liquid) formed by ordering. It was supposed that density maximum phenomena should be accomplished by the growth of LDL and HDL (high-density liquid) of low quality with a consistent population of HDL and a drop of LDL of high quality as temperature lowered.     

An analysis of molecular packing and chemical association in liquid water using quasichemical theory

We calculate the hydration free energy of liquid TIP3P water at 298 K and 1 bar using a quasi-chemical theory framework in which interactions between a distinguished water molecule and the surrounding water molecules are partitioned into chemical associations with proximal (inner-shell) waters and classical electrostatic-dispersion interactions with the remaining (outer-shell) waters. The calculated free energy is found to be independent of this partitioning, as expected, and in excellent agreement with values derived from the literature. An analysis of the spatial distribution of inner-shell water molecules as a function of the inner-shell volume reveals that water molecules are preferentially excluded from the interior of large volumes as the occupancy number decreases. The driving force for water exclusion is formulated in terms of a free energy for rearranging inner-shell water molecules under the influence of the field exerted by outer-shell waters in order to accommodate one water molecule at the center. The results indicate a balance between chemical association and molecular packing in liquid water that becomes increasingly important as the inner-shell volume grows in size.     

Molecular modelling of liquid crystal systems: An internal coordinate Monte Carlo approach

A Monte Carlo scheme is presented which is designed to provide a convenient mechanism to model accurately the internal molecular structure of liquid crystalline molecules. The technique stores atomic positions in terms of bond lengths, bond angles and dihedral angles within a Z-matrix, and the Monte Carlo scheme involves generating trial configurations from changes to the Z-matrix using the MM2 molecular mechanics potential to describe energy changes between different molecular conformations. The technique is applied to the liquid crystal molecule 4-n-pentyl-4'-cyanobiphenyl (5CB), and results are presented for the conformational populations and dihedral angle distributions of 5CB in the gas phase at 300 K. The effect of a nematic mean field on the distribution of molecular conformations is also examined via the addition of a conformation-dependent potential of mean torque to the internal energy.     

Average sequential water molecule binding enthalpies of M(H2O)(19-124)2+ (M = Co, Fe, Mn, and Cu) measured with ultraviolet photodissociation at 193 and 248 nm

The average sequential water molecule binding enthalpies to large water clusters (between 19 and 124 water molecules) containing divalent ions were obtained by measuring the average number of water molecules lost upon absorption of an UV photon (193 or 248 nm) and using a statistical model to account for the energy released into translations, rotations, and vibrations of the products. These values agree well with the trend established by more conventional methods for obtaining sequential binding enthalpies to much smaller hydrated divalent ions. The average binding enthalpies decrease to a value of ~10.4 kcal/mol for n > ~40 and are insensitive to the ion identity at large cluster size. This value is close to that of the bulk heat of vaporization of water (10.6 kcal/mol) and indicates that the structure of water in these clusters may more closely resemble that of bulk liquid water than ice, owing either to a freezing point depression or rapid evaporative cooling and kinetic trapping of the initial liquid droplet. A discrete implementation of the Thomson equation using parameters for liquid water at 0 °C generally fits the trend in these data but provides values that are ~0.5 kcal/mol too low.     

On the application of accelerated molecular dynamics to liquid water simulations

Our group recently proposed a robust bias potential function that can be used in an efficient all-atom accelerated molecular dynamics (MD) approach to simulate the transition of high energy barriers without any advance knowledge of the potential-energy landscape. The main idea is to modify the potential-energy surface by adding a bias, or boost, potential in regions close to the local minima, such that all transitions rates are increased. By applying the accelerated MD simulation method to liquid water, we observed that this new simulation technique accelerates the molecular motion without losing its microscopic structure and equilibrium properties. Our results showed that the application of a small boost energy on the potential-energy surface significantly reduces the statistical inefficiency of the simulation while keeping all the other calculated properties unchanged. On the other hand, although aggressive acceleration of the dynamics simulation increases the self-diffusion coefficient of water molecules greatly and dramatically reduces the correlation time of the simulation, configurations representative of the true structure of liquid water are poorly sampled. Our results also showed the strength and robustness of this simulation technique, which confirm this approach as a very useful and promising tool to extend the time scale of the all-atom simulations of biological system with explicit solvent models. However, we should keep in mind that there is a compromise between the strength of the boost applied in the simulation and the reproduction of the ensemble average properties.     

Spectroscopic characterization of microscopic hydrogen-bonding disparities in supercritical water

The local hydrogen-bonding environment in supercritical water (380 degrees C, 300 bars, density 0.54 gcm3) was studied by x-ray Raman scattering at the oxygen K edge. The spectra are compared to those of the gas phase, liquid surface, bulk liquid, and bulk ice, as well as to calculated spectra. The experimental model systems are used to assign spectral features and to quantify specific local hydrogen-bonding situations in supercritical water. The first coordination shell of the molecules is characterized in more detail with the aid of the calculations. Our analysis suggests that approximately 65% of the molecules in supercritical water are hydrogen bonded in configurations that are distinctly different from those in liquid water and ice. In contrast to liquid water the bonded molecules in supercritical water have four intact hydrogen bonds and in contrast to ice large variations of bond angles and distances are observed. The remaining approximately 35% of the molecules exhibit two free O-H bonds and are thus either not involved in hydrogen bonding at all or have one or two hydrogen bonds on the oxygen side. We determine an average O-O distance of 3.1+/-0.1 A in supercritical water for the H bonded molecules at the conditions studied here. This and the corresponding hydrogen bond lengths are shown to agree with neutron- and x-ray-diffraction data at similar conditions. Our results on the local hydrogen-bonding environment with mainly two disparate hydrogen-bonding configurations are consistent with an extended structural model of supercritical water as a heterogeneous system with small patches of bonded molecules in various tetrahedral configurations and surrounding nonbonded gas-phase-like molecules.     

Orientation of a benzene molecule inside a carbon nanotube

Benzene molecules confined in carbon nanotubes of varying radii are employed as semiconductors in electronic nanodevices, and their orientation determines the electrical properties of the system. In this paper, we investigate the interaction energy of all the possible configurations of a benzene molecule inside various carbon nanotubes and then we determine the equilibrium configuration. We adopt the continuous approach together with the semi-empirical Lennard-Jones potential function to model van der Waals interaction between a benzene molecule and a carbon nanotube. This approach results in an analytical expression, which accurately approximates the interaction energy and can be readily used to generate numerical data. We find that horizontal, tilted and perpendicular configurations on the axis of the carbon nanotube are all possible equilibrium configurations of the benzene molecule when the radius of the carbon nanotube is less than 5.580 Å. However, when the radius of the carbon nanotube is larger than 5.580 Å an offset horizontal orientation is the only possible equilibrium configuration of the benzene molecule. In the limiting case, the orientation of a benzene molecule on a graphene sheet can be derived simply by letting the radius of the carbon nanotube tend to infinity.     

Ab initio molecular dynamics study of carbon dioxide and bicarbonate hydration and the nucleophilic attack of hydroxide on CO2

We apply ab initio molecular dynamics (AIMD) to study the hydration structures of the carbon dioxide molecule and the bicarbonate and carbonate anions in liquid water. We also compute the free energy change associated with the nucleophilic attack of the hydroxide ion on carbon dioxide. CO2 behaves like a hydrophobic species and exhibits weak interactions with water molecules. The bicarbonate and carbonate ions are strongly hydrated and coordinate to an average of 6.9 and 8.7 water molecules, respectively. The energetics for the reaction in the gas phase are investigated using density functional theory and second-order M?ller-Plesset perturbation theory (MP2) in conjunction with high-quality basis sets. Using umbrella sampling techniques, we compute the standard state, aqueous phase free energy difference associated with the reaction CO2+OH--->HCO3- after correcting AIMD energies with MP2 results. Our predictions are in good agreement with experiments. The hydration structures along the reaction coordinate, which give rise to a predicted 9.7 kcal/mol standard state free energy barrier, are further analyzed.     

On the Debye-Waller factor of hexagonal ice: a computer simulation study

We investigate by molecular dynamics (MD) simulations the temperature dependence of the Debye-Waller (DW) factor of hexagonal ice with 25 different proton-disordered configurations. Each initial configuration is composed of 288 water molecules with no net dipole moment. The intermolecular interaction of water is described by TIP4P potential. Each production run of the simulation is 15 ns or longer. We observe a change in slope of the DW factor around 200 K, which cannot be explained within the framework of either classical or quantum harmonic approximation. Configurations generated by MD simulations are subjected to the steepest descent energy minimization. Analysis of the local energy minimum structures reveals that water molecules above 200 K jump to other lattice sites via some local energy minimum structures which contain some water molecules sitting on the locations other than the lattice sites. As time evolves, these defect molecules move back and forth to the lattice sites yielding defect-free structures. Those motions are responsible for the unusual increase in the DW factor at high temperatures. In making a transition from an energy-minimum structure to another one, a small number of water molecules are involved in a highly cooperative fashion. The larger DW factor at higher temperature arises from jump-like motions of water molecules among these locally stable configurations which may or may not be a family of the proton-disordered ice forms satisfying the "ice rule".     

Models of hydration and isomeric transitions of glucose molecules in aqueous solutions

The kinematic viscosity of aqueous glucose solutions is studied. It is found that the hydrodynamic radius of monosaccharide molecule in an aqueous solution depends on temperature in the range of 290–355 K. Using a bimodal model of the energy states of the volume in which the glucose molecule is located and local equilibrium is established, it is shown that the above-mentioned dependence can be attributed to disturbances in the equilibrium of isomeric transitions, induced by variations in temperature. The parameters of isomeric transitions for a glucose molecule in an aqueous solvent, the probability of “chair” and “boat” configurations occurring for glucose molecules, and the number of water molecules in the hydration shells of these configurations are calculated; the strain of the chemical bonds in the chair configuration of a glucose molecule is estimated.     

Adsorption of 1-octanol at the free water surface as studied by Monte Carlo simulation

The adsorption of 1-octanol at the free water surface has been investigated by Monte Carlo computer simulation. Six different systems, built up by an aqueous and a vapor phase, the latter also containing various number of octanol molecules, have been simulated. The number of the octanol molecules has been chosen in such a way that the octanol surface density varies in a broad range, between 0.27 and 7.83 micromol/m(2) in the six systems simulated. For reference, the interfacial system containing bulk liquid octanol in the apolar phase has also been simulated. The results have shown that the formation of hydrogen bonds between the interfacial water and adsorbed octanol molecules is of key importance in determining the properties of the adsorbed layer. At low octanol surface concentration values all the octanol molecules are strongly (i.e., by hydrogen bonds) bound to the aqueous phase, whereas their interaction with each other is negligibly small. Hence, they are preferentially oriented in such a way that their own binding energy (and thus their own free energy) is minimized. In this preferred orientation the O-H bond of the octanol molecule points flatly toward the aqueous phase, declining by about 30 degrees from the interfacial plane, irrespectively from whether the octanol molecule is the H-donor or the H-acceptor partner in the hydrogen bond. Hence, in its preferred orientation the octanol molecule can form at least two low energy hydrogen bonds with water: one as a H-donor and another one as a H-acceptor. Moreover, the preferred orientation of the hydrogen bonded water partners is close to one of the two preferred interfacial water alignments, in which the plane of the water molecule is parallel with the interface. When increasing the octanol surface density, the water surface gets saturated with hydrogen bonded octanols, and hence any further octanol molecule can just simply condense to the layer of the adsorbed octanols. The surface density value at which this saturation occurs is estimated to be about 1.7 micromol/m(2). Above this surface density value the hydrogen bonded octanols and their water partners are oriented in such a way that the number of the water-octanol hydrogen bonds is maximized. Hence, the preferred alignment of the O...O axes of these hydrogen bonds is perpendicular to the interface. This orientation is far from the optimal alignment of the individual octanol molecules, which is also reflected in the observed fact that, unlike in the case of many other adsorbents, the average molecular binding energy of the adsorbed octanol molecules increases (i.e., becomes less negative) with increasing octanol surface density.     

A Monte Carlo study of water clusters

The free energy, internal energy and entropy for water clusters of 8 and 64 molecules have been calculated using the Monte Carlo method and both the Lie-Clementi (C-XII) and Stillinger (ST2) potential functions. Detailed structural data for a cluster of 64 water molecules was also obtained by averaging over a large number of configurations. Average dipole orientations and radial density distributions indicated that individual water molecules were not particularly restricted to fixed orientations and that general ordering tendencies were small.     

Determination of the chemical potential using energy-biased sampling

An energy-biased method to evaluate ensemble averages requiring test-particle insertion is presented. The method is based on biasing the sampling within the subdomains of the test-particle configurational space with energies smaller than a given value freely assigned. These energy wells are located via unbiased random insertion over the whole configurational space and are sampled using the so-called Hit-and-Run algorithm, which uniformly samples compact regions of any shape immersed in a space of arbitrary dimensions. Because the bias is defined in terms of the energy landscape it can be exactly corrected to obtain the unbiased distribution. The test-particle energy distribution is then combined with the Bennett relation for the evaluation of the chemical potential. We apply this protocol to a system with relatively small probability of low-energy test-particle insertion, liquid argon at high density and low temperature, and show that the energy-biased Bennett method is around five times more efficient than the standard Bennett method. A similar performance gain is observed in the reconstruction of the energy distribution.     

Self-assembled monolayers and preorganization of organosilanes prior to surface grafting onto silica: a quantum mechanical study

Quantum chemical calculations have been carried out on the grafting of chain organosilane compounds on SiO(2)-hydroxylated solid surfaces. It is shown that a single molecule interacting with the surface lies flat to it, inhibiting further homogeneous film growth. This adsorption exhibits two molecule/surface interactions: a covalent bond on one side of the molecule and a hydrogen bond on the other side. We then investigate the possible preorganization of the molecules before grafting due to the presence of water molecules either in the gas/liquid phase or near the surface. This gives rise to the formation of dimerized chains. We then demonstrate that this preorganization process prevents subsequent lying flat of the molecules to the substrate after grafting. Energetics and associated configurations of the overall deposition process are discussed in detail and provide new insights on the understanding of the formation of self-assembled homogeneous organic films on microelectronics-type substrates.