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.

Combined chemical separation of Lu, Hf, Sm, Nd, and REEs from a single rock digest: Precise and accurate isotope Determinations of Lu-Hf and Sm-Nd using multicollector-ICPMS

A combined procedure for separating Lu, Hf, Sm, Nd, and rare earth elements (REEs) from a single sample digest is presented. The procedure consists of the following five steps: (1) sample dissolution via sodium peroxide sintering; (2) separation of the high field strength elements from the REEs and other matrix elements by a HF-free anion-exchange column procedure; (3) purification of Hf on a cation-exchange resin; (4) separation of REEs from other matrix elements by cation exchange; (5) Lu, Sm, and Nd separation from the other REEs by reversed-phase ion chromatography. Analytical reproducibilities of Sm-Nd and Lu-Hf isotope systematics are demonstrated for standard solutions and international rock reference materials. Results show overall good reproducibilities for Sm-Nd systematics independent of the rock type analyzed. For the Lu-Hf systematics, the reproducibility of the parent/daughter ratio is much better for JB-1 (basalt) than for two analyzed felsic crustal rocks (DR-N and an Archaean granitoid). It is demonstrated that this poorer reproducibility of the Lu/Hf ratio is truly caused by sample heterogeneity; thus, results are geologically reasonable.

Axial mixing in rotating drums using magnetic resonance imaging using bran as a model for solid state fermentations

Magnetic resonance imaging (MRI) is an easily automated, reliable technique to investigate axial mixing within rotating drums. Moist bran can be clearly differentiated from dry bran using MRI allowing a non-segregating tracer for axial mixing. For a 20-cm diameter drum, the axial dispersion coefficient in the particle bed was 0.51 cm s(-2). Axial dispersion is scale-dependent.

The strength of concentrated Mg-Zn solid solutions

Solid solution effects on the hardness and flow stress have been studied for zinc contents between 0.2 and 2.4 at% (0.5 and 6.9 wt%) in Mg. The alloys were grain refined with 0.6 wt% zirconium to ensure a similar grain size at all compositions. The hardness increases with the zinc content as Hv(10) (kg mm(-2)) = 9 Zn (at%) + 33. At low solute concentrations the (0.2%) proof strength does not change significantly with concentration. At concentrations above 0.7 at%, within the supersaturated solid solution region, the rate of solid solution hardening is high, following a c(2) rule, where c is the atom fraction of Zn. It is suggested that short-range order may account for most of the observed strengthening in concentrated Mg-Zn alloys.

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).

Modeling of gas adsorption equilibrium over a wide range of pressure: A thermodynamic approach based on equation of state

A thermodynamic approach based on the Bender equation of state is suggested for the analysis of supercritical gas adsorption on activated carbons at high pressure. The approach accounts for the equality of the chemical potential in the adsorbed phase and that in the corresponding bulk phase and the distribution of elements of the adsorption volume (EAV) over the potential energy for gas-solid interaction. This scheme is extended to subcritical fluid adsorption and takes into account the phase transition in EAV The method is adapted to gravimetric measurements of mass excess adsorption and has been applied to the adsorption of argon, nitrogen, methane, ethane, carbon dioxide, and helium on activated carbon Norit R I in the temperature range from 25 to 70 C. The distribution function of adsorption volume elements over potentials exhibits overlapping peaks and is consistently reproduced for different gases. It was found that the distribution function changes weakly with temperature, which was confirmed by its comparison with the distribution function obtained by the same method using nitrogen adsorption isotherm at 77 K. It was shown that parameters such as pore volume and skeleton density can be determined directly from adsorption measurements, while the conventional approach of helium expansion at room temperature can lead to erroneous results due to the adsorption of helium in small pores of activated carbon. The approach is a convenient tool for analysis and correlation of excess adsorption isotherms over a wide range of pressure and temperature. This approach can be readily extended to the analysis of multicomponent adsorption systems. (C) 2002 Elsevier Science (USA).

Stress correlation function evolution in lattice solid elasto-dynamic models of shear and fracture zones and earthquake prediction

