Equilibration of linear polyethylene melts with pre-defined molecular weight distributions employing united atom Monte Carlo simulations.

Possessing control over the molecular size (molecular weight/chain length/degree of polymerization) distribution of a polymeric material is extremely important in applications. This is manifested de facto by the extensive contemporary scientific literature on processes for controlling this distribution experimentally. Yet, the literature on computational techniques for achieving prescribed molecular size distributions in simulations and exploring their impact on properties is much less abundant than its experimental/technical counterpart. Here, we develop-on the basis of united atom melt simulations employing connectivity-altering Monte Carlo moves-a new Metropolis selection criterion that drives the multichain system to a prescribed but otherwise arbitrary distribution of molecular sizes. The new formulation is a generalization of that originally proposed [P. V. K. Pant and D. N. Theodorou, Macromolecules 28, 7224 (1995)], but simpler and more computationally efficient. It requires knowledge solely of the target distribution, which need not be normalized. We have implemented the new formulation on long-chain linear polyethylene melts, obtaining excellent results. The target molecular size distribution can be provided in tabulated form, allowing absolute freedom as to the types of chain size profiles that can be simulated. Distributions for which equilibration has been achieved here for linear polyethylene include a truncated most probable, a truncated Schulz-Zimm, an arbitrary one defined in tabulated form, a broad truncated Gaussian, and a bimodal Gaussian. The last two are comparable to those encountered in industrial applications. The impact of the molecular size distribution on the properties of the simulated melts, such as density, chain dimensions, and mixing thermodynamics, is explored.

A coarse-grained model for capturing the helical behavior of isotactic polypropylene.

Understanding the process-property relations of helical polymers using molecular simulations has been an attractive research field over the years. Specifically, isotactic polypropylene still remains a challenge for current computational experimentation, as it exhibits phenomena such as crystallization that emerge on large spatial and temporal scales. Coarse-graining is an efficient technique for approaching such phenomena, although previous coarse-grained models lack in preserving important atomistic and structural details. In this paper we develop a new coarse-grained model, based on the popular MARTINI force field, that is able to reproduce the helical behavior of isotactic polypropylene. To test the model, the predicted statistical and structural properties (characteristic ratio, density, entanglement molecular weight, solubility parameter in the melt) are compared with previous simulation results and available experimental data. For the development of the new coarse-grained force field, a single unperturbed chain Monte Carlo algorithm has been implemented: an efficient algorithm which samples conformations representative of a melt by simulating just a single chain.

Multiscale simulations of polyzwitterions in aqueous bulk solutions and brush array configurations.

Zwitterionic polymers are very promising candidates for antifouling materials that exhibit high chemical stability as compared to polyethylene glycol-based systems. A number of simulation and experimental studies have emerged over recent years for the investigation of sulfobetaine-based zwitterionic polymers. Investigating the structural and thermodynamic properties of such polymers requires access to broad time and length regimes, thus necessitating the development of multiscale simulation strategies. The present article advocates a mesoscopic dissipative particle dynamics (DPD) model capable of addressing a wide range of time and length scales. The mesoscopic force field was developed hand-in-hand with atomistic simulations based on the OPLS force field through a bottom-up parameterization procedure that matches the atomistically calculated strand-length, strand-angle and pair distribution functions. The DPD model is validated against atomistic simulations conducted in this work, and against relevant atomistic simulation studies, theoretical predictions and experimental correlations from the literature. Properties examined include the conformations of SPE polymers in dilute bulk aqueous solution, the density profile and thickness of brush arrays as functions of the grafting density and chain length. In addition, we compute the potential of mean force of an approaching hydrophilic or hydrophobic foulant via umbrella sampling as a function of its position relative to the poly-zwitterion-covered surface. The aforementioned observables lead to important insights regarding the conformational tendencies of grafted polyzwitterions and their antifouling properties.

Structure and thermodynamics of grafted silica/polystyrene dilute nanocomposites investigated through self-consistent field theory.

