Stochastic linear multistep methods for the simulation of chemical kinetics

In this paper, we introduce the Stochastic Adams-Bashforth (SAB) and Stochastic Adams-Moulton (SAM) methods as an extension of the tau-leaping framework to past information. Using the theta-trapezoidal tau-leap method of weak order two as a starting procedure, we show that the k-step SAB method with k >= 3 is order three in the mean and correlation, while a predictor-corrector implementation of the SAM method is weak order three in the mean but only order one in the correlation. These convergence results have been derived analytically for linear problems and successfully tested numerically for both linear and non-linear systems. A series of additional examples have been implemented in order to demonstrate the efficacy of this approach.

Transition to an unsteady flow induced by a fin on the sidewall of a differentially heated air-filled square cavity and heat transfer

The transition from a steady to an unsteady flow induced by an adiabatic fin on the sidewall of a differentially heated air-filled cavity is numerically investigated. Numerical simulations have been performed over the range of Rayleigh numbers from Ra = 105–109. The temporal development and spatial structures of natural convection flows in the cavity with a fin are described. It has been demonstrated that the fin may induce the transition to an unsteady flow and the critical Rayleigh number for the occurrence of the transition is between 3.72 × 106 and 3.73 × 106. Furthermore, the peak frequencies of the oscillations triggered by different mechanisms are obtained through spectral analysis. It has been found that the flow rate through the cavity with a fin is larger than that without a fin under the unsteady flow, indicating that the fin may improve the unsteady flow in the cavity.

Directional synthesis of tin oxide@graphene nanocomposites via a one-step up-scalable wet-mechanochemical route for lithium ion batteries

Directional synthesis of SnO2@graphene nanocomposites via a one-step, low-cost, and up-scalable wetmechanochemical method is achieved using graphene oxide and SnCl2 as precursors. The graphene oxides are reduced to graphene while the SnCl2 is oxidized to SnO2 nanoparticles that are in situ anchored onto the graphene sheets evenly and densely, resulting in uniform SnO2@graphene nanocomposites. The prepared nanocomposites possess excellent electrochemical performance and outstanding cycling in Li-ion batteries.

Efficient simulation of stochastic chemical kinetics with the Stochastic Bulirsch-Stoer extrapolation method

Background Biochemical systems with relatively low numbers of components must be simulated stochastically in order to capture their inherent noise. Although there has recently been considerable work on discrete stochastic solvers, there is still a need for numerical methods that are both fast and accurate. The Bulirsch-Stoer method is an established method for solving ordinary differential equations that possesses both of these qualities. Results In this paper, we present the Stochastic Bulirsch-Stoer method, a new numerical method for simulating discrete chemical reaction systems, inspired by its deterministic counterpart. It is able to achieve an excellent efficiency due to the fact that it is based on an approach with high deterministic order, allowing for larger stepsizes and leading to fast simulations. We compare it to the Euler τ-leap, as well as two more recent τ-leap methods, on a number of example problems, and find that as well as being very accurate, our method is the most robust, in terms of efficiency, of all the methods considered in this paper. The problems it is most suited for are those with increased populations that would be too slow to simulate using Gillespie’s stochastic simulation algorithm. For such problems, it is likely to achieve higher weak order in the moments. Conclusions The Stochastic Bulirsch-Stoer method is a novel stochastic solver that can be used for fast and accurate simulations. Crucially, compared to other similar methods, it better retains its high accuracy when the timesteps are increased. Thus the Stochastic Bulirsch-Stoer method is both computationally efficient and robust. These are key properties for any stochastic numerical method, as they must typically run many thousands of simulations.

Pre-lithiation of onion-like carbon/MoS2 nano-urchin anodes for high-performance rechargeable lithium ion batteries

