Branching Particle Systems and Superprocess

In this paper, we describe a system of particles that perform independent random motions in space and at the end of their lifetimes give birth to a random number of offspring. We show that the system in the large density, small mass, rapid branching or long time scale limit converges to a measure-valued diffusion called the superprocess.

Mathematical Model of Viral Kinetics In Vitro Estimates the Number of E2-CD81 Complexes Necessary for Hepatitis C Virus Entry

Interaction between the hepatitis C virus (HCV) envelope protein E2 and the host receptor CD81 is essential for HCV entry into target cells. The number of E2-CD81 complexes necessary for HCV entry has remained difficult to estimate experimentally. Using the recently developed cell culture systems that allow persistent HCV infection in vitro, the dependence of HCV entry and kinetics on CD81 expression has been measured. We reasoned that analysis of the latter experiments using a mathematical model of viral kinetics may yield estimates of the number of E2-CD81 complexes necessary for HCV entry. Here, we constructed a mathematical model of HCV viral kinetics in vitro, in which we accounted explicitly for the dependence of HCV entry on CD81 expression. Model predictions of viral kinetics are in quantitative agreement with experimental observations. Specifically, our model predicts triphasic viral kinetics in vitro, where the first phase is characterized by cell proliferation, the second by the infection of susceptible cells and the third by the growth of cells refractory to infection. By fitting model predictions to the above data, we were able to estimate the threshold number of E2-CD81 complexes necessary for HCV entry into human hepatoma-derived cells. We found that depending on the E2-CD81 binding affinity, between 1 and 13 E2-CD81 complexes are necessary for HCV entry. With this estimate, our model captured data from independent experiments that employed different HCV clones and cells with distinct CD81 expression levels, indicating that the estimate is robust. Our study thus quantifies the molecular requirements of HCV entry and suggests guidelines for intervention strategies that target the E2-CD81 interaction. Further, our model presents a framework for quantitative analyses of cell culture studies now extensively employed to investigate HCV infection.

The force distribution function of an oscillator model of polymer stretching at constant velocity

Recent simulations of the stretching of tethered biopolymers at a constant speed v (Ponmurugan and Vemparala, 2011 Phys. Rev. E 84 060101(R)) have suggested that for any time t, the distribution of the fluctuating forces f responsible for chain deformation is governed by a relation of the form P(+ f)/ P(- f) = expgamma f], gamma being a coefficient that is solely a function of v and the temperature T. This result, which is reminiscent of the fluctuation theorems applicable to stochastic trajectories involving thermodynamic variables, is derived in this paper from an analytical calculation based on a generalization of Mazonka and Jarzynski's classic model of dragged particle dynamics Mazonka and Jarzynski, 1999 arXiv:cond-\textbackslashmat/9912121v1]. However, the analytical calculations suggest that the result holds only if t >> 1 and the force fluctuations are driven by white rather than colored noise; they further suggest that the coefficient gamma in the purported theorem varies not as v(0.15)T-(0.7), as indicated by the simulations, but as vT(-1).

Piecewise Linearization Technique for Compact Charge Modeling of Independent DG MOSFET

Charge linearization techniques have been used over the years in advanced compact models for bulk and double-gate MOSFETs in order to approximate the position along the channel as a quadratic function of the surface potential (or inversion charge densities) so that the terminal charges can be expressed as a compact closed-form function of source and drain end surface potentials (or inversion charge densities). In this paper, in case of the independent double-gate MOSFETs, we show that the same technique could be used to model the terminal charges quite accurately only when the 1-D Poisson solution along the channel is fully hyperbolic in nature or the effective gate voltages are same. However, for other bias conditions, it leads to significant error in terminal charge computation. We further demonstrate that the amount of nonlinearity that prevails between the surface potentials along the channel actually dictates if the conventional charge linearization technique could be applied for a particular bias condition or not. Taking into account this nonlinearity, we propose a compact charge model, which is based on a novel piecewise linearization technique and shows excellent agreement with numerical and Technology Computer-Aided Design (TCAD) simulations for all bias conditions and also preserves the source/drain symmetry which is essential for Radio Frequency (RF) circuit design. The model is implemented in a professional circuit simulator through Verilog-A, and simulation examples for different circuits verify good model convergence.

