Thermal diffusion in molecular dynamics simulations of thin film diamond deposition

Molecular dynamics simulations of carbon atom depositions are used to investigate energy diffusion from the impact zone. A modified Stillinger-Weber potential models the carbon interactions for both sp2 and sp3 bonding. Simulations were performed on 50 eV carbon atom depositions onto the (111) surface of a 3.8 x 3.4 x 1.0 nm diamond slab containing 2816 atoms in 11 layers of 256 atoms each. The bottom layer was thermostated to 300 K. At every 100th simulation time step (27 fs), the average local kinetic energy, and hence local temperature, is calculated. To do this the substrate is divided into a set of 15 concentric hemispherical zones, each of thickness one atomic diameter (0.14 nm) and centered on the impact point. A 50-eV incident atom heats the local impact zone above 10 000 K. After the initial large transient (200 fs) the impact zone has cooled below 3000 K, then near 1000 K by 1 ps. Thereafter the temperature profile decays approximately as described by diffusion theory, perturbed by atomic scale fluctuations. A continuum model of classical energy transfer is provided by the traditional thermal diffusion equation. The results show that continuum diffusion theory describes well energy diffusion in low energy atomic deposition processes, at distance and time scales larger than 1.5 nm and 1-2 ps, beyond which the energy decays essentially exponentially. (C) 1998 Published by Elsevier Science S.A. All rights reserved.

From point defects to plate tectonic faults

Understanding and explaining emergent constitutive laws in the multi-scale evolution from point defects, dislocations and two-dimensional defects to plate tectonic scales is an arduous challenge in condensed matter physics. The Earth appears to be the only planet known to have developed stable plate tectonics as a means to get rid of its heat. The emergence of plate tectonics out of mantle convection appears to rely intrinsically on the capacity to form extremely weak faults in the top 100 km of the planet. These faults have a memory of at least several hundred millions of years, yet they appear to rely on the effects of water on line defects. This important phenomenon was first discovered in laboratory and dubbed ``hydrolytic weakening''. At the large scale it explains cycles of co-located resurgence of plate generation and consumption (the Wilson cycle), but the exact physics underlying the process itself and the enormous spanning of scales still remains unclear. We present an attempt to use the multi-scale non-equilibrium thermodynamic energy evolution inside the deforming lithosphere to move phenomenological laws to laws derived from basic scaling quantities, develop self-consistent weakening laws at lithospheric scale and give a fully coupled deformation-weakening constitutive framework. At meso- to plate scale we encounter in a stepwise manner three basic domains governed by the diffusion/reaction time scales of grain growth, thermal diffusion and finally water mobility through point defects in the crystalline lattice. The latter process governs the planetary scale and controls the stability of its heat transfer mode.

Growth of crystals in optical tweezers

We report here on the use of optical tweezers in the growth and manipulation of protein and inorganic crystals. Sodium chloride and hen egg-white lysozyme crystals were grown in a batch process, and then seeds from the solution were introduced into the optical tweezers. The regular and controllable shape and the known optical birefringence in these structures allowed a detailed study of the orientation effects in the beam due to both polarization and gradient forces. Additionally, we determined that the laser tweezers could be used to suspend a crystal for three-dimensional growth under varying conditions. Studies included increasing the protein concentration, thermal cycling, and a diffusion-induced increase in precipitant concentration. Preliminary studies on the use of the tweezers to create a localized seed for growth from polyethylene oxide solutions are also reported.

Thermal analysis, nuclear magnetic resonance spectroscopy, and impedance spectroscopy of N,N-dimethyl-pyrrolidinium iodide: An ionic solid exhibiting rotator phases

N,N-dimethyl-pyrrolidinium iodide has been investigated using differential scanning calorimetry, nuclear magnetic resonance (NMR) spectroscopy, second moment calculations, and impedance spectroscopy. This pyrrolidinium salt exhibits two solid-solid phase transitions, one at 373 K having an entropy change, Delta S, of 38 J mol(-1) K-1 and one at 478 K having Delta S of 5.7 J mol(-1) K-1. The second moment calculations relate the lower temperature transition to a homogenization of the sample in terms of the mobility of the cations, while the high temperature phase transition is within the temperature region of isotropic tumbling of the cations. At higher temperatures a further decrease in the H-1 NMR linewidth is observed which is suggested to be due to diffusion of the cations. (C) 2005 American Institute of Physics.

