Nonlinear theory of the collective Cherenkov interaction of a dense relativistic electron beam with a dense plasma in a waveguide

A nonlinear theory of the instability of a straight relativistic dense electron beam in a plasma waveguide is derived for conditions of the stimulated collective Cherenkov effect. A study is made of a waveguide with a dense plasma such that the plasma wave excited by the beam during the instability can be escribed, with a good degree of accuracy, as a potential wave. General relativistic nonlinear equations are btained that describe the temporal dynamics of beam-plasma instabilities with allowance for plasma nonlinearity and the generation of harmonics of the initial perturbation. Under the assumption that the resonant interaction between the beam waves and the plasma waves is weak, the general equations are reduced to relativistic equations with cubic nonlinearities by using the method of expansion in small perturbations of the trajectories and momenta of the beam and plasma electrons. The reduced equations are solved analytically, the time scales on which the instability saturates are determined, and the nonlinear saturation amplitudes are obtained. A comparison between analytical solutions to the reduced equations and numerical solutions to the general nonlinear equations shows them to be in good agreement. Nonlinear processes caused by the relativistic nature of the beam are found to prevent stochastization of the system in the nonlinear stage of the well-developed instability. In contrast, a nonrelativistic electron beam is found to be subject to significant anomalous nonlinear stochastization.     

Accurate measurements of coupling constants from two-dimensional nuclear magnetic resonance spectra of proteins and determination of phi-angles

A new and simple method to measure 3JHNH alpha coupling constants of proteins by adding and subtracting traces from corresponding two-dimensional nuclear Overhauser enhanced spectroscopy and two-dimensional correlated spectroscopy cross peaks after scaling is proposed. The optimal scaling for the addition and the subtraction of the two traces is obtained by minimizing an error function. The method was proven to give accurate and precise measurements of coupling constants when tested with a series of simulated spectra. The accuracy of the method was better than 0.1 Hz for all test cases including the limiting case of J = 2.0 Hz and line-width = 11.0 Hz. The accuracy of the method was better than 0.1 Hz for all test cases including The 3JHNH alpha coupling constants were measured in two-dimensional nuclear magnetic resonance spectra of the two proteins barley serine proteinase inhibitor (CI-2) and the bacterial ribonuclease (barnase) of Bacillus amyloliquefaciens. The experimentally measured coupling constants were used to calculate the constants in a Karplus equation to be: 3JHNH alpha = 6.7 cos2(phi-60) -1.3 cos(phi-60) +1.5. These constants are in good accordance with those obtained for basic pancreatic trypsin inhibitor (BPTI). In addition, special emphasis is given to the measurements of positive phi-angles, and to the contribution of molecular dynamics on the apparent coupling constants.     

A Global Numerical Analysis of the “Central Incisor / Local Maxillary Bone” System using a Meshless Method

In this work the maxillary central incisor is numerically analysed with an advance discretization technique – Natural Neighbour Radial Point Interpolation Method (NNRPIM). The NNRPIM permits to organically determine the nodal connectivity, which is essential to construct the interpolation functions. The NNRPIM procedure, based uniquely in the computational nodal mesh discretizing the problem domain, allows to obtain autonomously the required integration mesh, permitting to numerically integrate the differential equations ruling the studied physical phenomenon. A numerical analysis of a tooth structure using a meshless method is presented for the first time. A two-dimensional model of the maxillary central incisor, based on the clinical literature, is established and two distinct analyses are performed. First, a complete elasto-static analysis of the incisor/maxillary structure using the NNRPIM is evaluated and then a non-linear iterative bone tissue remodelling analysis of the maxillary bone, surrounding the central incisive, is performed. The obtained NNRPIM solutions are compared with other numerical methods solutions available in the literature and with clinical cases. The results show that the NNRPIM is a suitable numerical method to analyse numerically dental biomechanics problems.     

Computational Approaches to Solving Equations Arising from Wound Healing

In the wound healing process, the cell movement associated with chemotaxis generally outweighs the movement associated with random motion, leading to advection-dominated mathematical models of wound healing. The equations in these models must be solved with care, but often inappropriate approaches are adopted. Two one-dimensional test problems arising from advection-dominated models of wound healing are solved using four algorithms—MATLAB’s inbuilt routine pdepe.m, the Numerical Algorithms Group routine d03pcf.f, and two finite volume methods. The first finite volume method is based on a first-order upwinding treatment of chemotaxis terms and the second on a flux limiting approach. The first test problem admits an analytic solution which can be used to validate the numerical results by analyzing two measures of the error for each method: the average absolute difference and a mass balance error. These criteria as well as the visual comparison between the numerical methods and the exact solution lead us to conclude that flux limiting is the best approach to solving advection-dominated wound healing problems numerically in one dimension. The second test problem is a coupled nonlinear three species model of wound healing angiogenesis. Measurement of the mass balance error for this test problem further confirms our hypothesis that flux limiting is the most appropriate method for solving advection-dominated governing equations in wound healing models. We also consider two two-dimensional test problems arising from wound healing, one that admits an analytic solution and a more complicated problem of blood vessels growth into a devascularized wound bed. The results from the two-dimensional test problems also demonstrate that the flux limiting treatment of advective terms is ideal for an advection-dominated problem.     