Polymer/matrix nanocomposites (PNCs) are materials with exceptional properties. They offer a plethora of promising applications in key industrial sectors. In most cases, it is preferable to disperse the nanoparticles (NPs) homogeneously across the matrix phase. However, under certain conditions NPs might lump together and lead to a composite material with undesirable properties. A common strategy to stabilize the NPs is to graft on their surface polymer chains of the same chemical constitution as the matrix chains. There are several unresolved issues concerning the optimal molar mass and areal density of grafted chains that would ensure best dispersion, given the nanoparticles and the polymer matrix. We propose a model for the prediction of key structural and thermodynamic properties of PNC and apply it to a single spherical silica (SiO2) nanoparticle or planar surface grafted with polystyrene chains embedded at low concentration in a matrix phase of the same chemical constitution. Our model is based on self-consistent field theory, formulated in terms of the Edwards diffusion equation. The properties of the PNC are explored across a broad parameter space, spanning the mushroom regime (low grafting densities, small NPs and chain lengths), the dense brush regime, and the crowding regime (large grafting densities, NP diameters, and chain lengths). We extract several key quantities regarding the distributions and the configurations of the polymer chains, such as the radial density profiles and their decomposition into contributions of adsorbed and free chains, the chains/area profiles, and the tendency of end segments to segregate at the interfaces. Based on our predictions concerning the brush thickness, we revisit the scaling behaviors proposed in the literature and we compare our findings with experiment, relevant simulations, and analytic models, such as Alexander's model for incompressible brushes.

Molecular dynamics simulations of EPON-862/DETDA epoxy networks: structure, topology, elastic constants, and local dynamics.

Structural, topological, mechanical and dynamical properties of EPON-862/DETDA epoxy networks are investigated with Molecular Dynamics (MD) simulations. The epoxy networks are composed of the resin Diglycidyl Ether Bisphenol F (DGEBF), also known as EPON-862, and the hardener Diethyl Toluene Diamine (DETDA). Systems with four different crosslinking degrees are examined; the effect of the degree of crosslinking on studied properties is thus determined. The computed quantities are retrieved by employing several simulation strategies and numerical methods of statistical mechanics in order to gain a rigorous and solid understanding of the aforementioned properties as well as to assess the accuracy and applicability of the methods employed. We quantify and analyze the local structure of the EPON-862/DETDA epoxy networks through the partial pair distribution functions, the Faber-Ziman partial structure factors and through simulated X-ray diffraction patterns, demonstrating good agreement with an experimental spectrum from a similar epoxy resin. The topology of the networks is examined with the aim of assessing percolation of connectivity, the properties of network fragments (subnetworks), and the distribution of functionalities of the crosslinks. The elastic constants of the systems are retrieved by employing two equilibrium (analysis of volume fluctuations, Parrinello-Rahman strain fluctuation relation) and one nonequilibrium (uniaxial tension/compression deformations at prescribed rate) method. Finally, the glass temperatures of the systems are estimated by calculating the density as a function of temperature and by analyzing the reorientational dynamics of bond vectors which describe relaxation processes at the segment level.

SimNano: A Trust Region Strategy for Large-Scale Molecular Systems Energy Minimization Based on Exact Second-Order Derivative Information.

