Estudos experimentais sôbre a origem do milho

1) It may seem rather strange that, in spite of the efforts of a considerable number of scientists, the problem of the origin of indian corn or maize still has remained an open question. There are no fossil remains or archaeological relics except those which are quite identical with types still existing. (Fig. 1). The main difficulty in finding the wild ancestor- which may still exist - results from the fact that it has been somewhat difficult to decide what it should be like and also where to look for it. 2) There is no need to discuss the literature since an excellent review has recently been published by MANGELSDORF and REEVES (1939). It may be sufficient to state that there are basically two hypotheses, that of ST. HILAIRE (1829) who considered Brazilian pod corn as the nearest relative of wild corn still existing, and that of ASCHERSON (1875) who considered Euchlaena from Central America as the wild ancestor of corn. Later hypotheses represent or variants of these two hypotheses or of other concepts, howewer generally with neither disproving their predecessors nor showing why the new hypotheses were better than the older ones. Since nearly all possible combinations of ideas have thus been put forward, it har- dly seems possible to find something theoretically new, while it is essential first to produce new facts. 3) The studies about the origin of maize received a new impulse from MANGELSDORF and REEVES'S experimental work on both Zea-Tripsacum and Zea-Euchlaena hybrids. Independently I started experiments in 1937 with the hope that new results might be obtained when using South American material. Having lost priority in some respects I decided to withold publication untill now, when I can put forward more concise ideas about the origin of maize, based on a new experimental reconstruction of the "wild type". 4) The two main aspects of MANGELSDORF and REEVES hypothesis are discussed. We agree with the authors that ST. HILAIRE's theory is probably correct in so far as the tunicata gene is a wild type relic gene, but cannot accept the reconstruction of wild corn as a homozygous pod corn with a hermaphroditic tassel. As shown experimentally (Fig. 2-3) these tassels have their central spike transformed into a terminal, many rowed ear with a flexible rachis, while possessing at the same time the lateral ear. Thus no explanation is given of the origin of the corn ear, which is the main feature of cultivated corn (BRIEGER, 1943). The second part of the hypothesis referring to the origin of Euchlaena from corn, inverting thus ASCHERSON's theory, cannot be accepted for several reasons, stated in some detail. The data at hand justify only the conclusion that both genera, Euchlaena and Zea, are related, and there is as little proof for considering the former as ancestor of the latter as there is for the new inverse theory. 5) The analysis of indigenous corn, which will be published in detail by BRIEGER and CUTLER, showed several very primitive characters, but no type was found which was in all characters sufficiently primitive. A genetical analysis of Paulista Pod Corn showed that it contains the same gene as other tunicates, in the IV chromosome, the segregation being complicated by a new gametophyte factor Ga3. The full results of this analysis shall be published elsewhere. (BRIEGER). Selection experiments with Paulista Pod Corn showed that no approximation to a wild ancestor may be obtained when limiting the studies to pure corn. Thus it seemed necessary to substitute "domesticated" by "wild type" modifiers, and the only means for achieving this substitution are hybridizations with Euchlaena. These hybrids have now been analysed init fourth generation, including backcrosses, and, again, the full data will be published elsewhere, by BRIEGER and ADDISON. In one present publication three forms obtained will be described only, which represent an approximation to wild type corn. 6) Before entering howewer into detail, some arguments against ST. HILAIRE's theory must be mentioned. The premendelian argument, referring to the instability of this character, is explained by the fact that all fertile pod corn plants are heterozygous for the dominant Tu factor. But the sterility of the homozygous TuTu, which phenotypically cannot be identified, is still unexplained. The most important argument against the acceptance of the Tunicata faetor as wild type relic gene was removed recently by CUTLER (not yet published) who showed that this type has been preserved for centuries by the Bolivian indians as a mystical "medicine". 7) The main botanical requirements for transforming the corn ear into a wild type structure are stated, and alternative solutions given. One series of these characters are found in Tripsacum and Euchlaena : 2 rows on opposite sides of the rachis, protection of the grains by scales, fragility of the rachis. There remains the other alternative : 4 rows, possibly forming double rows of female and male spikelets, protection of kernels by their glumes, separation of grains at their base from the cob which is thin and flexible. 8) Three successive stages in the reconstruction of wild corn, obtained experimentally, are discussed and illustrated, all characterized by the presence of the Tu gene. a) The structure of the Fl hybrids has already been described in 1943. The main features of the Tunicata hybrids (Fig. -8), when compared with non-tunicate hybrids (Fig. 5-6), consist in the absence of scaly protections, the fragility of the rachis and finally the differentiation of the double rows into one male and one female spikelet. As has been pointed out, these characters represent new phenotypic effects of the tunicate factor which do not appear in the presence of pure maize modifiers. b) The next step was observed among the first backcross to teosinte (Fig. 9). As shown in the photography, Fig. 9D, the features are essencially those of the Fl plants, except that the rachis is more teosinte like, with longer internodes, irregular four-row-arrangement and a complete fragility on the nodes. c) In the next generation a completely new type appeared (Fig. 10) which resembles neither corn nor teosinte, mainly in consequence of one character: the rachis is thin and flexible and not fragile, while the grains have an abscission layer at the base, The medium sized, pointed, brownish and hard granis are protected by their well developed corneous glumes. This last form may not yet be the nearest approach to a wild grass, and I shall try in further experiments to introduce other changes such as an increase of fertile flowers per spikelet, the reduction of difference between terminal and lateral inflorescences, etc.. But the nature of the atavistic reversion is alveadwy such that it alters considerably our expectation when looking for a still existing wild ancestor of corn. 9) The next step in our deductions must now consist in an reversion of our question. We must now explain how we may obtain domesticated corn, starting from a hypothetical wild plant, similar to type c. Of the several changes which must have been necessary to attract the attention of the Indians, the following two seem to me the most important: the disappearance of all abscission layers and the reduction of the glumes. This may have been brought about by an accumulation of mutations. But it seems much more probable to assume that some crossing with a tripsacoid grass or even with Tripsacum australe may have been responsible. In such a cross, the two types of abscission layer would be counterbalanced as shown by the Flhybrids of corn, Tripsacum and Euchlaena. Furthermore in later generations a.tu-allele of Tripsacum may become homozygous and substitute the wild tunicate factor of corn. The hypothesis of a hybrid origin of cultivated corn is not completely new, but has been discussed already by HARSHBERGER and COLLINS. Our hypothesis differs from that of MANGELSDORF and REEVES who assume that crosses with Tripsacum are responsible only for some features of Central and North American corn. 10) The following arguments give indirects evidence in support of our hypothesis: a) Several characters have been observed in indigenous corn from the central region of South America, which may be interpreted as "tripsacoid". b) Equally "zeoid" characters seem to be present in Tripsacum australe of central South-America. c) A system of unbalanced factors, combined by the in-tergeneric cross, may be responsible for the sterility of the wild type tunicata factor when homozygous, a result of the action of modifiers, brought in from Tripsacum together with the tuallele. d) The hybrid theory may explain satisfactorily the presence of so many lethals and semilethals, responsible for the phenomenon of inbreeding in cultivated corn. It must be emphasized that corn does not possess any efficient mechanism to prevent crossing and which could explain the accumulation of these mutants during the evolutionary process. Teosinte which'has about the same mechanism of sexual reproduction has not accumulated such genes, nor self-sterile plants in spite of their pronounced preference for crossing. 11) The second most important step in domestication must have consisted in transforming a four rowed ear into an ear with many rows. The fusion theory, recently revived byLANGHAM is rejected. What happened evidently, just as in succulent pXants (Cactus) or in cones os Gymnosperms, is that there has been a change in phyllotaxy and a symmetry of longitudinal rows superimposed on the original spiral arrangement. 12) The geographical distribution of indigenous corn in South America has been discussed. So far, we may distinguish three zones. The most primitive corn appears in the central lowlands of what I call the Central Triangle of South America: east of the Andies, south of the Amazone-Basin, Northwest of a line formed by the rivers São Prancisco-Paraná and including the Paraguay-Basin. The uniformity of the types found in this extremely large zone is astonishing (BRIEGER and CUTLER). To the west, there is the well known Andian region, characterized by a large number of extremely diverse types from small pop corn to large Cuszco, from soft starch to modified sweet corn, from large cylindrical ears to small round ears, etc.. The third region extends along the atlantic coast in the east, from the Caribean Sea to the Argentine, and is characterized by Cateto, an orange hard flint corn. The Andean types must have been obtained very early, and undoubtedly are the result of the intense Inca agriculture. The Cateto type may be obtained easily by crosses, for instance, of "São Paulo Pointed Pop" to some orange soft corn of the central region. The relation of these three South American zones to Central and North America are not discussed, and it seems essential first to study the intermediate region of Ecuador, Colombia and Venezuela. The geograprical distribution of chromosome knobs is rapidly discussed; but it seems that no conclusions can be drawn before a large number of Tripsacum species has been analysed.

