XSophe, a computer simulation software suite for the analysis of electron paramagnetic resonance spectra

The XSophe computer simulation software suite consisting of a daemon, the XSophe interface and the computational program Sophe is a state of the art package for the simulation of electron paramagnetic resonance spectra. The Sophe program performs the computer simulation and includes a number of new technologies including; the SOPHE partition and interpolation schemes, a field segmentation algorithm, homotopy, parallelisation and spectral optimisation. The SOPHE partition and interpolation scheme along with a field segmentation algorithm greatly increases the speed of simulations for most systems. Multidimensional homotopy provides an efficient method for accurately tracing energy levels and hence tracing transitions in the presence of energy level anticrossings and looping transitions and allowing computer simulations in frequency space. Recent enhancements to Sophe include the generalised treatment of distributions of orientational parameters, termed the mosaic misorientation linewidth model and a faster more efficient algorithm for the calculation of resonant field positions and transition probabilities. For complex systems the parallelisation enables the simulation of these systems on a parallel computer and the optimisation algorithms in the suite provide the experimentalist with the possibility of finding the spin Hamiltonian parameters in a systematic manner rather than a trial-and-error process. The XSophe software suite has been used to simulate multifrequency EPR spectra (200 MHz to 6 00 GHz) from isolated spin systems (S > ~½) and coupled centres (Si, Sj _> I/2). Griffin, M.; Muys, A.; Noble, C.; Wang, D.; Eldershaw, C.; Gates, K.E.; Burrage, K.; Hanson, G.R."XSophe, a Computer Simulation Software Suite for the Analysis of Electron Paramagnetic Resonance Spectra", 1999, Mol. Phys. Rep., 26, 60-84.

Stress corrosion cracking and hydrogen embrittlement of an Al-Zn-Mg-Cu alloy

The age hardening, stress corrosion cracking (SCC) and hydrogen embrittlement (HE) of an Al-Zn-Mg-Cu 7175 alloy were investigated experimentally. There were two peak-aged states during ageing. For ageing at 413 K, the strength of the second peak-aged state was slightly higher than that of the first one, whereas the SCC susceptibility was lower, indicating that it is possible to heat treat 7175 to high strength and simultaneously to have high SCC resistance. The SCC susceptibility increased with increasing Mg segregation at the grain boundaries. Hydrogen embrittlement (HE) increased with increased hydrogen charging and decreased with increasing ageing time for the same hydrogen charging conditions. Computer simulations were carried out of (a) the Mg grain boundary segregation using the embedded atom method and (b) the effect of Mg and H segregation on the grain boundary strength using a quasi-chemical approach. The simulations showed that (a) Mg grain boundary segregation in Al-Zn-Mg-Cu alloys is spontaneous, (b) Mg segregation decreases the grain boundary strength, and (c) H embrittles the grain boundary more seriously than does Mg. Therefore, the SCC mechanism of Al-Zn-Mg Cu alloys is attributed to the combination of HE and Mg segregation induced grain boundary embrittlement. (C) 2004 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

PVT behaviour of fluids with potential models optimised for phase equilibria

NPT and NVT Monte Carlo simulations are applied to models for methane and water to predict the PVT behaviour of these fluids over a wide range of temperatures and pressures. The potential models examined in this paper have previously been presented in the literature with their specific parameters optimised to fit phase coexistence data. The exponential-6 potential for methane gives generally good prediction of PVT behaviour over the full range of temperature and pressures studied with the only significant deviation from experimental data seen at high temperatures and pressures. The NSPCE water model shows very poor prediction of PVT behaviour, particularly at dense conditions. To improve this. the charge separation in the NSPCE model is varied with density. Improvements for vapour and liquid phase PVT predictions are achieved with this variation. No improvement was found in the prediction of the oxygen-oxygen radial distribution by varying charge separation under dense phase conditions. (C) 2004 Elsevier B.V. All rights reserved.

Effect of surface-perturbed intermolecular interaction on adsorption of simple gases on a graphitized carbon surface

