Comparison of two representations of a random cut of identical sphere packing

We compare two-dimensional froths obtained by radical tessellation of random planar cuts made through disordered assemblies of monosize spheres at different packing fractions C from C=0 to C=0.64 with two-dimensional stereological cuts obtained through the three-dimensional froths made with the same packing. We have built numerically the packings using different algorithms. The study of both topological and metric properties shows significant differences between the two representations. Received 26 May 1999 and Received in final form 13 November 1999     

Statistical verification of crystallization in hard sphere packings under densification

The performance of various structure characteristics in the task of indicating structural peculiarities in packings of hard spheres is investigated. Various characteristics based on Voronoi polyhedra, spherical harmonics, and Delaunay simplices are considered together with the pair correlation function and the mean number of r-close triples. They are applied to a set of hard sphere packings of density φ from 0.62 to 0.72. It turns out that all used structure characteristics are able to indicate changes of order from non-crystalline to crystalline packings. However, not all of them are sensitive enough to indicate different stages of structure transformation under densification. The characteristics based on Delaunay simplices turn out to be the most sensitive for this purpose. For the models considered three principal structure classes are found: packings of densities lower than the known critical value 0.64 showing a non-crystalline behavior; packings with considerable crystalline regions for φ up to 0.66–0.67; rather complete crystals although with numerous defects for φ above 0.67.     

A new handshaking of Tight-Binding and Molecular Dynamics in multi-scale simulation

A new handshake scheme is presented for tight-binding (TB) and molecular dynamics (MD) for multi-scale simulation of covalent crystals. In the present scheme, when calculating the forces on MD atoms in the handshake region, the TB atoms in close proximity to the MD atoms are treated as MD atoms. The scheme is thus seamless for calculation of MD atoms. When determining the electronic states of the TB subsystem, instead of the four basic atomic orbitals, hybrid orbitals are employed as bases in TB method and also as representing the action of MD atoms on TB atoms. The present handshaking methodology has several advantages. Firstly, it avoids determining the physical parameters required by introducing a new orbital model. Secondly, the “seam” almost decreases by one order of magnitude compared to that of Silogen model. Thirdly, the whole scheme is stable for dynamic simulation.     

Interatomic potential effects on dynamical heterogeneities in liquid SiO2

Dynamical heterogeneities (DH) in low density liquid SiO₂ have been investigated by molecular dynamics (MD) method. Simulations were done in the basic cube under periodic boundary conditions containing 3000 particles with the pair interatomic potentials, which have a weak electrostatic interaction and a Morse type short range interaction (PMSI). We have evaluated the non-Gaussian parameter for the self part of the van Hove correlation function and we found a clear evidence of the existence of DH in low density liquid SiO₂. Moreover, the atomic displacement distribution (ADD) in a model has been obtained and it deviates from a Gaussian form. The results have been compared with those obtained in another liquid SiO₂ system with the Born-Mayer interatomic potentials (BMP) in order to observe the interatomic potential effects on the DH in the system and indeed, the effects are strong. Calculations showed that particles of extremely low or fast mobility have a tendency to form a cluster and mean cluster size of most mobile and immobile particles in PMSI models increases with decreasing temperature. In contrast, no systematic changes have been obtained for the most mobile and immobile particles in BMP models. Calculations show that there is no relation between local particle environment and particle mobility in the system.     

Can a difference in molecular weights cause an eruption in a driven flow of self-organizing immiscible system?

Driven flow of a non-equilibrium non-conservative (NENC) system with a mixture of immiscible particles (A,B of molecular weight M_A, M_B) exhibits self-organizing patterns (segregation, phase-separation, etc.) in steady-state. The flow response (v) of mass flux density (j) to bias (H), in steady-state is found to be sensitive to molecular weight ratio (α = M_B/M_A). While the flux density (j) responds linearly to bias for both components (A, B) at α = 1, onset of eruptive response occurs at extreme bias (H ↦ 1) at α > 1 where v ↦∞ for heavier (B) and v ↦- ∞ for lighter (A) constituents. Difference in molecular weights (M_A, M_B) is not only critical to eruptive flow but also in controlling the flow response prior to this crossover.     

Field-induced realignment of a smectic nanodroplet in an external magnetic field: A numerical investigation

