Zero-variance principle for Monte Carlo algorithms

We present a general approach to greatly increase at little cost the efficiency of Monte Carlo algorithms. To each observable to be computed we associate a renormalized observable (improved estimator) having the same average but a different variance. By writing down the zero-variance condition a fundamental equation determining the optimal choice for the renormalized observable is derived (zero-variance principle for each observable separately). We show, with several examples including classical and quantum Monte Carlo calculations, that the method can be very powerful.Comment: 9 pages, Latex, to appear in Phys. Rev. Let

Simulation of Potts models with real q and no critical slowing down

A Monte Carlo algorithm is proposed to simulate ferromagnetic q-state Potts model for any real q>0. A single update is a random sequence of disordering and deterministic moves, one for each link of the lattice. A disordering move attributes a random value to the link, regardless of the state of the system, while in a deterministic move this value is a state function. The relative frequency of these moves depends on the two parameters q and beta. The algorithm is not affected by critical slowing down and the dynamical critical exponent z is exactly vanishing. We simulate in this way a 3D Potts model in the range 2<q<3 for estimating the critical value q_c where the thermal transition changes from second-order to first-order, and find q_c=2.620(5).Comment: 5 pages, 3 figures slightly extended version, to appear in Phys. Rev.

Auger Recombination in Semiconductor Quantum Wells

The principal mechanisms of Auger recombination of nonequilibrium carriers in semiconductor heterostructures with quantum wells are investigated. It is shown for the first time that there exist three fundamentally different Auger recombination mechanisms of (i) thresholdless, (ii) quasi-threshold, and (iii) threshold types. The rate of the thresholdless Auger process depends on temperature only slightly. The rate of the quasi-threshold Auger process depends on temperature exponentially. However, its threshold energy essentially varies with quantum well width and is close to zero for narrow quantum wells. It is shown that the thresholdless and the quasi-threshold Auger processes dominate in narrow quantum wells, while the threshold and the quasi-threshold processes prevail in wide quantum wells. The limiting case of a three-dimensional (3D)Auger process is reached for infinitely wide quantum wells. The critical quantum well width is found at which the quasi-threshold and threshold Auger processes merge into a single 3D Auger process. Also studied is phonon-assisted Auger recombination in quantum wells. It is shown that for narrow quantum wells the act of phonon emission becomes resonant, which in turn increases substantially the coefficient of phonon-assisted Auger recombination. Conditions are found under which the direct Auger process dominates over the phonon-assisted Auger recombination at various temperatures and quantum well widths.Comment: 38 pages, 7 figure

Quantum Monte Carlo Loop Algorithm for the t-J Model

We propose a generalization of the Quantum Monte Carlo loop algorithm to the t-J model by a mapping to three coupled six-vertex models. The autocorrelation times are reduced by orders of magnitude compared to the conventional local algorithms. The method is completely ergodic and can be formulated directly in continuous time. We introduce improved estimators for simulations with a local sign problem. Some first results of finite temperature simulations are presented for a t-J chain, a frustrated Heisenberg chain, and t-J ladder models.Comment: 22 pages, including 12 figures. RevTex v3.0, uses psf.te

Hemodynamic Effects of Anthrax Toxins in the Rabbit Model and the Cardiac Pathology Induced by Lethal Toxin

Anthrax lethal toxin (LeTx) and edema toxin (EdTx) have been shown to alter hemodynamics in the rodent model, while LeTx primarily is reported to induce extensive tissue pathology. However, the rodent model has limitations when used for comparison to higher organisms such as humans. The rabbit model, on the other hand, has gained recognition as a useful model for studying anthrax infection and its pathophysiological effects. In this study, we assessed the hemodynamic effects of lethal toxin (LeTx) and edema toxin (EdTx) in the rabbit model using physiologically relevant amounts of the toxins. Moreover, we further examine the pathological effects of LeTx on cardiac tissue. We intravenously injected Dutch-belted rabbits with either low-dose and high-dose recombinant LeTx or a single dose of EdTx. The animals’ heart rate and mean arterial pressure were continuously monitored via telemetry until either 48 or 72 h post-challenge. Additional animals challenged with LeTx were used for cardiac troponin I (cTnI) quantitation, cardiac histopathology, and echocardiography. LeTx depressed heart rate at the lower dose and mean arterial pressure (MAP) at the higher dose. EdTx, on the other hand, temporarily intensified heart rate while lowering MAP. Both doses of LeTx caused cardiac pathology with the higher dose having a more profound effect. Lastly, left-ventricular dilation due to LeTx was not apparent at the given time-points. Our study demonstrates the hemodynamic effects of anthrax toxins, as well as the pathological effects of LeTx on the heart in the rabbit model, and it provides further evidence for the toxins’ direct impact on the heart

