Mesophases in Nearly 2D Room-Temperature Ionic Liquids

Computer simulations of (i) a [C(12)mim][Tf2N] film of nanometric thickness squeezed at kbar pressure by a piecewise parabolic confining potential reveal a mesoscopic in-plane density and composition modulation reminiscent of mesophases seen in 3D samples of the same room-temperature ionic liquid (RTIL). Near 2D confinement, enforced by a high normal load, as well as relatively long aliphatic chains are strictly required for the mesophase formation, as confirmed by computations for two related systems made of (ii) the same [C(12)mim][Tf2N] adsorbed at a neutral solid surface and (iii) a shorter-chain RTIL ([C(4)mim][Tf2N]) trapped in the potential well of part i. No in-plane modulation is seen for ii and iii. In case ii, the optimal arrangement of charge and neutral tails is achieved by layering parallel to the surface, while, in case iii, weaker dispersion and packing interactions are unable to bring aliphatic tails together into mesoscopic islands, against overwhelming entropy and Coulomb forces. The onset of in-plane mesophases could greatly affect the properties of long-chain RTILs used as lubricants.

Phase behavior of mixtures of colloidal platelets and nonadsorbing polymers

The phase behavior of a model system of colloidal platelets and nonadsorbing polymers is studied using computer simulations and perturbation theory. The equation of state for the pure platelet reference system is obtained by Monte Carlo simulations, and the free volume fraction accessible to polymers is measured by a trial insertion method. The free volume fraction is also calculated using scaled particle theory. Subsequently, the phase diagram for platelet-polymer mixtures is calculated. For a platelet aspect ratio L/D=0.1 and a polymer to platelet size ratio d/D>0.2, we observe coexistence between two isotropic phases with different densities. For smaller polymers d/D

A Monte Carlo model of DNA double-strand break clustering and rejoining kinetics for the analysis of pulsed-field gel electrophoresis data

In studies of radiation-induced DNA fragmentation and repair, analytical models may provide rapid and easy-to-use methods to test simple hypotheses regarding the breakage and rejoining mechanisms involved. The random breakage model, according to which lesions are distributed uniformly and independently of each other along the DNA, has been the model most used to describe spatial distribution of radiation-induced DNA damage. Recently several mechanistic approaches have been proposed that model clustered damage to DNA. In general, such approaches focus on the study of initial radiation-induced DNA damage and repair, without considering the effects of additional (unwanted and unavoidable) fragmentation that may take place during the experimental procedures. While most approaches, including measurement of total DNA mass below a specified value, allow for the occurrence of background experimental damage by means of simple subtractive procedures, a more detailed analysis of DNA fragmentation necessitates a more accurate treatment. We have developed a new, relatively simple model of DNA breakage and the resulting rejoining kinetics of broken fragments. Initial radiation-induced DNA damage is simulated using a clustered breakage approach, with three free parameters: the number of independently located clusters, each containing several DNA double-strand breaks (DSBs), the average number of DSBs within a cluster (multiplicity of the cluster), and the maximum allowed radius within which DSBs belonging to the same cluster are distributed. Random breakage is simulated as a special case of the DSB clustering procedure. When the model is applied to the analysis of DNA fragmentation as measured with pulsed-field gel electrophoresis (PFGE), the hypothesis that DSBs in proximity rejoin at a different rate from that of sparse isolated breaks can be tested, since the kinetics of rejoining of fragments of varying size may be followed by means of computer simulations. The problem of how to account for background damage from experimental handling is also carefully considered. We have shown that the conventional procedure of subtracting the background damage from the experimental data may lead to erroneous conclusions during the analysis of both initial fragmentation and DSB rejoining. Despite its relative simplicity, the method presented allows both the quantitative and qualitative description of radiation-induced DNA fragmentation and subsequent rejoining of double-stranded DNA fragments. (C) 2004 by Radiation Research Society.

Dynamical characteristics of solitary waves, shocks and envelope modes in kappa-distributed non-thermal plasmas: an overview

Space plasmas provide abundant evidence of highly energetic particle population, resulting in a long-tailed non-Maxwellian distribution. Furthermore, the first stages in the evolution of plasmas produced during laser-matter interaction are dominated by nonthermal electrons, as confirmed by experimental observation and computer simulations. This phenomenon is efficiently modelled via a kappa-type distribution. We present an overview, from first principles, of the effect of superthermality on the characteristics of electrostatic plasma waves. We rely on a fluid model for ion-acoustic excitations, employing a kappa distribution function to model excess superthermality of the electron distribution. Focusing on nonlinear excitations (solitons), in the form of solitary waves (pulses), shocks and envelope solitons, and employing standard methodological tools of nonlinear plasmadynamical analysis, we discuss the role of excess superthermality in their propagation dynamics (existence laws, stability profile), geometric characteristics and stability. Numerical simulations are employed to confirm theoretical predictions, namely in terms of the stability of electrostatic pulses, as well as the modulational stability profile of bright- and dark-type envelope solitons.

