A numerical approach to modeling the catalytic voltammetry of surface-confined redox enzymes

A finite difference method for simulating voltammograms of electrochemically driven enzyme catalysis is presented. The method enables any enzyme mechanism to be simulated. The finite difference equations can be represented as a matrix equation containing a nonlinear sparse matrix. This equation has been solved using the software package Mathematica. Our focus is on the use of cyclic voltammetry since this is the most commonly employed electrochemical method used to elucidate mechanisms. The use of cyclic voltammetry to obtain data from systems obeying Michaelis-Menten kinetics is discussed, and we then verify our observations on the Michaelis-Menten system using the finite difference simulation. Finally, we demonstrate how the method can be used to obtain mechanistic information on a real redox enzyme system, the complex bacterial molybdoenzyme xanthine dehydrogenase.

Bitwise bell-inequality violations for an entangled state involving 2N ions

Following on from previous work [J.-A. Larsson, Phys. Rev. A 67, 022108 (2003)], Bell inequalities based on correlations between binary digits are considered for a particular entangled state involving 2N trapped ions. These inequalities involve applying displacement operations to half of the ions and then measuring correlations between pairs of corresponding bits in the binary representations of the number of center-of-mass phonons of N particular ions. It is shown that the state violates the inequalities and thus displays nonclassical correlations. It is also demonstrated that it violates a Bell inequality when the displacements are replaced by squeezing operations.

Fock-state dynamics in Raman photoassociation of Bose-Einstein condensates

By stochastic modeling of the process of Raman photoassociation of Bose-Einstein condensates, we show that, the farther the initial quantum state is from a coherent state, the farther the one-dimensional predictions are from those of the commonly used zero-dimensional approach. We compare the dynamics of condensates, initially in different quantum states, finding that, even when the quantum prediction for an initial coherent state is relatively close to the Gross-Pitaevskii prediction, an initial Fock state gives qualitatively different predictions. We also show that this difference is not present in a single-mode type of model, but that the quantum statistics assume a more important role as the dimensionality of the model is increased. This contrasting behavior in different dimensions, well known with critical phenomena in statistical mechanics, makes itself plainly visible here in a mesoscopic system and is a strong demonstration of the need to consider physically realistic models of interacting condensates.

Effect of early marine diagenesis on coral reconstructions of surface-ocean 13C/12C and carbonate saturation state

Recent research suggests that future decreases in the carbonate saturation state of surface seawater associated with the projected build-up of atmospheric CO2 could cause a global decline in coral reef-building capacity. Whether significant reductions in coral calcification are underway is a matter of considerable debate. Multicentury records of skeletal calcification extracted from massive corals have the potential to reconstruct the progressive effect of anthropogenic changes in carbonate saturation on coral reefs. However, early marine aragonite cements are commonly precipitated from pore waters in the basal portions of massive coral skeletons and, if undetected, could result in apparent nonlinear reductions in coral calcification toward the present. To address this issue, we present records of coral skeletal density, extension rate, calcification rate, δ13C, and δ18O for well preserved and diagenetically altered coral cores spanning ∼1830-1994 A.D. at Ningaloo Reef Marine Park, Western Australia. The record for the pristine coral shows no significant decrease in skeletal density or δ13C indicative of anthropogenic changes in carbonate saturation state or δ13C of surface seawater (oceanic Suess effect). In contrast, progressive addition of early marine inorganic aragonite toward the base of the altered coral produces an apparent ∼25% decrease in skeletal density toward the present, which misleadingly matches the nonlinear twentieth century decrease in coral calcification predicted by recent modeling and experimental studies. In addition, the diagenetic aragonite is enriched in 13C, relative to coral aragonite, resulting in a nonlinear decrease in δ13C toward the present that mimics the decrease in δ13C expected from the oceanic Suess effect. Taken together, these diagenetic changes in skeletal density and δ13C could be misinterpreted to reflect changes in surface-ocean carbonate saturation state driven by the twentieth century build-up of atmospheric CO2. Copyright 2004 by the American Geophysical Union.

Modeling the development of acquired clinical immunity to Plasmodium falciparum malaria