It has been argued that power-law time-to-failure fits for cumulative Benioff strain and an evolution in size-frequency statistics in the lead-up to large earthquakes are evidence that the crust behaves as a Critical Point (CP) system. If so, intermediate-term earthquake prediction is possible. However, this hypothesis has not been proven. If the crust does behave as a CP system, stress correlation lengths should grow in the lead-up to large events through the action of small to moderate ruptures and drop sharply once a large event occurs. However this evolution in stress correlation lengths cannot be observed directly. Here we show, using the lattice solid model to describe discontinuous elasto-dynamic systems subjected to shear and compression, that it is for possible correlation lengths to exhibit CP-type evolution. In the case of a granular system subjected to shear, this evolution occurs in the lead-up to the largest event and is accompanied by an increasing rate of moderate-sized events and power-law acceleration of Benioff strain release. In the case of an intact sample system subjected to compression, the evolution occurs only after a mature fracture system has developed. The results support the existence of a physical mechanism for intermediate-term earthquake forecasting and suggest this mechanism is fault-system dependent. This offers an explanation of why accelerating Benioff strain release is not observed prior to all large earthquakes. The results prove the existence of an underlying evolution in discontinuous elasto-dynamic, systems which is capable of providing a basis for forecasting catastrophic failure and earthquakes.

Simulation of the load-unload response ratio and critical sensitivity in the lattice solid model

The Load-Unload Response Ratio (LURR) method is an intermediate-term earthquake prediction approach that has shown considerable promise. It involves calculating the ratio of a specified energy release measure during loading and unloading where loading and unloading periods are determined from the earth tide induced perturbations in the Coulomb Failure Stress on optimally oriented faults. In the lead-up to large earthquakes, high LURR values are frequently observed a few months or years prior to the event. These signals may have a similar origin to the observed accelerating seismic moment release (AMR) prior to many large earthquakes or may be due to critical sensitivity of the crust when a large earthquake is imminent. As a first step towards studying the underlying physical mechanism for the LURR observations, numerical studies are conducted using the particle based lattice solid model (LSM) to determine whether LURR observations can be reproduced. The model is initialized as a heterogeneous 2-D block made up of random-sized particles bonded by elastic-brittle links. The system is subjected to uniaxial compression from rigid driving plates on the upper and lower edges of the model. Experiments are conducted using both strain and stress control to load the plates. A sinusoidal stress perturbation is added to the gradual compressional loading to simulate loading and unloading cycles and LURR is calculated. The results reproduce signals similar to those observed in earthquake prediction practice with a high LURR value followed by a sudden drop prior to macroscopic failure of the sample. The results suggest that LURR provides a good predictor for catastrophic failure in elastic-brittle systems and motivate further research to study the underlying physical mechanisms and statistical properties of high LURR values. The results provide encouragement for earthquake prediction research and the use of advanced simulation models to probe the physics of earthquakes.

Simulation of the influence of rate- and state-dependent friction on the macroscopic behavior of complex fault zones with the lattice solid model

In order to understand the earthquake nucleation process, we need to understand the effective frictional behavior of faults with complex geometry and fault gouge zones. One important aspect of this is the interaction between the friction law governing the behavior of the fault on the microscopic level and the resulting macroscopic behavior of the fault zone. Numerical simulations offer a possibility to investigate the behavior of faults on many different scales and thus provide a means to gain insight into fault zone dynamics on scales which are not accessible to laboratory experiments. Numerical experiments have been performed to investigate the influence of the geometric configuration of faults with a rate- and state-dependent friction at the particle contacts on the effective frictional behavior of these faults. The numerical experiments are designed to be similar to laboratory experiments by DIETERICH and KILGORE (1994) in which a slide-hold-slide cycle was performed between two blocks of material and the resulting peak friction was plotted vs. holding time. Simulations with a flat fault without a fault gouge have been performed to verify the implementation. These have shown close agreement with comparable laboratory experiments. The simulations performed with a fault containing fault gouge have demonstrated a strong dependence of the critical slip distance D-c on the roughness of the fault surfaces and are in qualitative agreement with laboratory experiments.

Electron exchange coupling for single-donor solid-state spin qubits

Intervalley interference between degenerate conduction band minima has been shown to lead to oscillations in the exchange energy between neighboring phosphorus donor electron states in silicon [B. Koiller, X. Hu, and S. Das Sarma, Phys. Rev. Lett. 88, 027903 (2002); Phys. Rev. B 66, 115201 (2002)]. These same effects lead to an extreme sensitivity of the exchange energy on the relative orientation of the donor atoms, an issue of crucial importance in the construction of silicon-based spin quantum computers. In this article we calculate the donor electron exchange coupling as a function of donor position incorporating the full Bloch structure of the Kohn-Luttinger electron wave functions. It is found that due to the rapidly oscillating nature of the terms they produce, the periodic part of the Bloch functions can be safely ignored in the Heitler-London integrals as was done by Koiller, Hu, and Das Sarma, significantly reducing the complexity of calculations. We address issues of fabrication and calculate the expected exchange coupling between neighboring donors that have been implanted into the silicon substrate using an 15 keV ion beam in the so-called top down fabrication scheme for a Kane solid-state quantum computer. In addition, we calculate the exchange coupling as a function of the voltage bias on control gates used to manipulate the electron wave functions and implement quantum logic operations in the Kane proposal, and find that these gate biases can be used to both increase and decrease the magnitude of the exchange coupling between neighboring donor electrons. The zero-bias results reconfirm those previously obtained by Koiller, Hu, and Das Sarma.

