Mean-field equations for cigar- and disc-shaped Bose and Fermi superfluids

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Precipitation In Small Systems--II. Mean Field Equations More Effective Than Population Balance

Part I (Manjunath et al., 1994, Chem. Engng Sci. 49, 1451-1463) of this paper showed that the random particle numbers and size distributions in precipitation processes in very small drops obtained by stochastic simulation techniques deviate substantially from the predictions of conventional population balance. The foregoing problem is considered in this paper in terms of a mean field approximation obtained by applying a first-order closure to an unclosed set of mean field equations presented in Part I. The mean field approximation consists of two mutually coupled partial differential equations featuring (i) the probability distribution for residual supersaturation and (ii) the mean number density of particles for each size and supersaturation from which all average properties and fluctuations can be calculated. The mean field equations have been solved by finite difference methods for (i) crystallization and (ii) precipitation of a metal hydroxide both occurring in a single drop of specified initial supersaturation. The results for the average number of particles, average residual supersaturation, the average size distribution, and fluctuations about the average values have been compared with those obtained by stochastic simulation techniques and by population balance. This comparison shows that the mean field predictions are substantially superior to those of population balance as judged by the close proximity of results from the former to those from stochastic simulations. The agreement is excellent for broad initial supersaturations at short times but deteriorates progressively at larger times. For steep initial supersaturation distributions, predictions of the mean field theory are not satisfactory thus calling for higher-order approximations. The merit of the mean field approximation over stochastic simulation lies in its potential to reduce expensive computation times involved in simulation. More effective computational techniques could not only enhance this advantage of the mean field approximation but also make it possible to use higher-order approximations eliminating the constraints under which the stochastic dynamics of the process can be predicted accurately.

Modeling electrocortical activity through improved local approximations of integral neural field equations

Neural field models of firing rate activity typically take the form of integral equations with space-dependent axonal delays. Under natural assumptions on the synaptic connectivity we show how one can derive an equivalent partial differential equation (PDE) model that properly treats the axonal delay terms of the integral formulation. Our analysis avoids the so-called long-wavelength approximation that has previously been used to formulate PDE models for neural activity in two spatial dimensions. Direct numerical simulations of this PDE model show instabilities of the homogeneous steady state that are in full agreement with a Turing instability analysis of the original integral model. We discuss the benefits of such a local model and its usefulness in modeling electrocortical activity. In particular, we are able to treat “patchy” connections, whereby a homogeneous and isotropic system is modulated in a spatially periodic fashion. In this case the emergence of a “lattice-directed” traveling wave predicted by a linear instability analysis is confirmed by the numerical simulation of an appropriate set of coupled PDEs.

Conformally and gauge invariant spin-2 field equations

Using an approach based on the Casimir operators of the de Sitter group, conformally invariant equations for a fundamental spin-2 field are obtained, and their consistency is discussed. It is shown that only when the spin-2 field is interpreted as a 1-form assuming values in the Lie algebra of the translation group, rather than a symmetric second-rank tensor, the field equation is both conformally and gauge invariant. © 2013 Pleiades Publishing, Ltd.

Evans functions for integral neural field equations with Heaviside firing rate function

In this paper we show how to construct the Evans function for traveling wave solutions of integral neural field equations when the firing rate function is a Heaviside. This allows a discussion of wave stability and bifurcation as a function of system parameters, including the speed and strength of synaptic coupling and the speed of axonal signals. The theory is illustrated with the construction and stability analysis of front solutions to a scalar neural field model and a limiting case is shown to recover recent results of L. Zhang [On stability of traveling wave solutions in synaptically coupled neuronal networks, Differential and Integral Equations, 16, (2003), pp.513-536.]. Traveling fronts and pulses are considered in more general models possessing either a linear or piecewise constant recovery variable. We establish the stability of coexisting traveling fronts beyond a front bifurcation and consider parameter regimes that support two stable traveling fronts of different speed. Such fronts may be connected and depending on their relative speed the resulting region of activity can widen or contract. The conditions for the contracting case to lead to a pulse solution are established. The stability of pulses is obtained for a variety of examples, in each case confirming a previously conjectured stability result. Finally we show how this theory may be used to describe the dynamic instability of a standing pulse that arises in a model with slow recovery. Numerical simulations show that such an instability can lead to the shedding of a pair of traveling pulses.

Dynamic instabilities in scalar neural field equations with space-dependent delays

In this paper we consider a class of scalar integral equations with a form of space-dependent delay. These non-local models arise naturally when modelling neural tissue with active axons and passive dendrites. Such systems are known to support a dynamic (oscillatory) Turing instability of the homogeneous steady state. In this paper we develop a weakly nonlinear analysis of the travelling and standing waves that form beyond the point of instability. The appropriate amplitude equations are found to be the coupled mean-field Ginzburg-Landau equations describing a Turing-Hopf bifurcation with modulation group velocity of O(1). Importantly we are able to obtain the coefficients of terms in the amplitude equations in terms of integral transforms of the spatio-temporal kernels defining the neural field equation of interest. Indeed our results cover not only models with axonal or dendritic delays but those which are described by a more general distribution of delayed spatio-temporal interactions. We illustrate the predictive power of this form of analysis with comparison against direct numerical simulations, paying particular attention to the competition between standing and travelling waves and the onset of Benjamin-Feir instabilities.

