Comparative Monte Carlo Efficiency by Monte Carlo Analysis

We propose a modified power method for computing the subdominant eigenvalue $\lambda_2$ of a matrix or continuous operator. Here we focus on defining simple Monte Carlo methods for its application. The methods presented use random walkers of mixed signs to represent the subdominant eigenfuction. Accordingly, the methods must cancel these signs properly in order to sample this eigenfunction faithfully. We present a simple procedure to solve this sign problem and then test our Monte Carlo methods by computing the $\lambda_2$ of various Markov chain transition matrices. We first computed ${\lambda_2}$ for several one and two dimensional Ising models, which have a discrete phase space, and compared the relative efficiencies of the Metropolis and heat-bath algorithms as a function of temperature and applied magnetic field. Next, we computed $\lambda_2$ for a model of an interacting gas trapped by a harmonic potential, which has a mutidimensional continuous phase space, and studied the efficiency of the Metropolis algorithm as a function of temperature and the maximum allowable step size $\Delta$. Based on the $\lambda_2$ criterion, we found for the Ising models that small lattices appear to give an adequate picture of comparative efficiency and that the heat-bath algorithm is more efficient than the Metropolis algorithm only at low temperatures where both algorithms are inefficient. For the harmonic trap problem, we found that the traditional rule-of-thumb of adjusting $\Delta$ so the Metropolis acceptance rate is around 50% range is often sub-optimal. In general, as a function of temperature or $\Delta$, $\lambda_2$ for this model displayed trends defining optimal efficiency that the acceptance ratio does not. The cases studied also suggested that Monte Carlo simulations for a continuum model are likely more efficient than those for a discretized version of the model.Comment: 23 pages, 8 figure

Unfolding cross-linkers as rheology regulators in F-actin networks

We report on the nonlinear mechanical properties of a statistically homogeneous, isotropic semiflexible network cross-linked by polymers containing numerous small unfolding domains, such as the ubiquitous F-actin cross-linker Filamin. We show that the inclusion of such proteins has a dramatic effect on the large strain behavior of the network. Beyond a strain threshold, which depends on network density, the unfolding of protein domains leads to bulk shear softening. Past this critical strain, the network spontaneously organizes itself so that an appreciable fraction of the Filamin cross-linkers are at the threshold of domain unfolding. We discuss via a simple mean-field model the cause of this network organization and suggest that it may be the source of power-law relaxation observed in in vitro and in intracellular microrheology experiments. We present data which fully justifies our model for a simplified network architecture.Comment: 11 pages, 4 figures. to appear in Physical Review

Indeterminacy of Spatiotemporal Cardiac Alternans

Cardiac alternans, a beat-to-beat alternation in action potential duration (at the cellular level) or in ECG morphology (at the whole heart level), is a marker of ventricular fibrillation, a fatal heart rhythm that kills hundreds of thousands of people in the US each year. Investigating cardiac alternans may lead to a better understanding of the mechanisms of cardiac arrhythmias and eventually better algorithms for the prediction and prevention of such dreadful diseases. In paced cardiac tissue, alternans develops under increasingly shorter pacing period. Existing experimental and theoretical studies adopt the assumption that alternans in homogeneous cardiac tissue is exclusively determined by the pacing period. In contrast, we find that, when calcium-driven alternans develops in cardiac fibers, it may take different spatiotemporal patterns depending on the pacing history. Because there coexist multiple alternans solutions for a given pacing period, the alternans pattern on a fiber becomes unpredictable. Using numerical simulation and theoretical analysis, we show that the coexistence of multiple alternans patterns is induced by the interaction between electrotonic coupling and an instability in calcium cycling.Comment: 20 pages, 10 figures, to be published in Phys. Rev.

Child health in Syria: recognising the lasting effects of warfare on health

The war in Syria, now in its fourth year, is one of the bloodiest in recent times. The legacy of war includes damage to the health of children that can last for decades and affect future generations. In this article we discuss the effects of the war on Syria's children, highlighting the less documented longer-term effects. In addition to their present suffering, these children, and their own children, are likely to face further challenges as a result of the current conflict. This is essential to understand both for effective interventions and for ethical reasons

Are superflares on solar analogues caused by extra-solar planets?

Stellar flares with ${10^2-10^7}$ times more energy than the largest solar flare have been detected from 9 normal F and G main sequence stars (Schaefer, King & Deliyannis 1999). These superflares have durations of hours to days and are visible from at least x-ray to optical frequencies. The absence of world-spanning aurorae in historical records and of anomalous extinctions in the geological record indicate that our Sun likely does not suffer superflares. In seeking to explain this new phenomenon, we are struck by its similarity to large stellar flares on RS Canum Venaticorum binary systems, which are caused by magnetic reconnection events associated with the tangling of magnetic fields between the two stars. The superflare stars are certainly not of this class, although we propose a similar flare mechanism. That is, superflares are caused by magnetic reconnection between fields of the primary star and a close-in Jovian planet. Thus, by only invoking known planetary properties and reconnection scenarios, we can explain the energies, durations, and spectra of superflares, as well as explain why our Sun does not have such events.Comment: 13 pages, Accepted for publication in Ap

Multidimensional Geriatric Assessment: Back to the Future Multidimensional Preventive Home Visit Programs for Community-Dwelling Older Adults: A Systematic Review and Meta-Analysis of Randomized Controlled Trials

Background. Multidimensional preventive home visit programs aim at maintaining health and autonomy of older adults and preventing disability and subsequent nursing home admission, but results of randomized controlled trials (RCTs) have been inconsistent. Our objective was to systematically review RCTs examining the effect of home visit programs on mortality, nursing home admissions, and functional status decline. Methods. Data sources were MEDLINE, EMBASE, Cochrane CENTRAL database, and references. Studies were reviewed to identify RCTs that compared outcome data of older participants in preventive home visit programs with control group outcome data. Publications reporting 21 trials were included. Data on study population, intervention characteristics, outcomes, and trial quality were double-extracted. We conducted random effects meta-analyses. Results. Pooled effects estimates revealed statistically nonsignificant favorable, and heterogeneous effects on mortality (odds ratio [OR] 0.92, 95% confidence interval [CI], 0.80-1.05), functional status decline (OR 0.89, 95% CI, 0.77-1.03), and nursing home admission (OR 0.86, 95% CI, 0.68-1.10). A beneficial effect on mortality was seen in younger study populations (OR 0.74, 95% CI, 0.58-0.94) but not in older populations (OR 1.14, 95% CI, 0.90-1.43). Functional decline was reduced in programs including a clinical examination in the initial assessment (OR 0.64, 95% CI, 0.48-0.87) but not in other trials (OR 1.00, 95% CI, 0.88-1.14). There was no single factor explaining the heterogenous effects of trials on nursing home admissions. Conclusion. Multidimensional preventive home visits have the potential to reduce disability burden among older adults when based on multidimensional assessment with clinical examination. Effects on nursing home admissions are heterogeneous and likely depend on multiple factors including population factors, program characteristics, and health care settin