Erratum: Self-consistent formulations for stochastic nonlinear neuronal dynamics [Phys. Rev. E 101, 042124 (2020)].

This corrects the article DOI: 10.1103/PhysRevE.101.042124.

Gell-Mann-Low Criticality in Neural Networks.

Criticality is deeply related to optimal computational capacity. The lack of a renormalized theory of critical brain dynamics, however, so far limits insights into this form of biological information processing to mean-field results. These methods neglect a key feature of critical systems: the interaction between degrees of freedom across all length scales, required for complex nonlinear computation. We present a renormalized theory of a prototypical neural field theory, the stochastic Wilson-Cowan equation. We compute the flow of couplings, which parametrize interactions on increasing length scales. Despite similarities with the Kardar-Parisi-Zhang model, the theory is of a Gell-Mann-Low type, the archetypal form of a renormalizable quantum field theory. Here, nonlinear couplings vanish, flowing towards the Gaussian fixed point, but logarithmically slowly, thus remaining effective on most scales. We show this critical structure of interactions to implement a desirable trade-off between linearity, optimal for information storage, and nonlinearity, required for computation.

Self-consistent formulations for stochastic nonlinear neuronal dynamics.

Neural dynamics is often investigated with tools from bifurcation theory. However, many neuron models are stochastic, mimicking fluctuations in the input from unknown parts of the brain or the spiking nature of signals. Noise changes the dynamics with respect to the deterministic model; in particular classical bifurcation theory cannot be applied. We formulate the stochastic neuron dynamics in the Martin-Siggia-Rose de Dominicis-Janssen (MSRDJ) formalism and present the fluctuation expansion of the effective action and the functional renormalization group (fRG) as two systematic ways to incorporate corrections to the mean dynamics and time-dependent statistics due to fluctuations in the presence of nonlinear neuronal gain. To formulate self-consistency equations, we derive a fundamental link between the effective action in the Onsager-Machlup (OM) formalism, which allows the study of phase transitions, and the MSRDJ effective action, which is computationally advantageous. These results in particular allow the derivation of an OM effective action for systems with non-Gaussian noise. This approach naturally leads to effective deterministic equations for the first moment of the stochastic system; they explain how nonlinearities and noise cooperate to produce memory effects. Moreover, the MSRDJ formulation yields an effective linear system that has identical power spectra and linear response. Starting from the better known loopwise approximation, we then discuss the use of the fRG as a method to obtain self-consistency beyond the mean. We present a new efficient truncation scheme for the hierarchy of flow equations for the vertex functions by adapting the Blaizot, Méndez, and Wschebor approximation from the derivative expansion to the vertex expansion. The methods are presented by means of the simplest possible example of a stochastic differential equation that has generic features of neuronal dynamics.

Ab initio alpha-alpha scattering.

Processes such as the scattering of alpha particles ((4)He), the triple-alpha reaction, and alpha capture play a major role in stellar nucleosynthesis. In particular, alpha capture on carbon determines the ratio of carbon to oxygen during helium burning, and affects subsequent carbon, neon, oxygen, and silicon burning stages. It also substantially affects models of thermonuclear type Ia supernovae, owing to carbon detonation in accreting carbon-oxygen white-dwarf stars. In these reactions, the accurate calculation of the elastic scattering of alpha particles and alpha-like nuclei--nuclei with even and equal numbers of protons and neutrons--is important for understanding background and resonant scattering contributions. First-principles calculations of processes involving alpha particles and alpha-like nuclei have so far been impractical, owing to the exponential growth of the number of computational operations with the number of particles. Here we describe an ab initio calculation of alpha-alpha scattering that uses lattice Monte Carlo simulations. We use lattice effective field theory to describe the low-energy interactions of protons and neutrons, and apply a technique called the 'adiabatic projection method' to reduce the eight-body system to a two-cluster system. We take advantage of the computational efficiency and the more favourable scaling with system size of auxiliary-field Monte Carlo simulations to compute an ab initio effective Hamiltonian for the two clusters. We find promising agreement between lattice results and experimental phase shifts for s-wave and d-wave scattering. The approximately quadratic scaling of computational operations with particle number suggests that it should be possible to compute alpha scattering and capture on carbon and oxygen in the near future. The methods described here can be applied to ultracold atomic few-body systems as well as to hadronic systems using lattice quantum chromodynamics to describe the interactions of quarks and gluons.

Roles of the Wnt effector POP-1/TCF in the C. elegans endomesoderm specification gene network.

In C. elegans the 4-cell stage blastomere EMS is an endomesodermal precursor. Its anterior daughter, MS, makes primarily mesodermal cells, while its posterior daughter E generates the entire intestine. The gene regulatory network underlying specification of MS and E has been the subject of study for more than 15 years. A key component of the specification of the two cells is the involvement of the Wnt/beta-catenin asymmetry pathway, which through its nuclear effector POP-1, specifies MS and E as different from each other. Loss of pop-1 function results in the mis-specification of MS as an E-like cell, because POP-1 directly represses the end-1 and end-3 genes in MS, which would otherwise promote an endoderm fate. A long-standing question has been whether POP-1 plays a role in specifying MS fate beyond repression of endoderm fate. This question has been difficult to ask because the only chromosomal lesions that remove both end-1 and end-3 are large deletions removing hundreds of genes. Here, we report the construction of bona fide end-1 end-3 double mutants. In embryos lacking activity of end-1, end-3 and pop-1 together, we find that MS fate is partially restored, while E expresses early markers of MS fate and adopts characteristics of both MS and C. Our results suggest that POP-1 is not critical for MS specification beyond repression of endoderm specification, and reveal that Wnt-modified POP-1 and END-1/3 further reinforce E specification by repressing MS fate in E. By comparison, a previous work suggested that in the related nematode C. briggsae, Cb-POP-1 is not required to repress endoderm specification in MS, in direct contrast with Ce-POP-1, but is critical for repression of MS fate in E. The findings reported here shed new light on the flexibility of combinatorial control mechanisms in endomesoderm specification in Caenorhabditis.

Multipion systems in lattice QCD and the three-pion interaction.

The ground-state energies of 2, 3, 4, and 5 pi(+)'s in a spatial volume V approximately (2.5 fm)(3) are computed with lattice QCD. By eliminating the leading contribution from three-pi(+) interactions, particular combinations of these n-pi(+) ground-state energies provide precise extractions of the pi(+)pi(+) scattering length in agreement with that obtained from calculations involving only two pi(+)'s. The three-pi(+) interaction can be isolated by forming other combinations of the n-pi(+) ground-state energies. We find a result that is consistent with a repulsive three-pi(+) interaction for m_(pi) less, similar352 MeV.