Modulating neuronal activity produces specific and long-lasting changes in numerical competence

Around 20% of the population exhibits moderate to severe numerical disabilities [1-3], and a further percentage loses its numerical competence during the lifespan as a result of stroke or degenerative diseases [4]. In this work, we investigated the feasibility of using noninvasive stimulation to the parietal lobe during numerical learning to selectively improve numerical abilities. We used transcranial direct current stimulation (TDCS), a method that can selectively inhibit or excitate neuronal populations by modulating GABAergic (anodal stimulation) and glutamatergic (cathodal stimulation) activity [5, 6]. We trained subjects for 6 days with artificial numerical symbols, during which we applied concurrent TDCS to the parietal lobes. The polarity of the brain stimulation specifically enhanced or impaired the acquisition of automatic number processing and the mapping of number into space, both important indices of numerical proficiency [7-9]. The improvement was still present 6 months after the training. Control tasks revealed that the effect of brain stimulation was specific to the representation of artificial numerical symbols. The specificity and longevity of TDCS on numerical abilities establishes TDCS as a realistic tool for intervention in cases of atypical numerical development or loss of numerical abilities because of stroke or degenerative illnesses.     

The numerical response: rate of increase and food limitation in herbivores and predators

Two types of numerical response function have evolved since Solomon first introduced the term to generalize features of Lotka-Volterra predator-prey models: (i) the demographic numerical response, which links change in consumer demographic rates to food availability; and (ii) the isocline numerical response, which links consumer abundance per se to food availability. These numerical responses are interchangeable because both recognize negative feedback loops between consumer and food abundance resulting in population regulation. We review how demographic and isocline numerical responses have been used to enhance our understanding of population regulation of kangaroos and possums, and argue that their utility may be increased by explicitly accounting for non-equilibrium dynamics (due to environmental variability and/or biological interactions) and the existence of multiple limiting factors. Interferential numerical response functions may help bridge three major historical dichotomies in population ecology (equilibrium versus non-equilibrium dynamics, extrinsic versus intrinsic regulation and demographic versus isocline numerical responses).     

Relative numerical middle in rhesus monkeys

Animals show vast numerical competence in tasks that require both ordinal and cardinal numerical representations, but few studies have addressed whether animals can identify the numerical middle in a sequence. Two rhesus monkeys (Macaca mulatta) learned to select the middle dot in a horizontal sequence of three dots on a touchscreen. When subsequently presented with longer sequences composed of 5, 7 or 9 items, monkeys transferred the middle rule. Accuracy decreased as the length of the sequence increased. In a second test, we presented monkeys with asymmetrical sequences composed of nine items, where the numerical and spatial middle were distinct and both monkeys selected the numerical middle over the spatial middle. Our results demonstrate that rhesus macaques can extract an abstract numerical rule to bisect a discrete set of items.     

Effects of development and enculturation on number representation in the brain

A striking way in which humans differ from non-human primates is in their ability to represent numerical quantity using abstract symbols and to use these 'mental tools' to perform skills such as exact calculations. How do functional brain circuits for the symbolic representation of numerical magnitude emerge? Do neural representations of numerical magnitude change as a function of development and the learning of mental arithmetic? Current theories suggest that cultural number symbols acquire their meaning by being mapped onto non-symbolic representations of numerical magnitude. This Review provides an evaluation of this contention and proposes hypotheses to guide investigations into the neural mechanisms that constrain the acquisition of cultural representations of numerical magnitude.     

Functional imaging of numerical processing in adults and 4-y-old children