In this work, a new energy minimization strategy is presented that achieves better convergence properties than the standard algorithms employed in the field (fewer steps and usually a lower minimum) and is also computationally efficient; therefore, it becomes suitable for dealing with large-scale molecular systems. The proposed strategy is integrated into the SimNano energy minimization platform that is also described herein. SimNano relies on the analytical calculation of the molecular systems' gradient vectors and Hessian matrices using the computational modeling framework proposed by the authors ( Chatzieleftheriou , S. ; Adendorff , M. R. ; Lagaros , N. D. Generalized Potential Energy Finite Elements for Modeling Molecular Nanostructures . J. Chem. Inf. Model. 2016 , 56 ( 10 ), 1963 - 1978 ). The basis of the proposed minimization strategy is a trust region algorithm based on exact second-order derivative information. Taking advantage of the Hessian matrices' sparsity, a specialized treatment of the data structure is implemented. The latter is beneficial and often rather necessary, especially in the case of large-scale molecular systems, improving the speed and reducing the memory requirements. In order to demonstrate the efficiency of the proposed energy minimization strategy, several test examples are examined, and the results achieved are compared with those obtained by one of the most popular molecular simulation software packages, i.e., the Large-Scale Atomic/Molecular Massively Parallel Simulator (LAMMPS). The results indicate that the proposed minimization strategy exhibits superior convergence properties compared with the typical algorithms (i.e., nonlinear conjugate gradient algorithm, limited-memory Broyden-Fletcher-Goldfarb-Shanno (LBFGS) algorithm, etc.). The SimNano energy minimization platform can be downloaded from the site http://users.ntua.gr/nlagaros/simnano.html , enabling researchers in the field to build molecular systems and perform energy minimization runs using input files in LAMMPS format.

Exploring the interactions of irbesartan and irbesartan-2-hydroxypropyl-ß-cyclodextrin complex with model membranes.

The interactions of irbesartan (IRB) and irbesartan-2-hydroxypropyl-ß-cyclodextrin (HP-ß-CD) complex with dipalmitoyl phosphatidylcholine (DPPC) bilayers have been explored utilizing an array of biophysical techniques ranging from differential scanning calorimetry (DSC), small angle X-ray scattering (SAXS), ESI mass spectrometry (ESI-MS) and solid state nuclear magnetic resonance (ssNMR). Molecular dynamics (MD) calculations have been also conducted to complement the experimental results. Irbesartan was found to be embedded in the lipid membrane core and to affect the phase transition properties of the DPPC bilayers. SAXS studies revealed that irbesartan alone does not display perfect solvation since some coexisting irbesartan crystallites are present. In its complexed form IRB gets fully solvated in the membranes showing that encapsulation of IRB in HP-ß-CD may have beneficial effects in the ADME properties of this drug. MD experiments revealed the topological and orientational integration of irbesartan into the phospholipid bilayer being placed at about 1nm from the membrane centre.

Development of a Meshless Kernel-Based Scheme for Particle-Field Brownian Dynamics Simulations.

We develop a meshless discretization scheme for particle-field Brownian dynamics simulations. The density is assigned on the particle level using a weighting kernel with finite support. The system's free energy density is derived from an equation of state (EoS) and includes a square gradient term. The numerical stability of the scheme is evaluated in terms of reproducing the thermodynamics (equilibrium density and compressibility) and dynamics (diffusion coefficient) of homogeneous samples. Using a reduced description to simplify our analysis, we find that numerical stability depends strictly on reduced reference compressibility, kernel range, time step in relation to the friction factor, and reduced external pressure, the latter being relevant under isobaric conditions. Appropriate parametrization yields precise thermodynamics, further improved through a simple renormalization protocol. The dynamics can be restored exactly through a trivial manipulation of the time step and friction coefficient. A semiempirical formula for the upper bound on the time step is derived, which takes into account variations in compressibility, friction factor, and kernel range. We test the scheme on realistic mesoscopic models of fluids, involving both simple (Helfand) and more sophisticated (Sanchez-Lacombe) equations of state.

Interactions of a Novel Anthracycline with Oligonucleotide DNA and Cyclodextrins in an Aqueous Environment.

