Memory-induced anomalous dynamics in a minimal random walk model

A discrete-time dynamics of a non-Markovian random walker is analyzed using a minimal model where memory of the past drives the present dynamics. In recent work N. Kumar et al., Phys. Rev. E 82, 021101 (2010)] we proposed a model that exhibits asymptotic superdiffusion, normal diffusion, and subdiffusion with the sweep of a single parameter. Here we propose an even simpler model, with minimal options for the walker: either move forward or stay at rest. We show that this model can also give rise to diffusive, subdiffusive, and superdiffusive dynamics at long times as a single parameter is varied. We show that in order to have subdiffusive dynamics, the memory of the rest states must be perfectly correlated with the present dynamics. We show explicitly that if this condition is not satisfied in a unidirectional walk, the dynamics is only either diffusive or superdiffusive (but not subdiffusive) at long times.

Optimal Blood Glucose Regulation using Single Network Adaptive Critics

Diabetes is a serious disease during which the body's production and use of insulin is impaired, causing glucose concentration level toincrease in the bloodstream. Regulating blood glucose levels as close to normal as possible, leads to a substantial decrease in long term complications of diabetes. In this paper, an intelligent neural network on-line optimal feedback treatment strategy based on nonlinear optimal control theory is presented for the disease using subcutaneous treatment strategy. A simple mathematical model of the nonlinear dynamics of glucose and insulin interaction in the blood system is considered based on the Bergman's minimal model. A glucose infusion term representing the effect of glucose intake resulting from a meal is introduced into the model equations. The efficiency of the proposed controllers is shown taking random parameters and random initial conditions in presence of physical disturbances like food intake. A comparison study with linear quadratic regulator theory brings Out the advantages of the nonlinear control synthesis approach. Simulation results show that unlike linear optimal control, the proposed on-line continuous infusion strategy never leads to severe hypoglycemia problems.

A Triggering Mechanism for Enhanced Star-Formation in Colliding Galaxies

We propose a physical mechanism to explain the origin of the intense burst of massive-star formation seen in colliding/merging, gas-rich, field spiral galaxies. We explicitly take account of the different parameters for the two main mass components, H-2 and H I, of the interstellar medium within a galaxy and follow their consequent different evolution during a collision between two galaxies. We also note that, in a typical spiral galaxy-like our galaxy, the Giant Molecular Clouds (GMCs) are in a near-virial equilibrium and form the current sites of massive-star formation, but have a low star formation rate. We show that this star formation rate is increased following a collision between galaxies. During a typical collision between two field spiral galaxies, the H I clouds from the two galaxies undergo collisions at a relative velocity of approximately 300 km s-1. However, the GMCs, with their smaller volume filling factor, do not collide. The collisions among the H I clouds from the two galaxies lead to the formation of a hot, ionized, high-pressure remnant gas. The over-pressure due to this hot gas causes a radiative shock compression of the outer layers of a preexisting GMC in the overlapping wedge region. This makes these layers gravitationally unstable, thus triggering a burst of massive-star formation in the initially barely stable GMCs.The resulting value of the typical IR luminosity from the young, massive stars from a pair of colliding galaxies is estimated to be approximately 2 x 10(11) L., in agreement with the observed values. In our model, the massive-star formation occurs in situ in the overlapping regions of a pair of colliding galaxies. We can thus explain the origin of enhanced star formation over an extended, central area approximately several kiloparsecs in size, as seen in typical colliding galaxies, and also the origin of starbursts in extranuclear regions of disk overlap as seen in Arp 299 (NGC 3690/IC 694) and in Arp 244 (NGC 4038/39). Whether the IR emission from the central region or that from the surrounding extranuclear galactic disk dominates depends on the geometry and the epoch of the collision and on the initial radial gas distribution in the two galaxies. In general, the central starburst would be stronger than that in the disks, due to the higher preexisting gas densities in the central region. The burst of star formation is expected to last over a galactic gas disk crossing time approximately 4 x 10(7) yr. We can also explain the simultaneous existence of nearly normal CO galaxy luminosities and shocked H-2 gas, as seen in colliding field galaxies.This is a minimal model, in that the only necessary condition for it to work is that there should be a sufficient overlap between the spatial gas distributions of the colliding galaxy pair.

