Acceleration Schemes for Ab-Initio Molecular Dynamics and Electronic Structure Calculations

We study the convergence and the stability of fictitious dynamical methods for electrons. First, we show that a particular damped second-order dynamics has a much faster rate of convergence to the ground-state than first-order steepest descent algorithms while retaining their numerical cost per time step. Our damped dynamics has efficiency comparable to that of conjugate gradient methods in typical electronic minimization problems. Then, we analyse the factors that limit the size of the integration time step in approaches based on plane-wave expansions. The maximum allowed time step is dictated by the highest frequency components of the fictitious electronic dynamics. These can result either from the large wavevector components of the kinetic energy or from the small wavevector components of the Coulomb potential giving rise to the so called {\it charge sloshing} problem. We show how to eliminate large wavevector instabilities by adopting a preconditioning scheme that is implemented here for the first-time in the context of Car-Parrinello ab-initio molecular dynamics simulations of the ionic motion. We also show how to solve the charge-sloshing problem when this is present. We substantiate our theoretical analysis with numerical tests on a number of different silicon and carbon systems having both insulating and metallic character.Comment: RevTex, 9 figures available upon request, to appear in Phys. Rev.

Diagnosing human African trypanosomiasis in Angola using a card agglutination test: observational study of active and passive case finding strategies.

OBJECTIVE: To assess the operational feasibility of detecting human African trypanosomiasis by active and passive case finding using the card agglutination test with serial dilution of serum to guide treatment. SETTING: Trypanosomiasis control programme in the Negage focus, northern Angola, during a period of civil war. DESIGN: Observational study. PARTICIPANTS: 359 patients presenting themselves to health centres with symptoms (passive case finding) and 14,446 people actively screened in villages. MAIN OUTCOME MEASURES: Whole blood and serological tests at different dilutions using the card agglutination test, and detection of parasites by microscopy. RESULTS: Active case finding identified 251 people with a positive card agglutination test result, 10 of whom had confirmed parasites. In those presenting for investigation 34 of 51 with a positive card agglutination test result at the dilution of 1:8 or more used to guide treatment had parasites in blood, lymph node fluid, or cerebrospinal fluid, compared with 10 of 76 in those detected by active case finding: positive predictive values of 67% for passive case detection and 13% for active case detection. Only at a cut-off dilution more than 1:32 was the positive predictive value in active case detection reasonable (46%) and at this dilution 40% of microscopically proved cases were missed. CONCLUSIONS: The card agglutination test is useful for initial screening in active detection of cases with human African trypanosomiasis but, given the toxicity of the drugs, serology using the card agglutination test should be not used alone to guide treatment after active case finding. A second confirmatory test is needed

Dynamic Structure Factor of Liquid and Amorphous Ge From Ab Initio Simulations

We calculate the dynamic structure factor S(k,omega) of liquid Ge (l-Ge) at temperature T = 1250 K, and of amorphous Ge (a-Ge) at T = 300 K, using ab initio molecular dynamics. The electronic energy is computed using density-functional theory, primarily in the generalized gradient approximation, together with a plane wave representation of the wave functions and ultra-soft pseudopotentials. We use a 64-atom cell with periodic boundary conditions, and calculate averages over runs of up to 16 ps. The calculated liquid S(k,omega) agrees qualitatively with that obtained by Hosokawa et al, using inelastic X-ray scattering. In a-Ge, we find that the calculated S(k,omega) is in qualitative agreement with that obtained experimentally by Maley et al. Our results suggest that the ab initio approach is sufficient to allow approximate calculations of S(k,omega) in both liquid and amorphous materials.Comment: 31 pages and 8 figures. Accepted for Phys. Rev.

The Hessian Estimation Evolution Strategy

We present a novel black box optimization algorithm called Hessian Estimation Evolution Strategy. The algorithm updates the covariance matrix of its sampling distribution by directly estimating the curvature of the objective function. This algorithm design is targeted at twice continuously differentiable problems. For this, we extend the cumulative step-size adaptation algorithm of the CMA-ES to mirrored sampling. We demonstrate that our approach to covariance matrix adaptation is efficient by evaluation it on the BBOB/COCO testbed. We also show that the algorithm is surprisingly robust when its core assumption of a twice continuously differentiable objective function is violated. The approach yields a new evolution strategy with competitive performance, and at the same time it also offers an interesting alternative to the usual covariance matrix update mechanism

Thermodynamic Behavior of a Model Covalent Material Described by the Environment-Dependent Interatomic Potential

Using molecular dynamics simulations we study the thermodynamic behavior of a single-component covalent material described by the recently proposed Environment-Dependent Interatomic Potential (EDIP). The parameterization of EDIP for silicon exhibits a range of unusual properties typically found in more complex materials, such as the existence of two structurally distinct disordered phases, a density decrease upon melting of the low-temperature amorphous phase, and negative thermal expansion coefficients for both the crystal (at high temperatures) and the amorphous phase (at all temperatures). Structural differences between the two disordered phases also lead to a first-order transition between them, which suggests the existence of a second critical point, as is believed to exist for amorphous forms of frozen water. For EDIP-Si, however, the unusual behavior is associated not only with the open nature of tetrahedral bonding but also with a competition between four-fold (covalent) and five-fold (metallic) coordination. The unusual behavior of the model and its unique ability to simulation the liquid/amorphous transition on molecular-dynamics time scales make it a suitable prototype for fundamental studies of anomalous thermodynamics in disordeered systems.Comment: 48 pages (double-spaced), 13 figure

Self-organized stable pacemakers near the onset of birhythmicity

