Efficient Hybrid Particle-Field Coarse-Grained Model of Polymer Filler Interactions: Multiscale Hierarchical Structure of Carbon Black Particles in Contact with Polyethylene.

In the present study, we propose, validate, and give first applications for large-scale systems of coarse-grained models suitable for filler/polymer interfaces based on carbon black (CB) and polyethylene (PE). The computational efficiency of the proposed approach, based on hybrid particle-field models (hPF), allows large-scale simulations of CB primary particles of realistic size (∼20 nm) embedded in PE melts. The molecular detailed models, here introduced, allow a microscopic description of the bound layer, through the analysis of the conformational behavior of PE chains adsorbed on different surface sites of CB primary particles, where the conformational behavior of adsorbed chains is different from models based on flat infinite surfaces. On the basis of the features of the systems, an optimized version of OCCAM code for large-scale (up to more than 8 million of beads) parallel runs is proposed and benchmarked. The computational efficiency of the proposed approach opens the possibility of a computational screening of the bound layer, involving the optimal combination of surface chemistry, size, and shape of CB aggregates and the molecular weight distribution of the polymers achieving an important tool to address the polymer/fillers interface and interphase engineering in the polymer industry.

Molecular structure and multi-body potential of mean force in silica-polystyrene nanocomposites.

We perform a systematic application of the hybrid particle-field molecular dynamics technique [Milano, et al., J. Chem. Phys., 2009, 130, 214106] to study interfacial properties and potential of mean force (PMF) for separating nanoparticles (NPs) in a melt. Specifically, we consider Silica NPs bare or grafted with Polystyrene chains, aiming to shed light on the interactions among free and grafted chains affecting the dispersion of NPs in the nanocomposite. The proposed hybrid models show good performances in catching the local structure of the chains, and in particular their density profiles, documenting the existence of the "wet-brush-to-dry-brush" transition. By using these models, the PMF between pairs of ungrafted and grafted NPs in Polystyrene matrix are calculated. Moreover, we estimate the three-particle contribution to the total PMF and its role in regulating the phase separation on the nanometer scale. In particular, the multi-particle contribution to the PMF is able to give an explanation of the complex experimental morphologies observed at low grafting densities. More in general, we propose this approach and the models utilized here for a molecular understanding of specific systems and the impact of the chemical nature of the systems on the composite final properties.

Molecular Insights into the Eutectic Tripalmitin/Tristearin Binary System.

A molecular interpretation of the eutectic behavior of a binary mixture of tristearin (SSS) and tripalmitin (PPP) triglycerides was formulated using computer simulations and experimental techniques (calorimetry and X-ray scattering). A eutectic composition was identified using both experimental and computer simulation techniques at a composition of 70% PPP and 30% SSS, in agreement with previous findings in the literature. The decrease in the melting temperature at the eutectic composition can be ascribed to an interplay between enthalpic and entropic effects. In particular, a lower global melting enthalpy at the eutectic composition was detected here, caused by a less efficient packing of the triglycerides in the crystal. On the other hand, a higher crystalline disorder is reflected in a lower change in the entropy of melting. The simultaneous decrease in global enthalpy and entropy has a contrasting effect on the melting temperature, with a slight melting point depression found in both experiment and simulations, resulting from a combination of enthalpic and entropic factors. Computer simulations showed, in fact, that the eutectic effect can be ascribed to the reduction of crystalline order when SSS molecules are incorporated into the PPP crystal structure. This decrease of the crystalline order is due to the protrusion of SSS end-chains (last three carbons of each alkyl chain) into the interlamellar space between adjacent lamella. These end-chains disturb the orderly stacking of the lamella, as evidenced by low-density regions in the interlamellar space. Thus, the greater disorder of the last atoms of the SSS alkyl chains is consequently due to the greater conformational freedom. At molecular level, in fact, the conformational freedom of terminal atoms of SSS surrounded by shorter PPP molecules is larger than the conformational freedom of longer SSS in the neighborhood of short PPP.

On the calculation of the potential of mean force between atomistic nanoparticles.