Composition and quality of harvest residues and soil organic matter under windrow residue management in young hoop pine plantations as revealed by solid-state C-13 NMR spectroscopy

Solid-state C-13 nuclear magnetic resonance (NMR) with cross-polarisation (CP) and magic-angle-spinning (MAS) was used to: (a) examine the changes in carbon (C) composition of windrowed harvest residues during the first 3 years of hoop pine plantations in subtropical Australia; (b) assess the impacts of windrowed harvest residues on soil organic matter (SOM) composition and quality in the 0-10 cm soil layer. Harvest residues were collected from 0-, 1-, 2- and 3-year-old windrows of ca. 2.5 m width (15 m apart for 0-, 1- and 2-year-old sites and 10 m apart for 3-year-old site). Soils from the 0 to 10 cm soil layer were collected from the 1-, 2- and 3-year-old sites. The 13C NMR spectra of the harvest residues indicated the presence of lignin in the hoop pine wood, foliage and newly incorporated organic matter (NIOM). Condensed tannin structures were found in the decay-resistant bark, small wood and foliage, but were absent in other residue components and SOM. The NMR spectra of small wood samples contained condensed tannin structures because the outer layer of bark was not removed. NIOM showed a shift from foliage-like structures (celluloses) to lignin-type structures, indicating an incorporation of woody residues from the decomposing harvest residues. Suberins were also present in the small wood, foliage and bark. The 13C CP NMR spectra of SOM indicated that in areas where windrows were present, SOM did not show compositional changes. However, an increase in SOM quality under the windrows in the second year after their formation as characterised by the alkyl C/O-alkyl C (A/O-A) ratio was mainly due to inputs from the decomposition of the labile, readily available components of the windrowed harvest residues. (C) 2002 Published by Elsevier Science B.V.

Characterisation of grafted supports used for solid-phase synthesis

A series of 'pellicular' type supports were fabricated by direct gamma-radiation-mediated graft polymerisation of styrene onto polypropylene, followed by aminomethylation. Raman spectroscopy was used for measuring the level of penetration of polystyrene graft into polypropylene, and other structural features such as density of graft and depth of functionalisation. The kinetics of the coupling of fluorenylmethylcarbamate (Fmoc)-labelled amino acids, to the aminomethylated polystyrene grafts have been measured by UV absorption followed cleavage of the Fmoc chromophore. The Raman spectroscopy results showed that for this series of experiments the calculated rate coefficient for coupling of Fmoc-labelled amino acids was primarily dependent on graft thickness, but was also influenced by the proportion of polystyrene graft to polypropylene. In general, it was also shown that with increasing loading capacity of support the calculated rate coefficient for amino-acid coupling decreased correspondingly. In addition, a support that had both a high rate coefficient and a high loading capacity was prepared from polypropylene base material with a co-continuous porous structure (high surface area). (C) 2003 Society of Chemical Industry.

Use of confocal scanning laser microscopy to measure the concentrations of aerial and penetrative hyphae during growth of Rhizopus oligosporus on a solid surface

In order to develop a method for use in investigations of spatial biomass distribution in solid-state fermentation systems, confocal scanning laser microscopy was used to determine the concentrations of aerial and penetrative biomass against height and depth above and below the substrate surface, during growth of Rhizopus oligosporus on potato dextrose agar. Penetrative hyphae had penetrated to a depth of 0.445 cm by 64 h and showed rhizoid morphology, in which the maximum biomass concentration, of 4.45 mg dry wt cm(-3), occurred at a depth of 0.075 cm. For aerial biomass the maximum density of 39.54 mg dry wt(-3) occurred at the substrate surface. For both aerial and penetrative biomass, there were two distinct regions in which the biomass concentration decayed exponentially with distance from the surface. For aerial biomass, the first exponential decay region was up to 0.1 cm height. The second region above the height of 0.1 cm corresponded to that in which sporangiophores dominated. This work lays the foundation for deeper studies into what controls the growth of fungal hyphae above and below the surfaces of solid substrates. (C) Wiley Periodicals, Inc.