Individuals living in regions where malaria is endemic develop an acquired immunity to malaria which enables them to remain asymptomatic while still carrying parasites. Field studies indicate that cumulative exposure to a variety of diverse Plasmodium parasites is required for the transition from symptomatic to asymptomatic malaria. This study used a simulation model of the within-host dynamics of P. falciparum to investigate the development of acquired clinical immunity under different transmission conditions and levels of parasite diversity. Antibodies developed to P. falciparum erythrocyte membrane protein 1 (PfEMP1), a clonally variant molecule, were assumed to be a key human immunological response to P. falciparum infection, along with responses to clonally conserved but polymorphic antigens. The time to the development of clinical immunity was found to be proportional to parasite diversity and inversely proportional to transmission intensity. The effect of early termination of symptomatic infections by chemotherapy was investigated and found not to inhibit the host's ability to develop acquired immunity. However, the time required to achieve this state was approximately double that compared to when no treatment was administered. This study demonstrates that an immune response primarily targeted against PfEMP1 has the ability to reduce clinical symptoms of infections irrespective of whether treatment is administered, supporting its role in the development of acquired clinical immunity. The results also illustrate a novel use for simulation models of P. falciparum infections, investigation of the influence of intervention strategies on the development of naturally acquired clinical immunity.

Quantum-gate characterization in an extended Hilbert space

We describe an approach for characterizing the process performed by a quantum gate using quantum process tomography, by first modeling the gate in an extended Hilbert space, which includes nonqubit degrees of freedom. To prevent unphysical processes from being predicted, present quantum process tomography procedures incorporate mathematical constraints, which make no assumptions as to the actual physical nature of the system being described. By contrast, the procedure presented here assumes a particular class of physical processes, and enforces physicality by fitting the data to this model. This allows quantum process tomography to be performed using a smaller experimental data set, and produces parameters with a direct physical interpretation. The approach is demonstrated by example of mode matching in an all-optical controlled-NOT gate. The techniques described are general and could be applied to other optical circuits or quantum computing architectures.

Modeling diffraction and imaging of laser beams by the mode-expansion method

We present a new method of modeling imaging of laser beams in the presence of diffraction. Our method is based on the concept of first orthogonally expanding the resultant diffraction field (that would have otherwise been obtained by the laborious application of the Huygens diffraction principle) and then representing it by an effective multimodal laser beam with different beam parameters. We show not only that the process of obtaining the new beam parameters is straightforward but also that it permits a different interpretation of the diffraction-caused focal shift in laser beams. All of the criteria that we have used to determine the minimum number of higher-order modes needed to accurately represent the diffraction field show that the mode-expansion method is numerically efficient. Finally, the characteristics of the mode-expansion method are such that it allows modeling of a vast array of diffraction problems, regardless of the characteristics of the incident laser beam, the diffracting element, or the observation plane. (C) 2005 Optical Society of America.

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.

Quantum Information Processing with spin systems: from modeling to possible implementations

The physical implementation of quantum information processing is one of the major challenges of current research. In the last few years, several theoretical proposals and experimental demonstrations on a small number of qubits have been carried out, but a quantum computing architecture that is straightforwardly scalable, universal, and realizable with state-of-the-art technology is still lacking. In particular, a major ultimate objective is the construction of quantum simulators, yielding massively increased computational power in simulating quantum systems. Here we investigate promising routes towards the actual realization of a quantum computer, based on spin systems. The first one employs molecular nanomagnets with a doublet ground state to encode each qubit and exploits the wide chemical tunability of these systems to obtain the proper topology of inter-qubit interactions. Indeed, recent advances in coordination chemistry allow us to arrange these qubits in chains, with tailored interactions mediated by magnetic linkers. These act as switches of the effective qubit-qubit coupling, thus enabling the implementation of one- and two-qubit gates. Molecular qubits can be controlled either by uniform magnetic pulses, either by local electric fields. We introduce here two different schemes for quantum information processing with either global or local control of the inter-qubit interaction and demonstrate the high performance of these platforms by simulating the system time evolution with state-of-the-art parameters. The second architecture we propose is based on a hybrid spin-photon qubit encoding, which exploits the best characteristic of photons, whose mobility is exploited to efficiently establish long-range entanglement, and spin systems, which ensure long coherence times. The setup consists of spin ensembles coherently coupled to single photons within superconducting coplanar waveguide resonators. The tunability of the resonators frequency is exploited as the only manipulation tool to implement a universal set of quantum gates, by bringing the photons into/out of resonance with the spin transition. The time evolution of the system subject to the pulse sequence used to implement complex quantum algorithms has been simulated by numerically integrating the master equation for the system density matrix, thus including the harmful effects of decoherence. Finally a scheme to overcome the leakage of information due to inhomogeneous broadening of the spin ensemble is pointed out. Both the proposed setups are based on state-of-the-art technological achievements. By extensive numerical experiments we show that their performance is remarkably good, even for the implementation of long sequences of gates used to simulate interesting physical models. Therefore, the here examined systems are really promising buildingblocks of future scalable architectures and can be used for proof-of-principle experiments of quantum information processing and quantum simulation.

