Molecular Dynamics of Water Embedded Carbon Nanocones: Surface Waves Observation

All-Atom Molecular Dynamics Simulations of the Temperature Response of Poly(glycidyl ether)s with Oligooxyethylene Side Chains Terminated with Alkyl Groups

Recently, experimental investigations of a class of temperature-responsive polymers tethered to oligooxyethylene side chains terminated with alkyl groups have been conducted. In this study, aqueous solutions of poly(glycidyl ether)s (PGE) with varying numbers of oxyethylene units, poly(methyl(oligooxyethylene)_n glycidyl ether) (poly(Me(EO)_nGE)), and poly(ethyl(oligooxyethylene)_n glycidyl ether) (poly(Et(EO)_nGE) (n = 0, 1, and 2) were investigated by all-atom molecular dynamics simulations, focusing on the thermal responses of their chain extensions, the recombination of intrapolymer and polymer-water hydrogen bonds, and water-solvation shells around the alkyl groups. No clear relationship was established between the phase-transition temperature and the polymer-chain extensions unlike the case for the coil-globule transition of poly(N-isopropylacrylamide). However, the temperature response of the first water-solvation shell around the alkyl group exhibited a notable correlation with the phase-transition temperature. In addition, the temperature at which the hydrophobic hydration shell strength around the terminal alkyl group equals the bulk water density (TCRP) was slightly lower than the cloud point temperature (TCLP) for the methyl-terminated poly(Me(EO)_nGE) and slightly higher for the ethyl-terminated poly(Et(EO)_nGE). It was concluded that the polymer-chain fluctuation affects the relationship between TCRP and TCLP.

Umbrella Sampling Simulations of Carbon Nanoparticles Crossing Immiscible Solvents

We use molecular dynamics to compute the free energy of carbon nanoparticles crossing a hydrophobic–hydrophilic interface. The simulations are performed on a biphasic system consisting of immiscible solvents (i.e., cyclohexane and water). We solvate a carbon nanoparticle into the cyclohexane layer and use a pull force to drive the nanoparticle into water, passing over the interface. Next, we accumulate a series of umbrella sampling simulations along the path of the nanoparticle and compute the solvation free energy with respect to the two solvents. We apply the method on three carbon nanoparticles (i.e., a carbon nanocone, a nanotube, and a graphene nanosheet). In addition, we record the water-accessible surface area of the nanoparticles during the umbrella simulations. Although we detect complete wetting of the external surface of the nanoparticles, the internal surface of the nanotube becomes partially wet, whereas that of the nanocone remains dry. This is due to the nanoconfinement of the particular nanoparticles, which shields the hydrophobic interactions encountered inside the pores. We show that cyclohexane molecules remain attached on the concave surface of the nanotube or the nanocone without being disturbed by the water molecules entering the cavity.

On the Consistency of the Exfoliation Free Energy of Graphenes by Molecular Simulations

Monolayer graphene is now produced at significant yields, by liquid phase exfoliation of graphites in solvents. This has increased the interest in molecular simulation studies to give new insights in the field. We use decoupling simulations to compute the exfoliation free energy of graphenes in a liquid environment. Starting from a bilayer graphene configuration, we decouple the Van der Waals interactions of a graphene monolayer in the presence of saline water. Then, we introduce the monolayer back into water by coupling its interactions with water molecules and ions. A different approach to compute the graphene exfoliation free energy is to use umbrella sampling. We apply umbrella sampling after pulling the graphene monolayer on the shear direction up to a distance from a bilayer. We show that the decoupling and umbrella methods give highly consistent free energy results for three bilayer graphene samples with different size. This strongly suggests that the systems in both methods remain closely in equilibrium as we move between the states before and after the exfoliation. Therefore, the amount of nonequilibrium work needed to peel the two layers apart is minimized efficiently.

Binding Free Energy Calculations of Bilayer Graphenes Using Molecular Dynamics

Bilayer graphenes are dimeric assemblies of single graphene layers bound together by π-complexation interactions. Controlling these assemblies can be complicated, as the layered compounds disperse in solvents or aggregate into higher columnar configurations and clusters. One way to assess the interactions that contribute to the stability of the layered compounds is to use molecular simulation. We perform pulling molecular dynamics on bilayer graphenes with different sizes and obtain the normal and shear force profiles of dissociation. We generate pathways of dissociation along the two directions and calculate the binding free energies of the structures with umbrella sampling simulations. We show that the dissociation process is direction-dependent. Along the shear direction, we compute the same free energy for the different samples, which validates the consistency of our simulations. We notice that the dissociation is less adiabatic on the normal than the shear direction, having an entropic contribution to the Gibbs energy. This contribution is more enhanced for the larger bilayer graphenes.

Computational Approaches in Membrane Science and Engineering

Computational modelling and simulation form a consolidated branch in the multidisciplinary field of membrane science and technology [...]