Molecular dynamics simulation of organic-inorganic nanocomposites: Layering behavior and interlayer structure of organoclays

Isothermal-isobaric (NPT) molecular dynamics simulation has been performed to investigate the layering behavior and structure of nanoconfined quaternary alkylammoniums in organoclays. This work is focused on systems consisting of two clay layers and a number of alkylammoniums, and involves the use of modified Dreiding force field. The simulated basal spacings of organoclays agree satisfactorily with the experimental results in the literature. The atomic density profiles in the direction normal to the clay surface indicate that the alkyl chains within the interlayer space of montmorillonite exhibit an obvious layering behavior. The headgroups of long alkyl chains are distributed within two layers close to the clay surface, whereas the distributions of methyl and methylene groups are strongly dependent on the alkyl chain length and clay layer charge. Monolayer, bilayer, and pseudo-trilayer structures are found in organoclays modified with single long alkyl chains, which are identical to the structural models based on the measured basal spacings. A pseudo-quadrilayer structure, for the first time to our knowledge, is also identified in organoclays with double long alkyl chains. In the mixture structure of paraffin-type and multilayer, alkyl chains do not lie flat within a single layer but interlace, and also jump to the next layer in pseudo-trilayer as well as next nearest layer in pseudo-quadrilayer.

A novel Mo-W interlayer approach for CVD diamond deposition on steel

Steel is the most widely used material in engineering for its cost/performance ratio and coatings are routinely applied on its surface to further improve its properties. Diamond coated steel parts are an option for many demanding industrial applications through prolonging the lifetime of steel parts, enhancement of tool performance as well as the reduction of wear rates. Direct deposition of diamond on steel using conventional chemical vapour deposition (CVD) processes is known to give poor results due to the preferential formation of amorphous carbon on iron, nickel and other elements as well as stresses induced from the significant difference in the thermal expansion coefficients of those materials. This article reports a novel approach of deposition of nanocrystalline diamond coatings on high-speed steel (M42) substrates using a multi-structured molybdenum (Mo) - tungsten (W) interlayer to form steel/Mo/Mo-W/W/diamond sandwich structures which overcome the adhesion problem related to direct magnetron sputtering deposition of pure tungsten. Surface, interface and tribology properties were evaluated to understand the role of such an interlayer structure. The multi-structured Mo-W interlayer has been proven to improve the adhesion between diamond films and steel substrates by acting as an effective diffusion barrier during the CVD diamond deposition.

Nanocomposites: the contribution of multiscale modelling to improve the dispersion levels

Tese de Doutoramento em Ciência e Engenharia de Polímeros e Compósitos

Molecular dynamics simulation of the structural and dynamic properties of dioctadecyldimethyl ammoniums in organoclays

The structural and dynamic properties of dioctadecyldimethylammoniums (DODDMA) intercalated into 2:1 layered clays are investigated using isothermal-isobaric (NPT) molecular dynamics (MD) simulation. The simulated results are in reasonably good agreement with the available experimental measurements, such as X-ray diffraction (XRD), atom force microscopy (AFM), Fourier transform infrared (FTIR), and nuclear magnetic resonance (NMR) spectroscopies. The nitrogen atoms are found to be located mainly within two layers close to the clay surface whereas methylene groups form a pseudoquadrilayer structure. The results of tilt angle and order parameter show that interior two-bond segments of alkyl chains prefer an arrangement parallel to the clay surface, whereas the segments toward end groups adopt a random orientation. In addition, the alkyl chains within the layer structure lie almost parallel to the clay surface whereas those out of the layer structure are essentially perpendicular to the surface. The trans conformations are predominant in all cases although extensive gauche conformations are observed, which is in agreement with previous simulations on n-butane. Moreover, an odd-even effect in conformation distributions is observed mainly along the chains close to the head and tail groups. The diffusion constants of both nitrogen atoms and methylene groups in these nanoconfined alkyl chains increase with the temperature and methelene position toward the tail groups.

Synthesis, structure and magnetic characterisation of a new layered ammonium manganese(II) diphosphate hydrate, (NH4)(2)[Mn-3(P2O7)(2)(H2O)(2)]