Berubicin, a chemotherapy medication belonging to the class of anthracyclines, is simulated in double-stranded DNA sequences and cyclodextrins in an aqueous environment via full-atom molecular dynamics simulations on the time scale of microseconds. The drug is studied in both the neutral and protonated states so as to better comprehend the role of its charge in the formed complexes. The noncovalent berubicin-DNA and berubicin-cyclodextrin complexes are investigated in detail, paying special attention to their thermodynamic description by employing the double decoupling method, the solvent balance method, the weighted solvent accessible surface model, and the linear interaction energy method. A novel approach for extracting the desolvation thermodynamics of the binding process is also presented. Both the binding and desolvation Gibbs energies are decomposed into entropic and enthalpic contributions so as to elucidate the nature of complexation and its driving forces. Selected structural and geometrical properties of all the complexes, which are all stable, are analyzed. Both cyclodextrins under consideration are widely utilized for drug delivery purposes, and a comparative investigation between their bound states with berubicin is carried out.

Buckling kinetics of graphene membranes under uniaxial compression.

Despite past investigations of the buckling instability, the kinetics of the buckling process is not well understood. We develop a generic framework for determining the buckling kinetics of membranes under compressive stress (σ_{b}) via molecular dynamics simulations. The buckling time (t_{b}) is modeled by an extended Boltzmann-Arrhenius-Zhurkov equation accounting for temperature (T) and scale-dependent bending rigidity. We discern three regimes: (I) t_{b} decreases with T; (II) t_{b} increases with T; (III) t_{b} is T independent. Regime II coheres with the predictions of the theory of fluctuating sheets (TFS). Regime I is seen at small scales due to fluctuations about equilibrium and is not predicted by the TFS.

Measuring Oxygen Solubility in Amorphous and Semicrystalline Polyolefins Using Test Particle Insertion: A Comparative Study of Polyethylene and Isotactic Polypropylene.

The test particle insertion method is used to study the solubility of oxygen in two commonly used polymers: polyethylene (PE) and isotactic polypropylene (iPP). Amorphous samples for both polymers were prepared by means of Monte Carlo and molecular dynamics simulations, and the oxygen solubility was measured across different temperatures. The solubility-temperature dependence for iPP proved to be nonmonotonic due to the interplay between binding and reorganizational enthalpy, while for PE, it appeared to be monotonic based on the available data in the studied temperature range. A broad comparison to experiments and simulations is included. Further oxygen insertions were also performed in semicrystalline PE and iPP samples at ambient temperature, and the obtained values were compared to a linear relationship which correlates the solubility in the purely amorphous phase with the solubility in the crystalline phase. The solubility of PE closely follows the linear relationship, while iPP exhibits some divergence. All the semicrystalline samples were previously annealed at elevated temperatures for long periods (a few µs), and a strong effect of annealing was observed on the structure and the solubility of iPP. A well-developed iPP lamellar structure emerged at longer annealing times, while PE develops that structure already in the early crystallization stages. The solubility of semicrystalline iPP samples with lamellar morphology exhibited better agreement with extrapolated solubility values of the amorphous state─the extrapolation was made using a linear relationship connecting solubility in the purely amorphous phase and solubility in mixed phases (amorphous and crystalline). Results on the correlation of the solubility with the local structural ordering are also present.

Temporal disconnectivity of the energy landscape in glassy systems.

An alternative graphical representation of the potential energy landscape (PEL) has been developed and applied to a binary Lennard-Jones glassy system, providing insight into the unique topology of the system's potential energy hypersurface. With the help of this representation one is able to monitor the different explored basins of the PEL, as well as how--and mainly when--subsets of basins communicate with each other via transitions in such a way that details of the prior temporal history have been erased, i.e., local equilibration between the basins in each subset has been achieved. In this way, apart from detailed information about the structure of the PEL, the system's temporal evolution on the PEL is described. In order to gather all necessary information about the identities of two or more basins that are connected with each other, we consider two different approaches. The first one is based on consideration of the time needed for two basins to mutually equilibrate their populations according to the transition rate between them, in the absence of any effect induced by the rest of the landscape. The second approach is based on an analytical solution of the master equation that explicitly takes into account the entire explored landscape. It is shown that both approaches lead to the same result concerning the topology of the PEL and dynamical evolution on it. Moreover, a "temporal disconnectivity graph" is introduced to represent a lumped system stemming from the initial one. The lumped system is obtained via a specially designed algorithm [N. Lempesis, D. G. Tsalikis, G. C. Boulougouris, and D. N. Theodorou, J. Chem. Phys. 135, 204507 (2011)]. The temporal disconnectivity graph provides useful information about both the lumped and the initial systems, including the definition of "metabasins" as collections of basins that communicate with each other via transitions that are fast relative to the observation time. Finally, the two examined approaches are compared to an "on the fly" molecular dynamics-based algorithm [D. G. Tsalikis, N. Lempesis, G. C. Boulougouris, and D. N. Theodorou, J. Chem. Theory Comput. 6, 1307 (2010)].