Solution of a set of Maxwell-Lorentz equations for a ring relativistic electron beam

Nonlinear solutions to a set of Maxwell’s equations and the relativistic equations of electron motion are obtained that describe the equilibrium of a high-power ring relativistic electron beam against the background of immobile ions. By transforming the basic equations, a set of equations for a three-component vortex vector field is derived that describes ring beam configurations for plasma confinement. An example of a numerical calculation of the steady state of a compact beam torus of immobile ions and relativistic electrons is presented.     

Adaptive Finite Element Analysis of the Anisotropic Biphasic Theory of Tissue-Equivalent Mechanics

The nonlinear partial differential equations of the anisotropic biphasic theory of tissue-equivalent mechanics are solved with axial symmetry by an adaptive finite element system. The adaptive procedure operates within a method-of-lines framework using finite elements in space and backward difference software in time. Spatial meshes are automatically refined, coarsened, and relocated in response to error indications and material deformation. Problems with arbitrarily complex two-dimensional regions may be addressed. With meshes graded in high-error regions, the adaptive solutions have fewer degrees of freedom than solutions with comparable accuracy obtained on fixed quasi-uniform meshes. The adaptive software is used to address problems involving an isometric cell traction assay, where a cylindrical tissue equivalent is adhered at its end to fixed circular platens; a prototypical bioartificial artery; and a novel configuration that is intended as an initial step in a study to determine bioartificial arteries having optimal collagen and cell concentrations.     

Solution of large compartmental models using numerical transform inversion

Compartmental models of biological or physical systems are often described by a system of “stiff” differential equations. In this paper an algorithm for solving a system with linear coefficients is presented that employs numerical inversion of the Laplace transform of the model equations. The inversion algorithms and Gear's backward differentiation method are compared for two stiff test problems and a differential system governing a 27-compartment model of bile acid transport and metabolism. The inversion algorithm is reliable, requires modest computation time on a desktop computer and provides better accuracy than Gear's method, especially for the extremely stiff example.     

Boundary conditions for the electromagnetic field in a waveguide with a thin-walled annular plasma in the context of application to plasma microwave electronics

Effective boundary conditions for the electromagnetic field of the slow surface waves of a thinwalled annular plasma in a metal waveguide are derived and justified. With the boundary conditions obtained, there is no need to solve field equations in the plasma region of the waveguide, so that the dispersion properties of plasma waveguides can be investigated analytically for an arbitrary strength of the external magnetic field. Examples are given that show how to use the effective boundary conditions in order to describe surface waves with a normal and an anomalous dispersion. The boundary conditions are then employed to construct a theory of the radiative Cherenkov instabilities of a thin-walled annular electron beam in a waveguide with a thinwalled annular plasma. The single-particle and collective Cherenkov effects associated with low-and high-frequency surface waves in an arbitrary external magnetic field are studied analytically. The method of the effective boundary conditions is justified in the context of application to the problems of plasma relativistic microwave electronics.     

Nonstationary Model of a Low-Density Plasma Jet Ejected from a Stationary Plasma Thruster

A time-dependent problem on finding the parameters of a low-density plasma jet ejected from a stationary plasma thruster is analyzed. In contrast to earlier works in which a stationary problem was solved using a Crook-type kinetic model, a set of kinetic equations is derived at the accuracy level accepted in the kinetic theory with allowance for resonant charge-exchange interaction between ions and neutrals. A method to solve the obtained set of equations is developed and used to perform test simulations, some results of which are presented.     

An integrate-and-fire model for synchronized bursting in a network of cultured cortical neurons

It has been suggested that spontaneous synchronous neuronal activity is an essential step in the formation of functional networks in the central nervous system. The key features of this type of activity consist of bursts of action potentials with associated spikes of elevated cytoplasmic calcium. These features are also observed in networks of rat cortical neurons that have been formed in culture. Experimental studies of these cultured networks have led to several hypotheses for the mechanisms underlying the observed synchronized oscillations. In this paper, bursting integrate-and-fire type mathematical models for regular spiking (RS) and intrinsic bursting (IB) neurons are introduced and incorporated through a small-world connection scheme into a two-dimensional excitatory network similar to those in the cultured network. This computer model exhibits spontaneous synchronous activity through mechanisms similar to those hypothesized for the cultured experimental networks. Traces of the membrane potential and cytoplasmic calcium from the model closely match those obtained from experiments. We also consider the impact on network behavior of the IB neurons, the geometry and the small world connection scheme. Action Editor: David Golomb     

Nonlinear relativistic plasma resonance: Renormalization group approach