LOCAL TWO-DIMENSIONAL PARTICLE-IN-CELL SIMULATIONS OF THE COLLISIONLESS MAGNETOROTATIONAL INSTABILITY

The magnetorotational instability (MRI) is a crucial mechanism of angular momentum transport in a variety of astrophysical accretion disks. In systems accreting at well below the Eddington rate, such as the central black hole in the Milky Way (Sgr A*), the plasma in the disk is essentially collisionless. We present a nonlinear study of the collisionless MRI using first-principles particle-in-cell plasma simulations. We focus on local two-dimensional (axisymmetric) simulations, deferring more realistic three-dimensional simulations to future work. For simulations with net vertical magnetic flux, the MRI continuously amplifies the magnetic field, B, until the Alfven velocity, v(A), is comparable to the speed of light, c (independent of the initial value of v(A)/c). This is consistent with the lack of saturation of MRI channel modes in analogous axisymmetric MHD simulations. The amplification of the magnetic field by the MRI generates a significant pressure anisotropy in the plasma (with the pressure perpendicular to B being larger than the parallel pressure). We find that this pressure anisotropy in turn excites mirror modes and that the volume-averaged pressure anisotropy remains near the threshold for mirror mode excitation. Particle energization is due to both reconnection and viscous heating associated with the pressure anisotropy. Reconnection produces a distinctive power-law component in the energy distribution function of the particles, indicating the likelihood of non-thermal ion and electron acceleration in collisionless accretion disks. This has important implications for interpreting the observed emission-from the radio to the gamma-rays-of systems such as Sgr A*.

Probing ultrafast carrier dynamics, nonlinear absorption and refraction in core-shell silicon nanowires

We investigate the relaxation dynamics of photogenerated carriers in silicon nanowires consisting of a crystalline core and a surrounding amorphous shell, using femtosecond time-resolved differential reflectivity and transmission spectroscopy at 3.15 eV and 1.57 eV photon energies. The complex behaviour of the differential transmission and reflectivity transients is the mixed contributions from the crystalline core and the amorphous silicon on the nanowire surface and the substrate where competing effects of state-filling and photoinduced absorption govern the carrier dynamics. Faster relaxation rates are observed on increasing the photogenerated carrier density. Independent experimental results on crystalline silicon-on-sapphire (SOS) help us in separating the contributions from the carrier dynamics in crystalline core and the amorphous regions in the nanowire samples. Further, single-beam z-scan nonlinear transmission experiments at 1.57 eV in both open- and close-aperture configurations yield two-photon absorption coefficient beta (similar to 3 cm/GW) and nonlinear refraction coefficient gamma (-2.5 x 10 (-aEuro parts per thousand 4) cm(2)/GW).

Asymptotic expansions for the coupled wavenumbers in an infinite orthotropic flexible fluid-filled cylindrical shell

Analytical expressions are found for the coupled wavenumbers in flexible, fluid-filled, circular cylindrical orthotropic shells using the asymptotic methods. These expressions are valid for arbitrary circumferential orders. The Donnell-Mushtari shell theory is used to model the shell and the effect of the fluid is introduced through the fluid-loading parameter mu. The orthotropic problem is posed as a perturbation on the corresponding isotropic problem by defining a suitable orthotropy parameter epsilon, which is a measure of the degree of orthotropy. For the first study, an isotropic shell is considered (by setting epsilon = 0) and expansions are found for the coupled wavenumbers using a regular perturbation approach. In the second study, asymptotic expansions are found for the coupled wavenumbers in the limit of small orthotropy (epsilon << 1). For each study, isotropy and orthotropy, expansions are found for small and large values of the fluid-loading parameter mu. All the asymptotic solutions are compared with numerical solutions to the coupled dispersion relation and the match is seen to be good. The differences between the isotropic and orthotropic solutions are discussed. The main contribution of this work lies in extending the existing literature beyond in vacuo studies to the case of fluid-filled shells (isotropic and orthotropic).

