Lattice damage produced in GaN by swift heavy ions

Wurtzite GaN epilayers bombarded at 300 K with 200 MeV Au-197(16+) ions are studied by a combination of transmission electron microscopy (TEM) and Rutherford backscattering/channeling spectrometry (RBS/C). Results reveal the formation of near-continuous tracks propagating throughout the entire similar to1.5-mum-thick GaN film. These tracks, similar to100 Angstrom in diameter, exhibit a large degree of structural disordering but do not appear to be amorphous. Throughout the bombarded epilayer, high-resolution TEM reveals planar defects which are parallel to the basal plane of the GaN film. The gross level of lattice disorder, as measured by RBS/C, gradually increases with increasing ion fluence up to similar to10(13) cm(-2). For larger fluences, delamination of the nitride film from the sapphire substrate occurs. Based on these results, physical mechanisms of the formation of lattice disorder in GaN in such a high electronic stopping power regime are discussed. (C) 2004 American Institute of Physics.

Study of hexane adsorption in nanoporous MCM-41 silica

We study here the adsorption of hexane on nanoporous MCM-41 silica at 303, 313, and 323 K, for various pore diameters between 2.40 and 4.24 nm. Adsorption equilibria, measured thermogravimetrically, show that all the isotherms, that are somewhat akin to those of type V, exhibit remarkably sharp capillary adsorption phase transition steps and are reversible. The position of the phase transition step gradually shifts from low to high relative pressure with an increase in the temperature as well as the pore sizes. The isosteric heats of adsorption derived from the equilibrium information using the Clapeyron equation reveal a gradual decrease with increasing adsorbed amount because of the surface heterogeneity but approach a constant value near the phase transition. A decrease in the pore size results in an increase in the isosteric heat of adsorption because of the increased dispersion forces. A simple strategy, based on the Broekhoff and De Boer adsorption theory, successfully interprets the hexane adsorption isotherms for the different pore size MCM-41 samples. The parameters of an empirical expression, used to represent the potential of interaction between the adsorbate and adsorbent, are obtained by fitting the monolayer region prior to capillary condensation and the experimental phase transition simultaneously, for some pore sizes. Subsequently, the parameters are used to predict the adsorption isotherm on other pore size samples, which showed good agreement with experimental data.

Statistical tests of load-unload response ratio signals by lattice solid model: Implication to tidal triggering and earthquake prediction

Statistical tests of Load-Unload Response Ratio (LURR) signals are carried in order to verify statistical robustness of the previous studies using the Lattice Solid Model (MORA et al., 2002b). In each case 24 groups of samples with the same macroscopic parameters (tidal perturbation amplitude A, period T and tectonic loading rate k) but different particle arrangements are employed. Results of uni-axial compression experiments show that before the normalized time of catastrophic failure, the ensemble average LURR value rises significantly, in agreement with the observations of high LURR prior to the large earthquakes. In shearing tests, two parameters are found to control the correlation between earthquake occurrence and tidal stress. One is, A/(kT) controlling the phase shift between the peak seismicity rate and the peak amplitude of the perturbation stress. With an increase of this parameter, the phase shift is found to decrease. Another parameter, AT/k, controls the height of the probability density function (Pdf) of modeled seismicity. As this parameter increases, the Pdf becomes sharper and narrower, indicating a strong triggering. Statistical studies of LURR signals in shearing tests also suggest that except in strong triggering cases, where LURR cannot be calculated due to poor data in unloading cycles, the larger events are more likely to occur in higher LURR periods than the smaller ones, supporting the LURR hypothesis.

A parallel implementation of the lattice solid model for the simulation of rock mechanics and earthquake dynamics

The Lattice Solid Model has been used successfully as a virtual laboratory to simulate fracturing of rocks, the dynamics of faults, earthquakes and gouge processes. However, results from those simulations show that in order to make the next step towards more realistic experiments it will be necessary to use models containing a significantly larger number of particles than current models. Thus, those simulations will require a greatly increased amount of computational resources. Whereas the computing power provided by single processors can be expected to increase according to Moore's law, i.e., to double every 18-24 months, parallel computers can provide significantly larger computing power today. In order to make this computing power available for the simulation of the microphysics of earthquakes, a parallel version of the Lattice Solid Model has been implemented. Benchmarks using large models with several millions of particles have shown that the parallel implementation of the Lattice Solid Model can achieve a high parallel-efficiency of about 80% for large numbers of processors on different computer architectures.

