Molecular dynamics simulation of observed c(4x4) and c(4x3) C60 alignments on the Si(100) reconstructed surface.

We have implemented a large-scale classical molecular dynamics simulation at constant temperature to provide a theoretical insight into the results of a recently performed experiment on the monolayer and multi-layer formations of molecular films on the Si(100) reconstructed dimerized surface. Our simulation has successfully reproduced all of the morphologies observed on the monolayer film by this experiment. We have obtained the formation of both c(4 4) and c(4 3) structures of the molecules and have also obtained phase transitions of the former into the latter.

Molecular dynamics simulation of adsorption of Ag particles on a graphite substrate

We have performed for the first time a molecular dynamics simulation of the adsorption of gas-phase Ag particles on a graphite substrate to provide an insight into the results of a comprehensive STM-based experiment on this system. Both pair-wise and many-body interatomic potentials have been employed, and a Morse-type Ag–C potential was specifically constructed to describe the interactions at the interface. Our simulation has successfully reproduced a significant portion of the experimental findings. We have also observed the intercalation of silver in graphite.

Molecular dynamics simulation of fractures using an N-body potential

Molecular dynamics has been employed to model the fracture of a twodimensional triangular atomic lattice. The N-body Sutton-Chen potential developed for fcc metals and its extended version (Rafii-Tabar and Sutton) for fcc random binary alloys were used for the interatomic interactions. It is shown that at low temperatures cleavage fractures can occur in both an elemental metal and an alloy. At elevated temperatures the nucleation of dislocations is shown to cause a brittle-to-ductile transition. For the brittle crack propagation in the elemental metal, crack propagation speeds have been computed for different stress rates, and a crack instability found to exist as the speed reaches a critical value of about 32% of the Rayleigh wave speed. For the random alloy, we find that the dislocation movement can be affected by the distorted lattice.

Molecular dynamics simulation of fractures using an N-body potential

Molecular dynamics has been employed to model the fracture of a two dimensional triangular atomic lattice. The N-body Sutton-Chen potential developed for fcc metals and its extended version (Rafii-Tabar and Sutton) for fcc random binary alloys were used for the interatomic interactions. It is shown that at low temperatures cleavage fractures can occur in both an elemental metal and an alloy. At elevated temperatures the nucleation of dislocations is shown to cause a brittle-to-ductile transition. For the brittle crack propagation in the elemental metal, crack propagation speeds have been computed for different stress rates, and a crack instability found to exist as the speed reaches a critical value of about 32% of the Rayleigh wave speed. For the random alloy, we find that the dislocation movement can be affected by the distorted lattice.

Numerical modelling of adsorption of metallic particles on graphite substrate via molecular dynamics simulation

A computer-based numerical modelling of the adsorption process of gas phase metallic particles on the surface of a graphite substrate has been performed via the application of molecular dynamics simulation method. The simulation relates to an extensive STM-based experiment performed in this field, and reproduces part of the experimental results. Both two-body and many-body inter-atomic potentials have been employed. A Morse-type potential describing the metal-carbon interactions at the interface was specifically formulated for this modelling. Intercalation of silver in graphite has been observed as well as the correct alignments of monomers, dimers and two-dimensional islands on the surface. PACS numbers: 02.60.Cb, 07.05.Tp, 68.55.-a, 81.05.Tp

Multiscale numerical modelling of crack propagation in two-dimensional metal plate

A novel multiscale model of brittle crack propagation in an Ag plate with macroscopic dimensions has been developed. The model represents crack propagation as stochastic drift-diffusion motion of the crack tip atom through the material, and couples the dynamics across three different length scales. It integrates the nanomechanics of bond rupture at the crack tip with the displacement and stress field equations of continuum based fracture theories. The finite element method is employed to obtain the continuum based displacement and stress fields over the macroscopic plate, and these are then used to drive the crack tip forward at the atomic level using the molecular dynamics simulation method based on many-body interatomic potentials. The linkage from the nanoscopic scale back to the macroscopic scale is established via the Ito stochastic calculus, the stochastic differential equation of which advances the tip to a new position on the macroscopic scale using the crack velocity and diffusion constant obtained on the nanoscale. Well known crack characteristics, such as the roughening transitions of the crack surfaces, crack velocity oscillations, as well as the macroscopic crack trajectories, are obtained.

Modelling the adsorption and imaging of C60 molecules on a graphite substrate

The adsorption of a C60 monolayer on a graphite substrate was modelled via molecular dynamics simulation covering a significant period of 160 picoseconds. The final configuration of C60s agrees closely with that observed in a scanning tunnelling microscopy (STM) experiment. Clusters of adsorbed molecules were then selected and their STM-like images were computed via the Keldysh Green function method.

Collaborating components in mesh-based electronic packaging

From the model geometry creation to the model analysis, the stages in between such as mesh generation are the most manpower intensive phase in a mesh-based computational mechanics simulation process. On the other hand the model analysis is the most computing intensive phase. Advanced computational hardware and software have significantly reduced the computing time - and more importantly the trend is downward. With the kind of models envisaged coming, which are larger, more complex in geometry and modelling, and multiphysics, there is no clear trend that the manpower intensive phase is to decrease significantly in time - in the present way of operation it is more likely to increase with model complexity. In this paper we address this dilemma in collaborating components for models in electronic packaging application.

Novel inclusion complex of ibuprofen tromethamine with cyclodextrins: Physico-chemical characterization

