Order-of-magnitude speedup for steady states and traveling waves via Stokes preconditioning in Channelflow and Openpipeflow

Steady states and traveling waves play a fundamental role in understanding hydrodynamic problems. Even when unstable, these states provide the bifurcation-theoretic explanation for the origin of the observed states. In turbulent wall-bounded shear flows, these states have been hypothesized to be saddle points organizing the trajectories within a chaotic attractor. These states must be computed with Newton's method or one of its generalizations, since time-integration cannot converge to unstable equilibria. The bottleneck is the solution of linear systems involving the Jacobian of the Navier-Stokes or Boussinesq equations. Originally such computations were carried out by constructing and directly inverting the Jacobian, but this is unfeasible for the matrices arising from three-dimensional hydrodynamic configurations in large domains. A popular method is to seek states that are invariant under numerical time integration. Surprisingly, equilibria may also be found by seeking flows that are invariant under a single very large Backwards-Euler Forwards-Euler timestep. We show that this method, called Stokes preconditioning, is 10 to 50 times faster at computing steady states in plane Couette flow and traveling waves in pipe flow. Moreover, it can be carried out using Channelflow (by Gibson) and Openpipeflow (by Willis) without any changes to these popular spectral codes. We explain the convergence rate as a function of the integration period and Reynolds number by computing the full spectra of the operators corresponding to the Jacobians of both methods.Comment: in Computational Modelling of Bifurcations and Instabilities in Fluid Dynamics, ed. Alexander Gelfgat (Springer, 2018

Order by disorder and spiral spin liquid in frustrated diamond lattice antiferromagnets

Frustration refers to competition between different interactions that cannot be simultaneously satisfied, a familiar feature in many magnetic solids. Strong frustration results in highly degenerate ground states, and a large suppression of ordering by fluctuations. Key challenges in frustrated magnetism are characterizing the fluctuating spin-liquid regime and determining the mechanism of eventual order at lower temperature. Here, we study a model of a diamond lattice antiferromagnet appropriate for numerous spinel materials. With sufficiently strong frustration a massive ground state degeneracy develops amongst spirals whose propagation wavevectors reside on a continuous two-dimensional ``spiral surface'' in momentum space. We argue that an important ordering mechanism is entropic splitting of the degenerate ground states, an elusive phenomena called order-by-disorder. A broad ``spiral spin-liquid'' regime emerges at higher temperatures, where the underlying spiral surface can be directly revealed via spin correlations. We discuss the agreement between these predictions and the well characterized spinel MnSc2S4

Coxiella burnetii Phagocytosis Is Regulated by GTPases of the Rho Family and the RhoA Effectors mDia1 and ROCK

The GTPases belonging to the Rho family control the actin cytoskeleton rearrangements needed for particle internalization during phagocytosis. ROCK and mDia1 are downstream effectors of RhoA, a GTPase involved in that process. Coxiella burnetii, the etiologic agent of Q fever, is internalized by the host´s cells in an actin-dependent manner. Nevertheless, the molecular mechanism involved in this process has been poorly characterized. This work analyzes the role of different GTPases of the Rho family and some downstream effectors in the internalization of C. burnetii by phagocytic and non-phagocytic cells. The internalization of C. burnetii into HeLa and RAW cells was significantly inhibited when the cells were treated with Clostridium difficile Toxin B which irreversibly inactivates members of the Rho family. In addition, the internalization was reduced in HeLa cells that overexpressed the dominant negative mutants of RhoA, Rac1 or Cdc42 or that were knocked down for the Rho GTPases. The pharmacological inhibition or the knocking down of ROCK diminished bacterium internalization. Moreover, C. burnetii was less efficiently internalized in HeLa cells overexpressing mDia1-N1, a dominant negative mutant of mDia1, while the overexpression of the constitutively active mutant mDia1-ΔN3 increased bacteria uptake. Interestingly, when HeLa and RAW cells were infected, RhoA, Rac1 and mDia1 were recruited to membrane cell fractions. Our results suggest that the GTPases of the Rho family play an important role in C. burnetii phagocytosis in both HeLa and RAW cells. Additionally, we present evidence that ROCK and mDia1, which are downstream effectors of RhoA, are involved in that processFil: Salinas Ojeda, Romina Paola. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Mendoza. Instituto de Histología y Embriología de Mendoza Dr. Mario H. Burgos. Universidad Nacional de Cuyo. Facultad de Cienicas Médicas. Instituto de Histología y Embriología de Mendoza Dr. Mario H. Burgos; ArgentinaFil: Ortiz Flores, Rodolfo Matias. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Mendoza. Instituto de Histología y Embriología de Mendoza Dr. Mario H. Burgos. Universidad Nacional de Cuyo. Facultad de Cienicas Médicas. Instituto de Histología y Embriología de Mendoza Dr. Mario H. Burgos; ArgentinaFil: Distel, Jesús Sebastián. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Mendoza. Instituto de Histología y Embriología de Mendoza Dr. Mario H. Burgos. Universidad Nacional de Cuyo. Facultad de Cienicas Médicas. Instituto de Histología y Embriología de Mendoza Dr. Mario H. Burgos; ArgentinaFil: Aguilera, Milton Osmar. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Mendoza. Instituto de Histología y Embriología de Mendoza Dr. Mario H. Burgos. Universidad Nacional de Cuyo. Facultad de Cienicas Médicas. Instituto de Histología y Embriología de Mendoza Dr. Mario H. Burgos; ArgentinaFil: Colombo, Maria Isabel. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Mendoza. Instituto de Histología y Embriología de Mendoza Dr. Mario H. Burgos. Universidad Nacional de Cuyo. Facultad de Cienicas Médicas. Instituto de Histología y Embriología de Mendoza Dr. Mario H. Burgos; ArgentinaFil: Beron, Walter. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Mendoza. Instituto de Histología y Embriología de Mendoza Dr. Mario H. Burgos. Universidad Nacional de Cuyo. Facultad de Cienicas Médicas. Instituto de Histología y Embriología de Mendoza Dr. Mario H. Burgos; Argentin