On solving the master equation in spatially periodic systems.

We present a new method for solving the master equation for a system evolving on a spatially periodic network of states. The network contains 2(ν) images of a "unit cell" of n states, arranged along one direction with periodic boundary conditions at the ends. We analyze the structure of the symmetrized (2(ν)n) × (2(ν)n) rate constant matrix for this system and derive a recursive scheme for determining its eigenvalues and eigenvectors, and therefore analytically expressing the time-dependent probabilities of all states in the network, based on diagonalizations of n × n matrices formed by consideration of a single unit cell. We apply our new method to the problem of low-temperature, low-occupancy diffusion of xenon in the zeolite silicalite-1 using the states, interstate transitions, and transition state theory-based rate constants previously derived by June et al. [J. Phys. Chem. 95, 8866 (1991)]. The new method yields a diffusion tensor for this system which differs by less than 3% from the values derived previously via kinetic Monte Carlo (KMC) simulations and confirmed by new KMC simulations conducted in the present work. The computational requirements of the new method are compared against those of KMC, numerical solution of the master equation by the Euler method, and direct molecular dynamics. In the problem of diffusion of xenon in silicalite-1, the new method is shown to be faster than these alternative methods by factors of about 3.177 × 10(4), 4.237 × 10(3), and 1.75 × 10(7), respectively. The computational savings and ease of setting up calculations afforded by the new method of master equation solution by recursive reduction of dimensionality in diagonalizing the rate constant matrix make it attractive as a means of predicting long-time dynamical phenomena in spatially periodic systems from atomic-level information.

Molecular simulations of fluoxetine in hydrated lipid bilayers, as well as in aqueous solutions containing ß-cyclodextrin.

Fluoxetine, which is a well-known antidepressant drug, is studied in hydrated cholesterol-free and cholesterol-containing lipid bilayers through unbiased and biased atomistic molecular dynamics simulations. The latter are conducted for the calculation of the potential of mean force (PMF) of fluoxetine along an axis perpendicular to the two leaflets of the bilayer. The PMF indicates that the drug prefers to reside inside the lipid phase and allows us to calculate important thermodynamic properties, such as the Gibbs energy difference of partitioning from the water to the lipid phase and the Gibbs energy barrier for hopping events between the two leaflets of the bilayer. The results from the biased simulations are in accord with the mass density profiles calculated from the unbiased simulations. Moreover, we estimate the effect of fluoxetine mole fraction on the order parameters of the lipid alkyl chains and on the area per lipid. It is also found that fluoxetine forms a hydrogen bond network with lipids and water molecules penetrating into the lipid phase. In addition, fluoxoetine is studied in detail in aqueous solutions containing ß-cyclodextrin. It is observed from unbiased molecular dynamics simulations that the two aforementioned molecules form a noncovalent complex spontaneously and the calculated binding free energy is in agreement with the literature.

Solvation Free Energy of Dilute Grafted (Nano)Particles in Polymer Melts via the Self-Consistent Field Theory.