Gas-surface interactions in the thermal and sub-thermal regime: A molecular dynamics study

Molecular dynamics simulations are used to study the interaction of low-energy Ar atoms with the Ni(001) surface, Angular scattering distributions, in and out of the plane of incidence, are investigated as a function of incident energy, angles of incidence, crystallographic orientation of the incident beam and surface temperature. The results show a clear transition to the structure scattering regime at around 2 eV. However, at lower energies, two sub-regimes are revealed by the simulations, Far energies up to 250 meV, scattering is mainly diffuse, and significant trapping on the surface is observed, At energies above this level, lobular patterns start to form and trapping decreases with the increase in energy, Generally, there is a weak temperature dependence, but variations in the angle of incidence and/or changes in the crystallographic direction, generate significant changes in the scattering patterns.

Non-equilibrium energy and momentum accommodation coefficients of Ar atoms scattered from Ni(001) in the thermal regime: A molecular dynamics study

Molecular dynamics simulations are used to study energy and momentum transfer of low-energy Ar atoms scattered from the Ni(001) surface. The investigation concentrates on the dependence of these processes on incident energy, angles of incidence and surface temperature. Energy transfer exhibits a strong dependence on the surface temperature, at incident energies below 500 meV, and incident angles close to specular incidence. Above 500 meV, the surface temperature dependence vanishes, and a limiting value in the amount of energy transferred to the surface is attained. Momentum exchange is investigated in terms of tangential and normal components. Both components exhibit a weak surface temperature dependence, but they have opposite behaviours at all incidence angles. In each component, momentum can be lost or gained following the interaction with the surface. (C) 1997 Elsevier Science B.V.

Fast, scalable master equation solution algorithms. IV. Lanczos iteration with diffusion approximation preconditioned iterative inversion

In this paper we propose a second linearly scalable method for solving large master equations arising in the context of gas-phase reactive systems. The new method is based on the well-known shift-invert Lanczos iteration using the GMRES iteration preconditioned using the diffusion approximation to the master equation to provide the inverse of the master equation matrix. In this way we avoid the cubic scaling of traditional master equation solution methods while maintaining the speed of a partial spectral decomposition. The method is tested using a master equation modeling the formation of propargyl from the reaction of singlet methylene with acetylene, proceeding through long-lived isomerizing intermediates. (C) 2003 American Institute of Physics.

Fast, scalable master equation solution algorithms. III. Direct time propagation accelerated by a diffusion approximation preconditioned iterative solver

In this paper we propose a novel fast and linearly scalable method for solving master equations arising in the context of gas-phase reactive systems, based on an existent stiff ordinary differential equation integrator. The required solution of a linear system involving the Jacobian matrix is achieved using the GMRES iteration preconditioned using the diffusion approximation to the master equation. In this way we avoid the cubic scaling of traditional master equation solution methods and maintain the low temperature robustness of numerical integration. The method is tested using a master equation modelling the formation of propargyl from the reaction of singlet methylene with acetylene, proceeding through long lived isomerizing intermediates. (C) 2003 American Institute of Physics.

Thermal Cycling Stability of Silica Membranes for Gas Separation

Hydrogen is being seen as an alternative energy carrier to conventional hydrocarbons to reduce greenhouse gas emissions. High efficiency separation technologies to remove hydrogen from the greenhouse gas, carbon dioxide, are therefore in growing demand. Traditional thermodynamic separation systems utilise distillation, absorption and adsorption, but are limited in efficiency at compact scales. Molecular sieve silica (MSS) membranes can perform this separation as they have high permselectivity of hydrogen to carbon dioxide, but their stability under thermal cycling is not well reported. In this work we exposed a standard MSS membrane and a carbonised template MSS (CTMSS) membrane to thermal cycling from 100 to 450°C. The standard MSS and carbonised template CTMSS membranes both showed permselectivity of helium to nitrogen dropping from around 10 to 6 in the first set of cycles, remaining stable until the last test. The permselectivity drop was due to small micropore collapse, which occurred via structure movement during cycling. Simulating single stage membrane separation with a 50:50 molar feed of H2:CO2, H2 exiting the permeate stream would start at 79% and stabilise at 67%. Higher selectivity membranes showed less of a purity drop, indicating the margin at which to design a stable membrane separation unit for CO2 capture.