In this paper, we investigate the effect of the solid surface on the fluid-fluid intermolecular potential energy. This modified fluid-fluid interaction energy due to the inducement of a solid surface is used in the grand canonical Monte Carlo (GCMC) simulation of various noble gases, nitrogen, and methane on graphitized thermal carbon black. This effect is such that the effective interaction potential energy between two particles close to surface is less than the potential energy if the solid substrate is not present. With this modification the GCMC simulation results agree extremely well with the experimental data over a wide range of pressures while the simulation results with the unmodified potential energy give rise to a shoulder near the neighborhood of monolayer coverage and the significant overprediction of the second and higher layer coverages. In particular the unmodified GCMC results exhibit very sharp change in those higher layers while the experimental data have a much gradual change in the uptake. We will illustrate this theory with adsorption data of argon, xenon, neon, nitrogen, and methane on graphitized thermal carbon black.

Adsorption of quadrupolar, diatomic nitrogen onto graphitized thermal carbon black and in slit-shaped carbon pores. Effects of surface mediation

In this paper, we study the effect of solid surface mediation on the intermolecular potential energy of nitrogen, and its impact on the adsorption of nitrogen on a graphitized carbon black surface and in carbon slit-shaped pores. This effect arises from the lower effective interaction potential energy between two particles close to the surface compared to the potential energy of the same two particles when they are far away from the surface. A simple equation is proposed to calculate the reduction factor and this is used in the Grand Canonical Monte Carlo (GCMC) simulation of nitrogen adsorption on graphitized thermal carbon black. With this modification, the GCMC simulation results agree extremely well with the experimental data over a wide range of pressure; the simulation results with the original potential energy (i.e. no surface mediation) give rise to a shoulder in the neighbourhood of monolayer coverage and a significant over-prediction of the second and higher layer coverages. The influence of this surface mediation on the dependence of the pore-filling pressure on the pore width is also studied. It is shown that such surface mediation has a significant effect on the pore-filling pressure. This implies that the use of the local isotherms obtained from the potential model without surface mediation could give rise to a serious error in the determination of the pore-size distribution.

Effects of surface heterogeneity on the adsorption of nitrogen on graphitized thermal carbon black

In this paper, we study the surface heterogeneity and the surface mediation on the intermolecular potential energy for nitrogen adsorption on graphitized thermal carbon black (GTCB). The surface heterogeneity is modeled as the random distribution of effective carbonyl functional groups on the graphite surface. The molecular parameters and the discrete charges of this carbonyl group are taken from Jorgensen, et al. (J. Am. Chem. Soc., (1984) 106, 6638) while those for nitrogen (dispersive parameters and discrete charges) are taken from Murthy et al. (Mol. Phys., (1983) 50, 531) in our Grand Canonical Monte Carlo (GCMC) simulation. The solid surface mediation in the reduction of intermolecular potential energy between two fluid molecules was taken from a recent work by Do et al. (Langmuir, (2004) 20, 7623). Our simulation results accounting for the surface heterogeneity and surface mediation on intermolecular potential energy were compared with the experimental data of nitrogen at 77 and 90 K. The solid-fluid dispersive parameters are determined from the Lorentz-Berthelot (LB) rule. The fraction of the graphite surface covered with carbonyl functional groups was then derived from the consideration of the Henry constant, and for the data of Kruk et al. (Langmuir, (1999) 15, 1435) we have found that 1% of their GTCB surface is covered with effective carbonyl functional groups. The damping constant, due to surface mediation, was determined from the consideration of the portion of the adsorption isotherm where the first layer is being completed, and it was found to take a value of 0.0075. With these parameters, we have found that the GCMC simulation results describe the data over the complete range of pressure substantially better than any other MC models in the literature. The implication of this work is demonstrated with local adsorption isotherms of 10 and 20 A slit pores. One was obtained without allowance for surface mediation, while the other correctly accounts for these factors. The two local isotherms differ substantially, and the implication is that if we used incorrect local isotherms (i.e. without the surface mediation) the pore size distribution would be incorrectly derived.

Cooperative and competitive adsorption of ethylene, ethane, nitrogen and argon on graphitized carbon black and in slit pores

In this paper we investigate the mixture adsorption of ethylene, ethane, nitrogen and argon on graphitized thermal carbon black and in slit pores by means of the Grand Canonical Monte Carlo simulations. Pure component adsorption isotherms on graphitized thermal carbon black are first characterized with the GCMC method, and then mixture simulations are carried out over a wide range of pore width, temperature, pressure and composition to investigate the cooperative and competitive adsorption of all species in the mixture. Results of mixture simulations are compared with the experimental data of ethylene and ethane (Friederich and Mullins, 1972) on Sterling FTG-D5 (homogeneous carbon black having a BET surface area of 13 m(2)/g) at 298 K and a pressure range of 1.3-93 kPa. Because of the co-operative effect, the Henry constant determined by the traditional chromatography method is always greater than that obtained from the volumetric method.