Guest-host interactions of ibuprofen tromethamine salt (Ibu.T) with native and modified cyclodextrins (CyDs) have been investigated using several techniques, namely phase solubility diagrams (PSDs), proton nuclear magnetic resonance (H-1 NMR), differential scanning calorimetry (DSC), Fourier transform infrared spectroscopy (FT-IR), X-ray powder diffractometry (XRPD). scanning-electron microscopy (SEM) and molecular mechanics (MM). From the analysis of PSD data (A(L)-type) it is concluded that the anionic tromethamine salt of ibuprofen (pK(a) = 4.55) forms 1: 1 soluble complexes with all CyDs investigated in buffered water at pH 7.0, while the neutral form of Ibu forms an insoluble complex with beta-CyD (B-S-type) in buffered water at pH 2.0. Ibu.T has a lower tendency to complex with beta-CyD (K-11 = 58 M-1 at pH 7.0) compared with the neutral Ibu (K-11 = 4200 M (1)) in water. Complex formation of Ibu.T with beta-CyD (Delta G degrees = -20.4 kJ/mol) is enthalpy driven (Delta H degrees = -22.9 kJ/mol) and is accompanied by a small unfavorable entropy (Delta S degrees = -8.4 J/mol K) change. H-1 NMR studies and MM computations revealed that, on complexation, the hydrophobic central benzene ring of lbu.T and part of the isobutyl group reside within the beta-CyD cavity leaving the peripheral groups (carboxylate, tromethamine and methyl groups) located near the hydroxyl group networks at either rim of beta-CyD. PSD, H-1 NMR, DSC, FT-IR, XRPD, SEM and MM studies confirmed the formation of Ibu.T/beta-CyD inclusion complex in solution and the solid state. (C) 2009 Elsevier B.V. All rights reserved.

Surface physics at the nano-scale via scanning probe microscopy and molecular dynamics simulations

Probe-based scanning microscopes, such as the STM and the AFM, are used to obtain the topographical and electronic structure maps of material surfaces, and to modify their morphologies on nanoscopic scales. They have generated new areas of research in condensed matter physics and materials science. We will review some examples from the fields of experimental nano-mechanics, nano-electronics and nano-magnetism. These now form the basis of the emerging field of Nano-technology. A parallel development has been brought about in the field of Computational Nano-science, using quantum-mechanical techniques and computer-based numerical modelling, such as the Molecular Dynamics (MD) simulation method. We will report on the simulation of nucleation and growth of nano-phase films on supporting substrates. Furthermore, a theoretical modelling of the formation of STM images of metallic clusters on metallic substrates will also be discussed within the non-equilibrium Keldysh Green function method to study the effects of coherent tunnelling through different atomic orbitals in a tip-sample geometry.

A multi-scale numerical modelling of crack propagation in a 2D metallic plate

A new multi-scale model of brittle fracture growth in an Ag plate with macroscopic dimensions is proposed in which the crack propagation is identified with the stochastic drift-diffusion motion of the crack-tip atom through the material. The model couples molecular dynamics simulations, based on many-body interatomic potentials, with the continuum-based theories of fracture mechanics. The Ito stochastic differential equation is used to advance the tip position on a macroscopic scale before each nano-scale simulation is performed. Well-known crack characteristics, such as the roughening transitions of the crack surfaces, as well as the macroscopic crack trajectories are obtained.

Parallel performance in multi-physics simulation

A comprehensive simulation of solidification/melting processes requires the simultaneous representation of free surface fluid flow, heat transfer, phase change, non-linear solid mechanics and, possibly, electromagnetics together with their interactions in what is now referred to as "multi-physics" simulation. A 3D computational procedure and software tool, PHYSICA, embedding the above multi-physics models using finite volume methods on unstructured meshes (FV-UM) has been developed. Multi-physics simulations are extremely compute intensive and a strategy to parallelise such codes has, therefore, been developed. This strategy has been applied to PHYSICA and evaluated on a range of challenging multi-physics problems drawn from actual industrial cases.

Numerical simulation of microwave thawing using a coupled solver approach

Thawing of a frozen food product in a domestic microwave oven is numerically simulated using a coupled solver approach. The approach consists of a dedicated electromagnetic FDTD solver and a closely coupled UFVM multi-physics package. Two overlapping numerical meshes are defined; the food material and container were meshed for heat transfer and phase change solution, whilst the microwave oven cavity and waveguide were meshed for the microwave irradiation. The two solution domains were linked using a cross-mapping routine. This approach allowed the rotation of the food load to be captured. Power densities obtained on the structured FDTD mesh were interpolated onto the UFVM mesh for each timestep/turntable position. The UFVM solver utilised the power density data to advance the temperature and phase distribution solution. The temperature-dependant dielectric and thermo-physical properties of the food load were updated prior to revising the electromagnetic solution. Changes in thermal/electric properties associated with the phase transition were fully accounted for as well as heat losses from product to cavity. Two scenarios were investigated: a centric and eccentric placement on the turntable. Developing temperature fields predicted by the numerical solution are validated against experimentally obtained data. Presented results indicate the feasibility of fully coupled simulations of the microwave heating of a frozen product. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)

Multiphysics simulation of microwave processing using a multidomain coupled solver approach

Microwave processing of materials is numerically simulated using a coupled solver approach. Microwave heating is a complex coupled process due to the variation in dielectric properties during heating. The effects of heating an object in a electromagnetic field directly influence the manner in which it interacts with the field. Simplifying assumptions and empirical solutions do not capture the fundamental physics involved and, in general, do not provide usefully accurate solutions in a number of practical problems. In order to capture the underlying processes involved in microwave heating, the problem must be looked at in a holistic manner rather than a number of discrete processes. This contribution outlines a coupled-solver multiphysics analysis approach to the solution of practical microwave heating problems.