Nanocomposite Nafion-Silica membranes for direct methanol fuel cells

Commercially available proton exchange membranes such as Nafion do not meet the requirements for high power density direct methanol fuel cells, partly due to their high methanol permeability. The aim of this work is to develop a new class of high-proton conductivity membranes, with thermal and mechanical stability similar to Nafion and reduced methanol permeability. Nanocomposite membranes were produced by the in-situ sol-gel synthesis of silicon dioxide particles in preformed Nafion membranes. Microstructural modification of Nafion membranes with silica nanoparticles was shown in this work to reduce methanol crossover from 7.48x10-6 cm2s^-1 for pure Nafion® to 2.86 x10-6 cm2s^-1 for nanocomposite nafion membranes (Methanol 50% (v/v) solution, 75 degrees C). Best results were achieved with a silica composition of 2.6% (w/w). We propose that silica inhibits the conduction of methanol through Nafion by blocking sites necessary for methanol diffusion through the polymer electrolyte membrane. Effects of surface chemistry, nanoparticle formation and interactions with Nafion matrix are further addressed.

Thermal effects in the evolution of initially layered mantle material

A simplified model for anisotropic mantle convection based on a novel class of rheologies, originally developed for folding instabilities in multilayered rock (MUHLHAUS et al., 2002), is extended ¨ through the introduction of a thermal anisotropy dependent on the local layering. To examine the effect of the thermal anisotropy on the evolution of mantle material, a parallel implementation of this model was undertaken using the Escript modelling toolkit and the Finley finite-element computational kernel (DAVIES et al., 2004). For the cases studied, there appears too little if any effect. For comparative purposes, the effects of anisotropic shear viscosity and the introduced thermal anisotropy are also presented. These results contribute to the characterization of viscous anisotropic mantle convection subject to variation in thermal conductivities and shear viscosities.

Growth of Bose-Einstein condensates from thermal vapor

We report on a quantitative study of the growth process of 87Rb Bose-Einstein condensates. By continuous evaporative cooling we directly control the thermal cloud from which the condensate grows. We compare the experimental data with the results of a theoretical model based on quantum kinetic theory. We find quantitative agreement with theory for the situation of strong cooling, whereas in the weak cooling regime a distinctly different behavior is found in the experiment.

A finite integral transform technique for solving the diffusion-reaction equation with Michaelis-Menten kinetics

An approximate analytical technique employing a finite integral transform is developed to solve the reaction diffusion problem with Michaelis-Menten kinetics in a solid of general shape. A simple infinite series solution for the substrate concentration is obtained as a function of the Thiele modulus, modified Sherwood number, and Michaelis constant. An iteration scheme is developed to bring the approximate solution closer to the exact solution. Comparison with the known exact solutions for slab geometry (quadrature) and numerically exact solutions for spherical geometry (orthogonal collocation) shows excellent agreement for all values of the Thiele modulus and Michaelis constant.

Transient response of diffusion cell containing a porous solid

Transient response of an adsorbing or non-adsorbing tracer injected as step or square pulse input in a diffusion cell with two flowing streams across the pellet is theoretically investigated in this paper. Exact solutions and the asymptotic solutions in the time domain and in three different limits are obtained by using an integral transform technique and a singular perturbation technique, respectively. Parametric dependence of the concentrations in the top and bottom chambers can be revealed by investigating the asymptotic solutions, which are far simpler than their exact counterpart. In the time domain investigation, it is found that the bottom-chamber concentration is very sensitive to the value of the macropore effective diffusivity. Therefore this concentration could be used to extract diffusivity by fitting in the time domain. The bottom-chamber concentration is also sensitive to flow rate, pellet length chamber volume and the type of input (step and square input).