Hybrid urchin-like nanostructures composed of a spherical onion-like carbon (OLC) core and MoS2 nanoleaves were synthesized by a simple solvothermal method followed by thermal annealing treatment. Compared to commercial MoS2 powder, MoS2/OLC nanocomposites exhibit enhanced electrochemical performance as anode materials of lithium-ion batteries (LIBs) with a specific capacity of 853 mA h g−1 at a current density of 50 mA g−1 after 60 cycles, and a moderate initial coulombic efficiency of 71.1%. Furthermore, a simple pre-lithiation method based on direct contact of lithium foil with MoS2/OLC nano-urchins was used to achieve a very high coulombic efficiency of 97.6% in the first discharge/charge cycle, which is at least 26% higher compared to that of pristine MoS2/OLC nano-urchins. This pre-lithiation method can be generalized to develop other carbon-metal sulfide nanohybrids for LIB anode materials. These results may open up a new avenue for the development of the next-generation high-performance LIBs.

SnO2 decorated graphene nanocomposite anode materials prepared via an up-scalable wet-mechanochemical process for sodium ion batteries

A facile and up-scalable wet-mechanochemical process is designed for fabricating ultra-fine SnO2 nanoparticles anchored on graphene networks for use as anode materials for sodium ion batteries. A hierarchical structure of the SnO2@graphene composite is obtained from the process. The resultant rechargeable SIBs achieved high rate capability and good cycling stability.

Multifunctional three-dimensional T-junction graphene micro-wells: Energy-efficient, plasma-enabled growth and instant water-based transfer for flexible device applications

The “third-generation” 3D graphene structures, T-junction graphene micro-wells (T-GMWs) are produced on cheap polycrystalline Cu foils in a single-step, low-temperature (270 °C), energy-efficient, and environment-friendly dry plasma-enabled process. T-GMWs comprise vertical graphene (VG) petal-like sheets that seemlessly integrate with each other and the underlying horizontal graphene sheets by forming T-junctions. The microwells have the pico-to-femto-liter storage capacity and precipitate compartmentalized PBS crystals. The T-GMW films are transferred from the Cu substrates, without damage to the both, in de-ionized or tap water, at room temperature, and without commonly used sacrificial materials or hazardous chemicals. The Cu substrates are then re-used to produce similar-quality T-GMWs after a simple plasma conditioning. The isolated T-GMW films are transferred to diverse substrates and devices and show remarkable recovery of their electrical, optical, and hazardous NO2 gas sensing properties upon repeated bending (down to 1 mm radius) and release of flexible trasparent display plastic substrates. The plasma-enabled mechanism of T-GMW isolation in water is proposed and supported by the Cu plasma surface modification analysis. Our GMWs are suitable for various optoelectronic, sesning, energy, and biomedical applications while the growth approach is potentially scalable for future pilot-scale industrial production.

UWB channel measurement and data transfer analysis for multiuser Infostation applications

In this paper, the results of the time dispersion parameters obtained from a set of channel measurements conducted in various environments that are typical of multiuser Infostation application scenarios are presented. The measurement procedure takes into account the practical scenarios typical of the positions and movements of the users in the particular Infostation network. To provide one with the knowledge of how much data can be downloaded by users over a given time and mobile speed, data transfer analysis for multiband orthogonal frequency division multiplexing (MB-OFDM) is presented. As expected, the rough estimate of simultaneous data transfer in a multiuser Infostation scenario indicates dependency of the percentage of download on the data size, number and speed of the users, and the elapse time.

Catalyst engineering for lithium ion batteries: The catalytic role of Ge in enhancing the electrochemical performance of SnO2(GeO2)0.13/G anodes