Fundamental characteristics of space-charge waves in photorefractive sillenite crystals

We investigate experimentally the fundamental characteristics of space-charge waves excited in a photorefractive crystal of Bi12SiO20. Features such as their transient rise and decay as well as their steady-state frequency response are investigated. Based on this, we find the dependence of the space-charge waves' quality factor on spatial frequency and electric-field biasing. The experimental findings are compared with the linear space-charge wave theory developed previously by Sturman et al. [J. Opt. Sec. Am. B 10, 1919 (1993)].

Hydrodynamic modeling of mineral wool fiber suspensions in a two-dimensional channel flow

A consequence of a loss of coolant accident is the damage of adjacent insulation materials (IM). IM may then be transported to the containment sump strainers where water is drawn into the ECCS (emergency core cooling system). Blockage of the strainers by IM lead to an increased pressure drop acting on the operating ECCS pumps. IM can also penetrate the strainers, enter the reactor coolant system and then accumulate in the reactor pressure vessel. An experimental and theoretical study that concentrates on mineral wool fiber transport in the containment sump and the ECCS is being performed. The study entails fiber generation and the assessment of fiber transport in single and multi-effect experiments. The experiments include measurement of the terminal settling velocity, the strainer pressure drop, fiber sedimentation and resuspension in a channel flow and jet flow in a rectangular tank. An integrated test facility is also operated to assess the compounded effects. Each experimental facility is used to provide data for the validation of equivalent computational fluid dynamic models. The channel flow facility allows the determination of the steady state distribution of the fibers at different flow velocities. The fibers are modeled in the Eulerian-Eulerian reference frame as spherical wetted agglomerates. The fiber agglomerate size, density, the relative viscosity of the fluid-fiber mixture and the turbulent dispersion of the fibers all affect the steady state accumulation of fibers at the channel base. In the current simulations, two fiber phases are separately considered. The particle size is kept constant while the density is modified, which affects both the terminal velocity and volume fraction. The relative viscosity is only significant at higher concentrations. The numerical model finds that the fibers accumulate at the channel base even at high velocities; therefore, modifications to the drag and turbulent dispersion forces can be made to reduce fiber accumulation.

Fundamental characteristics of space-charge waves in photorefractive sillenite crystals

We investigate experimentally the fundamental characteristics of space-charge waves excited in a photorefractive crystal of Bi12SiO20. Features such as their transient rise and decay as well as their steady-state frequency response are investigated. Based on this, we find the dependence of the space-charge waves' quality factor on spatial frequency and electric-field biasing. The experimental findings are compared with the linear space-charge wave theory developed previously by Sturman et al. [J. Opt. Sec. Am. B 10, 1919 (1993)].

Study of neutral Fe(III) complexes of pyridoxal-N-substituted thiosemicarbazone with desolvation induced spin-state transformation above room temperature