We study the potential of mean force (PMF) between atomistic silica and gold nanoparticles in the vacuum by using molecular dynamics simulations. Such an investigation is devised in order to fully characterize the effective interactions between atomistic nanoparticles, a crucial step to describe the PMF in high-density coarse-grained polymer nanocomposites. In our study, we first investigate the behavior of silica nanoparticles, considering cases corresponding to different particle sizes and assessing results against an analytic theory developed by Hamaker for a system of Lennard-Jones interacting particles (H.C. Hamaker, Physica A 4, 1058 (1937)). Once validated the procedure, we calculate effective interactions between gold nanoparticles, which are considered both bare and coated with polyethylene chains, in order to investigate the effects of the grafting density [Formula: see text] on the PMF. Upon performing atomistic molecular dynamics simulations, it turns out that silica nanoparticles experience similar interactions regardless of the particle size, the most remarkable difference being a peak in the PMF due to surface interactions, clearly apparent for the larger size. As for bare gold nanoparticles, they are slightly interacting, the strength of the effective force increasing for the coated cases. The profile of the resulting PMF resembles a Lennard-Jones potential for intermediate [Formula: see text], becoming progressively more repulsive for high [Formula: see text] and low interparticle separations.

Biomembrane solubilization mechanism by Triton X-100: a computational study of the three stage model.

The solubilization mechanism of lipid membranes in the presence of Triton X-100 (TX-100) is investigated at molecular resolution using molecular dynamics (MD) simulations. Thanks to the large time and length scales accessible by the hybrid particle-field formulation of the models employed here, the complex process of membrane solubilization has been studied, with the goal of verifying the three stage model reported in the literature. DPPC lipid bilayers and vesicles have been studied at different concentrations of the TX-100 detergent employing coarse grained (CG) models. Systems up to ∼600.000 beads, corresponding to more than 2 millions heavy atoms, have been simulated. Moreover, in order to clarify several experimental pieces of evidence, both slow and fast detergent partition scenarios have been investigated. Flat and curved (vesicles) lipid bilayer surfaces, interacting with TX-100, have been considered to study the curvature effects on the detergent partition rate in the membrane. Shape and conformational changes of mixed DPPC/TX-100 vesicles, as a function of TX-100 content, have also been studied. In particular, high curvature surfaces, corresponding to a higher local TX-100 content, promote a membrane rupture. In flat lipid surfaces, on the time scale simulated the detergent partition is almost absent, following a different pathway of the solubilization membrane mechanism.

MARTINI Coarse-Grained Model of Triton TX-100 in Pure DPPC Monolayer and Bilayer Interfaces.

The coarse-grained MARTINI model of Triton TX-100 has been validated by direct comparison of the experimental and calculated area increase in pure DPPC lipid bilayers and monolayers at water/air interfaces in the presence of surfactant and by comparison of electron density profiles calculated with more detailed atomistic models based on the CHARMM force field. Bilayer simulations have been performed and compared with monolayers and with atomistic models. The validated CG model has been employed to study the phase separation of TX-100 molecules in lipid bilayers and the effect of the lipid bilayer curvature.

Communication: molecular dynamics and (1)H NMR of n-hexane in liquid crystals.

The NMR spectrum of n-hexane orientationally ordered in the nematic liquid crystal ZLI-1132 is analysed using covariance matrix adaptation evolution strategy (CMA-ES). The spectrum contains over 150 000 transitions, with many sharp features appearing above a broad, underlying background signal that results from the plethora of overlapping transitions from the n-hexane as well as from the liquid crystal. The CMA-ES requires initial search ranges for NMR spectral parameters, notably the direct dipolar couplings. Several sets of such ranges were utilized, including three from MD simulations and others from the modified chord model that is specifically designed to predict hydrocarbon-chain dipolar couplings. In the end, only inaccurate dipolar couplings from an earlier study utilizing proton-proton double quantum 2D-NMR techniques on partially deuterated n-hexane provided the necessary estimates. The precise set of dipolar couplings obtained can now be used to investigate conformational averaging of n-hexane in a nematic environment.