General amplitude equations for reaction-diffusion systems near to the soft onset of birhythmicity described by a supercritical pitchfork-Hopf bifurcation are derived. Using these equations and applying singular perturbation theory, we show that stable autonomous pacemakers represent a generic kind of spatiotemporal patterns in such systems. This is verified by numerical simulations, which also show the existence of breathing and swinging pacemaker solutions. The drift of self-organized pacemakers in media with spatial parameter gradients is analytically and numerically investigated.Comment: 4 pages, 4 figure

Cognitive Computation sans Representation

The Computational Theory of Mind (CTM) holds that cognitive processes are essentially computational, and hence computation provides the scientific key to explaining mentality. The Representational Theory of Mind (RTM) holds that representational content is the key feature in distinguishing mental from non-mental systems. I argue that there is a deep incompatibility between these two theoretical frameworks, and that the acceptance of CTM provides strong grounds for rejecting RTM. The focal point of the incompatibility is the fact that representational content is extrinsic to formal procedures as such, and the intended interpretation of syntax makes no difference to the execution of an algorithm. So the unique 'content' postulated by RTM is superfluous to the formal procedures of CTM. And once these procedures are implemented in a physical mechanism, it is exclusively the causal properties of the physical mechanism that are responsible for all aspects of the system's behaviour. So once again, postulated content is rendered superfluous. To the extent that semantic content may appear to play a role in behaviour, it must be syntactically encoded within the system, and just as in a standard computational artefact, so too with the human mind/brain - it's pure syntax all the way down to the level of physical implementation. Hence 'content' is at most a convenient meta-level gloss, projected from the outside by human theorists, which itself can play no role in cognitive processing

Factors associated with Anaplasma spp. seroprevalence among dogs in the United States

Background Dogs in the United States are hosts to a diverse range of ticks and tick-borne pathogens, including A. phagocytophilum, an important emerging canine and human pathogen. Previously, a Companion Animal Parasite Council (CAPC)-sponsored workshop proposed factors purported to be associated with the infection risk for tick-transmitted pathogens in dogs in the United States, including climate conditions, socioeconomic characteristics, local topography, and vector distribution. Methods Approximately four million test results from routine veterinary diagnostic tests from 2011–2013, which were collected on a county level across the contiguous United States, are statistically analyzed with the proposed factors via logistic regression and generalized estimating equations. Spatial prevalence maps of baseline Anaplasma spp. prevalence are constructed from Kriging and head-banging smoothing methods. Results All of the examined factors, with the exception of surface water coverage, were significantly associated with Anaplasma spp. prevalence. Overall, Anaplasma spp. prevalence increases with increasing precipitation and forestation coverage and decreases with increasing temperature, population density, relative humidity, and elevation. Interestingly, socioeconomic status and deer/vehicle collisions were positively and negatively correlated with canine Anaplasma seroprevalence, respectively. A spatial map of the canine Anaplasma hazard is an auxiliary product of the analysis. Anaplasma spp. prevalence is highest in New England and the Upper Midwest. Conclusions The results from the two posited statistical models (one that contains an endemic areas assumption and one that does not) are in general agreement, with the major difference being that the endemic areas model estimates a larger prevalence in Western Texas, New Mexico, and Colorado. As A. phagocytophilum is zoonotic, the results of this analysis could also help predict areas of high risk for human exposure to this pathogen

Pumpless Extracorporeal Hemadsorption Technique (pEHAT) : A Proof-of-Concept Animal Study

Background: Extracorporeal hemadsorption eliminates proinflammatory mediators in critically ill patients with hyperinflammation. The use of a pumpless extracorporeal hemadsorption technique allows its early usage prior to organ failure and the need for an additional medical device. In our animal model, we investigated the feasibility of pumpless extracorporeal hemadsorption over a wide range of mean arterial pressures (MAP). Methods: An arteriovenous shunt between the femoral artery and femoral vein was established in eight pigs. The hemadsorption devices were inserted into the shunt circulation; four pigs received CytoSorb® and four Oxiris® hemadsorbers. Extracorporeal blood flow was measured in a range between mean arterial pressures of 45–85 mmHg. Mean arterial pressures were preset using intravenous infusions of noradrenaline, urapidil, or increased sedatives. Results: Extracorporeal blood flows remained well above the minimum flows recommended by the manufacturers throughout all MAP steps for both devices. Linear regression resulted in CytoSorb® blood flow [mL/min] = 4.226 × MAP [mmHg] − 3.496 (R-square 0.8133) and Oxiris® blood flow [mL/min] = 3.267 × MAP [mmHg] + 57.63 (R-square 0.8708), respectively. Conclusion: Arteriovenous pumpless extracorporeal hemadsorption resulted in sufficient blood flows through both the CytoSorb® and Oxiris® devices over a wide range of mean arterial blood pressures and is likely an intriguing therapeutic option in the early phase of septic shock or hyperinflammatory syndromes

In silico selection of RNA aptamers

In vitro selection of RNA aptamers that bind to a specific ligand usually begins with a random pool of RNA sequences. We propose a computational approach for designing a starting pool of RNA sequences for the selection of RNA aptamers for specific analyte binding. Our approach consists of three steps: (i) selection of RNA sequences based on their secondary structure, (ii) generating a library of three-dimensional (3D) structures of RNA molecules and (iii) high-throughput virtual screening of this library to select aptamers with binding affinity to a desired small molecule. We developed a set of criteria that allows one to select a sequence with potential binding affinity from a pool of random sequences and developed a protocol for RNA 3D structure prediction. As verification, we tested the performance of in silico selection on a set of six known aptamer–ligand complexes. The structures of the native sequences for the ligands in the testing set were among the top 5% of the selected structures. The proposed approach reduces the RNA sequences search space by four to five orders of magnitude—significantly accelerating the experimental screening and selection of high-affinity aptamers