On the Coupling Time of the Heat-Bath Process for the Fortuin–Kasteleyn Random–Cluster Model

We consider the coupling from the past implementation of the random-cluster heat-bath process, and study its random running time, or coupling time. We focus on hypercubic lattices embedded on tori, in dimensions one to three, with cluster fugacity at least one. We make a number of conjectures regarding the asymptotic behaviour of the coupling time, motivated by rigorous results in one dimension and Monte Carlo simulations in dimensions two and three. Amongst our findings, we observe that, for generic parameter values, the distribution of the appropriately standardized coupling time converges to a Gumbel distribution, and that the standard deviation of the coupling time is asymptotic to an explicit universal constant multiple of the relaxation time. Perhaps surprisingly, we observe these results to hold both off criticality, where the coupling time closely mimics the coupon collector's problem, and also at the critical point, provided the cluster fugacity is below the value at which the transition becomes discontinuous. Finally, we consider analogous questions for the single-spin Ising heat-bath process

The reference frame for encoding and retention of motion depends on stimulus set size

YesThe goal of this study was to investigate the reference frames used in perceptual encoding and storage of visual motion information. In our experiments, observers viewed multiple moving objects and reported the direction of motion of a randomly selected item. Using a vector-decomposition technique, we computed performance during smooth pursuit with respect to a spatiotopic (nonretinotopic) and to a retinotopic component and compared them with performance during fixation, which served as the baseline. For the stimulus encoding stage, which precedes memory, we found that the reference frame depends on the stimulus set size. For a single moving target, the spatiotopic reference frame had the most significant contribution with some additional contribution from the retinotopic reference frame. When the number of items increased (Set Sizes 3 to 7), the spatiotopic reference frame was able to account for the performance. Finally, when the number of items became larger than 7, the distinction between reference frames vanished. We interpret this finding as a switch to a more abstract nonmetric encoding of motion direction. We found that the retinotopic reference frame was not used in memory. Taken together with other studies, our results suggest that, whereas a retinotopic reference frame may be employed for controlling eye movements, perception and memory use primarily nonretinotopic reference frames. Furthermore, the use of nonretinotopic reference frames appears to be capacity limited. In the case of complex stimuli, the visual system may use perceptual grouping in order to simplify the complexity of stimuli or resort to a nonmetric abstract coding of motion information

Small Heat Shock Proteins Potentiate Amyloid Dissolution by Protein Disaggregases from Yeast and Humans

The authors define how small heat-shock proteins synergize to regulate the assembly and disassembly of a beneficial prion, and then they exploit this knowledge to identify the human amyloid depolymerase

The Mammalian Disaggregase Machinery: Hsp110 Synergizes with Hsp70 and Hsp40 to Catalyze Protein Disaggregation and Reactivation in a Cell-Free System

Bacteria, fungi, protozoa, chromista and plants all harbor homologues of Hsp104, a AAA+ ATPase that collaborates with Hsp70 and Hsp40 to promote protein disaggregation and reactivation. Curiously, however, metazoa do not possess an Hsp104 homologue. Thus, whether animal cells renature large protein aggregates has long remained unclear. Here, it is established that mammalian cytosol prepared from different sources possesses a potent, ATP-dependent protein disaggregase and reactivation activity, which can be accelerated and stimulated by Hsp104. This activity did not require the AAA+ ATPase, p97. Rather, mammalian Hsp110 (Apg-2), Hsp70 (Hsc70 or Hsp70) and Hsp40 (Hdj1) were necessary and sufficient to slowly dissolve large disordered aggregates and recover natively folded protein. This slow disaggregase activity was conserved to yeast Hsp110 (Sse1), Hsp70 (Ssa1) and Hsp40 (Sis1 or Ydj1). Hsp110 must engage substrate, engage Hsp70, promote nucleotide exchange on Hsp70, and hydrolyze ATP to promote disaggregation of disordered aggregates. Similarly, Hsp70 must engage substrate and Hsp110, and hydrolyze ATP for protein disaggregation. Hsp40 must harbor a functional J domain to promote protein disaggregation, but the J domain alone is insufficient. Optimal disaggregase activity is achieved when the Hsp40 can stimulate the ATPase activity of Hsp110 and Hsp70. Finally, Hsp110, Hsp70 and Hsp40 fail to rapidly remodel amyloid forms of the yeast prion protein, Sup35, or the Parkinson's disease protein, alpha-synuclein. However, Hsp110, Hsp70 and Hsp40 enhanced the activity of Hsp104 against these amyloid substrates. Taken together, these findings suggest that Hsp110 fulfils a subset of Hsp104 activities in mammals. Moreover, they suggest that Hsp104 can collaborate with the mammalian disaggregase machinery to rapidly remodel amyloid conformers