Validation of a lattice Boltzmann model for gas-solid reactions with experiments

A lattice Boltzmann method is used to model gas-solid reactions where the composition of both the gas and solid phase changes with time, while the boundary between phases remains fixed. The flow of the bulk gas phase is treated using a multiple relaxation time MRT D3Q19 model; the dilute reactant is treated as a passive scalar using a single relaxation time BGK D3Q7 model with distinct inter- and intraparticle diffusivities. A first-order reaction is incorporated by modifying the method of Sullivan et al. [13] to include the conversion of a solid reactant. The detailed computational model is able to capture the multiscale physics encountered in reactor systems. Specifically, the model reproduced steady state analytical solutions for the reaction of a porous catalyst sphere (pore scale) and empirical solutions for mass transfer to the surface of a sphere at Re=10 (particle scale). Excellent quantitative agreement between the model and experiments for the transient reduction of a single, porous sphere of Fe 2O 3 to Fe 3O 4 in CO at 1023K and 10 5Pa is demonstrated. Model solutions for the reduction of a packed bed of Fe 2O 3 (reactor scale) at identical conditions approached those of experiments after 25 s, but required prohibitively long processor times. The presented lattice Boltzmann model resolved successfully mass transport at the pore, particle and reactor scales and highlights the relevance of LB methods for modelling convection, diffusion and reaction physics. © 2012 Elsevier Inc.

A simple bioconjugate attachment protocol for use in single molecule force spectroscopy experiments based on mixed self-assembled monolayers.

Single molecule force spectroscopy is a technique that can be used to probe the interaction force between individual biomolecular species. We focus our attention on the tip and sample coupling chemistry, which is crucial to these experiments. We utilised a novel approach of mixed self-assembled monolayers of alkanethiols in conjunction with a heterobifunctional crosslinker. The effectiveness of the protocol is demonstrated by probing the biotin-avidin interaction. We measured unbinding forces comparable to previously reported values measured at similar loading rates. Specificity tests also demonstrated a significant decrease in recognition after blocking with free avidin.

Topologically invariant reaction coordinates for simulating multistate chemical reactions.

Evaluating free energy profiles of chemical reactions in complex environments such as solvents and enzymes requires extensive sampling, which is usually performed by potential of mean force (PMF) techniques. The reliability of the sampling depends not only on the applied PMF method but also the reaction coordinate space within the dynamics is biased. In contrast to simple geometrical collective variables that depend only on the positions of the atomic coordinates of the reactants, the E(gap) reaction coordinate (the energy difference obtained by evaluating a suitable force field using reactant and product state topologies) has the unique property that it is able to take environmental effects into account leading to better convergence, a more faithful description of the transition state ensemble and therefore more accurate free energy profiles. However, E(gap) requires predefined topologies and is therefore inapplicable for multistate reactions, in which the barrier between the chemically equivalent topologies is comparable to the reaction activation barrier, because undesired "side reactions" occur. In this article, we introduce a new energy-based collective variable by generalizing the E(gap) reaction coordinate such that it becomes invariant to equivalent topologies and show that it yields more well behaved free energy profiles than simpler geometrical reaction coordinates.

Tests of an adaptive QM/MM calculation on free energy profiles of chemical reactions in solution.

We present reaction free energy calculations using the adaptive buffered force mixing quantum mechanics/molecular mechanics (bf-QM/MM) method. The bf-QM/MM method combines nonadaptive electrostatic embedding QM/MM calculations with extended and reduced QM regions to calculate accurate forces on all atoms, which can be used in free energy calculation methods that require only the forces and not the energy. We calculate the free energy profiles of two reactions in aqueous solution: the nucleophilic substitution reaction of methyl chloride with a chloride anion and the deprotonation reaction of the tyrosine side chain. We validate the bf-QM/MM method against a full QM simulation, and show that it correctly reproduces both geometrical properties and free energy profiles of the QM model, while the electrostatic embedding QM/MM method using a static QM region comprising only the solute is unable to do so. The bf-QM/MM method is not explicitly dependent on the details of the QM and MM methods, so long as it is possible to compute QM forces in a small region and MM forces in the rest of the system, as in a conventional QM/MM calculation. It is simple, with only a few parameters needed to control the QM calculation sizes, and allows (but does not require) a varying and adapting QM region which is necessary for simulating solutions.

Single event molecular signalling for estimation and control

Cell biology is characterised by low molecule numbers and coupled stochastic chemical reactions with intrinsic noise permeating and dominating the interactions between molecules. Recent work [9] has shown that in such environments there are hard limits on the accuracy with which molecular populations can be controlled and estimated. These limits are predicated on a continuous diffusion approximation of the target molecule (although the remainder of the system is non-linear and discrete). The principal result of [9] assumes that the birth rate of the signalling species is linearly dependent on the target molecule population size. In this paper, we investigate the situation when the entire system is kept discrete, and arbitrary non-linear coupling is allowed between the target molecule and downstream signalling molecules. In this case it is possible, by relying solely on the event triggered nature of control and signalling reactions, to define non-linear reaction rate modulation schemes that achieve improved performance in certain parameter regimes. These schemes would not appear to be biologically relevant, raising the question of what are an appropriate set of assumptions for obtaining biologically meaningful results. © 2013 EUCA.

Kinetic modelling of growth and storage molecule production in microalgae under mixotrophic and autotrophic conditions.

