Population-based local search for protein folding simulation in the MJ energy model and cubic lattices

We present experimental results on benchmark problems in 3D cubic lattice structures with the Miyazawa-Jernigan energy function for two local search procedures that utilise the pull-move set: (i) population-based local search (PLS) that traverses the energy landscape with greedy steps towards (potential) local minima followed by upward steps up to a certain level of the objective function; (ii) simulated annealing with a logarithmic cooling schedule (LSA). The parameter settings for PLS are derived from short LSA-runs executed in pre-processing and the procedure utilises tabu lists generated for each member of the population. In terms of the total number of energy function evaluations both methods perform equally well, however. PLS has the potential of being parallelised with an expected speed-up in the region of the population size. Furthermore, both methods require a significant smaller number of function evaluations when compared to Monte Carlo simulations with kink-jump moves. (C) 2009 Elsevier Ltd. All rights reserved.

Molecular dynamics simulation model for the quantitative assessment of tool wear during single point diamond turning of cubic silicon carbide

Silicon carbide (SiC) is a material of great technological interest for engineering applications concerning hostile environments where silicon-based components cannot work (beyond 623 K). Single point diamond turning (SPDT) has remained a superior and viable method to harness process efficiency and freeform shapes on this harder material. However, it is extremely difficult to machine this ceramic consistently in the ductile regime due to sudden and rapid tool wear. It thus becomes non trivial to develop an accurate understanding of tool wear mechanism during SPDT of SiC in order to identify measures to suppress wear to minimize operational cost.

In this paper, molecular dynamics (MD) simulation has been deployed with a realistic analytical bond order potential (ABOP) formalism based potential energy function to understand tool wear mechanism during single point diamond turning of SiC. The most significant result was obtained using the radial distribution function which suggests graphitization of diamond tool during the machining process. This phenomenon occurs due to the abrasive processes between these two ultra hard materials. The abrasive action results in locally high temperature which compounds with the massive cutting forces leading to sp3–sp2 order–disorder transition of diamond tool. This represents the root cause of tool wear during SPDT operation of cubic SiC. Further testing led to the development of a novel method for quantitative assessment of the progression of diamond tool wear from MD simulations.

Molecular dynamics simulation of nanoindentation of Fe3C and tetrahedral Fe4C

Study of nanomechanical response of iron carbides is important because presence of iron carbides greatly influences the performance and longevity of steel components. This work contributes to the literature by exploring nanoindentation of nanocrystalline Fe3C and tetrahedral-Fe4C using molecular dynamics simulation. The chemical interactions of iron and carbon were described through an analytical bond order inter-atomic potential (ABOP) energy function. The indentations were performed at an indentation speed of 50 m/sec and a repeat trial was performed at 5 m/sec. Load-displacement (P-h) curve for both these carbides showed residual indentation depth and maximum indentation depth (hf/hmax) ratio to be higher than 0.7 i.e. a circumstance where Oliver and Pharr method was not appropriate to be applied to evaluate the material properties. Alternate evaluation revealed Fe3C to be much harder than Fe4C. Gibbs free energy of formation and radial distribution function, coupled with state of the average local temperature and von Mises stresses indicate the formation of a new phase of iron-carbide. Formation of this newer phase was found to be due to deviatoric strain rather than the high temperature induced in the substrate during nanoindentation

Cyclic nanoindentation and nano-impact fatigue characterisation of functionally graded TiN/TiNi shape memory film

This paper investigates the mechanism of nanoscale fatigue of functionally graded TiN/TiNi films using nano-impact and multiple-loading-cycle nanoindentation tests. The functionally graded films were deposited on silicon substrate, in which TiNi films maintain shape memory and pseudo elastic behavior, while a modified TiN surface layer provides tribological and anti-corrosion properties. Nanomechanical tests were performed to comprehend the localized film performance and failure modes of the functionally graded film using NanoTestTM equipped with Berkovich and conical indenter between 100 μN to 500 mN loads. The loading mechanism and load history are critical to define film failure modes (i.e. backward depth deviation) including the shape memory effect of the functionally graded layer. The results are sensitive to the applied load, loading type (e.g. semi-static, dynamic) and probe geometry. Based on indentation force-depth profiles, depth-time data and post-test surface observations of films, it is concluded that the shape of the nanoindenter is critical in inducing the localized indentation stress and film failure, including shape recovery at the lower load range. Elastic-plastic finite element (FE) simulation during nanoindentation loading indicated that the location of subsurface maximum stress near the interface influences the backward depth deviation type of film failure. A standalone, molecular dynamics simulation was performed with the help of a long range potential energy function to simulate the tensile test of TiN nanowire with two different aspect ratios to investigate the theory of its failure mechanism.