DNA content of a functioning chicken kinetochore

© The Author(s) 2014. In order to understand the three-dimensional structure of the functional kinetochore in vertebrates, we require a complete list and stoichiometry for the protein components of the kinetochore, which can be provided by genetic and proteomic experiments. We also need to know how the chromatin-containing CENP-A, which makes up the structural foundation for the kinetochore, is folded, and how much of that DNA is involved in assembling the kinetochore. In this MS, we demonstrate that functioning metaphase kinetochores in chicken DT40 cells contain roughly 50 kb of DNA, an amount that corresponds extremely closely to the length of chromosomal DNA associated with CENP-A in ChIP-seq experiments. Thus, during kinetochore assembly, CENP-A chromatin is compacted into the inner kinetochore plate without including significant amounts of flanking pericentromeric heterochromatin. © 2014 The Author(s).Wellcome Trust [grant number 073915]; Wellcome Trust Centre for Cell Biology (core grant numbers 077707 and 092076); Darwin Trust of Edinburg

Using self-organizing maps to investigate environmental factors regulating colony size and breeding success of the White Stork (Ciconia ciconia)

We studied variations in the size of breeding colonies and in breeding performance of White Storks Ciconia ciconia in 2006–2008 in north-east Algeria. Each colony site was characterized using 12 environmental variables describing the physical environment, land-cover categories, and human activities, and by three demographic parameters: the number of breeding pairs, the number of pairs with chicks, and the number of fledged chicks per pair. Generalized linear mixed models and the self-organizing map algorithm (SOM, neural network) were used to investigate effects of biotic, abiotic, and anthropogenic factors on demographic parameters and on their relationships. Numbers of breeding pairs and of pairs with chicks were affected by the same environmental factors, mainly anthropogenic, which differed from those affecting the number of fledged chicks per pair. Numbers of fledged chicks per pair was not affected by colony size or by the number of nests with chicks. The categorization of the environmental variables into natural and anthropogenic, in connection with demographic parameters, was relevant to detect factors explaining variation in colony size and breeding parameters. The SOM proved a relevant tool to help determine actual dynamics in White Stork colonies, and thus to support effective conservation decisions at a regional scale

Search for new phenomena in final states with an energetic jet and large missing transverse momentum in pp collisions at √ s = 8 TeV with the ATLAS detector

Results of a search for new phenomena in final states with an energetic jet and large missing transverse momentum are reported. The search uses 20.3 fb−1 of √ s = 8 TeV data collected in 2012 with the ATLAS detector at the LHC. Events are required to have at least one jet with pT > 120 GeV and no leptons. Nine signal regions are considered with increasing missing transverse momentum requirements between Emiss T > 150 GeV and Emiss T > 700 GeV. Good agreement is observed between the number of events in data and Standard Model expectations. The results are translated into exclusion limits on models with either large extra spatial dimensions, pair production of weakly interacting dark matter candidates, or production of very light gravitinos in a gauge-mediated supersymmetric model. In addition, limits on the production of an invisibly decaying Higgs-like boson leading to similar topologies in the final state are presente