Order and conformation of biphenyl in cyanobiphenyl liquid crystals: a combined atomistic molecular dynamics and 1H NMR study.

The alignment of biphenyl (2P) in the liquid-crystal phases of 4-n-pentyl-4'-cyanobiphenyl (5CB) and 4-n-octyl-4'-cyanobiphenyl (8CB) is investigated by using a combination of predictive atomistic molecular dynamics (MD) simulations and (1)H liquid-crystal nuclear magnetic resonance (LXNMR) residual dipolar coupling measurements. A detailed comparison and validation of the MD results with LXNMR is provided, showing a good agreement between the simulated and experimental dipolar couplings at the same reduced temperature. MD is then used to examine the location of 2P in the smectic phase, which is unavailable to LXNMR, and 2P is found to be rather uniformly distributed. The combination of MD and NMR spectroscopy provides detailed information about the order, interconnection between orientation and conformation, local positional order, and interactions with the liquid-crystalline solvent.

Supramolecular organization of functional organic materials in the bulk and at organic/organic interfaces: a modeling and computer simulation approach.

The molecular organization of functional organic materials is one of the research areas where the combination of theoretical modeling and experimental determinations is most fruitful. Here we present a brief summary of the simulation approaches used to investigate the inner structure of organic materials with semiconducting behavior, paying special attention to applications in organic photovoltaics and clarifying the often obscure jargon hindering the access of newcomers to the literature of the field. Special attention is paid to the choice of the computational "engine" (Monte Carlo or Molecular Dynamics) used to generate equilibrium configurations of the molecular system under investigation and, more importantly, to the choice of the chemical details in describing the molecular interactions. Recent literature dealing with the simulation of organic semiconductors is critically reviewed in order of increasing complexity of the system studied, from low molecular weight molecules to semiflexible polymers, including the challenging problem of determining the morphology of heterojunctions between two different materials.

An atomistic description of the nematic and smectic phases of 4-n-octyl-4' cyanobiphenyl (8CB).

We report the results of atomistic molecular dynamics simulations of 4-n-octyl-4' cyanobiphenyl (8CB) on samples of 750 and 3000 molecules showing the spontaneous formation of the nematic phase and then of smectic layers by gradually cooling down from the isotropic phase. Orientational, positional, and mixed order parameters, layer spacing, translational diffusion tensor components and their temperature dependence are reported. A detailed comparison with available experimental data validates the model and force field employed and clarifies the molecular organization of this important liquid crystal often used as reference smectic material.

Predicting the anchoring of liquid crystals at a solid surface: 5-cyanobiphenyl on cristobalite and glassy silica surfaces of increasing roughness.

We employ atomistic molecular dynamics simulations to predict the alignment and anchoring strength of a typical nematic liquid crystal, 4-n-pentyl-4'-cyano biphenyl (5CB), on different forms of silica. In particular, we study a thin (~20 nm) film of 5CB supported on surfaces of crystalline (cristobalite) and amorphous silica of different roughness. We find that the orientational order at the surface and the anchoring strength depend on the morphology of the silica surface and its roughness. Cristobalite yields a uniform planar orientation and increases the order at the surface with respect to the bulk whereas amorphous glass has a disordering effect. Despite the low order at the amorphous surfaces, a planar orientation is established with a persistence length into the film higher than the one obtained for cristobalite.

Efficient analysis of highly complex nuclear magnetic resonance spectra of flexible solutes in ordered liquids by using molecular dynamics.

The NMR spectra of n-pentane as solute in the liquid crystal 5CB are measured at several temperatures in the nematic phase. Atomistic molecular dynamics simulations of this system are carried out to predict the dipolar couplings of the orientationally ordered pentane, and the spectra predicted from these simulations are compared with the NMR experimental ones. The simulation predictions provide an excellent starting point for analysis of the experimental NMR spectra using the covariance matrix adaptation evolutionary strategy. This shows both the power of atomistic simulations for aiding spectral analysis and the success of atomistic molecular dynamics in modeling these anisotropic systems.