TIME-RESOLVED OPTICAL-EMISSION AND ELECTRON-ENERGY DISTRIBUTION FUNCTION MEASUREMENTS IN RF PLASMAS

Comparisons between experimentally measured time-dependent electron energy distribution functions and optical emission intensities are reported for low-frequency (100 and 400 kHz) radio-frequency driven discharges in argon. The electron energy distribution functions were measured with a time-resolved Langmuir probe system. Time-resolved optical emissions of argon resonance lines at 687.1 and 750.4 nm were determined by photon-counting methods. Known ground-state and metastable-state excitation cross sections were used along with the measured electron energy distribution functions to calculate the time dependence of the optical emission intensity. It was found that a calculation using only the ground-state cross sections gave the best agreement with the time dependence of the measured optical emission. Time-dependent electron density, electron temperature, and plasma potential measurements are also reported.

TIME-RESOLVED ELECTRON-ENERGY DISTRIBUTION FUNCTION MEASUREMENTS IN A PULSED MAGNETIC MULTIPOLE HYDROGEN DISCHARGE

The time evolution of measured plasma parameters, including the electron energy distribution function (EEDF), in the discharge and post-discharge regime of a pulsed hydrogen magnetic multipole plasma is presented. The time necessary for the plasma to reach equilibrium has been established as 160-mu-s. The present results clarify the mechanisms which initiate the discharge. The decay rates of the charged-particle density and energy in the post-discharge have been measured. These measurements indicate that particle transport to the wall is the dominant loss mechanism for both charged-particle density and energy. The time-resolved EEDF is found to be non-Maxwellian in the discharge and Maxwellian in the late post-discharge.

Exchange potential from the common energy denominator approximation for the Kohn-Sham Green's function: Application to (hyper)polarizabilities of molecular chains

An approximate Kohn-Sham (KS) exchange potential v(xsigma)(CEDA) is developed, based on the common energy denominator approximation (CEDA) for the static orbital Green's function, which preserves the essential structure of the density response function. v(xsigma)(CEDA) is an explicit functional of the occupied KS orbitals, which has the Slater v(Ssigma) and response v(respsigma)(CEDA) potentials as its components. The latter exhibits the characteristic step structure with "diagonal" contributions from the orbital densities \psi(isigma)\(2), as well as "off-diagonal" ones from the occupied-occupied orbital products psi(isigma)psi(j(not equal1)sigma). Comparison of the results of atomic and molecular ground-state CEDA calculations with those of the Krieger-Li-Iafrate (KLI), exact exchange (EXX), and Hartree-Fock (HF) methods show, that both KLI and CEDA potentials can be considered as very good analytical "closure approximations" to the exact KS exchange potential. The total CEDA and KLI energies nearly coincide with the EXX ones and the corresponding orbital energies epsilon(isigma) are rather close to each other for the light atoms and small molecules considered. The CEDA, KLI, EXX-epsilon(isigma) values provide the qualitatively correct order of ionizations and they give an estimate of VIPs comparable to that of the HF Koopmans' theorem. However, the additional off-diagonal orbital structure of v(xsigma)(CEDA) appears to be essential for the calculated response properties of molecular chains. KLI already considerably improves the calculated (hyper)polarizabilities of the prototype hydrogen chains H-n over local density approximation (LDA) and standard generalized gradient approximations (GGAs), while the CEDA results are definitely an improvement over the KLI ones. The reasons of this success are the specific orbital structures of the CEDA and KLI response potentials, which produce in an external field an ultranonlocal field-counteracting exchange potential. (C) 2002 American Institute of Physics.

Performance analysis of an improved power-saving medium access protocol for IEEE 802.11 point coordination function WLAN

Wireless enabled portable devices must operate with the highest possible energy efficiency while still maintaining a minimum level and quality of service to meet the user's expectations. The authors analyse the performance of a new pointer-based medium access control protocol that was designed to significantly improve the energy efficiency of user terminals in wireless local area networks. The new protocol, pointer controlled slot allocation and resynchronisation protocol (PCSAR), is based on the existing IEEE 802.11 point coordination function (PCF) standard. PCSAR reduces energy consumption by removing the need for power saving stations to remain awake and listen to the channel. Using OPNET, simulations were performed under symmetric channel loading conditions to compare the performance of PCSAR with the infrastructure power saving mode of IEEE 802.11, PCF-PS. The simulation results demonstrate a significant improvement in energy efficiency without significant reduction in performance when using PCSAR. For a wireless network consisting of an access point and 8 stations in power saving mode, the energy saving was up to 31% while using PCSAR instead of PCF-PS, depending upon frame error rate and load. The results also show that PCSAR offers significantly reduced uplink access delay over PCF-PS while modestly improving uplink throughput.