Electron transfer reaction in the marcus inverted region: role of high frequency vibrational modes

A theoretical study of the dynamics of photo-electron transfer reactions in the Marcus inverted regime is presented. This study is motivated partly by the recent proposal of Barbara et al. (J. Phys. Chem. 96, 3728, 1991) that a minimal model of an electron transfer reaction should consist of a polar solvent mode (X), a low frequency vibrational mode (Q) and one high frequency mode (q). Interplay between these modes may be responsible for the crossover observed in the dynamics from a solvent controlled to a vibrational controlled electron transfer. The following results have been obtained. (i) In the case of slowly relaxing solvents, the proximity of the point of excitation to an effective sink on the excited surface is critical in determining the decay of the reactant population. This is because the Franck-Condon overlap between the reactant ground and the product excited states decreases rapidly with increase in the quantum number of the product vibrational state. (ii) Non-exponential solvation dynamics has an important effect in determining the rates of electron transfer. Especially, a biphasic solvation and a large coupling between the reactant and the product states both may be needed to explain the experimental results. ©1996 American Institute of Physics

Band Structures of Bilayer Graphene Superlattices

We formulate a low energy effective Hamiltonian to study superlattices in bilayer graphene (BLG) using a minimal model which supports quadratic band touching points. We show that a one dimensional (1D) periodic modulation of the chemical potential or the electric field perpendicular to the layers leads to the generation of zero-energy anisotropic massless Dirac fermions and finite energy Dirac points with tunable velocities. The electric field superlattice maps onto a coupled chain model comprised of ``topological'' edge modes. 2D superlattice modulations are shown to lead to gaps on the mini-Brillouin zone boundary but do not, for certain symmetries, gap out the quadratic band touching point. Such potential variations, induced by impurities and rippling in biased BLG, could lead to subgap modes which are argued to be relevant to understanding transport measurements.

Dynamics of intermittent force fluctuations in vesicular nanotubulation

Irregular force fluctuations are seen in most nanotubulation experiments. The dynamics behind their presence has, however, been neither commented upon nor modeled. A simple estimate of the mean energy dissipated in force drops turns out to be several times the thermal energy. This coupled with the rate dependent nature of the deformation reported in several experiments point to a dynamical origin of the serrations. We simplify the whole process of tether formation through a three-stage model of successive deformations of sphere to ellipsoid, neck-formation, and tubule birth and extension. Based on this, we envisage a rate-softening frictional force at the neck that must be overcome before a nanotube can be pulled out. Our minimal model includes elastic and visco-elastic deformation of the vesicle, and has built-in dependence on pull velocity, vesicle radius, and other material parameters, enabling us to capture various kinds of serrated force-extension curves for different parameter choices. Serrations are predicted in the nanotubulation region. Other features of force-extension plots reported in the literature such as a plateauing serrated region beyond a force drop, serrated flow region with a small positive slope, an increase in the elastic threshold with pull velocity, force-extension curves for vesicles with larger radius lying lower than those for smaller radius, are all also predicted by the model. A toy model is introduced to demonstrate that the role of the friction law is limited to inducing stick-slip oscillations in the force, and all other qualitative and quantitative features emerging from the model can only be attributed to other physical mechanisms included in the deformation dynamics of the vesicle. (C) 2014 AIP Publishing LLC.

Active Viscoelastic Matter: From Bacterial Drag Reduction to Turbulent Solids

A paradigm for internally driven matter is the active nematic liquid crystal, whereby the equations of a conventional nematic are supplemented by a minimal active stress that violates time-reversal symmetry. In practice, active fluids may have not only liquid-crystalline but also viscoelastic polymer degrees of freedom. Here we explore the resulting interplay by coupling an active nematic to a minimal model of polymer rheology. We find that adding a polymer can greatly increase the complexity of spontaneous flow, but can also have calming effects, thereby increasing the net throughput of spontaneous flow along a pipe (a ``drag-reduction'' effect). Remarkably, active turbulence can also arise after switching on activity in a sufficiently soft elastomeric solid.