The field-induced realignment of a smectic-A phase is in principle a complicated process involving the director rotation via the interaction with the field and the layer rotation via the molecular interactions. Time-resolved X-ray scattering experiments have revealed major phenomena concerning the maintenance of the integrity of the smectic-A layer structure during the alignment process. In order to obtain a deeper insight into this process, we have carried out a dissipative particle dynamics study of the realignment kinetics of a nanodroplet of a smectic-A liquid crystal suspended in an isotropic fluid following a switch in the direction of an applied magnetic field. The strength of the mesogen-field interaction is small compared to the inter-molecular interactions. The reaction of the smectic configuration to the field switch was found to depend on the balance between the inter-molecular interactions stabilising the formation of the smectic layering and the interaction of the mesogens with the external field. It is found that the rotational behaviour of the smectic layers under the influence of an external magnetic field arises from a combination of stochastic translational displacements and rotational motions of the centres of mass of the mesogens in the nanodroplets. The simulations indicate that X-ray scattering and NMR experiments monitoring the orientational order are sensitive to different aspects of the realignment process.     

A dissipative particle dynamics study of the realignment of a nanodroplet of a nematic in a weak external magnetic field

We present a dissipative particle dynamics (DPD) approach for simulating the realignment of a nematic nanodroplet suspended in an isotropic fluid following a switch in the direction of an applied external magnetic field. The interaction of the mesogens with the external field is weak relative to the inter-molecular interactions. The simulations were used to investigate the way orientational equilibrium is re-established. The results reveal that the realignment process of the nanodroplet is consistent with its fluid structure. The reorientation of the nanodroplet as a whole is found to be caused by an internal structural rearrangement rather than a coherent rotation of the centres of mass of the mesogens about the centre of the nanodroplet. The switch in the field direction furthermore is found to induce a transient spatial variation in the orientational order of the long axes of the mesogens: the orientational order parameters decreases on moving from the core of the nanodroplet to the surface in contact with the isotropic environment. The results highlight differences between the time evolution of the orientation of the long molecular axes in the field and the rotations of the centres of mass of the mesogens about the centre of the nanodroplet.     

Living polymers in random media: a 2D Monte-Carlo investigation on a square lattice

The evolution of the mean chain length and mean end to end square radius of a two dimensional system of living polymers at constant monomer concentration is studied as a function of the obstacle density . The fact that the system adapts the mean chain length in order to reduce the entropic constraint does not lead to a different asymptotic dependence of on than what is observed for dead polymers. The change of the molecular weight distribution form in the presence of obstacles suggests that a Levy flight could appear in system of wormlike micelles in a porous medium. Received: 1 September 1997 / Received in final form: 22 January 1998 / Accepted: 30 January 1998     

Ab initio studies of the atomic and electronic structure of pure and hydrogenated a-Si

We propose a method to simulate a-Si and a-Si:H using an ab initio approach based on the Harris functional and thermally-amorphisized periodically-continued cells with at least 64 atoms. Hydrogen incorporation was achieved via diffusive addition. In preparing samples that may simulate the distributions of atoms in the amorphous materials, simulated annealing calculations were carried out from given starting conditions using short and long time steps. The different time-steps led to samples having distinctly different topological disorder. The radial distribution functions (partial and total) of the resulting samples were calculated and compared with measured distributions; the agreement is very good. These comparisons allowed some tentative conclusions to be made as regards the kind of disorder prevailing in the real samples. In addition, we studied the effect of the topological disorder on the electronic densities of states of the samples; the passivating effect of hydrogen can be observed. Received 17 May 2001     

Structure and magnetism on iron oxide clusters Fe nO m ( n = 1-5): Calculation from first principles

We have studied structural and magnetic properties in small iron oxide clusters, Fe_nO_m (n = 1-5), by means of the first-principles calculation based on the density functional theory. We have used not only the usual spin polarized scheme, but also the scheme for noncollinear magnetism to carry out efficient optimization in magnetic structure. The result of FeO_m (m = 1-4) is in good agreement with the previous work. We found the stable adduct clusters in FeO₅ and FeO₆. The bridge site of oxygen atom is more favorable in energy than any other site for the clusters of Fe_nO (n = 2-5). As increasing the number of oxygen atoms, the alignment of Fe magnetic moments changes from ferromagnetic configuration to antiferromagnetic one at Fe_nO_n (n = 2-4). Received 10 September 2002 Published online 3 July 2003     

Structure and magnetism in carbon nanotubes including magnetic wire

