An atomistic simulation of solid state sintering using Monte Carlo methods

This paper presents results on the simulation of the solid state sintering of copper wires using Monte Carlo techniques based on elements of lattice theory and cellular automata. The initial structure is superimposed onto a triangular, two-dimensional lattice, where each lattice site corresponds to either an atom or vacancy. The number of vacancies varies with the simulation temperature, while a cluster of vacancies is a pore. To simulate sintering, lattice sites are picked at random and reoriented in terms of an atomistic model governing mass transport. The probability that an atom has sufficient energy to jump to a vacant lattice site is related to the jump frequency, and hence the diffusion coefficient, while the probability that an atomic jump will be accepted is related to the change in energy of the system as a result of the jump, as determined by the change in the number of nearest neighbours. The jump frequency is also used to relate model time, measured in Monte Carlo Steps, to the actual sintering time. The model incorporates bulk, grain boundary and surface diffusion terms and includes vacancy annihilation on the grain boundaries. The predictions of the model were found to be consistent with experimental data, both in terms of the microstructural evolution and in terms of the sintering time. (C) 2002 Elsevier Science B.V. All rights reserved.

Fluid-mobile trace element constraints on the role of slab melting and implications for Archaean crustal growth models

We use published and new trace element data to identify element ratios which discriminate between arc magmas from the supra-subduction zone mantle wedge and those formed by direct melting of subducted crust (i.e. adakites). The clearest distinction is obtained with those element ratios which are strongly fractionated during refertilisation of the depleted mantle wedge, ultimately reflecting slab dehydration. Hence, adakites have significantly lower Pb/Nd and B/Be but higher Nb/Ta than typical arc magmas and continental crust as a whole. Although Li and Be are also overenriched in continental crust, behaviour of Li/Yb and Be/Nd is more complex and these ratios do not provide unique signatures of slab melting. Archaean tonalite-trondhjemite-granodiorites (TTGs) strongly resemble ordinary mantle wedge-derived arc magmas in terms of fluid-mobile trace element content, implying that they-did not form by slab melting but that they originated from mantle which was hydrated and enriched in elements lost from slabs during prograde dehydration. We suggest that Archaean TTGs formed by extensive fractional crystallisation from a mafic precursor. It is widely claimed that the time between the creation and subduction of oceanic lithosphere was significantly shorter in the Archaean (i.e. 20 Ma) than it is today. This difference was seen as an attractive explanation for the presumed preponderance of adakitic magmas during the first half of Earth's history. However, when we consider the effects of a higher potential mantle temperature on the thickness of oceanic crust, it follows that the mean age of oceanic lithosphere has remained virtually constant. Formation of adakites has therefore always depended on local plate geometry and not on potential mantle temperature.

The martensitic transformation in ceramics - its role in transformation toughening

This paper reviews the current knowledge and understanding of martensitic transformations in ceramics - the tetragonal to monoclinic transformation in zirconia in particular. This martensitic transformation is the key to transformation toughening in zirconia ceramics. A very considerable body of experimental data on the characteristics of this transformation is now available. In addition, theoretical predictions can be made using the phenomenological theory of martensitic transformations. As the paper will illustrate, the phenomenological theory is capable of explaining all the reported microstructural and crystallographic features of the transformation in zirconia and in some other ceramic systems. Hence the theory, supported by experiment, can be used with considerable confidence to provide the quantitative data that is essential for developing a credible, comprehensive understanding of the transformation toughening process. A critical feature in transformation toughening is the shape strain that accompanies the transformation. This shape strain, or nucleation strain, determines whether or not the stress-induced martensitic transformation can occur at the tip of a potentially dangerous crack. If transformation does take place, then it is the net transformation strain left behind in the transformed region that provides toughening by hindering crack growth. The fracture mechanics based models for transformation toughening, therefore, depend on having a full understanding of the characteristics of the martensitic transformation and, in particular, on being able to specify both these strains. A review of the development of the models for transformation toughening shows that their refinement and improvement over the last couple of decades has been largely a result of the inclusion of more of the characteristics of the stress-induced martensitic transformation. The paper advances an improved model for the stress-induced martensitic transformation and the strains resulting from the transformation. This model, which separates the nucleation strain from the subsequent net transformation strain, is shown to be superior to any of the constitutive models currently available. (C) 2002 Elsevier Science Ltd. All rights reserved.