Adult humans, infants, pre-school children, and non-human animals appear to share a system of approximate numerical processing for non-symbolic stimuli such as arrays of dots or sequences of tones. Behavioral studies of adult humans implicate a link between these non-symbolic numerical abilities and symbolic numerical processing (e.g., similar distance effects in accuracy and reaction-time for arrays of dots and Arabic numerals). However, neuroimaging studies have remained inconclusive on the neural basis of this link. The intraparietal sulcus (IPS) is known to respond selectively to symbolic numerical stimuli such as Arabic numerals. Recent studies, however, have arrived at conflicting conclusions regarding the role of the IPS in processing non-symbolic, numerosity arrays in adulthood, and very little is known about the brain basis of numerical processing early in development. Addressing the question of whether there is an early-developing neural basis for abstract numerical processing is essential for understanding the cognitive origins of our uniquely human capacity for math and science. Using functional magnetic resonance imaging (fMRI) at 4-Tesla and an event-related fMRI adaptation paradigm, we found that adults showed a greater IPS response to visual arrays that deviated from standard stimuli in their number of elements, than to stimuli that deviated in local element shape. These results support previous claims that there is a neurophysiological link between non-symbolic and symbolic numerical processing in adulthood. In parallel, we tested 4-y-old children with the same fMRI adaptation paradigm as adults to determine whether the neural locus of non-symbolic numerical activity in adults shows continuity in function over development. We found that the IPS responded to numerical deviants similarly in 4-y-old children and adults. To our knowledge, this is the first evidence that the neural locus of adult numerical cognition takes form early in development, prior to sophisticated symbolic numerical experience. More broadly, this is also, to our knowledge, the first cognitive fMRI study to test healthy children as young as 4 y, providing new insights into the neurophysiology of human cognitive development.     

DNA数进制与酶催化的信息传递

通过对DNA中四种碱基进行排列组合的计算,发现DNA是一个多种的数进制系统;20进制信息约占全部数进制信息的4.515%;每种数进制信息由它的产生方式而表现出不同的调控功能;相同的DNA片段用不同的读取方式,可以得到不同的数进制信息.从能量守恒定理和小分子RNA具有酶催化功能推断,酶的催化机制是一个信息传递的过程.     

An open 3D CFD model for the investigation of flow environments experienced by freshwater fish

Computational fluid dynamics (CFD) provides a powerful numerical tool to simulate and study many of the complex fluid-body interactions experienced by freshwater fish. However, major gaps remain in the application of CFD to study the fluid-body interactions of fish, including the absence of an openly available reference body geometry, the lack of a detailed study on suitable numerical methods and a deficit of available velocity laboratory measurements for model calibration and validation. To address these gaps, we provide a set of numerical models based on the open-source CFD toolkit OpenFOAM. The contributions of this work are two-fold: First, to provide a validated openly available numerical setup using a realistic fish model geometry including laboratory velocity measurements. Second, to determine the best-performing turbulence models and near-wall treatments using Reynolds-Averaged Navier-Stokes (RANS) numerical simulations. Finally, we conclude with a critical evaluation of the effects and trade-offs of resolving or modelling the boundary layer (BL) in numerical studies of fish-shaped bodies.     

我国植被数量分析方法的研究概况和发展趋势

植被数量分析是现代植被研究的重要手段,数量分类和排序是现代植被生态学研究最重要的,也是应用最广泛的生态学技术。数量分析方法在20世纪50年代引入植被生态学研究领域,我国学者在70年代后期开始研究植被的数量分类和排序。本文主要从相关资料、应用研究、新方法研究三个方面论述了我国植被数量分析方法的发展,重点阐述了从20世纪70年代后期以来出现的并且被广泛应用的新方法及其应用研究概况,并在此基础上分析了植被数量分析方法未来的发展趋势。     

Numerical taxonomy: criticisms and critiques: Presidential address to the Systematics Association, November 1970

Some recent criticisms and critiques of numerical taxonomy are reviewed, together with some of its present shortcomings. It is pointed out that most of the problems are equally severe for orthodox taxonomy, and many of them can only be investigated by numerical techniques. The reasons for the general success of numerical methods in bacterial classification are discussed. Besides bringing deeper insights into taxonomy as a whole, numerical taxonomy is entering a new and heuristic phase, which includes potential applications to the study of evolution.     

Pseudo-Lyapunov exponents and predictability of Hodgkin-Huxley neuronal network dynamics