Diffusion of n-decane in mesoporous MCM-41 silicas

We investigate here the diffusion of n-decane in nanoporous MCM-41 silicas with pore diameters between 3.0 and 4.3 nm, and at various temperatures and purge flow rates, by the Zero Length Column method. A complete-time-range analysis of desorption curves is used to derive the diffusion coefficient, and the effect of pore size, purge flow rate and temperature on the diffusion character is systematically studied. The results show that the calculated low-coverage diffusivity values are strongly dependent on temperature but only weakly dependent on pore size. The study reveals that transport is controlled by intracrystalline diffusion and dominated by sorbate-sorbent interaction, with the experimental isosteric heat matching the potential energy of flat-lying n-decane molecules on the surface, determined using a united atom model. The diffusion activation energy and adsorption isosteric heat at zero loading for the different pore size MCM-41 samples vary in a narrow range respectively, and their ratio is essentially constant over the pore size range studied. The study shows that the ZLC method is an effective tool to investigate the diffusion kinetics of hydrocarbons in mesoporous MCM-41 materials. (c) 2005 Elsevier Inc. All rights reserved.

Adsorption of argon and nitrogen in cylindrical pores of MCM-41 materials: application of density functional theory

Adsorption of argon and nitrogen at their respective boiling points in cylindrical pores of MCM-41 type silica-like adsorbents is studied by means of a non-local density functional theory (NLDFT), which is modified to deal with amorphous solids. By matching the theoretical results of the pore filling pressure versus pore diameter against the experimental data, we arrive at a conclusion that the adsorption branch (rather than desorption) corresponds to the true thermodynamic equilibrium. If this is accepted, we derive the optimal values for the solid–fluid molecular parameters for the system amorphous silica–Ar and amorphous silica–N2, and at the same time we could derive reliably the specific surface area of non-porous and mesoporous silica-like adsorbents, without a recourse to the BET method. This method is then logically extended to describe the local adsorption isotherms of argon and nitrogen in silica-like pores, which are then used as the bases (kernel) to determine the pore size distribution. We test this with a number of adsorption isotherms on the MCM-41 samples, and the results are quite realistic and in excellent agreement with the XRD results, justifying the approach adopted in this paper.

Adsorption of pure fluids in activated carbon with a modified lattice gas model

In this article we study the effects of adsorbed phase compression, lattice structure, and pore size distribution on the analysis of adsorption in microporous activated carbon. The lattice gas approach of Ono-Kondo is modified to account for the above effects. Data of nitrogen adsorption at 77 K onto a number of activated carbon samples are analyzed to investigate the pore filling pressure versus pore width, the packing effect, and the compression of the adsorbed phase. It is found that the PSDs obtained from this analysis are comparable to those obtained by the DFT method. The discrete nature of the PSDs derived from the modified lattice gas theory is due to the inherent assumption of discrete layers of molecules. Nevertheless, it does provide interesting information on the evolution of micropores during the activation process.

Improvement of the Derjaguin-Broekhoff-de Boer theory for capillary condensation/evaporation of nitrogen in mesoporous systems and its implications for pore size analysis of MCM-41 silicas and related materials

In this work, we propose an improvement of the classical Derjaguin-Broekhoff-de Boer (DBdB) theory for capillary condensation/evaporation in mesoporous systems. The primary idea of this improvement is to employ the Gibbs-Tolman-Koenig-Buff equation to predict the surface tension changes in mesopores. In addition, the statistical film thickness (so-called t-curve) evaluated accurately on the basis of the adsorption isotherms measured for the MCM-41 materials is used instead of the originally proposed t-curve (to take into account the excess of the chemical potential due to the surface forces). It is shown that the aforementioned modifications of the original DBdB theory have significant implications for the pore size analysis of mesoporous solids. To verify our improvement of the DBdB pore size analysis method (IDBdB), a series of the calcined MCM-41 samples, which are well-defined materials with hexagonally ordered cylindrical mesopores, were used for the evaluation of the pore size distributions. The correlation of the IDBdB method with the empirically calibrated Kruk-Jaroniec-Sayari (KJS) relationship is very good in the range of small mesopores. So, a major advantage of the IDBdB method is its applicability for small mesopores as well as for the mesopore range beyond that established by the KJS calibration, i.e., for mesopore radii greater than similar to4.5 nm. The comparison of the IDBdB results with experimental data reported by Kruk and Jaroniec for capillary condensation/evaporation as well as with the results from nonlocal density functional theory developed by Neimark et al. clearly justifies our approach. Note that the proposed improvement of the classical DBdB method preserves its original simplicity and simultaneously ensures a significant improvement of the pore size analysis, which is confirmed by the independent estimation of the mean pore size by the powder X-ray diffraction method.