Modeling electrocortical activity through improved local approximations of integral neural field equations

Neural field models of firing rate activity typically take the form of integral equations with space-dependent axonal delays. Under natural assumptions on the synaptic connectivity we show how one can derive an equivalent partial differential equation (PDE) model that properly treats the axonal delay terms of the integral formulation. Our analysis avoids the so-called long-wavelength approximation that has previously been used to formulate PDE models for neural activity in two spatial dimensions. Direct numerical simulations of this PDE model show instabilities of the homogeneous steady state that are in full agreement with a Turing instability analysis of the original integral model. We discuss the benefits of such a local model and its usefulness in modeling electrocortical activity. In particular we are able to treat "patchy'" connections, whereby a homogeneous and isotropic system is modulated in a spatially periodic fashion. In this case the emergence of a "lattice-directed" traveling wave predicted by a linear instability analysis is confirmed by the numerical simulation of an appropriate set of coupled PDEs. Article published and (c) American Physical Society 2007

f Gravity and gravitational singularities

TRAUTMAN has postulated1 that the usual space−time singularity occurring in classical cosmological models and in the gravitational collapse of massive objects could be averted if intrinsic spin effects are incorporated into general relativity by adding torsion terms to the usual Einstein field equations, that is through the Einstein−Cartan theory. Invoking a primordial magnetic field for aligning all the individual nuclear spins he shows that his universe consisting of 1080 aligned neutrons collapses to a minimum radius of the order of 1 cm with a corresponding matter density of 1055 g cm-3.

f Gravity and gravitational singularities

TRAUTMAN has postulated1 that the usual space−time singularity occurring in classical cosmological models and in the gravitational collapse of massive objects could be averted if intrinsic spin effects are incorporated into general relativity by adding torsion terms to the usual Einstein field equations, that is through the Einstein−Cartan theory. Invoking a primordial magnetic field for aligning all the individual nuclear spins he shows that his universe consisting of 1080 aligned neutrons collapses to a minimum radius of the order of 1 cm with a corresponding matter density of 1055 g cm-3.

Purely magnetic spacetimes in general relativity

Given the deep links between Maxwell's equations and General Relativity, it is natural to ask whether there exist solutions of the field equations which are purely magnetic, like Dirac's magnetic monopole. The thesis succeeds in producing two entirely new families of purely magnetic solutions.

Shear-free perfect fluids with a solenoidal magnetic curvature

We investigate shear-free perfect fluid solutions of the Einstein field equations where the fluid pressure satisfies a barotropic equation of state and the spatial divergence of the magnetic part of the Weyl tensor is zero. We prove, with the exception of certain quite restricted special cases within the class of solutions in which there exists a Killing vector aligned with the vorticity and for which the magnitude of the vorticity ω is not a function of the matter density μ alone, that such a fluid is either non-rotating or non-expanding. In the restricted cases the equation of state must satisfy an over-determined differential system.

Shear-free perfect fluids with a solenoidal electric curvature

We prove that the vorticity or the expansion vanishes for any shear-free perfect fluid solution of the Einstein field equations where the pressure satisfies a barotropic equation of state and the spatial divergence of the electric part of the Weyl tensor is zero.

Criação da matéria em cosmologias com velocidade da luz variável e eletrônico não-linear

In this work we obtain the cosmological solutions and investigate the thermodynamics of matter creation in two diferent contexts. In the first we propose a cosmological model with a time varying speed of light c. We consider two diferent time dependence of c for a at Friedmann-Robertson- Walker (FRW) universe. We write the energy conservation law arising from Einstein equations and study how particles are created as c decreases with cosmic epoch. The variation of c is coupled to a cosmological Λ term and both singular and non-singular solutions are possible. We calculate the "adiabatic" particle creation rate and the total number of particles as a function of time and find the constrains imposed by the second law of thermodynamics upon the models. In the second scenario, we study the nonlinearity of the electrodynamics as a source of matter creation in the cosmological models with at FRW geometry. We write the energy conservation law arising from Einstein field equations with cosmological term Λ, solve the field equations and study how particles are created as the magnetic field B changes with cosmic epoch. We obtain solutions for the adiabatic particle creation rate, the total number of particles and the scale factor as a function of time in three cases: Λ = 0, Λ = constant and Λ α H2 (cosmological term proportional to the Hubble parameter). In all cases, the second law of thermodynamics demands that the universe is not contracting (H ≥ 0). The first two solutions are non-singular and exhibit in ationary periods. The third case studied allows an always in ationary universe for a suficiently large cosmological term

Black holes in realistic branes: Black string-like objects?

A realistic model describing a black string-like object in an expanding Universe is analyzed in the context of the McVittie's solution of the Einstein field equations. The bulk metric near the brane is provided analogously to previous solutions for black strings. In particular, we show that at least when the Hubble parameter on the brane is positive, a black string-like object seems to play a fundamental role in the braneworld scenario, generalizing the standard black strings in the context of a dynamical brane. © 2013 Elsevier B.V.