The phase diagram of water at negative pressures: Virtual ices

The phase diagram of water at negative pressures as obtained from computer simulations for two models of water, TIP4P/2005 and TIP5P is presented. Several solid structures with lower densities than ice Ih, so-called virtual ices, were considered as possible candidates to occupy the negative pressure region of the phase diagram of water. In particular the empty hydrate structures sI, sII, and sH and another, recently proposed, low-density ice structure. The relative stabilities of these structures at 0 K was determined using empirical water potentials and density functional theory calculations. By performing free energy calculations and Gibbs-Duhem integration the phase diagram of TIP4P/2005 was determined at negative pressures. The empty hydrates sII and sH appear to be the stable solid phases of water at negative pressures. The phase boundary between ice Ih and sII clathrate occurs at moderate negative pressures, while at large negative pressures sH becomes the most stable phase. This behavior is in reasonable agreement with what is observed in density functional theory calculations.

Designing and understanding permanent microporosity in liquids

Standard microporous materials are typically crystalline solids that exhibit a regular array of cavities of uniform size and shape. Packing and directional bonding between molecular building blocks give rise to interstitial pores that confer size and shape-specific sorption properties to the material. In the liquid state interstitial cavities are transient. However, permanent and intrinsic "pores'' can potentially be built into the structure of the molecules that constitute the liquid. With the aid of computer simulations we have designed, synthesised and characterised a series of liquids composed of hollow cage-like molecules, which are functionalised with hydrocarbon chains to make them liquid at accessible temperatures. Experiments and simulations demonstrate that chain length and size of terminal chain substituents can be used to tune, within certain margins, the permanence of intramolecular cavities in such neat liquids. Simulations identify a candidate "porous liquid'' in which 30% of the cages remain empty in the liquid state. Absorbed methane molecules selectively occupy these empty cavities.

Photo-ionized neon plasmas induced by radiation pulses of a laser-plasma EUV source and a free electron laser FLASH

In this work, a laser-produced plasma extreme ultraviolet source and a free electron laser were used to create Ne photo-ionized plasmas. In both cases, a radiation beam was focused onto a gas stream injected into a vacuum chamber synchronously with the radiation pulse. Extreme ultraviolet radiation from the plasma spanned a wide spectral range with pronounced maximum centered at lambda = 11 +/- 1 nm while the free electron laser pulses were emitted at a wavelength of 32 nm. The power density of the focused plasma radiation was approximately 2 x 10(7) W/cm(2) and was seven orders of magnitude lower compared with the focused free electron laser beam. Radiation fluences in both experimental conditions were comparable. Despite quite different spectral characteristics and extremely different power densities, emission spectra of both photo-ionized plasmas consist of the same spectral lines within a wavelength range of 20 to 50 nm, however, with different relative intensities of the corresponding lines. The dominating spectral lines originated from singly charged ions (Ne II); however, Ne III lines were also detected. Additionally, computer simulations of the emission spectra, obtained for photo-ionized plasmas, driven by the plasma extreme ultraviolet source, were performed. The corresponding measured and calculated spectra are presented. An electron temperature and ionic composition were estimated. Differences between the experimental spectra, obtained for both irradiation conditions, were analyzed. The differences were attributed mainly to different energies of driving photons.

Dual-Hop Cognitive Amplify-and-Forward Relaying Networks over η− μ Fading Channels

This paper presents a thorough performance analysis of dual-hop cognitive amplify-and-forward (AF) relaying networks under spectrum-sharing mechanism over independent non-identically distributed (i.n.i.d.) 􀀀 fading channels. In order to guarantee the quality-of-service (QoS) of primary networks, both maximum tolerable peak interference power Q at the primary users (PUs) and maximum allowable transmit power P at secondary users (SUs) are considered to constrain transmit power at the cognitive transmitters. For integer-valued fading parameters, a closed-form lower bound for the outage probability (OP) of the considered networks is obtained. Moreover, assuming arbitrary-valued fading parameters, the lower bound in integral form for the OP is derived. In order to obtain further insights on the OP performance, asymptotic expressions for the OP at high SNRs are derived, from which the diversity/coding gains and the diversity-multiplexing gain tradeoff (DMT) of the secondary network can be readily deduced. It is shown that the diversity gain and also the DMT are solely determined by the fading parameters of the secondary network whereas the primary network only affects the coding gain. The derived results include several others available in previously published works as special cases, such as those for Nakagami-m fading channels. In addition, performance evaluation results have been obtained by Monte Carlo computer simulations which have verified the accuracy of the theoretical analysis.