A determinação dos números de indivíduos mínimos necessários na experimentação genética

The main object of the present paper consists in giving formulas and methods which enable us to determine the minimum number of repetitions or of individuals necessary to garantee some extent the success of an experiment. The theoretical basis of all processes consists essentially in the following. Knowing the frequency of the desired p and of the non desired ovents q we may calculate the frequency of all possi- ble combinations, to be expected in n repetitions, by expanding the binomium (p-+q)n. Determining which of these combinations we want to avoid we calculate their total frequency, selecting the value of the exponent n of the binomium in such a way that this total frequency is equal or smaller than the accepted limit of precision n/pª{ 1/n1 (q/p)n + 1/(n-1)| (q/p)n-1 + 1/ 2!(n-2)| (q/p)n-2 + 1/3(n-3) (q/p)n-3... < Plim - -(1b) There does not exist an absolute limit of precision since its value depends not only upon psychological factors in our judgement, but is at the same sime a function of the number of repetitions For this reasen y have proposed (1,56) two relative values, one equal to 1-5n as the lowest value of probability and the other equal to 1-10n as the highest value of improbability, leaving between them what may be called the "region of doubt However these formulas cannot be applied in our case since this number n is just the unknown quantity. Thus we have to use, instead of the more exact values of these two formulas, the conventional limits of P.lim equal to 0,05 (Precision 5%), equal to 0,01 (Precision 1%, and to 0,001 (Precision P, 1%). The binominal formula as explained above (cf. formula 1, pg. 85), however is of rather limited applicability owing to the excessive calculus necessary, and we have thus to procure approximations as substitutes. We may use, without loss of precision, the following approximations: a) The normal or Gaussean distribution when the expected frequency p has any value between 0,1 and 0,9, and when n is at least superior to ten. b) The Poisson distribution when the expected frequecy p is smaller than 0,1. Tables V to VII show for some special cases that these approximations are very satisfactory. The praticai solution of the following problems, stated in the introduction can now be given: A) What is the minimum number of repititions necessary in order to avoid that any one of a treatments, varieties etc. may be accidentally always the best, on the best and second best, or the first, second, and third best or finally one of the n beat treatments, varieties etc. Using the first term of the binomium, we have the following equation for n: n = log Riim / log (m:) = log Riim / log.m - log a --------------(5) B) What is the minimun number of individuals necessary in 01der that a ceratin type, expected with the frequency p, may appaer at least in one, two, three or a=m+1 individuals. 1) For p between 0,1 and 0,9 and using the Gaussean approximation we have: on - ó. p (1-p) n - a -1.m b= δ. 1-p /p e c = m/p } -------------------(7) n = b + b² + 4 c/ 2 n´ = 1/p n cor = n + n' ---------- (8) We have to use the correction n' when p has a value between 0,25 and 0,75. The greek letters delta represents in the present esse the unilateral limits of the Gaussean distribution for the three conventional limits of precision : 1,64; 2,33; and 3,09 respectively. h we are only interested in having at least one individual, and m becomes equal to zero, the formula reduces to : c= m/p o para a = 1 a = { b + b²}² = b² = δ2 1- p /p }-----------------(9) n = 1/p n (cor) = n + n´ 2) If p is smaller than 0,1 we may use table 1 in order to find the mean m of a Poisson distribution and determine. n = m: p C) Which is the minimun number of individuals necessary for distinguishing two frequencies p1 and p2? 1) When pl and p2 are values between 0,1 and 0,9 we have: n = { δ p1 ( 1-pi) + p2) / p2 (1 - p2) n= 1/p1-p2 }------------ (13) n (cor) We have again to use the unilateral limits of the Gaussean distribution. The correction n' should be used if at least one of the valors pl or p2 has a value between 0,25 and 0,75. A more complicated formula may be used in cases where whe want to increase the precision : n (p1 - p2) δ { p1 (1- p2 ) / n= m δ = δ p1 ( 1 - p1) + p2 ( 1 - p2) c= m / p1 - p2 n = { b2 + 4 4 c }2 }--------- (14) n = 1/ p1 - p2 2) When both pl and p2 are smaller than 0,1 we determine the quocient (pl-r-p2) and procure the corresponding number m2 of a Poisson distribution in table 2. The value n is found by the equation : n = mg /p2 ------------- (15) D) What is the minimun number necessary for distinguishing three or more frequencies, p2 p1 p3. If the frequecies pl p2 p3 are values between 0,1 e 0,9 we have to solve the individual equations and sue the higest value of n thus determined : n 1.2 = {δ p1 (1 - p1) / p1 - p2 }² = Fiim n 1.2 = { δ p1 ( 1 - p1) + p1 ( 1 - p1) }² } -- (16) Delta represents now the bilateral limits of the : Gaussean distrioution : 1,96-2,58-3,29. 2) No table was prepared for the relatively rare cases of a comparison of threes or more frequencies below 0,1 and in such cases extremely high numbers would be required. E) A process is given which serves to solve two problemr of informatory nature : a) if a special type appears in n individuals with a frequency p(obs), what may be the corresponding ideal value of p(esp), or; b) if we study samples of n in diviuals and expect a certain type with a frequency p(esp) what may be the extreme limits of p(obs) in individual farmlies ? I.) If we are dealing with values between 0,1 and 0,9 we may use table 3. To solve the first question we select the respective horizontal line for p(obs) and determine which column corresponds to our value of n and find the respective value of p(esp) by interpolating between columns. In order to solve the second problem we start with the respective column for p(esp) and find the horizontal line for the given value of n either diretly or by approximation and by interpolation. 2) For frequencies smaller than 0,1 we have to use table 4 and transform the fractions p(esp) and p(obs) in numbers of Poisson series by multiplication with n. Tn order to solve the first broblem, we verify in which line the lower Poisson limit is equal to m(obs) and transform the corresponding value of m into frequecy p(esp) by dividing through n. The observed frequency may thus be a chance deviate of any value between 0,0... and the values given by dividing the value of m in the table by n. In the second case we transform first the expectation p(esp) into a value of m and procure in the horizontal line, corresponding to m(esp) the extreme values om m which than must be transformed, by dividing through n into values of p(obs). F) Partial and progressive tests may be recomended in all cases where there is lack of material or where the loss of time is less importent than the cost of large scale experiments since in many cases the minimun number necessary to garantee the results within the limits of precision is rather large. One should not forget that the minimun number really represents at the same time a maximun number, necessary only if one takes into consideration essentially the disfavorable variations, but smaller numbers may frequently already satisfactory results. For instance, by definition, we know that a frequecy of p means that we expect one individual in every total o(f1-p). If there were no chance variations, this number (1- p) will be suficient. and if there were favorable variations a smaller number still may yield one individual of the desired type. r.nus trusting to luck, one may start the experiment with numbers, smaller than the minimun calculated according to the formulas given above, and increase the total untill the desired result is obtained and this may well b ebefore the "minimum number" is reached. Some concrete examples of this partial or progressive procedure are given from our genetical experiments with maize.

Diet of Procyon cancrivorus (Carnivora, Procyonidae) in restinga and estuarine environments of southern Brazil

Despite its wide range and abundance on certain habitats, the crab-eating raccoon Procyon cancrivorus (G. Cuvier, 1798) is considered one of the less known Neotropical carnivore species. In the present study we analyzed the diet of P. cancrivorus in a peat forest and in an estuarine island in southernmost Brazil. Fruits of the gerivá palm tree Syagrus romanzoffiana were the most consumed item in the peat forest, followed by insects and mollusks. Small mammals, followed by Bromelia antiacantha (Bromeliaceae) fruits and brachyuran crustaceans were the most frequent items in the estuarine island. Other items found in lower frequencies were Solanum sp., Psidium sp., Smilax sp. and Dyospiros sp. fruits, diplopods, scorpions, fishes, anuran amphibians, reptiles (black tegu lizard and snakes), birds and medium-sized mammals (white-eared opossum, armadillo and coypu). Levin’s index values (peat forest: 0.38; estuarine island: 0.45) indicate an approximation to a median position between a specialist and a well distributed diet. Pianka’s index (0.80) showed a considerable diet similarity between the two systems. Procyon cancrivorus presented a varied diet in the studied areas and may play an important role as seed disperser on coastal environments in southernmost Brazil.

Approximating subtree distances between Phylogenies

We give a 5-approximation algorithm to the rooted Subtree-Prune-and-Regraft (rSPR) distance between two phylogenies, which was recently shown to be NP-complete by Bordewich and Semple [5]. This paper presents the first approximation result for this important tree distance. The algorithm follows a standard format for tree distances such as Rodrigues et al. [24] and Hein et al. [13]. The novel ideas are in the analysis. In the analysis, the cost of the algorithm uses a \cascading" scheme that accounts for possible wrong moves. This accounting is missing from previous analysis of tree distance approximation algorithms. Further, we show how all algorithms of this type can be implemented in linear time and give experimental results.

Convergence of the mass-transport steepest descent scheme for the sub-critical Patlak-Keller-Segel model

Variational steepest descent approximation schemes for the modified Patlak-Keller-Segel equation with a logarithmic interaction kernel in any dimension are considered. We prove the convergence of the suitably interpolated in time implicit Euler scheme, defined in terms of the Euclidean Wasserstein distance, associated to this equation for sub-critical masses. As a consequence, we recover the recent result about the global in time existence of weak-solutions to the modified Patlak-Keller-Segel equation for the logarithmic interaction kernel in any dimension in the sub-critical case. Moreover, we show how this method performs numerically in one dimension. In this particular case, this numerical scheme corresponds to a standard implicit Euler method for the pseudo-inverse of the cumulative distribution function. We demonstrate its capabilities to reproduce easily without the need of mesh-refinement the blow-up of solutions for super-critical masses.

Diseño de filtros con topología apilada mediante resonadores FBAR

Els continguts de la memòria es divideixen en dues parts fonamentals, precedides d’una breu explicació on s’introdueix al lector en el tema; primer s’exposen les bases d’un ressonador simple FBAR, d’on s’obtenen equacions de control vitals, i seguidament es plantegen les bases per a un ressonador apilat SCR, derivades de l’estudi previ d’un ressonador simple. La primera part fonamental del treball es centra en l’anàlisi d’un ressonador SCR. Aquest anàlisi es recolza sobre el punt de vista teòric dels paràmetres imatge, el nostre punt de sortida. Aquesta primera part és la més teòrica: obtenció i aplicació dels paràmetres imatge i obtenció dels elements discrets que conformen el circuit equivalent SCR, una xarxa de dos ports. Posteriorment, s’analitza què succeeix modificant els valors dels elements discrets sense variar determinats paràmetres imatge i finalment, es proposa una aproximació per a controlar determinades especificacions de disseny, com són l’ample de banda de transmissió de la xarxa i el factor de qualitat, en funció de les modificacions dels elements discrets. En la segona part s’analitzen, de forma qualitativa, xarxes compostes per diferents ressonadors apilats connectats en cascada. Aquest segon estudi es divideix en dues subparts. En la primera connectem N ressonadors idèntics, plantegem algunes equacions de control i analitzem les respostes. En la segona, es planteja la connexió de dos (N=2) ressonadors diferents amb freqüències de ressonància properes. Aquest segon anàlisi es també totalment qualitatiu, però ens aporta informació que amb la unió de N ressonadors iguals no aconseguíem. Finalitzant aquesta segona part, es planteja l’optimització dels resultats obtinguts per a N=2, mitjançant estructures N=3 ressonadors.

Ordering Policy Rules with an Unconditional Welfare Measure

The unconditional expectation of social welfare is often used to assess alternative macroeconomic policy rules in applied quantitative research. It is shown that it is generally possible to derive a linear - quadratic problem that approximates the exact non-linear problem where the unconditional expectation of the objective is maximised and the steady-state is distorted. Thus, the measure of pol icy performance is a linear combinat ion of second moments of economic variables which is relatively easy to compute numerically, and can be used to rank alternative policy rules. The approach is applied to a simple Calvo-type model under various monetary policy rules.

Tractable Consumer Choice

We derive a rational model of separable consumer choice which can also serve as a behavioral model. The central construct is [lambda] , the marginal utility of money, derived from the consumer's rest-of-life problem. We present a robust approximation of [lambda], and show how to incorporate liquidity constraints, indivisibilities and adaptation to a changing environment. We fi nd connections with numerous historical and recent constructs, both behavioral and neoclassical, and draw contrasts with standard partial equilibrium analysis. The result is a better grounded, more flexible and more intuitive description of consumer choice.

The Inflation Bias under Calvo and Rotemberg Pricing.

New Keynesian models rely heavily on two workhorse models of nominal inertia - price contracts of random duration (Calvo, 1983) and price adjustment costs (Rotemberg, 1982) - to generate a meaningful role for monetary policy. These alternative descriptions of price stickiness are often used interchangeably since, to a first order of approximation they imply an isomorphic Phillips curve and, if the steady-state is efficient, identical objectives for the policy maker and as a result in an LQ framework, the same policy conclusions. In this paper we compute time-consistent optimal monetary policy in bench-mark New Keynesian models containing each form of price stickiness. Using global solution techniques we find that the inflation bias problem under Calvo contracts is significantly greater than under Rotemberg pricing, despite the fact that the former typically significant exhibits far greater welfare costs of inflation. The rates of inflation observed under this policy are non-trivial and suggest that the model can comfortably generate the rates of inflation at which the problematic issues highlighted in the trend inflation literature emerge, as well as the movements in trend inflation emphasized in empirical studies of the evolution of inflation. Finally, we consider the response to cost push shocks across both models and find these can also be significantly different. The choice of which form of nominal inertia to adopt is not innocuous.

Bilateral oligopoly and quantity competition

Bilateral oligopoly is a strategic market game with two commodities, allowing strategic behavior on both sides of the market. When the number of buyers is large, such a game approximates a game of quantity competition played by sellers. We present examples which show that this is not typically a Cournot game. Rather, we introduce an alternative game of quantity competition (the market share game) and, appealing to results in the literature on contests, show that this yields the same equilibria as the many-buyer limit of bilateral oligopoly, under standard assumptions on costs and preferences. We also show that the market share and Cournot games have the same equilibria if and only if the price elasticity of the latter is one. These results lead to necessary and sufficient conditions for the Cournot game to be a good approximation to bilateral oligopoly with many buyers and to an ordering of total output when they are not satisfied.

The non-linear Dirichlet problem in Hadamard manifolds

We prove existence theorems for the Dirichlet problem for hypersurfaces of constant special Lagrangian curvature in Hadamard manifolds. The first results are obtained using the continuity method and approximation and then refined using two iterations of the Perron method. The a-priori estimates used in the continuity method are valid in any ambient manifold.

Quantification of the neurochemical profile using simulated macromolecule resonances at 3 T.

The broad resonances underlying the entire (1) H NMR spectrum of the brain, ascribed to macromolecules, can influence metabolite quantification. At the intermediate field strength of 3 T, distinct approaches for the determination of the macromolecule signal, previously used at either 1.5 or 7 T and higher, may become equivalent. The aim of this study was to evaluate, at 3 T for healthy subjects using LCModel, the impact on the metabolite quantification of two different macromolecule approaches: (i) experimentally measured macromolecules; and (ii) mathematically estimated macromolecules. Although small, but significant, differences in metabolite quantification (up to 23% for glutamate) were noted for some metabolites, 10 metabolites were quantified reproducibly with both approaches with a Cramer-Rao lower bound below 20%, and the neurochemical profiles were therefore similar. We conclude that the mathematical approximation can provide sufficiently accurate and reproducible estimation of the macromolecule contribution to the (1) H spectrum at 3 T. Copyright © 2013 John Wiley & Sons, Ltd.

Discontinuous Galerkin methods for the one-dimensional Vlasov-Poisson system

We construct a new family of semi-discrete numerical schemes for the approximation of the one-dimensional periodic Vlasov-Poisson system. The methods are based on the coupling of discontinuous Galerkin approximation to the Vlasov equation and several finite element (conforming, non-conforming and mixed) approximations for the Poisson problem. We show optimal error estimates for the all proposed methods in the case of smooth compactly supported initial data. The issue of energy conservation is also analyzed for some of the methods.

Application of standard and refined heat balance integral methods to one-dimensional Stefan problems

The work in this paper concerns the study of conventional and refined heat balance integral methods for a number of phase change problems. These include standard test problems, both with one and two phase changes, which have exact solutions to enable us to test the accuracy of the approximate solutions. We also consider situations where no analytical solution is available and compare these to numerical solutions. It is popular to use a quadratic profile as an approximation of the temperature, but we show that a cubic profile, seldom considered in the literature, is far more accurate in most circumstances. In addition, the refined integral method can give greater improvement still and we develop a variation on this method which turns out to be optimal in some cases. We assess which integral method is better for various problems, showing that it is largely dependent on the specified boundary conditions.

A numerical solver for a nonlinear Fokker-Planck equation representation of neuronal network dynamics

To describe the collective behavior of large ensembles of neurons in neuronal network, a kinetic theory description was developed in [13, 12], where a macroscopic representation of the network dynamics was directly derived from the microscopic dynamics of individual neurons, which are modeled by conductance-based, linear, integrate-and-fire point neurons. A diffusion approximation then led to a nonlinear Fokker-Planck equation for the probability density function of neuronal membrane potentials and synaptic conductances. In this work, we propose a deterministic numerical scheme for a Fokker-Planck model of an excitatory-only network. Our numerical solver allows us to obtain the time evolution of probability distribution functions, and thus, the evolution of all possible macroscopic quantities that are given by suitable moments of the probability density function. We show that this deterministic scheme is capable of capturing the bistability of stationary states observed in Monte Carlo simulations. Moreover, the transient behavior of the firing rates computed from the Fokker-Planck equation is analyzed in this bistable situation, where a bifurcation scenario, of asynchronous convergence towards stationary states, periodic synchronous solutions or damped oscillatory convergence towards stationary states, can be uncovered by increasing the strength of the excitatory coupling. Finally, the computation of moments of the probability distribution allows us to validate the applicability of a moment closure assumption used in [13] to further simplify the kinetic theory.