Projecting low-dimensional chaos from spatiotemporal dynamics in a model for plastic instability

We investigate the possibility of projecting low-dimensional chaos from spatiotemporal dynamics of a model for a kind of plastic instability observed under constant strain rate deformation conditions. We first discuss the relationship between the spatiotemporal patterns of the model reflected in the nature of dislocation bands and the nature of stress serrations. We show that at low applied strain rates, there is a one-to-one correspondence with the randomly nucleated isolated bursts of mobile dislocation density and the stress drops. We then show that the model equations are spatiotemporally chaotic by demonstrating the number of positive Lyapunov exponents and Lyapunov dimension scale with the system size at low and high strain rates. Using a modified algorithm for calculating correlation dimension density, we show that the stress-strain signals at low applied strain rates corresponding to spatially uncorrelated dislocation bands exhibit features of low-dimensional chaos. This is made quantitative by demonstrating that the model equations can be approximately reduced to space-independent model equations for the average dislocation densities, which is known to be low-dimensionally chaotic. However, the scaling regime for the correlation dimension shrinks with increasing applied strain rate due to increasing propensity for propagation of the dislocation bands.

Application of Artificial Neural Networks and Particle Swarm Optimization for the Management of Groundwater Resources

Ground management problems are typically solved by the simulation-optimization approach where complex numerical models are used to simulate the groundwater flow and/or contamination transport. These numerical models take a lot of time to solve the management problems and hence become computationally expensive. In this study, Artificial Neural Network (ANN) and Particle Swarm Optimization (PSO) models were developed and coupled for the management of groundwater of Dore river basin in France. The Analytic Element Method (AEM) based flow model was developed and used to generate the dataset for the training and testing of the ANN model. This developed ANN-PSO model was applied to minimize the pumping cost of the wells, including cost of the pipe line. The discharge and location of the pumping wells were taken as the decision variable and the ANN-PSO model was applied to find out the optimal location of the wells. The results of the ANN-PSO model are found similar to the results obtained by AEM-PSO model. The results show that the ANN model can reduce the computational burden significantly as it is able to analyze different scenarios, and the ANN-PSO model is capable of identifying the optimal location of wells efficiently.

Variance-Reduced Particle Filters for Structural System Identification Problems

A few variance reduction schemes are proposed within the broad framework of a particle filter as applied to the problem of structural system identification. Whereas the first scheme uses a directional descent step, possibly of the Newton or quasi-Newton type, within the prediction stage of the filter, the second relies on replacing the more conventional Monte Carlo simulation involving pseudorandom sequence with one using quasi-random sequences along with a Brownian bridge discretization while representing the process noise terms. As evidenced through the derivations and subsequent numerical work on the identification of a shear frame, the combined effect of the proposed approaches in yielding variance-reduced estimates of the model parameters appears to be quite noticeable. DOI: 10.1061/(ASCE)EM.1943-7889.0000480. (C) 2013 American Society of Civil Engineers.

Time variant reliability model updating in instrumented dynamical systems based on Girsanov's transformation

The problem of updating the reliability of instrumented structures based on measured response under random dynamic loading is considered. A solution strategy within the framework of Monte Carlo simulation based dynamic state estimation method and Girsanov's transformation for variance reduction is developed. For linear Gaussian state space models, the solution is developed based on continuous version of the Kalman filter, while, for non-linear and (or) non-Gaussian state space models, bootstrap particle filters are adopted. The controls to implement the Girsanov transformation are developed by solving a constrained non-linear optimization problem. Numerical illustrations include studies on a multi degree of freedom linear system and non-linear systems with geometric and (or) hereditary non-linearities and non-stationary random excitations.

Bimetallic core-shell nanocomposites using weak reducing agent and their transformation to alloy nanostructures