Temperature dependence of the crystal structure and g-values of [(HC(Ph2PO)(3))(2)Cu] (ClO4)(2)center dot 2H(2)O: Influence of dynamic Jahn-Teller coupling and lattice strain interactions

The temperature dependence of the X-ray crystal structure and powder EPR spectrum of [(HC(Ph2PO)(3))(2)CU]-(ClO4)(2)center dot 2H(2)O is reported, and the structure at room temperature confirms that reported previously. Below similar to 100 K, the data imply a geometry with near elongated tetragonal symmetry for the [(HC(Ph2PO)(3))(2)Cu](2+) complex, but on warming the two higher Cu-O bond lengths and g-values progressively converge, and by 340 K the bond lengths correspond to a compressed tetragonal geometry. The data may be interpreted satisfactorily assuming an equilibrium among the energy levels of a Cu-O-6 polyhedron subjected to Jahn-Teller vibronic coupling and a lattice strain. However, agreement with the experiment is obtained only if the orthorhombic component of the lattice strain decreases to a negligible value as the temperature approaches 340 K.

Comparison of nitrogen adsorption at 77 K on non-porous silica and pore wall of MCM-41 materials by means of density functional theory

Nitrogen adsorption on a surface of a non-porous reference material is widely used in the characterization. Traditionally, the enhancement of solid-fluid potential in a porous solid is accounted for by incorporating the surface curvature into the solid-fluid Potential of the flat reference surface. However, this calculation procedure has not been justified experimentally. In this paper, we derive the solid-fluid potential of mesoporous MCM-41 solid by using solely the adsorption isotherm of that solid. This solid-fluid potential is then compared with that of the non-porous reference surface. In derivation of the solid-fluid potential for both reference surface and mesoporous MCM-41 silica (diameter ranging front 3 to 6.5 nm) we employ the nonlocal density functional theory developed for amorphous solids. It is found that, to out, surprise, the solid-fluid potential of a porous solid is practically the same as that for the reference surface, indicating that there is no enhancement due to Surface curvature. This requires further investigations to explain this unusual departure from our conventional wisdom of curvature-induced enhancement. Accepting the curvature-independent solid-fluid potential derived from the non-porous reference surface, we analyze the hysteresis features of a series of MCM-41 samples. (c) 2005 Elsevier Inc. All rights reserved.

Influence of lattice interactions on the jahn-teller distortion of the [Cu(H2O)(6)](2+) ion: Dependence of the crystal structure of K-2[Cu(H2O)(6)](SO4)(2x)(SeO4)(2-2x) upon the sulfate/selenate ratio

The temperature dependence of the structure of the mixed-anion Tutton salt K-2[Cu(H2O)(6)](SO4)(2x)(SeO4)(2-2x) has been determined for crystals with 0, 17, 25, 68, 78, and 100% sulfate over the temperature range of 85-320 K. In every case, the [Cu(H2O)(6)](2+) ion adopts a tetragonally elongated coordination geometry with an orthorhombic distortion. However, for the compounds with 0, 17, and 25% sulfate, the long and intermediate bonds occur on a different pair of water molecules from those with 68, 78, and 100% sulfate. A thermal equilibrium between the two forms is observed for each crystal, with this developing more readily as the proportions of the two counterions become more similar. Attempts to prepare a crystal with approximately equal amounts of sulfate and selenate were unsuccessful. The temperature dependence of the bond lengths has been analyzed using a model in which the Jahn-Teller potential surface of the [Cu(H2O)(6)](2+) ion is perturbed by a lattice-strain interaction. The magnitude and sign of the orthorhombic component of this strain interaction depends on the proportion of sulfate to selenate. Significant deviations from Boltzmann statistics are observed for those crystals exhibiting a large temperature dependence of the average bond lengths, and this may be explained by cooperative interactions between neighboring complexes.

Testing predictions of macroscopic binary diffusion coefficients using lattice models with site heterogeneity