Fibroblast growth factor receptor 4 (FGFR4) expression in newborn murine calvaria and primary osteoblast cultures

Fibroblast growth factor receptor (FGFR) signalling is important in the initiation and regulation of osteogenesis. Although mutations in FGFR1, 2 and 3 genes are known to cause skeletal deformities, the expression of FGFR4 in bony tissue remains unclear. We have investigated the expression pattern of FGFR4 in the neonatal mouse calvaria and compared it to the expression pattern in cultures of primary osteoblasts. Immunohistochemistry demonstrated that FGFR4 was highly expressed in rudimentary membranous bone and strictly localised to the cellular components (osteoblasts) between the periosteal and endosteal layers. Cells in close proximity to the newly formed osteoid (preosteoblasts) also expressed FGFR4 on both the endosteal and periosteal surfaces. Immunocytochemical analysis of primary osteoblast cultures taken from the same cranial region also revealed high levels of FGFR4 expression, suggesting a similar pattern of cellular expression in vivo and in vitro. RT-PCR and Western blotting for FGFR4 confirmed its presence in primary osteoblast cultures. These results suggest that FGFR4 may be an important regulator of osteogenesis with involvement in preosteoblast proliferation and differentiation as well as osteoblast functioning during intramembranous ossification. The consistent expression of FGFR4 in vivo and in vitro supports the use of primary osteoblast cultures for elucidating the role of FGFR4 during osteogenesis.

Molecular orbital theory calculations of the H2O-carbon reaction

Carbon gasification with steam to produce H-2 and CO is an important reaction widely used in industry for hydrogen generation. Although the literature is vast, the. mechanism for the formation of H-2 is still unclear. In particular, little has, been done to investigate the potential of molecular orbital theory to distinguish different mechanism possibilities. In this work, we used molecular orbital theory to demonstrate a favorable energetic pathway where H2O is first physically adsorbed on the virgin graphite surface with negligible change in molecular structure. Chemisorption occurs via O approaching the carbon edge site with one H atom stretching away from the O in the transition state. This is followed by a local minimum. state in which the stretching H is further disconnected from the O atoms and the remaining OH group is still on the carbon edge site. The disconnected H then pivot around the OH group to bond with the H of the OH group and forms H-2. The O atom remaining on the carbon edge site is subsequently desorbed as CO. The reverse occurs when H-2 reacts with the surface oxygen to produce H2O.

Close packed transitions in slit-shaped pores: Density functional theory study of methane adsorption capacity in carbon

An important feature of improving lattice gas models and classical isotherms is the incorporation of a pore size dependent capacity, which has hitherto been overlooked. In this paper, we develop a model for predicting the temperature dependent variation in capacity with pore size. The model is based on the analysis of a lattice gas model using a density functional theory approach at the close packed limit. Fluid-fluid and solid-fluid interactions are modeled by the Lennard-Jones 12-6 potential and Steele's 10-4-3, potential respectively. The capacity of methane in a slit-shaped carbon pore is calculated from the characteristic parameters of the unit cell, which are extracted by minimizing the grand potential of the unit cell. The capacities predicted by the proposed model are in good agreement with those obtained from grand canonical Monte Carlo simulation, for pores that can accommodate up to three adsorbed layers. Single particle and pair distributions exhibit characteristic features that correspond to the sequence of buckling and rhombic transitions that occur as the slit pore width is increased. The model provides a useful tool to model continuous variation in the microstructure of an adsorbed phase, namely buckling and rhombic transitions, with increasing pore width. (C) 2002 American Institute of Physics.

Kinetics of adsorption on activated carbon: application of heterogeneous vacancy solution theory