We have studied the electronic structure of the carbon nanotubes which include Fe atomic wire with using the density functional theory. As the stable geometries, we obtained the straight and zigzag wires, which have ferromagnetic and antiferromagnetic alignments, respectively. The antiferromagnets consists of the two ferromagnetic dimers which couple in antiparallel alignment. We presents the band dispersions and the density of states for the magnetic nanotubes. The electronic structure at the Fermi level consists of the Fe 3d and C 2pπ states, which shows a strong hybridization between them.     

Structural properties of amorphous TiO<Subscript>2</Subscript> nanoparticles

Structural properties of amorphous TiO₂ spherical nanoparticles have been studied in models with different sizes of 2 nm, 3 nm, 4 nm and 5 nm under non-periodic boundary conditions. We use the pairwise interatomic potentials proposed by Matsui and Akaogi. Models have been obtained by cooling from the melt via molecular dynamics (MD) simulation. Structural properties of an amorphous nanoparticle obtained at 350 K have been analyzed in detail through the partial radial distribution functions (PRDFs), coordination number distributions, bond-angle distributions and interatomic distances. Moreover, we show the radial density profile in a nanoparticle. Calculations show that size effects on structure of a model are significant and that if the size is larger than 3 nm, amorphous TiO₂ nanoparticles have a distorted octahedral network structure with the mean coordination number Z_Ti-O ≈6.0 and Z_O-Ti ≈3.0 like those observed in the bulk. Surface structure and surface energy of nanoparticles have been obtained and presented.     

A model of binary assemblies of spheres

We propose a theoretical model of random binary assemblies of spheres at any packing fraction. We use the notion of geometrical neighborhood between grains that is defined through two generalizations of the Vorono? tessellation: the radical (or Laguerre) tessellation and the navigation map. The model is tested on different numerical packings. We find a weak local segregation for high packing fraction. We also find that the higher the size ratio of the particles, the more important the segregation. Received 19 February 2001 and Received in final form 27 June 2001     

Desiccation of a clay film: Cracking versus peeling

We report a simulation study on competition between cracking and peeling, in a layer of clay on desiccation and how this is affected by the rate of drying, as well as the roughness of the substrate. The system is based on a simple 2-dimensional spring model. A vertical section through the layer with finite thickness is represented by a rectangular array of nodes connected by linear springs on a square lattice. The effect of reduction of the natural length of the springs, which mimics the drying is studied. Varying the strength of adhesion between sample and substrate and the rate of penetration of the drying front produces an interesting phase diagram, showing cross-over from peeling to cracking behavior. Changes in the number and width of cracks on varying the layer thickness is observed to reproduce experimental reports.     

Homogeneous nucleation and growth in iron-platinum vapour investigated by molecular dynamics simulation

Homogeneous nucleation and growth from binary metal vapour is investigated by molecular dynamics simulation. It is focused here mainly on the iron-platinum system with a mole fraction of 0.5. The simulations are started in the highly supersaturated vapour phase. Argon is added as carrier gas removing the heat of condensation from the forming clusters. The embedded atom method is employed for modelling of the force field of iron and platinum. The simulation runs are evaluated with respect to the nucleation rate, monomer temperature, monomer amount, and with respect to the size of the largest cluster in the system including possible pure metal clusters. It turns out that depending on the composition of the complete system pure platinum clusters with sizes up to 10 to 15 atoms are formed in addition to binary clusters. Due to the high temperature of these clusters iron atoms less likely condense at the beginning of the particle formation simulation. This leads to temporary difference in the temperatures of the platinum and the iron subsystems, which eventually approach each other when only binary clusters are present. In summary, the results obtained from the cluster statistics show that pure platinum nucleation and growth can take place to some extent within the binary system.     

Structural properties of amorphous Fe2O3 nanoparticles

We have investigated the microstructure of amorphous Fe₂O₃ nanoparticles by using molecular dynamics (MD) simulations. Non-periodic boundary conditions with Born-Mayer type pair potentials were used to simulate a spherical model of different diameters of 2, 3, 4 and 5 nm. Structural properties of an amorphous model obtained at 350 K have been analyzed in detail through the partial radial distribution functions (PRPFs), coordination number distributions, bond-angle distributions and interatomic distances. Calculations showed that structural characteristics of the model are in qualitative agreement with the experimental data. The observation of a large amount of structural defects as the particle size is decreased suggested that surface structure strongly depends on the size of nanoparticles. In addition, surface structure of amorphous Fe₂O₃ nanoparticles have been studied and compared with that observed in the core and in the bulk counterpart. Radial density profiles and stoichiometry in morphous Fe₂O₃ nanoparticles were also found and discussed.     