An analytical solution to the nonlinear set of equations describing the electron dynamics and electric field structure in the vicinity of the critical density in a nonuniform plasma is constructed using the renormalization group approach with allowance for relativistic effects of electron motion. It is demonstrated that the obtained solution describes two regimes of plasma oscillations in the vicinity of the plasma resonance— stationary and nonstationary. For the stationary regime, the spatiotemporal and spectral characteristics of the resonantly enhanced electric field are investigated in detail and the effect of the relativistic nonlinearity on the spatial localization of the energy of the plasma relativistic field is considered. The applicability limits of the obtained solution, which are determined by the conditions of plasma wave breaking in the vicinity of the resonance, are established and analyzed in detail for typical laser and plasma parameters. The applicability limits of the earlier developed nonrelativistic theories are refined.     

A second-order high resolution finite difference scheme for a structured erythropoiesis model subject to malaria infection

We develop a second-order high-resolution finite difference scheme to approximate the solution of a mathematical model describing the within-host dynamics of malaria infection. The model consists of two nonlinear partial differential equations coupled with three nonlinear ordinary differential equations. Convergence of the numerical method to the unique weak solution with bounded total variation is proved. Numerical simulations demonstrating the achievement of the designed accuracy are presented.     

Adaptive Finite Element Analysis of the Anisotropic Biphasic Theory of Tissue-Equivalent Mechanics

Abstract

The nonlinear partial differential equations of the anisotropic biphasic theory of tissue-equivalent mechanics are solved with axial symmetry by an adaptive finite element system. The adaptive procedure operates within a method-of-lines framework using finite elements in space and backward difference software in time. Spatial meshes are automatically refined, coarsened, and relocated in response to error indications and material deformation. Problems with arbitrarily complex two-dimensional regions may be addressed. With meshes graded in high-error regions, the adaptive solutions have fewer degrees of freedom than solutions with comparable accuracy obtained on fixed quasi-uniform meshes. The adaptive software is used to address problems involving an isometric cell traction assay, where a cylindrical tissue equivalent is adhered at its end to fixed circular platens; a prototypical bioartificial artery; and a novel configuration that is intended as an initial step in a study to determine bioartificial arteries having optimal collagen and cell concentrations.     

Nonlineart theory of relativistic beam-plasma instabilities in the regime of the collective Cherenkov effect

A general mathematical model is proposed that is based on the Vlasov kinetic equation with a self-consistent field and describes the nonlinear dynamics of the electromagnetic instabilities of a relativistic electron beam in a spatially bounded plasma. Two limiting cases are analyzed, namely, high-frequency (HF) and low-frequency (LF) instabilities of a relativistic electron beam, of which the LF instability is a qualitatively new phenomenon in comparison with the known Cherenkov resonance effects. For instabilities in the regime of the collective Cherenkov effect, the equations containing cubic nonlinearities and describing the nonlinear saturation of the instabilities of a relativistic beam in a plasma are derived by using the methods of expansion in small perturbations of the trajectories and momenta of the beam electrons. Analytic expressions for the amplitudes of the interacting beam and plasma waves are obtained. The analytical results are shown to agree well with the exact solutions obtained numerically from the basic general mathematical model of the instabilities in question. The general mathematical model is also used to discuss the effects associated with variation in the constant component of the electron current in a beam-plasma system.     

Nonlinear extraordinary wave in dense plasma

Conditions for the propagation of a slow extraordinary wave in dense magnetized plasma are found. A solution to the set of relativistic hydrodynamic equations and Maxwell’s equations under the plasma resonance conditions, when the phase velocity of the nonlinear wave is equal to the speed of light, is obtained. The deviation of the wave frequency from the resonance frequency is accompanied by nonlinear longitudinal-transverse oscillations. It is shown that, in this case, the solution to the set of self-consistent equations obtained by averaging the initial equations over the period of high-frequency oscillations has the form of an envelope soliton. The possibility of excitation of a nonlinear wave in plasma by an external electromagnetic pulse is confirmed by numerical simulations.     

Accurate cancer classification using expressions of very few genes

We aim at finding the smallest set of genes that can ensure highly accurate classification of cancers from microarray data by using supervised machine learning algorithms. The significance of finding the minimum gene subsets is three-fold: 1) it greatly reduces the computational burden and "noise" arising from irrelevant genes. In the examples studied in this paper, finding the minimum gene subsets even allows for extraction of simple diagnostic rules which lead to accurate diagnosis without the need for any classifiers, 2) it simplifies gene expression tests to include only a very small number of genes rather than thousands of genes, which can bring down the cost for cancer testing significantly, 3) it calls for further investigation into the possible biological relationship between these small numbers of genes and cancer development and treatment. Our simple yet very effective method involves two steps. In the first step, we choose some important genes using a feature importance ranking scheme. In the second step, we test the classification capability of all simple combinations of those important genes by using a good classifier. For three "small" and "simple" data sets with two, three, and four cancer (sub)types, our approach obtained very high accuracy with only two or three genes. For a "large" and "complex" data set with 14 cancer types, we divided the whole problem into a group of binary classification problems and applied the 2-step approach to each of these binary classification problems. Through this "divide-and-conquer" approach, we obtained accuracy comparable to previously reported results but with only 28 genes rather than 16,063 genes. In general, our method can significantly reduce the number of genes required for highly reliable diagnosis