The kinetics of single component adsorption on activated carbon is investigated here using a heterogeneous vacancy solution theory (VST) of adsorption. The adsorption isotherm is developed to account for the adsorbate non-ideality due to the size difference between the adsorbate molecule and the vacant site, while incorporating adsorbent heterogeneity through a pore-width-related potential energy. The transport process in the bidisperse carbon considers coupled mass transfer in both macropore and micropore phases simultaneously. Adsorbate diffusion in the micropore network is modeled through effective medium theory, thus considering pore network connectivity in the adsorbent, with the activation energy for adsorbate diffusion related to the adsorption energy, represented by the Steele 10-4-3 potential for carbons. Experimental data of five hydrocarbons, CO2 and SO2 on Ajax carbon at multiple temperatures, as well as three hydrocarbons on Norit carbon at three temperatures are first fitted by the heterogeneous VST model to obtain the isotherm parameters, followed by application of the kinetic model to uptake data on carbon particles of different sizes and geometry at various temperatures. For the hydrocarbons studied, the model can successfully correlate the experimental data for both adsorption equilibrium and kinetics. However, there is some deviation in the fit of the desorption kinetics for polar compounds such as CO2 and SO2, due to the inadequacy of the L-J potential model in this case. The significance of viscous transport in the micropores is also considered here and found to be negligible, consistent with recent molecular simulation studies. (C) 2002 Elsevier Science Ltd. All rights reserved.

Vacancy solution theory for binary adsorption equilibria in heterogeneous carbon

A heterogeneous modified vacancy solution model of adsorption developed is evaluated. The new model considers the adsorption process through a mass-action law and is thermodynamically consistent, while maintaining the simplicity in calculation of multicomponent adsorption equilibria, as in the original vacancy solution theory. It incorporates the adsorbent heterogeneity through a pore-width-related potential energy, represented by Steele's 10-4-3 potential expression. The experimental data of various hydrocarbons, CO2 and SO2 on four different activated carbons - Ajax, Norit, Nuxit, and BPL - at multiple temperatures over a wide range of pressures were studied by the heterogeneous modified VST model to obtain the isotherm parameters and micropore-size distribution of carbons. The model successfully correlates the single-component adsorption equilibrium data for all compounds studied on various carbons. The fitting results for the vacancy occupancy parameter are consistent with the pressure change on different carbons, and the effect of pore heterogeneity is important in adsorption at elevated pressure. It predicts binary adsorption equilibria better than the IAST scheme, reflecting the significance of molecular size nonideality.

Density functional theory analysis of the influence of pore wall heterogeneity on adsorption in carbons

Density functional theory for adsorption in carbons is adapted here to incorporate a random distribution of pore wall thickness in the solid, and it is shown that the mean pore wall thickness is intimately related to the pore size distribution characteristics. For typical carbons the pore walls are estimated to comprise only about two graphene layers, and application of the modified density functional theory approach shows that the commonly used assumption of infinitely thick walls can severely affect the results for adsorption in small pores under both supercritical and subcritical conditions. Under supercritical conditions the Henry's law coefficient is overpredicted by as much as a factor of 2, while under subcritical conditions pore wall heterogeneity appears to modify transitions in small pores into a sequence of smaller ones corresponding to pores with different wall thicknesses. The results suggest the need to improve current pore size distrubution analysis methods to allow for pore wall heterogeneity. The density functional theory is further extended here to allow for interpore adsorbate interactions, and it appears that these interaction are negligible for small molecules such as nitrogen but significant for more strongly interacting heavier molecules such as butane, for which the traditional independent pore model may not be adequate.

Modeling of phase equilibrium and vapor adsorption on carbon black based on a combination of a lattice theory and equation of state

A thermodynamic approach is developed in this paper to describe the behavior of a subcritical fluid in the neighborhood of vapor-liquid interface and close to a graphite surface. The fluid is modeled as a system of parallel molecular layers. The Helmholtz free energy of the fluid is expressed as the sum of the intrinsic Helmholtz free energies of separate layers and the potential energy of their mutual interactions calculated by the 10-4 potential. This Helmholtz free energy is described by an equation of state (such as the Bender or Peng-Robinson equation), which allows us a convenient means to obtain the intrinsic Helmholtz free energy of each molecular layer as a function of its two-dimensional density. All molecular layers of the bulk fluid are in mechanical equilibrium corresponding to the minimum of the total potential energy. In the case of adsorption the external potential exerted by the graphite layers is added to the free energy. The state of the interface zone between the liquid and the vapor phases or the state of the adsorbed phase is determined by the minimum of the grand potential. In the case of phase equilibrium the approach leads to the distribution of density and pressure over the transition zone. The interrelation between the collision diameter and the potential well depth was determined by the surface tension. It was shown that the distance between neighboring molecular layers substantially changes in the vapor-liquid transition zone and in the adsorbed phase with loading. The approach is considered in this paper for the case of adsorption of argon and nitrogen on carbon black. In both cases an excellent agreement with the experimental data was achieved without additional assumptions and fitting parameters, except for the fluid-solid potential well depth. The approach has far-reaching consequences and can be readily extended to the model of adsorption in slit pores of carbonaceous materials and to the analysis of multicomponent adsorption systems. (C) 2002 Elsevier Science (USA).