An in situ seeding growth methodology towards the preparation of core-shell nanoparticles composed of noble metals has been developed by employing trimethylamine borane (TMAB) as the reducing agent. Being a weak reducing agent, TMAB is able to distinguish the smallest reduction potential window of any two metals which renders selective reduction of metal ions thus affording a core-shell architecture of the nanoparticles. A dramatic effect of solvent was noted during the reduction of Ag+ ions: an immediate reduction took place at room temperature when dry THF was used as solvent however, usage of wet THF (THF used directly from the bottle) brings out the reduction only at reflux conditions. In the case of Au and Pd nanoparticles, preparation was found to be independent of the quality of solvent used. Au nanoparticles are realized at room temperature whereas reflux conditions are required in the case of Pd nanoparticles. This difference in behavior of the monometallic nanoparticles was successfully exploited to construct different noble metal nanoparticles with core-shell architectures such as Au@Ag, Ag@Au, and Ag@Pd. Transformation of these core-shell nanoparticles to their thermodynamically stable alloy counterparts is also demonstrated under very mild conditions reported to date.

Implications of Higgs boson to diphoton decay rate in the bilinear R-parity violating supersymmetric model

The Large Hadron Collider has recently discovered a Higgs-like particle having a mass around 125 GeVand also indicated that there is an enhancement in the Higgs to diphoton decay rate as compared to that in the standard model. We have studied implications of these discoveries in the bilinear R-parity violating supersymmetric model, whose main motivation is to explain the nonzero masses for neutrinos. The R-parity violating parameters in this model are epsilon and b(epsilon), and these parameters determine the scale of neutrino masses. If the enhancement in the Higgs to diphoton decay rate is true, then we have found epsilon greater than or similar to 0.01 GeV and b epsilon similar to 1 GeV2 in order to be compatible with the neutrino oscillation data. Also, in the above mentioned analysis, we can determine the soft masses of sleptons (m(L)) and CP-odd Higgs boson mass (mA). We have estimated that m(L) greater than or similar to 300 GeV and m(A) greater than or similar to 700 GeV. We have also commented on the allowed values of epsilon and b(epsilon), in case there is no enhancement in the Higgs to diphoton decay rate. Finally, we present a model to explain the smallness of epsilon and b(epsilon).

EFFECT OF CORE-SHELL STRUCTURE OF HYDROGEL BEADS ON THE THRESHOLD CONCENTRATION OF WATER FOR SWELLING AND ITS PH SENSITIVITY

In this paper we investigate the effect of core-shell structure of Sodium Alginate based hydrogel beads and their size on certain activation threshold concentration of water for applications in swelling and pH sensing. This type of hydrogel experiences diffusive pressure due to transport of certain free charges across its interface with a solvent or electrolyte. This process is essentially a dynamic equilibrium of the electric force field, stress in the polymeric network with cage like structure and molecular diffusion including phase transformation due to pressure imbalance between the hydrogel and its surroundings. The effect of pH of the solvant on the swelling rate of these beads has been studied experimentally. A mathematical model of the swelling process has been developed by considering Nernst-Planck equation representing the migration of mobile ions and Er ions, Poisson equation representing the equilibrium of the electric field and mechanical field equation representing swelling of the gel. An attempt has been made to predict the experimentally observed phenomena using these numerical simulations. It is observed experimentally that certain minimum concentration called activation threshold concentration of the water molecules must be present in the hydrogel in order to activate the swelling process. For the required activation threshold concentration of water in the beads, the pH induced change in the rate of swelling is also investigated. This effect is analyzed for various different core-shell structures of the beads.

Bilinear tension softening diagrams of concrete mixes corresponding to their size-independent specific fracture energy

The study of the fracture behaviour of concrete structures using the fictitious crack model requires two fracture properties of the concrete mix, namely, the size-independent specific fracture energy G(F). and the corresponding tension softening relation sigma(w) between the residual stress carrying capacity sigma and the crack opening w in the fracture process zone ahead of a real crack. In this paper, bi-linear tension softening diagrams of three different concrete mixes, ranging in compressive strength from 57 to 122 MPa whose size-independent specific fracture energy has been previously determined, have been constructed in an inverse manner based on the concept of a non-linear hinge from the load-crack mouth opening plots of notched three-point bend beams. (C) 2013 Elsevier Ltd. All rights reserved.