Predicting the distribution of a chemical species across multiple phases is of critical importance to environmental protection, pharmaceuticals, and high added-value chemicals. Computationally, this problem is addressed by determining the free energy of solvation of the species in the different phases using a well-established thermodynamic formulation. Following recent developments in sterically stabilized colloids and nanocomposite materials, the solvation of polymer-grafted nanoparticles in different solvents or polymer melts has become relevant. We develop a Self-Consistent Field theoretical framework to determine the solvation free energy of grafted particles inside a molten polymer matrix phase at low concentrations. The solvation free energy is calculated based on the notion of a pseudochemical potential introduced by Ben-Naim. Grafted and matrix chains are taken to be of the same chemical constitution, but their lengths are varied systematically, as are the particle radius and the areal density of grafted chains. In addition, different affinities between the nanoparticle core and the polymer (contact angles) are considered. At very low or very high amounts of grafted material, solvation depends on the adhesion tension between the bare particle and the matrix or on the surface tension of the grafted polymer, respectively. The dependence of the solvation free energy on molecular characteristics is more complicated at intermediate grafting densities and high curvatures, where the contribution of the entropy of grafted chains becomes significant. In general, solvation is less favored in cases where the matrix chains are much shorter than the grafted ones. The former tend to penetrate and swell the brush, thus generating conformational and translational entropy penalties. This effect becomes more pronounced when considering large particles since the grafted chains have less available space and extend more. For extremely low amounts of grafted material, we observe the opposite trend, albeit weak. Based on our calculations, we propose a generic model for estimating the solvation free energies of grafted nanoparticles in polymer melts from their molecular characteristics. The model and associated SCF formulation, illustrated here for chemically identical grafted and matrix chains, can be extended to obtain partition coefficients of grafted nanoparticles between different polymer melts.

Lumping analysis for the prediction of long-time dynamics: from monomolecular reaction systems to inherent structure dynamics of glassy materials.

In this work we develop, test, and implement a methodology that is able to perform, in an automated manner, "lumping" of a high-dimensional, discrete dynamical system onto a lower-dimensional space. Our aim is to develop an algorithm which, without any assumption about the nature of the system's slow dynamics, is able to reproduce accurately the long-time dynamics with minimal loss of information. Both the original and the lumped systems conform to master equations, related via the "lumping" analysis introduced by Wei and Kuo [Ind. Eng. Chem. Fundam. 8, 114 (1969)], and have the same limiting equilibrium probability distribution. The proposed method can be used in a variety of processes that can be modeled via a first order kinetic reaction scheme. Lumping affords great savings in the computational cost and reveals the characteristic times governing the slow dynamics of the system. Our goal is to approach the best lumping scheme with respect to three criteria, in order for the lumped system to be able to fully describe the long-time dynamics of the original system. The criteria used are: (a) the lumping error arising from the reduction process; (b) a measure of the magnitude of singular values associated with long-time evolution of the lumped system; and (c) the size of the lumped system. The search for the optimum lumping proceeds via Monte Carlo simulation based on the Wang-Landau scheme, which enables us to overcome entrapment in local minima in the above criteria and therefore increases the probability of encountering the global optimum. The developed algorithm is implemented to reproduce the long-time dynamics of a glassy binary Lennard-Jones mixture based on the idea of "inherent structures," where the rate constants for transitions between inherent structures have been evaluated via hazard plot analysis of a properly designed ensemble of molecular dynamics trajectories.

In silico study of levodopa in hydrated lipid bilayers at the atomistic level.

This article presents atomistic molecular dynamics and umbrella sampling simulations of levodopa at various concentrations in hydrated cholesterol-free 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) and cholesterol-containing 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayers. Levodopa is the standard medication for Parkinson's disease and is marketed under various trade names; in the context of this article, the levodopa molecule is mostly studied in its zwitterionic form but some results concerning the neutral levodopa are presented as well for comparison purposes. The motivation is to study in detail how levodopa behaves in different hydrated lipid membranes, primarily from the thermodynamic point of view, and reveal aspects of mechanism of its permeation through them. Dependencies of properties on the levodopa concentration are also investigated. Special attention is paid to the calculation of mass density profiles, order parameters and self-diffusion coefficients. Levodopa zwitterions, which form a hydrogen bond network with water and phospholipid molecules, are found to be preferentially located at the water/lipid interface, as well as in the aqueous phase surrounding the cholesterol-free and cholesterol-containing bilayers. This is concluded from the potentials of mean force calculated by umbrella sampling simulations as levodopa is transferred from the lipid to the aqueous phase along an axis perpendicular to the two leaflets of the membranes.