A new layered ammonium manganese(II) diphosphate, (NH4)(2)[Mn-3(P2O7)(2)(H2O)(2)], has been synthesised under solvothermal conditions at 433 K in ethylene glycol and the structure determined at 293 K using single-crystal X-ray diffraction data (M-r = 584.82, monoclinic, space group P2(1)/a, a = 9.4610( 8), b = 8.3565( 7), c = 9.477(1) Angstrom, beta = 99.908(9) degrees, V = 738.07 Angstrom(3), Z = 2, R = 0.0351 and R-w = 0.0411 for 1262 observed data (I > 3(sigma(I))). The structure consists of chains of cis- and trans-edge sharing MnO6 octahedra linked via P2O7 units to form layers of formula [Mn3P4O14(H2O)(2)](2-) in the ab plane. Ammonium ions lie between the manganese-diphosphate layers. A network of interlayer and ammonium-layer based hydrogen bonding holds the structure together. Magnetic measurements indicate Curie - Weiss behaviour above 30 K with mu(eff) = 5.74(1) mu(B) and theta = -23(1) K, consistent with the presence of high-spin Mn2+ ions and antiferromagnetic interactions. However, the magnetic data reveal a spontaneous magnetisation at 5 K, indicating a canting of Mn2+ moments in the antiferromagnetic ground state. On heating (NH4)(2)[Mn-3(P2O7)(2)(H2O)(2)] in water at 433 K under hydrothermal conditions, Mn-5(HPO4)(2)(PO4)(2).4H(2)O, synthetic hureaulite, is formed.

The Structure of the chiral Pt{531} surface: a combined LEED and DFT study

The structure of the chiral kinked Pt{531} surface has been determined by low-energy electron diffraction intensity-versus-energy (LEED-IV) analysis and density functional theory (DFT). Large contractions and expansions of the vertical interlayer distances with respect to the bulk-terminated surface geometry were found for the first six layers (LEED: d(12) = 0.44 angstrom, d(23) = 0.69 angstrom, d(34) = 0.49 angstrom, d(45) = 0.95 angstrom, d(56) = 0.56 angstrom; DFT: d(12) = 0.51 angstrom, d(23) = 0.55 angstrom, d(34) = 0.74 angstrom, d(45) = 0.78 angstrom, d(56) = 0.63 angstrom; d(bulk) = 0.66 angstrom). Energy-dependent cancellations of LEED spots over unusually large energy ranges, up to 100 eV, can be explained by surface roughness and reproduced by applying a model involving 0.25 ML of vacancies and adatoms in the scattering calculations. The agreement between the results from LEED and DFT is not as good as in other cases, which could be due to this roughness of the real surface.

The surface geometries of the medium and high coverage carbon monoxide structures c(2 × 4)–(2CO) and p(root7- x root7)R19°–(4CO) structure

The surface geometries of the p (root7- x root7)R19degrees-(4CO) and c(2 x 4)-(2CO) layers on Ni {111} and the clean Ni {111} surface were determined by low energy electron diffraction structure analysis. For the clean surface small but significant contractions of d(12) and d(23) (both 2.02 Angstrom) were found with respect to the bulk interlayer distance (2.03 Angstrom). In the c(2 x 4)-(2CO) structure these distances are expanded, with values of d(12) = 2.08 Angstrom and d(23) = 2.06 Angstrom and buckling of 0.08 and 0.02 Angstrom, respectively, in the first and second layer. CO resides near hcp and fcc hollow sites with relatively large lateral shifts away from the ideal positions leading to unequal C-Ni bond lengths between 1.76 and 1.99 Angstrom. For the p(root7- x root7-)R19'-(4CO) layer two best fit geometries were found, which agree in most of their atomic positions, except for one out of four CO molecules, which is either near atop or between bridge and atop. The remaining three molecules reside near hcp and fcc sites, again with large lateral deviations from their ideal positions. The average C Ni bond length for these molecules is, however, the same as for CO on hollow sites at low coverage. The average CNi bond length at hollow sites, the interlayer distances, and buckling in the first Ni layer are similar to the c(2 x 4)(2CO) geometry, only the buckling in the second layer (0.08 Angstrom) is significantly larger. Lateral and vertical shifts of the Ni atoms in the first layer lead to unsymmetric environments for the CO molecules, which can be regarded as an imprint of the chiral p(root7- x root7-)R19degrees lattice geometry onto the substrate.

Electronic structure of Pd multi-layers on Re(0001): the role of charge transfer

Understanding the origin of the properties of metal-supported metal thin films is important for the rational design of bimetallic catalysts and other applications, but it is generally difficult to separate effects related to strain from those arising from interface interactions. Here we use density functional (DFT) theory to examine the structure and electronic behavior of few-layer palladium films on the rhenium (0001) surface, where there is negligible interfacial strain and therefore other effects can be isolated. Our DFT calculations predict stacking sequences and interlayer separations in excellent agreement with quantitative low-energy electron diffraction experiments. By theoretically simulating the Pd core-level X-ray photoemission spectra (XPS) of the films, we are able to interpret and assign the basic features of both low-resolution and high-resolution XPS measurements. The core levels at the interface shift to more negative energies, rigidly following the shifts in the same direction of the valence d-band center. We demonstrate that the valence band shift at the interface is caused by charge transfer from Re to Pd, which occurs mainly to valence states of hybridized s-p character rather than to the Pd d-band. Since the d-band filling is roughly constant, there is a correlation between the d-band center shift and its bandwidth. The resulting effect of this charge transfer on the valence d-band is thus analogous to the application of a lateral compressive strain on the adlayers. Our analysis suggests that charge transfer should be considered when describing the origin of core and valence band shifts in other metal / metal adlayer systems.

Pore structure characterization of kaolin, metakaolin, and their acid-treated products using small-angle X-ray scattering

Pore structure of dealuminated kaolin and metakaolin was studied by small-angle X-ray scattering (SAXS). Both parent kaolin and metakaolin have about 10% of the total pore volume provided by globular pores with 105 Å mean pore size. Their surface area is about 14 m2/g. Acid dealumination of kaolin causes an increase of its globular pore volume without an appreciable change in the mean pore size, its surface area increasing up to about 90 m2/g. Acid dealumination of metakaolin enhances the globular pore volume, although there is generation of slit-shaped pores with a narrow thickness distribution whose mean value is 14 Å. This interlayer spacing causes an increase in surface area of about 190 m2/g by SAXS. © 1994.

Synthesis, structure, and thermal stability of poly(methyl methacrylate)-co- poly(3-tri(methoxysilyil)propyl methacrylate)/ montmorillonite nanocomposites

The structure and the thermodegradation behavior of both poly(methyl methacrylate)-co-poly(3-tri(methoxysilyil)propyl methacrylate) polymer modified with silyl groups and of intercalated poly(methyl methacrylate)-co-poly(3- tri(methoxysilyil)propyl methacrylate)/Cloisite 15A™ nanocomposite have been in situ probed. The structural feature were comparatively studied by Fourier transform infrared spectroscopy (FTIR), 13C and 29Si nuclear magnetic resonance (NMR), and small angle X-ray scattering (SAXS) measurements. The intercalation of polymer in the interlayer galleries was evidenced by the increment of the basal distance from 31 to 45 Å. The variation of this interlayer distance as function of temperature was followed by in situ SAXS. Pristine polymer decomposition pathway depends on the atmosphere, presenting two steps under air and three under N2. The nanocomposites are more stable than polymer, and this thermal improvement is proportional to the clay loading. The experimental results indicate that clay nanoparticles play several different roles in polymer stabilization, among them, diffusion barrier, charring, and suppression of degradation steps by chemical reactions between polymer and clay. Charring is atmosphere dependent, occurring more pronounced under air. © 2012 Society of Plastics Engineers.

Multi-scale micro-structure generation strategy for up-scaling transport in clays

Clays and claystones are used as backfill and barrier materials in the design of waste repositories, because they act as hydraulic barriers and retain contaminants. Transport through such barriers occurs mainly by molecular diffusion. There is thus an interest to relate the diffusion properties of clays to their structural properties. In previous work, we have developed a concept for up-scaling pore-scale molecular diffusion coefficients using a grid-based model for the sample pore structure. Here we present an operational algorithm which can generate such model pore structures of polymineral materials. The obtained pore maps match the rock’s mineralogical components and its macroscopic properties such as porosity, grain and pore size distributions. Representative ensembles of grains in 2D or 3D are created by a lattice Monte Carlo (MC) method, which minimizes the interfacial energy of grains starting from an initial grain distribution. Pores are generated at grain boundaries and/or within grains. The method is general and allows to generate anisotropic structures with grains of approximately predetermined shapes, or with mixtures of different grain types. A specific focus of this study was on the simulation of clay-like materials. The generated clay pore maps were then used to derive upscaled effective diffusion coefficients for non-sorbing tracers using a homogenization technique. The large number of generated maps allowed to check the relations between micro-structural features of clays and their effective transport parameters, as is required to explain and extrapolate experimental diffusion results. As examples, we present a set of 2D and 3D simulations and investigated the effects of nanopores within particles (interlayer pores) and micropores between particles. Archie’s simple power law is followed in systems with only micropores. When nanopores are present, additional parameters are required; the data reveal that effective diffusion coefficients could be described by a sum of two power functions, related to the micro- and nanoporosity. We further used the model to investigate the relationships between particle orientation and effective transport properties of the sample.

A new bidimensional compound containing (μ-thiocyanate)(bpy)copper(II) molecules: synthesis, crystal structure and magnetic properties of [Cu(bpy)(NCS)2]n (bpy = 2,2'-bipyridyl)

The crystal structure of the first bidimensional copper(II) compound containing only thiocyanate as bridging ligands [Cu(bpy)(NCS) 2 ] n , where bpy=2,2'-bipyridyl, has been determined by X-ray diffraction on single-crystals. Two different environments for both types of copper(II) ions in the unit cell are apparent: a distorted octahedron and a square pyramid. A bidimensional structure with a deformed honeycomb-layer motif is formed, the bipyridyl ligands filling the interlayer space. The magnetic susceptibility data of the compound have been investigated between 280 and 1.8 K. The compound presents a very weak antiferromagnetic interaction that has been fitted by using the Bleaney-Bowers expression for a dimeric unit, whereby a J value of -1.01(1) cm - 1 (H=-JS 1 .S 2 ) and a g value of 2.08(1) have been obtained.

The interlayer swelling and molecular packing in organoclays

Understanding the interlayer swelling and molecular packing in organoclays is important to the formation and design of polymer nanocomposites. This paper presents recent experimental and molecular simulation studies on a variety of organoclays that show a linear relationship between the increase of d-spacing and the mass ratio between organic and clay. A denser molecular packing is observed in organoclays containing surfactants with hydroxyl-ethyl units. Moreover, our simulation results show that the head (nitrogen) groups are essentially tethered to the clay surface while the long hydrocarbon chains tend to adopt a layering structure with disordered conformation, which contrasts with the previous assumptions of either the chains lying parallel to the clay surface or being tilted at rather precise angles. (c) 2005 Elsevier Inc. All rights reserved.

Adsorption in slit pores and pore-size distribution: A molecular layer structure theory

A new approach is developed to analyze the thermodynamic properties of a sub-critical fluid adsorbed in a slit pore of activated carbon. The approach is based on a representation that an adsorbed fluid forms an ordered structure close to a smoothed solid surface. This ordered structure is modelled as a collection of parallel molecular layers. Such a structure allows us to express the Helmholtz free energy of a molecular layer as the sum of the intrinsic Helmholtz free energy specific to that layer and the potential energy of interaction of that layer with all other layers and the solid surface. The intrinsic Helmholtz free energy of a molecular layer is a function (at given temperature) of its two-dimensional density and it can be readily obtained from bulk-phase properties, while the interlayer potential energy interaction is determined by using the 10-4 Lennard-Jones potential. The positions of all layers close to the graphite surface or in a slit pore are considered to correspond to the minimum of the potential energy of the system. This model has led to accurate predictions of nitrogen and argon adsorption on carbon black at their normal boiling points. In the case of adsorption in slit pores, local isotherms are determined from the minimization of the grand potential. The model provides a reasonable description of the 0-1 monolayer transition, phase transition and packing effect. The adsorption of nitrogen at 77.35 K and argon at 87.29 K on activated carbons is analyzed to illustrate the potential of this theory, and the derived pore-size distribution is compared favourably with that obtained by the Density Functional Theory (DFT). The model is less time-consuming than methods such as the DFT and Monte-Carlo simulation, and most importantly it can be readily extended to the adsorption of mixtures and capillary condensation phenomena.

Hydrotalcite intercalated siRNA:computational characterization of the interlayer environment

Using molecular dynamics (MD) simulations, we explore the structural and dynamical properties of siRNA within the intercalated environment of a Mg:Al 2:1 Layered Double Hydroxide (LDH) nanoparticle. An ab initio force field (Condensed-phase Optimized Molecular Potentials for Atomistic Simulation Studies: COMPASS) is used for the MD simulations of the hybrid organic-inorganic systems. The structure, arrangement, mobility, close contacts and hydrogen bonds associated with the intercalated RNA are examined and contrasted with those of the isolated RNA. Computed powder X-ray diffraction patterns are also compared with related LDH-DNA experiments. As a method of probing whether the intercalated environment approximates the crystalline or rather the aqueous state, we explore the stability of the principle parameters (e.g., the major groove width) that differentiate both A- and A'- crystalline forms of siRNA and contrast this with recent findings for the same siRNA simulated in water. We find the crystalline forms remain structurally distinct when intercalated, whereas this is not the case in water. Implications for the stability of hybrid LDH-RNA systems are discussed.