Finite, three-dimensional adsorbed phase growth in activated carbon mesopores

The role of PACs (primary adsorption centers) in the mesopore (i.e., transport) region of activated carbons during adsorption of polar species, such as water, is unclear. A classical model of three-dimensional adsorption on finite PACs is presented. The model is a preliminary, theoretical investigation into adsorption on mesopore PACs and is intended to give some insight into the energetic and physical processes at work. Work processes are developed to obtain isotherms and three-dimensional sorbate growth on PACs of varying size and energetic characteristics. The work processes allow two forms of adsorbed phase growth: densification at constant boundary and boundary growth at constant density. Relatively strong sorbate-sorbent interactions and strong surface tension favor adsorbed phase densification over boundary growth. Conversely, relatively weak sorbate-sorbent interactions and weak surface tension favor boundary growth over densification. If sorbate-sorbate interactions are strong compared to sorbate-sorbent interactions, condensation with hysteresis occurs. This can also give rise to delayed boundary growth, where all initial adsorption occurs in the monolayer only. The results indicate that adsorbed phase growth on PACs may be quite complex.

Control of ovulation with a GnRH agonist after superstimulation of follicular growth in buffalo: fertilization and embryo recovery

The potential to use a GnRH agonist bioimplant and injection of exogenous LH to control the time of ovulation in a multiple ovulation and embryo transfer (MOET) protocol was examined in buffalo. Mixed-parity buffalo (Bubalus bubalis; 4-15-year-old; 529 13 kg LW) were randomly assigned to one of five groups (n = 6): Group 1, conventional MOET protocol; Group 2, conventional MOET with 12 It delay in injection of PGF(2alpha); Group 3, implanted with GnRH agonist to block the pre-ovulatory surge release of LH; Group 4, implanted with GnRH agonist and injected with exogenous LH (Lutropin(R), 25 mg) 24 h after 4 days of superstimulation with FSH; Group 5, implanted with GnRH agonist and injected with LH 36 h after superstimulation with FSH. Ovarian follicular growth in all buffaloes was stimulated by treatment with FSH (Folltropin-V(R), 200 mg) administered over 4 days, and was monitored by ovarian ultrasonography. At the time of estrus, the number of follicles greater than or equal to8 mm. was greater (P < 0.05) for buffaloes in Group 2 (12.8) than for buffaloes in Groups 1 (8.5), 3 (7.3), 4 (6.1) and 5 (6.8), which did not differ. All buffaloes were mated by AI after spontaneous (Groups 1-3) or induced (Groups 4 and 5) ovulation. The respective number of buffalo that ovulated, number of corpora lutea, ovulation rate (%), and embryos + oocytes recovered were: Group 1 (2, 1.8 +/- 1.6, 18.0 +/- 13.6, 0.2 +/- 0.2); Group 2 (4, 6.1 +/- 2.9, 40.5 +/- 17.5, 3.7 +/- 2.1); Group 3 (0, 0, 0, 0); Group 4 (6, 4.3 +/- 1.2, 69.3 +/- 14.2, 2.0 +/- 0.9); and Group 5 (1, 2.5 +/- 2.5, 15.5 +/- 15.5, 2.1 +/- 2.1). All buffaloes in Group 4 ovulated after injection of LH and had a relatively high ovulation rate (69%) and embryo recovery (46%). It has been shown that the GnRH agonist-LH protocol can be used to improve the efficiency of MOET in buffalo. (C) 2002 Elsevier Science Inc. All rights reserved.