Adsorption of flexible n-alkane on graphitized thermal carbon black: analysis of adsorption isotherm by means of GCMC simulation

In this paper, we investigate the suitability of the grand canonical Monte Carlo in the description of adsorption equilibria of flexible n-alkane (butane, pentane and hexane) on graphitized thermal carbon black. Potential model of n-alkane of Martin and Siepmann (J. Phys. Chem. 102 (1998) 2569) is employed in the simulation, and we consider the flexibility of molecule in the simulation. By this we study two models, one is the fully flexible molecular model in which n-alkane is subject to bending and torsion, while the other is the rigid molecular model in which all carbon atoms reside on the same plane. It is found that (i) the adsorption isotherm results of these two models are close to each other, suggesting that n-alkane model behaves mostly as rigid molecules with respect to adsorption although the isotherm for longer chain n-hexane is better described by the flexible molecular model (ii) the isotherms agree very well with the experimental data at least up to two layers on the surface.

Analysis of cardiac ventricular wall motion based on a three-dimensional electromechanical biventricular model

This paper describes a biventricular model, which couples the electrical and mechanical properties of the heart, and computer simulations of ventricular wall motion and deformation by means of a biventricular model. In the constructed electromechanical model, the mechanical analysis was based on composite material theory and the finite-element method; the propagation of electrical excitation was simulated using an electrical heart model, and the resulting active forces were used to calculate ventricular wall motion. Regional deformation and Lagrangian strain tensors were calculated during the systole phase. Displacements, minimum principal strains and torsion angle were used to describe the motion of the two ventricles. The simulations showed that during the period of systole, (1) the right ventricular free wall moves towards the septum, and at the same time, the base and middle of the free wall move towards the apex, which reduces the volume of the right ventricle; the minimum principle strain (E3) is largest at the apex, then at the middle of the free wall and its direction is in the approximate direction of the epicardial muscle fibres; (2) the base and middle of the left ventricular free wall move towards the apex and the apex remains almost static; the torsion angle is largest at the apex; the minimum principle strain E3 is largest at the apex and its direction on the surface of the middle wall of the left ventricle is roughly in the fibre orientation. These results are in good accordance with results obtained from MR tagging images reported in the literature. This study suggests that such an electromechanical biventricular model has the potential to be used to assess the mechanical function of the two ventricles, and also could improve the accuracy ECG simulation when it is used in heart torso model-based body surface potential simulation studies.

Adsorption of benzene on graphitized thermal carbon black: Reduction of the quadrupole moment in the adsorbed phase

The performance of intermolecular potential models on the adsorption of benzene on graphitized thermal carbon black at various temperatures is investigated. Two models contain only dispersive sites, whereas the other two models account explicitly for the dispersive and electrostatic sites. Using numerous data in the literature on benzene adsorption on graphitized thermal carbon black at various temperatures, we have found that the effect of surface mediation on interaction between adsorbed benzene molecules must be accounted for to describe correctly the adsorption isotherm as well as the isosteric heat. Among the two models with partial charges tested, the WSKS model of Wick et at. I that has only six dispersive sites and three discrete partial charges is better than the very expensive all-atom model of Jorgensen and Severance.(2) Adsorbed benzene molecules on graphitized thermal carbon black have a complex orientation with respect to distance from the surface and also with respect to loading. At low loadings, they adopt the parallel configuration relative to the graphene surface, whereas at higher loadings (still less than monolayer coverage) some molecules adopt a slant orientation to maximize the fluid-fluid interaction. For loadings in the multilayer region, the orientation of molecules in the first layer is influenced by the presence of molecules in the second layer. The data that are used in this article come from the work of Isirikyan and Kiselev,(3) Pierotti and Smallwood,(4) Pierce and Ewing,(5) Belyakova, Kiselev, and Kovaleva,(6) and Carrott et al.(7)

Studying phenotypic evolution using multivariate quantitative genetics