Potential of Mean Force between Bare or Grafted Silica/Polystyrene Surfaces from Self-Consistent Field Theory.

We investigate single and opposing silica plates, either bare of grafted, in contact with vacuum or melt phases, using self-consistent field theory. Solid-polymer and solid-solid nonbonded interactions are described by means of a Hamaker potential, in conjunction with a ramp potential. The cohesive nonbonded interactions are described by the Sanchez-Lacombe or the Helfand free energy densities. We first build our thermodynamic reference by examining single surfaces, either bare or grafted, under various wetting conditions in terms of the corresponding contact angles, the macroscopic wetting functions (i.e., the work of cohesion, adhesion, spreading and immersion), the interfacial free energies and brush thickness. Subsequently, we derive the potential of mean force (PMF) of two approaching bare plates with melt between them, each time varying the wetting conditions. We then determine the PMF between two grafted silica plates separated by a molten polystyrene film. Allowing the grafting density and the molecular weight of grafted chains to vary between the two plates, we test how asymmetries existing in a real system could affect steric stabilization induced by the grafted chains. Additionally, we derive the PMF between two grafted surfaces in vacuum and determine how the equilibrium distance between the two grafted plates is influenced by their grafting density and the molecular weight of grafted chains. Finally, we provide design rules for the steric stabilization of opposing grafted surfaces (or fine nanoparticles) by taking account of the grafting density, the chain length of the grafted and matrix chains, and the asymmetry among the opposing surfaces.

Diffusion via space discretization method to study the concentration dependence of self-diffusivity under confinement.

The concentration dependence of self-diffusivity is investigated by means of a novel method, extending our previously developed second-order Markov process model to periodic media. Introducing the concept of minimum-crossing surface, we obtain a unique decomposition of the self-diffusion coefficient into two parameters with specific physical meanings. Two case studies showing a maximum in self-diffusivity as a function of concentration are investigated, along with two cases where such a maximum cannot be present. Subsequently, the method is applied to the large cavity pore network of the ITQ-1 (Mobil tWenty tWo, MWW) zeolite for methane (displaying a maximum in self-diffusivity) and carbon dioxide (no maximum), explaining the diffusivity trend on the basis of the evolution of the model parameters as a function of concentration.

Kinetic concepts and local failure in the interfacial shear strength of epoxy-graphene nanocomposites.

Interfacial shear strength (IFSS) is a key property in the design of composites and nanocomposites. Many simulation studies quantify the interfacial characteristics of sandwichlike specimens in terms of the IFSS and pullout force; a common feature of these studies is that they employ finite model systems and are therefore subject to strong finite size effects. We propose an alternative approach which is applicable to both aperiodic and periodic computational specimens. The interfaces are subjected to multiple shear deformation simulations over a wide range of temperatures (T) and shear stresses (σ_{zx}). From these simulations we collect the failure times (t_{f}); by analyzing them in the framework of an extended Boltzmann-Arrhenius-Zhurkov kinetic equation we derive the IFSS, the limiting stress for barrierless transitions, the activation energy, the activation volume for failure, the sliding velocities, and a local elastic shear modulus for the interface. We test our methodology on epoxy diglycidyl ether bisphenol F-diethyl toluene diamine interfaces in contact with (i) pristine graphene, (ii) graphene with single-atom vacancies, and (iii) graphene with hydroxyl-ΟΗ groups. Differences in the mechanism of interfacial failure among these three systems are discussed.