SrTiO3(001)(2x1) reconstructions: First-principles calculations of surface energy and atomic structure compared with scanning tunneling microscopy images

(1x1) and (2x1) reconstructions of the (001) SrTiO3 surface were studied using the first-principles full-potential linear muffin-tin orbital method. Surface energies were calculated as a function of TiO2 chemical potential, oxygen partial pressure and temperature. The (1x1) unreconstructed surfaces were found to be energetically stable for many of the conditions considered. Under conditions of very low oxygen partial pressure the (2x1) Ti2O3 reconstruction [Martin R. Castell, Surf. Sci. 505, 1 (2002)] is stable. The question as to why STM images of the (1x1) surfaces have not been obtained was addressed by calculating charge densities for each surface. These suggest that the (2x1) reconstructions would be easier to image than the (1x1) surfaces. The possibility that the presence of oxygen vacancies would destabilise the (1x1) surfaces was also investigated. If the (1x1) surfaces are unstable then there exists the further possibility that the (2x1) DL-TiO2 reconstruction [Natasha Erdman Nature (London) 419, 55 (2002)] is stable in a TiO2-rich environment and for p(O2)>10(-18) atm.

Spectroscopic measurements of phase-resolved electron energy distribution functions in RF-excited discharges

The reliable measurement of the electron energy distribution function (EEDF) of plasmas is one of the most important subjects of plasma diagnostics, because this piece of information is the key to understand basic discharge mechanisms. Specific problems arise in the case of RF-excited plasmas, since the properties of electrons are subject to changes on a nanosecond time scale and show pronounced spatial anisotropy. We report on a novel spectroscopic method for phase- and space-resolved measurements of the electron energy distribution function of energetic (> 12 eV) electrons in RF discharges. These electrons dominate excitation and ionization processes and are therefore of particular interest. The technique is based on time-dependent measurements during the RF cycle of excited-state populations of rare gases admixed in small fractions. These measurements yield � in combination with an analytical model � detailed information on the excitation processes. Phase-resolved optical emission spectroscopy allows us to overcome the difficulties connected with the very low densities (107�109 cm�3) and the transient character of the electrons in the sheath region. The EEDF of electrons accelerated in the sheath region can be described by a shifted Maxwellian with a drift velocity component in direction of the electric field. The method yields the high-energy tail of the EEDF on an absolute scale. The applicability of the method is demonstrated at a capacitively coupled RF discharge in hydrogen.

Damage to plasmid DNA induced by low energy carbon ions

The damage induced in supercoiled plasmid DNA molecules by 1-6 keV carbon ions has been investigated as a function of ion exposure, energy and charge state. The production of short linear fragments through multiple double strand breaks has been demonstrated and exponential exposure responses for each of the topoisomers have been found. The cross section for the loss of supercoiling was calculated to be (2.2 +/- 0.5) x 10(-14) cm(2) for 2 keVC(+) ions. For singly charged carbon ions, increased damage was observed with increasing ion energy. In the case of 2 keV doubly charged ions, the damage was greater than for singly charged ions of the same energy. These observations demonstrate that ion induced damage is a function of both the kinetic and potential energies of the ion.

Energy efficiency in gait, activity, participation, and health status in children with cerebral palsy

The aim of the study was to establish if a relationship exists between the energy efficiency of gait, and measures of activity limitation, participation restriction, and health status in a representative sample of children with cerebral palsy (CP). Secondary aims were to investigate potential differences between clinical subtypes and gross motor classification, and to explore other relationships between the measures under investigation. A longitudinal study of a representative sample of 184 children with ambulant CP was conducted (112 males, 72 females; 94 had unilateral spastic C P, 84 had bilateral spastic C P, and six had non-spastic forms; age range 4-17y; Gross Motor Function Classification System Level I, n=57; Level II, n=91; Level III, n=22; and Level IV, n=14); energy efficiency (oxygen cost) during gait, activity limitation, participation restriction, and health status were recorded. Energy efficiency during gait was shown to correlate significantly with activity limitations; no relationship between energy efficiency during gait was found with either participation restriction or health status. With the exception of psychosocial health, all other measures showed significant differences by clinical subtype and gross motor classification. The energy efficiency of walking is not reflective of participation restriction or health status. Thus, therapies leading to improved energy efficiency may not necessarily lead to improved participation or general health.