Using tau polarization to probe the stau co-annihilation region of the minimal supergravity model at the LHC

The minimal supergravity model predicts the polarization of the tau coming from the stau to bino decay in the co-annihilation region to +1. This can be exploited to extract this soft tau signal at LHC and also to measure the tiny mass differences between the stau and the bi lightest superparticle. Moreover, this strategy will be applicable for a wider class of bino lightest superparticle models, where the lighter stau has a right component at least of similar size as the left.

Invisible decays of the heavier Higgs boson in the minimal supersymmetric standard model

We consider the possibility that the heavier CP-even Higgs boson (H-0) in the minimal supersymmetric standard model (MSSM) decays invisibly into neutralinos in the light of the recent discovery of the 126 GeV resonance at the CERN Large Hadron Collider (LHC). For this purpose we consider the minimal supersymmetric standard model with universal, nonuniversal and arbitrary boundary conditions on the supersymmetry breaking gaugino mass parameters at the grand unified scale. Typically, scenarios with universal and nonuniversal gaugino masses do not allow invisible decays of the lightest Higgs boson (h(0)), which is identified with the 126 GeV resonance, into the lightest neutralinos in the MSSM. With arbitrary gaugino masses at the grand unified scale, such an invisible decay is possible. The second lightest Higgs boson can decay into various invisible final states for a considerable region of the MSSM parameter space with arbitrary gaugino masses as well as with the gaugino masses restricted by universal and nonuniversal boundary conditions at the grand unified scale. The possibility of the second lightest Higgs boson of the MSSM decaying into invisible channels is more likely for arbitrary gaugino masses at the grand unified scale. The heavier Higgs boson decay into lighter particles leads to the intriguing possibility that the entire Higgs boson spectrum of the MSSM may be visible at the LHC even if it decays invisibly, during the searches for an extended Higgs boson sector at the LHC. In such a scenario the nonobservation of the extended Higgs sector of the MSSM may carefully be used to rule out regions of the MSSM parameter space at the LHC.

Lepton masses and flavor violation in the Randall-Sundrum model

Lepton masses and mixing angles via localization of 5-dimensional fields in the bulk are revisited in the context of Randall-Sundrum models. The Higgs is assumed to be localized on the IR brane. Three cases for neutrino masses are considered: (a) The higher-dimensional neutrino mass operator (LH.LH), (b) Dirac masses, and (c) Type I seesaw with bulk Majorana mass terms. Neutrino masses and mixing as well as charged lepton masses are fit in the first two cases using chi(2) minimization for the bulk mass parameters, while varying the O(1) Yukawa couplings between 0.1 and 4. Lepton flavor violation is studied for all the three cases. It is shown that large negative bulk mass parameters are required for the right-handed fields to fit the data in the LH.LH case. This case is characterized by a very large Kaluza-Klein (KK) spectrum and relatively weak flavor-violating constraints at leading order. The zero modes for the charged singlets are composite in this case, and their corresponding effective 4-dimensional Yukawa couplings to the KK modes could be large. For the Dirac case, good fits can be obtained for the bulk mass parameters, c(i), lying between 0 and 1. However, most of the ``best-fit regions'' are ruled out from flavor-violating constraints. In the bulk Majorana terms case, we have solved the profile equations numerically. We give example points for inverted hierarchy and normal hierarchy of neutrino masses. Lepton flavor violating rates are large for these points. We then discuss various minimal flavor violation schemes for Dirac and bulk Majorana cases. In the Dirac case with minimal-flavor-violation hypothesis, it is possible to simultaneously fit leptonic masses and mixing angles and alleviate lepton flavor violating constraints for KK modes with masses of around 3 TeV. Similar examples are also provided in the Majorana case.

Survey of vector-like fermion extensions of the Standard Model and their phenomenological implications