Measurement of the inclusive and dijet cross-sections of b-jets in pp collisions at sqrt(s) = 7 TeV with the ATLAS detector

The inclusive and dijet production cross-sections have been measured for jets containing b-hadrons (b-jets) in proton-proton collisions at a centre-of-mass energy of sqrt(s) = 7 TeV, using the ATLAS detector at the LHC. The measurements use data corresponding to an integrated luminosity of 34 pb^-1. The b-jets are identified using either a lifetime-based method, where secondary decay vertices of b-hadrons in jets are reconstructed using information from the tracking detectors, or a muon-based method where the presence of a muon is used to identify semileptonic decays of b-hadrons inside jets. The inclusive b-jet cross-section is measured as a function of transverse momentum in the range 20 < pT < 400 GeV and rapidity in the range |y| < 2.1. The bbbar-dijet cross-section is measured as a function of the dijet invariant mass in the range 110 < m_jj < 760 GeV, the azimuthal angle difference between the two jets and the angular variable chi in two dijet mass regions. The results are compared with next-to-leading-order QCD predictions. Good agreement is observed between the measured cross-sections and the predictions obtained using POWHEG + Pythia. MC@NLO + Herwig shows good agreement with the measured bbbar-dijet cross-section. However, it does not reproduce the measured inclusive cross-section well, particularly for central b-jets with large transverse momenta.Comment: 10 pages plus author list (21 pages total), 8 figures, 1 table, final version published in European Physical Journal

Environmental and Parental Influences on Offspring Health and Growth in Great Tits (Parus major)

PMCID: PMC3728352This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited

Search for direct pair production of the top squark in all-hadronic final states in proton-proton collisions at s√=8 TeV with the ATLAS detector

The results of a search for direct pair production of the scalar partner to the top quark using an integrated luminosity of 20.1fb−1 of proton–proton collision data at √s = 8 TeV recorded with the ATLAS detector at the LHC are reported. The top squark is assumed to decay via t˜→tχ˜01 or t˜→ bχ˜±1 →bW(∗)χ˜01 , where χ˜01 (χ˜±1 ) denotes the lightest neutralino (chargino) in supersymmetric models. The search targets a fully-hadronic final state in events with four or more jets and large missing transverse momentum. No significant excess over the Standard Model background prediction is observed, and exclusion limits are reported in terms of the top squark and neutralino masses and as a function of the branching fraction of t˜ → tχ˜01 . For a branching fraction of 100%, top squark masses in the range 270–645 GeV are excluded for χ˜01 masses below 30 GeV. For a branching fraction of 50% to either t˜ → tχ˜01 or t˜ → bχ˜±1 , and assuming the χ˜±1 mass to be twice the χ˜01 mass, top squark masses in the range 250–550 GeV are excluded for χ˜01 masses below 60 GeV

A simulation study on the effects of neuronal ensemble properties on decoding algorithms for intracortical brain-machine interfaces

Background: Intracortical brain-machine interfaces (BMIs) harness movement information by sensing neuronal activities using chronic microelectrode implants to restore lost functions to patients with paralysis. However, neuronal signals often vary over time, even within a day, forcing one to rebuild a BMI every time they operate it. The term &quot;rebuild&quot; means overall procedures for operating a BMI, such as decoder selection, decoder training, and decoder testing. It gives rise to a practical issue of what decoder should be built for a given neuronal ensemble. This study aims to address it by exploring how decoders&apos; performance varies with the neuronal properties. To extensively explore a range of neuronal properties, we conduct a simulation study. Methods: Focusing on movement direction, we examine several basic neuronal properties, including the signal-to-noise ratio of neurons, the proportion of well-tuned neurons, the uniformity of their preferred directions (PDs), and the non-stationarity of PDs. We investigate the performance of three popular BMI decoders: Kalman filter, optimal linear estimator, and population vector algorithm. Results: Our simulation results showed that decoding performance of all the decoders was affected more by the proportion of well-tuned neurons that their uniformity. Conclusions: Our study suggests a simulated scenario of how to choose a decoder for intracortical BMIs in various neuronal conditions