The catalytic role of germanium (Ge) was investigated to improve the electrochemical performance of tin dioxide grown on graphene (SnO(2)/G) nanocomposites as an anode material of lithium ion batteries (LIBs). Germanium dioxide (GeO(20) and SnO(2) nanoparticles (<10 nm) were uniformly anchored on the graphene sheets via a simple single-step hydrothermal method. The synthesized SnO(2)(GeO(2))0.13/G nanocomposites can deliver a capacity of 1200 mA h g(-1) at a current density of 100 mA g(-1), which is much higher than the traditional theoretical specific capacity of such nanocomposites (∼ 702 mA h g(-1)). More importantly, the SnO(2)(GeO(2))0.13/G nanocomposites exhibited an improved rate, large current capability (885 mA h g(-1) at a discharge current of 2000 mA g(-1)) and excellent long cycling stability (almost 100% retention after 600 cycles). The enhanced electrochemical performance was attributed to the catalytic effect of Ge, which enabled the reversible reaction of metals (Sn and Ge) to metals oxide (SnO(2) and GeO(2)) during the charge/discharge processes. Our demonstrated approach towards nanocomposite catalyst engineering opens new avenues for next-generation high-performance rechargeable Li-ion batteries anode materials.

The significance of controlled conditions in lentiviral vector titration and in the use of multiplicity of infection (MOI) for predicting gene transfer events

Background: Although lentiviral vectors have been widely used for in vitro and in vivo gene therapy researches, there have been few studies systematically examining various conditions that may affect the determination of the number of viable vector particles in a vector preparation and the use of Multiplicity of Infection (MOI) as a parameter for the prediction of gene transfer events. Methods: Lentiviral vectors encoding a marker gene were packaged and supernatants concentrated. The number of viable vector particles was determined by in vitro transduction and fluorescent microscopy and FACs analyses. Various factors that may affect the transduction process, such as vector inoculum volume, target cell number and type, vector decay, variable vector - target cell contact and adsorption periods were studied. MOI between 0-32 was assessed on commonly used cell lines as well as a new cell line. Results: We demonstrated that the resulting values of lentiviral vector titre varied with changes of conditions in the transduction process, including inoculum volume of the vector, the type and number of target cells, vector stability and the length of period of the vector adsorption to target cells. Vector inoculum and the number of target cells determine the frequencies of gene transfer event, although not proportionally. Vector exposure time to target cells also influenced transduction results. Varying these parameters resulted in a greater than 50-fold differences in the vector titre from the same vector stock. Commonly used cell lines in vector titration were less sensitive to lentiviral vector-mediated gene transfer than a new cell line, FRL 19. Within 0-32 of MOI used transducing four different cell lines, the higher the MOI applied, the higher the efficiency of gene transfer obtained. Conclusion: Several variables in the transduction process affected in in vitro vector titration and resulted in vastly different values from the same vector stock, thus complicating the use of MOI for predicting gene transfer events. Commonly used target cell lines underestimated vector titre. However, within a certain range of MOI, it is possible that, if strictly controlled conditions are observed in the vector titration process, including the use of a sensitive cell line, such as FRL 19 for vector titration, lentivector-mediated gene transfer events could be predicted. © 2004 Zhang et al; licensee BioMed Central Ltd.

Unsteady flow and heat transfer between rotating coaxial disks

A study is made on the flow and heat transfer of a viscous fluid confined between two parallel disks. The disks are allowed to rotate with different time dependent angular velocities, and the upper disk is made to approach the lower one with a constant speed. Numerical solutions of the governing parabolic partial differential equations are obtained through a fourth-order accurate compact finite difference scheme. The normal forces and torques that the fluid exerts on the rotating surfaces are obtained at different nondimensional times for different values of the rate of squeezing and disk angular velocities. The temperature distribution and heat transfer are also investigated in the present analysis.

Nonlinear Enzyme Kinetics Can Lead to High Metabolic Flux Control Coefficients: Implications for the Evolution of Dominance

In a classic study, Kacser & Burns (1981, Genetics 97, 639-666) demonstrated that given certain plausible assumptions, the flux in a metabolic pathway was more or less indifferent to the activity of any of the enzymes in the pathway taken singly. It was inferred from this that the observed dominance of most wild-type alleles with respect to loss-of-function mutations did not require an adaptive, meaning selectionist, explanation. Cornish-Bowden (1987, J. theor. Biol. 125, 333-338) showed that the Kacser-Burns inference was not valid when substrate concentrations were large relative to the relevant Michaelis constants. We find that in a randomly constructed functional pathway, even when substrate levels are small, one can expect high values of control coefficients for metabolic flux in the presence of significant nonlinearities as exemplified by enzymes with Hill coefficients ranging from two to six, or by the existence of oscillatory loops. Under these conditions the flux can be quite sensitive to changes in enzyme activity as might be caused by inactivating one of the two alleles in a diploid. Therefore, the phenomenon of dominance cannot be a trivial ''default'' consequence of physiology but must be intimately linked to the manner in which metabolic networks have been moulded by natural selection.

Structural and lithium ion transport studies in sol-gel-prepared lithium silicophosphate glasses

Lithium silicophosphate glasses have been prepared by a sol-gel route over a wide range of compositions. Their structural and electrical properties have been investigated. Infrared spectroscopic studies show the presence of hydroxyl groups attached to Si and P. MAS NMR investigations provide evidence for the presence of different phosphatic units in the structure. The variations of de conductivities at 423 K and activation energies have been studied as a function of composition, and both exhibit an increasing trend with the ratio of nonbridging oxygen to bridging oxygen in the structure. Ac conductivity behavior shows that the power law exponent, s, is temperature dependent and exhibits a minimum. Relaxation behavior has been examined in detail using an electrical modulus formalism, and modulus data were fitted to Kohlraush-William-Watts stretched exponential function. A structural model has been proposed and the unusual properties exhibited by this unique system of glasses have been rationalized using this model. Ion transport in these glasses appears to be confined to unidimensional conduits defined by modified phosphate chains and interspersed with unmodified silica units.

Synthesis and vesicle formation from hybrid bolaphile/amphiphile ion-pairs. Evidence of membrane property modulation by molecular design

Four new hybrid (bolaphile/amphiphile) ion-pairs were synthesized. Electron microscopy indicated that each of these forms bilayer membranes upon dispersion in aqueous media. Membrane properties have also been examined by differential scanning calorimetry, microcalorimetry, temperature-dependent fluorescence anisotropy measurements, and UV-vis spectroscopy. The T-m values for the vesicular 1, 2, 3, 4, and 5 were 38, 12, 85, 31.3, and 41.6 degrees C, respectively. Interestingly the T-m values for 1 and 3 were found to depend on their concentration. The entrapment of small solute and the release capability have also been examined to demonstrate that these bilayers form enclosed vesicles. X-ray diffraction of the cast films has been performed to understand the nature and the thickness of these membrane organizations. The membrane widths ranged from 33 to 47 Angstrom. Finally, the above observations have been analyzed in light of the results obtained from molecular modeling studies. Thus we have demonstrated that membrane properties can be modulated by simple structural changes at the amphiphile level. It was shown that by judicious incorporation of central, isomeric, disubstituted aromatic units as structural anchors into different bolaphiles, one can modulate the properties of the resulting vesicles.

Void Geometry Driven Spin Crossover in Zeolite-Encapsulated Cobalt Tris(bipyridyl) Complex Ion

The cobalt(II) tris(bipyridyl) complex ion encapsulated in zeolite-Y supercages exhibits a thermally driven interconversion between a low-spin and a high-spin state-a phenomenon not observed for this ion either in solid state or in solution. From a comparative study of the magnetism and optical spectroscopy of the encapsulated and unencapsulated complex ion, supported by molecular modeling, such spin behavior is shown to be intramolecular in origin. In the unencapsulated or free state, the [Co(bipy)(3)](2+) ion exhibits a marked trigonal prismatic distortion, but on encapsulation, the topology of the supercage forces it to adopt a near-octahedral geometry. An analysis using the angular overlap ligand field model with spectroscopically derived parameters shows that the geometry does indeed give rise to a low-spin ground state, and suggests a possible scenario for the spin state interconversion.