In order to improve algal biofuel production on a commercial-scale, an understanding of algal growth and fuel molecule accumulation is essential. A mathematical model is presented that describes biomass growth and storage molecule (TAG lipid and starch) accumulation in the freshwater microalga Chlorella vulgaris, under mixotrophic and autotrophic conditions. Biomass growth was formulated based on the Droop model, while the storage molecule production was calculated based on the carbon balance within the algal cells incorporating carbon fixation via photosynthesis, organic carbon uptake and functional biomass growth. The model was validated with experimental growth data of C. vulgaris and was found to fit the data well. Sensitivity analysis showed that the model performance was highly sensitive to variations in parameters associated with nutrient factors, photosynthesis and light intensity. The maximum productivity and biomass concentration were achieved under mixotrophic nitrogen sufficient conditions, while the maximum storage content was obtained under mixotrophic nitrogen deficient conditions.

Biofunctionalization of PET/SiO2 surfaces for single molecule experiments and medical applications

PET/SiO2 layers were chemically modified to maintain immobilization of functional single molecules. GFP molecules provide an ideal system due to their stability and intrinsic fluorescence. GFP in vivo biotinylated within its NH2-terminal region and attached on the substrate via the biotinstreptavidin bond was further investigated with confocal microscopy, atomic force microscopy (AFM) and spectroscopic ellipsometry (SE). AFM revealed monolayered donut-like structures representing assemblies of biotinstreptavidinbiotinGFP immobilized onto PET/SiO2 surfaces via mPEG. In particular, regions with an approximate height of 12 nm, which approaches the molecular dimensions of the above complex given by molecular modeling, could be detected. The dimensions of the donut-like structures suggest a close-to-each-other positioning of the GFP molecules - which, however, retain their functionality, as evidenced by confocal microscopy. © 2011 World Scientific Publishing Company.

Kinetic modelling of growth and storage molecule production in microalgae under mixotrophic and autotrophic conditions

In order to improve algal biofuel production on a commercial-scale, an understanding of algal growth and fuel molecule accumulation is essential. A mathematical model is presented that describes biomass growth and storage molecule (TAG lipid and starch) accumulation in the freshwater microalga Chlorella vulgaris, under mixotrophic and autotrophic conditions. Biomass growth was formulated based on the Droop model, while the storage molecule production was calculated based on the carbon balance within the algal cells incorporating carbon fixation via photosynthesis, organic carbon uptake and functional biomass growth. The model was validated with experimental growth data of C. vulgaris and was found to fit the data well. Sensitivity analysis showed that the model performance was highly sensitive to variations in parameters associated with nutrient factors, photosynthesis and light intensity. The maximum productivity and biomass concentration were achieved under mixotrophic nitrogen sufficient conditions, while the maximum storage content was obtained under mixotrophic nitrogen deficient conditions. © 2014 Elsevier Ltd.

Genomic organization, gene duplication, and expression analysis of interleukin-1 beta in channel catfish (Ictalurus punctatus)

Interleukin-1 beta (IL-1 beta) is one of the pivotal early response pro-inflammatory cytokines that enables organisms to respond to infection and induces a cascade of reactions leading to inflammation. In spite of its importance and two decades of studies in the mammalian species, genes encoding IL-1 beta were not identified from non-mammalian species until recently. Recent research, particularly with genomic approaches, has led to sequencing of IL-1 beta from many species. Clinical studies also Suggested IL-1 beta as an immunoreagulatory molecule potentially useful for enhancing vaccination. However, no IL-1 beta genes have been identified from channel catfish, the primary aquaculture species from the United States. In this study, we identified two distinct cDNAs encoding catfish IL-1 beta. Their encoding genes were identified, sequenced, and characterized. The catfish IL-1 beta genes were assigned to bacterial artificial chromosome (BAC) clones. Genomic studies indicated that the IL-1 beta genes were tandemly duplicated on the same chromosome. Phylogenetic analysis of various IL-1 beta genes indicated the possibility of recent species-specific gene duplications in channel catfish, and perhaps also in swine and carp. Expression analysis indicated that both IL-1 beta genes were expressed, but exhibited distinct expression profiles in various catfish tissues, and after bacterial infection with Edwardsiella ictaluri. (c) 2005 Elsevier Ltd. All rights reserved.

Role of interactions in electronic structure of a two-electron quantum dot molecule

We have studied a two-electron quantum dot molecule in a magnetic field. The electron interaction is treated accurately by the direct diagonalization of the Hamiltonian matrix. We calculate two lowest energy levels of the two-electron quantum dot molecule in a magnetic field. Our results show that the electron interactions are significant, as they can change the total spin of the two-electron ground state of the system by adjusting the magnetic field between S = 0 and S = 1. The energy difference DeltaE between the lowest S = 0 and S = 1 states is shown as a function of the axial magnetic field. We found that the energy difference between the lowest S = 0 and S = 1 states in the strong-B S = 0 state varies linearly. Our results provide a possible realization for a qubit to be fabricated by current growth techniques.

Electronic structures of N quantum dot molecule

The electronic structures of N quantum dot molecules (QDMs) are investigated theoretically in the framework of effective-mass envelope function theory. The electron and hole energy levels are calculated. In the calculations, the effects of finite offset and valence-band mixing are taken into account. The theoretical method can be used to calculate the electronic structures of any QDM. The results show that (1) electronic energy levels decrease monotonically and the energy difference between the N QDMs decreases as the quantum dot (QD) radius increases; (2) the electron energy level is lower and quantum confinement is smaller for the larger N QDM; (3) the hole ground state energy level is lower for the one dot QDM than N (greater 1) QDMs if the QD radius is larger than about 5 nm due to the valence-band mixing. The results are useful for the application of the N QDM to photoelectric devices.