Topics in the theory of amorphous materials

In this Colloquium, I describe some current frontiers in the physics of semiconducting amorphous materials and glasses, including a short, but self-contained discussion of techniques for creating computer models, among them the quench from the melt method, the Activation-Relaxation Technique, the decorate and relax method, and the experimentally constrained molecular relaxation scheme. A representative study of an interesting and important glass (amorphous GeSe₃:Ag) is provided. This material is a fast-ion conductor and a serious candidate to replace current FLASH memory. Next, I discuss the effects of topological disorder on electronic states. By computing the decay of the density matrix in real space, and also computing well-localized Wannier functions, we close with a quantitative discussion of Kohn’s Principle of Nearsightedness in amorphous silicon.     

Ab-initio molecular dynamical study of a single transition metal atom on fullerene C<Subscript>60</Subscript>: the case of Ta

We report the first-principles Car-Parrinello molecular dynamics study of the behaviour of a single transition metal Ta atom on fullerene C₆₀, at different temperatures, and for both neutral and charged clusters. We seek to characterise the motion of the lone Ta metal atom on the C₆₀ surface, contrasting its behaviour both with that of three Ta atoms, as well as with a single alkali metal atom on the cage surface. Our earlier simulations on C₆₀Ta₃ had revealed that the Ta atoms on the surface of the fullerene are affected by a rather high mobility, and that the motion of these atoms is highly correlated due to Ta-atom-Ta-atom attraction. Earlier, experimental studies of a single metal atom (K, Rb) on the surface of a C₆₀ molecule had led to the inference that at room temperature the metal atom skates freely over the surface, the first direct evidence for which was presented by us in earlier first principles molecular dynamical simulations.     

Temperature and size dependence of the optical response of small Ag clusters

The absorption spectra of small Ag ⁺ n clusters are calculated at finite vibrational temperature by using a microscopic tight-binding RPA method. We consider free clusters with sizes between n = 3 and n =13 and take into account explicitly the degrees of freedom corresponding to the 4 d-electrons. We analyze the optical absorption as a function of the cluster size. We show that the contribution of the d-electrons has an important influence on the size dependence of the energy of the Mie plasmon. We also perform ensemble averages to obtain the absorption spectra for different vibrational temperatures. We obtain relatively good agreement with experiment for a temperature . The dynamics of the 4 d-electrons, which shows in small clusters an incipient delocalized character for n >7, yields an important contribution to the absorption spectrum already for n =13. We find that the strength of this contribution can be controlled by varying the vibrational temperature. Received: 4 January 1999 / Received in final form: 12 May 1999     

First-principles study of ground- and metastable-state properties of XO (X = Be,Mg, Ca,Sr, Ba,Zn and Cd)

The evolution of the local atomic order of an amorphous Ni₄₆Ti₅₄ alloy produced by mechanical alloying as a function of temperature was studied by synchrotron X-ray diffraction (XRD) and differential scanning calorimetry (DSC) techniques. XRD measurements at several temperatures (25 °C, 350 °C, 412 °C, 430 °C, 450 °C and 515 °C) were performed and analyzed using the reverse Monte Carlo (RMC) simulations method or the Rietveld refinement procedure. The experimental total structure factor for samples at 25 °C and 350 °C, which are amorphous in nature, were simulated by using the RMC method, and the local structures of the alloy at both temperatures were determined, indicating a decrease in its density as the temperature increases. At 412 °C, the XRD pattern shows a partially crystalline sample, indicating that the crystallization process is in progress. At 430 °C, 450 °C and 515 °C, the XRD measurements indicate the presence of two crystalline phases, NiTi and NiTi₂, whose structural parameters (lattice parameters, coherently diffracting domains (CDD) sizes, microstrains and relative amount of phases) were determined using the Rietveld refinement procedure. DSC measurements at different heating rates furnished the crystallization temperature, enthalpy and activation energy of the crystallization process, and these values are similar to those found in other amorphous alloys of the Ni-Ti system. They also showed the existence of a second exothermic process, which was related to diffusive processes in the crystalline phases, which could be associated with the changes in the CDD sizes happening from 450 °C to 515 °C.