Internal cooling in rotational molding-A review

Rotational molding suffers from a relatively long cycle time, which hampers more widespread growth of the process. During each cycle, both the polymer and mold must be heated from room temperature to above polymer melting temperature and subsequently cooled to room temperature. The cooling time in this process is relatively long due to the poor thermal conductivity of plastics. Although rapid external cooling is possible, internal cooling rates are the major limitation. This causes the process to be uneconomical for large production runs of small parts. Various researchers have strived to minimize cycle times by applying various internal cooling procedures. This article presents a review of these methods, including computer simulations and practical investigations published to date. The effects of cooling rate on the morphology, shrinkage, warpage, and impact properties of rotationally molded polyolefins are also highlighted. In general, rapid and symmetrical cooling across the mold results in smaller spherulite size, increased mechanical properties and less potential warpage or distortion in moldings. POLYM. ENG. SCI., 2011. ©2011 Society of Plastics Engineers.

Surface alloying and mixing at the Mn/Fe(001) interface: Real-time photoelectron spectroscopy and modified embedded atom simulations

Structural and magnetic properties of thin Mn films on the Fe(001) surface have been investigated by a combination of photoelectron spectroscopy and computer simulation in the temperature range 300 Kless than or equal toTless than or equal to750 K. Room-temperature as deposited Mn overlayers are found to be ferromagnetic up to 2.5-monolayer (ML) coverage, with a magnetic moment parallel to that of the iron substrate. The Mn atomic moment decreases with increasing coverage, and thicker samples (4-ML and 4.5-ML coverage) are antiferromagnetic. Photoemission measurements performed while the system temperature is rising at constant rate (dT/dtsimilar to0.5 K/s) detect the first signs of Mn-Fe interdiffusion at T=450 K, and reveal a broad temperature range (610 Kless than or equal toTless than or equal to680 K) in which the interface appears to be stable. Interdiffusion resumes at Tgreater than or equal to680 K. Molecular dynamics and Monte Carlo simulations allow us to attribute the stability plateau at 610 Kless than or equal toTless than or equal to680 K to the formation of a single-layer MnFe surface alloy with a 2x2 unit cell and a checkerboard distribution of Mn and Fe atoms. X-ray-absorption spectroscopy and analysis of the dichroic signal show that the alloy has a ferromagnetic spin structure, collinear with that of the substrate. The magnetic moments of Mn and Fe atoms in the alloy are estimated to be 0.8mu(B) and 1.1mu(B), respectively.

Computer simulating a trial of a load-bearing implant: example of an intramedullary prosthesis

Computational modelling is becoming ever more important for obtaining regulatory approval for new medical devices. An accepted approach is to infer performance in a population from an analysis conducted for an idealised or ‘average’ patient; we present here a method for predicting the performance of an orthopaedic implant when released into a population—effectively simulating a clinical trial. Specifically we hypothesise that an analysis based on a method for predicting the performance in a population will lead to different conclusions than an analysis based on an idealised or ‘average’ patient. To test this hypothesis we use a finite element model of an intramedullary implant in a bone whose size and remodelling activity is different for each individual in the population. We compare the performance of a low Young’s modulus implant (View the MathML source) to one with a higher Young’s modulus (200 GPa). Cyclic loading is applied and failure is assumed when the migration of the implant relative to the bone exceeds a threshold magnitude. The analysis for an idealised of ‘average’ patient predicts that the lower modulus device survives longer whereas the analysis simulating a clinical trial predicts no statistically-significant tendency (p=0.77) for the low modulus device to perform better. It is concluded that population-based simulations of implant performance–simulating a clinical trial–present a very valuable opportunity for more realistic computational pre-clinical testing of medical devices.

The Influence of Friction Stir Welding Process Idealization on Residual Stress and Distortion Predictions for Future Airframe Assembly Simulations

The ability to accurately predict residual stresses and resultant distortions is a key product from process assembly simulations. Assembly processes necessarily consider large structural components potentially making simulations computationally expensive. The objective herein is to develop greater understanding of the influence of friction stir welding process idealization on the prediction of residual stress and distortion and thus determine the minimum required modeling fidelity for future airframe assembly simulations. The combined computational and experimental results highlight the importance of accurately representing the welding forging force and process speed. In addition, the results emphasize that increased CPU simulation times are associated with representing the tool torque, while there is potentially only local increase in prediction fidelity.