The preparation and characterization of two new neutral ferric complexes with desolvation-induced discontinuous spin-state transformation above room temperature are reported. The compounds, Fe(Hthpy)(thpy).CH3OH.3H2O (1) and Fe(Hmthpy)(mthpy).2H2O (2), are low-spin (LS) at room temperature and below, whereas their nonsolvated forms are high-spin (HS), exhibiting zero-field splitting. In these complexes, Hthpy, Hmthpy, and thpy, mthpy are the deprotonated forms of pyridoxal thiosemicarbazone and pyridoxal methylthiosemicarbazone, respectively; each is an O,N,S-tridentate ligand. The molecular structures have been determined at 100(1) K using single-crystal X-ray diffraction techniques and resulted in a triclinic system (space group P1) and monoclinic unit cell (space group P21/c) for 1 and 2, respectively. Structures were refined to the final error indices, where RF = 0.0560 for 1 and RF = 0.0522 for 2. The chemical inequivalence of the ligands was clearly established, for the "extra" hydrogen atom on the monodeprotonated ligands (Hthpy, Hmthpy) was found to be bound to the nitrogen of the pyridine ring. The ligands are all of the thiol form; the doubly deprotonated chelates (thpy, mthpy) have C-S bond lengths slightly longer than those of the singly deprotonated forms. There is a three-dimensional network of hydrogen bonds in both compounds. The discontinuous spin-state transformation is accompanied with liberation of solvate molecules. This is evidenced also from DSC analysis. Heat capacity data for the LS and HS phases are tabulated at selected temperatures, the values of the enthalpy and entropy changes connected with the change of spin state were reckoned at DeltaH = 12.5 0.3 kJ mol-1 and DeltaS = 33.3 0.8 J mol-1 K-1, respectively, for 1 and DeltaH = 6.5 0.3 kJ mol-1 and DeltaS = 17.6 0.8 J mol-1 K-1, respectively, for 2

Bandwidth management of WiMAX systems and performance modeling

WiMAX has been introduced as a competitive alternative for metropolitan broadband wireless access technologies. It is connection oriented and it can provide very high data rates, large service coverage, and flexible quality of services (QoS). Due to the large number of connections and flexible QoS supported by WiMAX, the uplink access in WiMAX networks is very challenging since the medium access control (MAC) protocol must efficiently manage the bandwidth and related channel allocations. In this paper, we propose and investigate a cost-effective WiMAX bandwidth management scheme, named the WiMAX partial sharing scheme (WPSS), in order to provide good QoS while achieving better bandwidth utilization and network throughput. The proposed bandwidth management scheme is compared with a simple but inefficient scheme, named the WiMAX complete sharing scheme (WCPS). A maximum entropy (ME) based analytical model (MEAM) is proposed for the performance evaluation of the two bandwidth management schemes. The reason for using MEAM for the performance evaluation is that MEAM can efficiently model a large-scale system in which the number of stations or connections is generally very high, while the traditional simulation and analytical (e.g., Markov models) approaches cannot perform well due to the high computation complexity. We model the bandwidth management scheme as a queuing network model (QNM) that consists of interacting multiclass queues for different service classes. Closed form expressions for the state and blocking probability distributions are derived for those schemes. Simulation results verify the MEAM numerical results and show that WPSS can significantly improve the network's performance compared to WCPS.

eHabitat, a multi-purpose Web Processing Service for ecological modeling

The number of interoperable research infrastructures has increased significantly with the growing awareness of the efforts made by the Global Earth Observation System of Systems (GEOSS). One of the Societal Benefit Areas (SBA) that is benefiting most from GEOSS is biodiversity, given the costs of monitoring the environment and managing complex information, from space observations to species records including their genetic characteristics. But GEOSS goes beyond simple data sharing to encourage the publishing and combination of models, an approach which can ease the handling of complex multi-disciplinary questions. It is the purpose of this paper to illustrate these concepts by presenting eHabitat, a basic Web Processing Service (WPS) for computing the likelihood of finding ecosystems with equal properties to those specified by a user. When chained with other services providing data on climate change, eHabitat can be used for ecological forecasting and becomes a useful tool for decision-makers assessing different strategies when selecting new areas to protect. eHabitat can use virtually any kind of thematic data that can be considered as useful when defining ecosystems and their future persistence under different climatic or development scenarios. The paper will present the architecture and illustrate the concepts through case studies which forecast the impact of climate change on protected areas or on the ecological niche of an African bird.