Quantitatively predicting mass transport rates for chemical mixtures in porous materials is important in applications of materials such as adsorbents, membranes, and catalysts. Because directly assessing mixture transport experimentally is challenging, theoretical models that can predict mixture diffusion coefficients using Only single-component information would have many uses. One such model was proposed by Skoulidas, Sholl, and Krishna (Langmuir, 2003, 19, 7977), and applications of this model to a variety of chemical mixtures in nanoporous materials have yielded promising results. In this paper, the accuracy of this model for predicting mixture diffusion coefficients in materials that exhibit a heterogeneous distribution of local binding energies is examined. To examine this issue, single-component and binary mixture diffusion coefficients are computed using kinetic Monte Carlo for a two-dimensional lattice model over a wide range of lattice occupancies and compositions. The approach suggested by Skoulidas, Sholl, and Krishna is found to be accurate in situations where the spatial distribution of binding site energies is relatively homogeneous, but is considerably less accurate for strongly heterogeneous energy distributions.

Excitation spectra of the spin-1/2 triangular-lattice Heisenberg antiferromagnet

We use series expansion methods to calculate the dispersion relation of the one-magnon excitations for the spin-(1)/(2) triangular-lattice nearest-neighbor Heisenberg antiferromagnet above a three-sublattice ordered ground state. Several striking features are observed compared to the classical (large-S) spin-wave spectra. Whereas, at low energies the dispersion is only weakly renormalized by quantum fluctuations, significant anomalies are observed at high energies. In particular, we find rotonlike minima at special wave vectors and strong downward renormalization in large parts of the Brillouin zone, leading to very flat or dispersionless modes. We present detailed comparison of our calculated excitation energies in the Brillouin zone with the spin-wave dispersion to order 1/S calculated recently by Starykh, Chubukov, and Abanov [Phys. Rev. B74, 180403(R) (2006)]. We find many common features but also some quantitative and qualitative differences. We show that at temperatures as low as 0.1J the thermally excited rotons make a significant contribution to the entropy. Consequently, unlike for the square lattice model, a nonlinear sigma model description of the finite-temperature properties is only applicable at temperatures < 0.1J. Finally, we review recent NMR measurements on the organic compound kappa-(BEDT-TTF)(2)Cu-2(CN)(3). We argue that these are inconsistent with long-range order and a description of the low-energy excitations in terms of interacting magnons, and that therefore a Heisenberg model with only nearest-neighbor exchange does not offer an adequate description of this material.

Implementation of particle-scale rotation in the 3-d Lattice Solid Model

In this study, 3-D Lattice Solid Model (LSMearth or LSM) was extended by introducing particle-scale rotation. In the new model, for each 3-D particle, we introduce six degrees of freedom: Three for translational motion, and three for orientation. Six kinds of relative motions are permitted between two neighboring particles, and six interactions are transferred, i.e., radial, two shearing forces, twisting and two bending torques. By using quaternion algebra, relative rotation between two particles is decomposed into two sequence-independent rotations such that all interactions due to the relative motions between interactive rigid bodies can be uniquely decided. After incorporating this mechanism and introducing bond breaking under torsion and bending into the LSM, several tests on 2-D and 3-D rock failure under uni-axial compression are carried out. Compared with the simulations without the single particle rotational mechanism, the new simulation results match more closely experimental results of rock fracture and hence, are encouraging. Since more parameters are introduced, an approach for choosing the new parameters is presented.

Parallel 3D Simulation of a Fault Gouge using the Lattice Solid Model

Despite the insight gained from 2-D particle models, and given that the dynamics of crustal faults occur in 3-D space, the question remains, how do the 3-D fault gouge dynamics differ from those in 2-D? Traditionally, 2-D modeling has been preferred over 3-D simulations because of the computational cost of solving 3-D problems. However, modern high performance computing architectures, combined with a parallel implementation of the Lattice Solid Model (LSM), provide the opportunity to explore 3-D fault micro-mechanics and to advance understanding of effective constitutive relations of fault gouge layers. In this paper, macroscopic friction values from 2-D and 3-D LSM simulations, performed on an SGI Altix 3700 super-cluster, are compared. Two rectangular elastic blocks of bonded particles, with a rough fault plane and separated by a region of randomly sized non-bonded gouge particles, are sheared in opposite directions by normally-loaded driving plates. The results demonstrate that the gouge particles in the 3-D models undergo significant out-of-plane motion during shear. The 3-D models also exhibit a higher mean macroscopic friction than the 2-D models for varying values of interparticle friction. 2-D LSM gouge models have previously been shown to exhibit accelerating energy release in simulated earthquake cycles, supporting the Critical Point hypothesis. The 3-D models are shown to also display accelerating energy release, and good fits of power law time-to-failure functions to the cumulative energy release are obtained.