With the renewed interest in vector-like fermion extensions of the Standard Model, we present here a study of multiple vector-like theories and their phenomenological implications. Our focus is mostly on minimal flavor conserving theories that couple the vector-like fermions to the SM gauge fields and mix only weakly with SM fermions so as to avoid flavor problems. We present calculations for precision electroweak and vector-like state decays, which are needed to investigate compatibility with currently known data. We investigate the impact of vector-like fermions on Higgs boson production and decay, including loop contributions, in a wide variety of vector-like extensions and their parameter spaces.

Turbulent electrical activity at sharp-edged inexcitable obstacles in a model for human cardiac tissue

Wave propagation around various geometric expansions, structures, and obstacles in cardiac tissue may result in the formation of unidirectional block of wave propagation and the onset of reentrant arrhythmias in the heart. Therefore, we investigated the conditions under which reentrant spiral waves can be generated by high-frequency stimulation at sharp-edged obstacles in the ten Tusscher-Noble-Noble-Panfilov (TNNP) ionic model for human cardiac tissue. We show that, in a large range of parameters that account for the conductance of major inward and outward ionic currents of the model fast inward Na+ current (INa), L-type slow inward Ca2+ current (I-CaL), slow delayed-rectifier current (I-Ks), rapid delayed-rectifier current (I-Kr), inward rectifier K+ current (I-K1)], the critical period necessary for spiral formation is close to the period of a spiral wave rotating in the same tissue. We also show that there is a minimal size of the obstacle for which formation of spirals is possible; this size is similar to 2.5 cm and decreases with a decrease in the excitability of cardiac tissue. We show that other factors, such as the obstacle thickness and direction of wave propagation in relation to the obstacle, are of secondary importance and affect the conditions for spiral wave initiation only slightly. We also perform studies for obstacle shapes derived from experimental measurements of infarction scars and show that the formation of spiral waves there is facilitated by tissue remodeling around it. Overall, we demonstrate that the formation of reentrant sources around inexcitable obstacles is a potential mechanism for the onset of cardiac arrhythmias in the presence of a fast heart rate.

Interplay of electrostatics and lipid packing determines the binding of charged polymer coated nanoparticles to model membranes

Understanding of nanoparticle-membrane interactions is useful for various applications of nanoparticles like drug delivery and imaging. Here we report on the studies of interaction between hydrophilic charged polymer coated semiconductor quantum dot nanoparticles with model lipid membranes. Atomic force microscopy and X-ray reflectivity measurements suggest that cationic nanoparticles bind and penetrate bilayers of zwitterionic lipids. Penetration and binding depend on the extent of lipid packing and result in the disruption of the lipid bilayer accompanied by enhanced lipid diffusion. On the other hand, anionic nanoparticles show minimal membrane binding although, curiously, their interaction leads to reduction in lipid diffusivity. It is suggested that the enhanced binding of cationic QDs at higher lipid packing can be understood in terms of the effective surface potential of the bilayers which is tunable through membrane lipid packing. Our results bring forth the subtle interplay of membrane lipid packing and electrostatics which determine nanoparticle binding and penetration of model membranes with further implications for real cell membranes.

Gravitational rescue of minimal gauge mediation

Gravity mediated supersymmetry breaking becomes comparable to gauge mediated supersymmetry breaking contributions when messenger masses are close to the GUT scale. By suitably arranging the gravity contributions, one can modify the soft supersymmetry breaking sector to generate a large stop mixing parameter and a light Higgs mass of 125 GeV. In this kind of hybrid models, however, the nice features of gauge mediation like flavor conservation, etc. are lost. To preserve the nice features, gravitational contributions should become important for lighter messenger masses and should be important only for certain fields. This is possible when the hidden sector contains multiple (at least two) spurions with hierarchical vacuum expectation values. In this case, the gravitational contributions can be organized to be ``just right.'' We present a complete model with two spurion hidden sector where the gravitational contribution is from a warped flavor model in a Randall-Sundrum setting. Along the way, we present simple expressions to handle renormalization group equations when supersymmetry is broken by two different sectors at two different scales.