We present a numerical analysis of the dynamics of all-to-all coupled Hodgkin-Huxley (HH) neuronal networks with Poisson spike inputs. It is important to point out that, since the dynamical vector of the system contains discontinuous variables, we propose a so-called pseudo-Lyapunov exponent adapted from the classical definition using only continuous dynamical variables, and apply it in our numerical investigation. The numerical results of the largest Lyapunov exponent using this new definition are consistent with the dynamical regimes of the network. Three typical dynamical regimes—asynchronous, chaotic and synchronous, are found as the synaptic coupling strength increases from weak to strong. We use the pseudo-Lyapunov exponent and the power spectrum analysis of voltage traces to characterize the types of the network behavior. In the nonchaotic (asynchronous or synchronous) dynamical regimes, i.e., the weak or strong coupling limits, the pseudo-Lyapunov exponent is negative and there is a good numerical convergence of the solution in the trajectory-wise sense by using our numerical methods. Consequently, in these regimes the evolution of neuronal networks is reliable. For the chaotic dynamical regime with an intermediate strong coupling, the pseudo-Lyapunov exponent is positive, and there is no numerical convergence of the solution and only statistical quantifications of the numerical results are reliable. Finally, we present numerical evidence that the value of pseudo-Lyapunov exponent coincides with that of the standard Lyapunov exponent for systems we have been able to examine.     

随机时滞Lotka-Volterra模型数值解的收敛性

通常情况下,随机时滞Lotka-Volterra模型没有解析解,因而数值逼近方法是研究其性质的有效工具.本文根据Euler数值方法,利用鞅不等式和Ito公式讨论了一类随机时滞Lotka-Volterra模型数值解的收敛性,给出了数值解收敛于解析解的条件.最后通过数值算例对数值计算方法进行了验证.     

Verification and comparison of four numerical schemes for a 1D viscoelastic blood flow model

A reliable and fast numerical scheme is crucial for the 1D simulation of blood flow in compliant vessels. In this paper, a 1D blood flow model is incorporated with a Kelvin–Voigt viscoelastic arterial wall. This leads to a nonlinear hyperbolic–parabolic system, which is then solved with four numerical schemes, namely: MacCormack, Taylor–Galerkin, monotonic upwind scheme for conservation law and local discontinuous Galerkin. The numerical schemes are tested on a single vessel, a simple bifurcation and a network with 55 arteries. The numerical solutions are checked favorably against analytical, semi-analytical solutions or clinical observations. Among the numerical schemes, comparisons are made in four important aspects: accuracy, ability to capture shock-like phenomena, computational speed and implementation complexity. The suitable conditions for the application of each scheme are discussed.     

An unconditionally stable numerical method for the Luo-Rudy 1 model used in simulations of defibrillation

Numerical simulations of defibrillation using the Bidomain model coupled to a model of membrane kinetics represent a serious numerical challenge. This is because very high voltages close to defibrillation electrodes demand that extreme time step restrictions be placed on standard numerical schemes, e.g. the forward Euler scheme. A common solution to this problem is to modify the cell model by simple if-tests applied to several equations and rate functions. These changes are motivated by numerical problems rather than physiology, and should therefore be avoided whenever possible. The purpose of this paper is to present a numerical scheme that handles the original model without modifications and which is unconditionally stable for the Luo-Rudy phase 1 model. This also shows that the cell model is mathematically well-behaved, even in the presence of very high voltages. Our theoretical results are illustrated by numerical computations.     

Advanced material modelling in numerical simulation of primary acetabular press-fit cup stability

Primary stability of artificial acetabular cups, used for total hip arthroplasty, is required for the subsequent osteointegration and good long-term clinical results of the implant. Although closed-cell polymer foams represent an adequate bone substitute in experimental studies investigating primary stability, correct numerical modelling of this material depends on the parameter selection. Material parameters necessary for crushable foam plasticity behaviour were originated from numerical simulations matched with experimental tests of the polymethacrylimide raw material. Experimental primary stability tests of acetabular press-fit cups consisting of static shell assembly with consecutively pull-out and lever-out testing were subsequently simulated using finite element analysis. Identified and optimised parameters allowed the accurate numerical reproduction of the raw material tests. Correlation between experimental tests and the numerical simulation of primary implant stability depended on the value of interference fit. However, the validated material model provides the opportunity for subsequent parametric numerical studies.     

Advanced material modelling in numerical simulation of primary acetabular press-fit cup stability

Primary stability of artificial acetabular cups, used for total hip arthroplasty, is required for the subsequent osteointegration and good long-term clinical results of the implant. Although closed-cell polymer foams represent an adequate bone substitute in experimental studies investigating primary stability, correct numerical modelling of this material depends on the parameter selection.

Material parameters necessary for crushable foam plasticity behaviour were originated from numerical simulations matched with experimental tests of the polymethacrylimide raw material. Experimental primary stability tests of acetabular press-fit cups consisting of static shell assembly with consecutively pull-out and lever-out testing were subsequently simulated using finite element analysis.