Quantitative genetics provides a powerful framework for studying phenotypic evolution and the evolution of adaptive genetic variation. Central to the approach is G, the matrix of additive genetic variances and covariances. G summarizes the genetic basis of the traits and can be used to predict the phenotypic response to multivariate selection or to drift. Recent analytical and computational advances have improved both the power and the accessibility of the necessary multivariate statistics. It is now possible to study the relationships between G and other evolutionary parameters, such as those describing the mutational input, the shape and orientation of the adaptive landscape, and the phenotypic divergence among populations. At the same time, we are moving towards a greater understanding of how the genetic variation summarized by G evolves. Computer simulations of the evolution of G, innovations in matrix comparison methods, and rapid development of powerful molecular genetic tools have all opened the way for dissecting the interaction between allelic variation and evolutionary process. Here I discuss some current uses of G, problems with the application of these approaches, and identify avenues for future research.

Transpiration efficiency in a segregating population of sunflower: Inheritance, correlation with other traits and association with hybrid grain yield

Plants incorporate isotopes of carbon into their tissue at different rates because of discrimination against 13C relative to 12C during photosynthesis. This difference in discrimination has been negatively correlated with transpiration efficiency (TE) in many C3 species and so, carbon isotope discrimination (Δ) of leaf tissues has been proposed as a potential tool for selecting genotypes with improved performance under water limited conditions. The relationship between Δ and TE in sunflower has been described previously using diverse genotypes, but this relationship has not been investigated with material selected from a segregating population. In this study, the TE of twenty recombinant inbred lines from a population (HAR4 x SA52) segregating for Δ was evaluated in a rainout shelter experiment. A strong negative genetic correlation between TE and Δ was observed (rg = -0.58), confirming previous studies of sunflower with unrelated lines. In addition, TE was strongly correlated to plant height at the final harvest (rg = 0.64) and TDW (rg = 0.58), and moderately correlated to SLW (rg = 0.46) and SPAD (rg = 0.21) but not leaf number (rg = 0.02). Estimates of narrow sense heritability of TE and Δ were very high (0.82 and 0.77, respectively) suggesting that selection for these traits could occur in early generations of segregating populations. Grain yield evaluations under field conditions of hybrids contrasting for Δ showed that low Δ (high TE) hybrids had a yield advantage between 22-35% in dry environments where the yield was less than 2t/ha. While this level of yield advantage may not be realized in commercial breeding programs, computer simulations suggest that 10-15% yield improvements may be possible. Low Δ material selected from the population HAR4 x SA52 has been distributed to private seed companies for further evaluation.

Effect of line of sight propagation on capacity of an indoor MIMO system

In this paper the performance of a multiple input multiple output (MIMO) wireless communication system operating in an indoor environment, featuring both line of sight (LOS) and non-line of sight (NLOS) signal propagation, is assessed. In the model the scattering objects are assumed to be uniformly distributed in an area surrounding the transmitting and receiving array antennas. Mutual coupling effects in the arrays are treated in an exact manner. However interactions with scattering objects are taken into account via a single bounce approach. Computer simulations are carried out for the system capacity for varying inter-element spacing in the receiving array for assumed values of LOS/NLOS power fraction and signal to noise ratio (SNR).

Tracking the states of a nonlinear system in the weight-space of a feed-forward neural network

Nonlinear, non-stationary signals are commonly found in a variety of disciplines such as biology, medicine, geology and financial modeling. The complexity (e.g. nonlinearity and non-stationarity) of such signals and their low signal to noise ratios often make it a challenging task to use them in critical applications. In this paper we propose a new neural network based technique to address those problems. We show that a feed forward, multi-layered neural network can conveniently capture the states of a nonlinear system in its connection weight-space, after a process of supervised training. The performance of the proposed method is investigated via computer simulations.

Assessing an FPGA implemented MIMO testbed with the use of channel emulator

In this paper, a channel emulator for assessing the performance of MIMO testbed implemented in a field programmable gate array technology is described. The FPGA based MIMO system includes a signal generator, modulation/demodulation and space time coding/decoding modules. The emulator uses information about a wireless channel from computer simulations or actual measurements. In simulations, a single bounce scattering model for an indoor environment is applied. The generated data is stored in the FPGA board. The tests are performed for a 2times2 MIMO system that uses Alamouti scheme for space coding/decoding. The performed tests show proper operation of the FPGA implemented MIMO testbed. Good agreement between the results using measured and simulated channel data is obtained.