Identified and optimised parameters allowed the accurate numerical reproduction of the raw material tests. Correlation between experimental tests and the numerical simulation of primary implant stability depended on the value of interference fit. However, the validated material model provides the opportunity for subsequent parametric numerical studies.     

A simple and highly accurate numerical differentiation method for sensitivity analysis of large-scale metabolic reaction systems

Numerical differentiation is known to be one of the most difficult numerical calculation methods to obtain reliable calculated values at all times. A simple numerical differentiation method using a combination of finite-difference formulas, derived by approximation of Taylor-series equations, is investigated in order to efficiently perform the sensitivity analysis of large-scale metabolic reaction systems. A result of the application to four basic mathematical functions reveals that the use of the eight-point differentiation formula with a non-dimensionalized stepsize close to 0.01 mostly provides more than 14 digits of accuracy in double precision for the numerical derivatives. Moreover, a result of the application to the modified TCA cycle model indicates that the numerical differentiation method gives the calculated values of steady-state metabolite concentrations within a range of round-off error and also makes it possible to transform the Michaelis-Menten equations into the S-system equations having the kinetic orders whose accuracies are mostly more than 14 significant digits. Because of the simple structure of the numerical differentiation formula and its promising high accuracy, it is evident that the present numerical differentiation method is useful for the analysis of large-scale metabolic reaction systems according to the systematic procedure of BST.     

Benchmark solution for the prediction of temperature distributions during radiofrequency ablation of cardiac tissue

Several studies on radiofrequency (RF) ablation are aimed at accurately predicting tissue temperature distributions by numerical solution of the bioheat equation. This paper describes the development of a solution that can serve as a benchmark for subsequent numerical solutions. The solution was obtained using integral transforms and evaluated using a C program. Temperature profiles were generated at various times and for different convection coefficients. In addition, a numerical model was developed using the same assumptions made in obtaining the benchmark solution. Comparison of surface and axial temperature profiles shows that the two solutions match very closely, cross validating the numerical methods used in evaluating both solutions.     

Distinction of microcytic disorders: comparison of expert, numerical-discriminant, and microcomputer analysis

Presumptive distinction between iron deficiency and heterozygous thalassemia by analysis of the automated blood count and differential continues to be a challenge. We compared two proposed numerical discriminants (MCV2 x MCH, and MCV2 x RDW/100 x Hb) with an analytic microcomputer program (BCDE2 Lea & Febiger). In 7114 subjects, the numerical discriminants and the BCDE2 program correctly identified greater than 90% of thalassemia. In subjects with iron deficiency, the BCDE program was greater than 90% sensitive and specific for positive identification, while the numerical discriminants were less than 75% sensitive and specific at inferential identification. The numerical discriminants, BCDE2, and 17 experts in blood counting were asked to interpret the blood-count data in 7 fully-defined actual cases. The mean experts' score was 5.65 cases correct out of 7. The BCDE program was correct in all cases. The numerical discriminants could not analyze all cases, and both were incorrect in at least one case. We conclude that for the task of analyzing blood counts for microcytic disorders, microcomputer analysis by BCDE outperforms both numerical discriminant functions and analysis by expert hematologists.     

Non-standard numerical methods applied to subsurface biobarrier formation models in porous media

Biofilm forming microbes have complex effects on the flow properties of natural porous media. Subsurface biofilms have the potential for the formation of biobarriers to inhibit contaminant migration in groundwater. Another example of beneficial microbial effects is the biotransformation of organic contaminants to less harmful forms, thereby providing an in situ method for treatment of contaminated groundwater supplies. Mathematical models that describe contaminant transport with biodegradation involve a set of coupled convection-dispersion equations with non-linear reactions. The reactive solute transport equation is one for which numerical solution procedures continue to exhibit significant limitations for certain problems of groundwater hydrology interest. Accurate numerical simulations are crucial to the development of contaminant remediation strategies. A new numerical method is developed for simulation of reactive bacterial transport in porous media. The non-standard numerical approach is based on the ideas of the ‘exact’ time-stepping scheme. It leads to solutions free from the numerical instabilities that arise from incorrect modeling of derivatives and reaction terms. Applications to different biofilm models are examined and numerical results are presented to demonstrate the performance of the proposed new method.