Science.gov

Sample records for relativistic mean-field models

  1. Parametrization of light clusters within relativistic mean field models

    SciTech Connect

    Ferreira, Marcio; Providencia, Constanca

    2013-06-10

    Light clusters are included in the equation of state of nuclearmatter within the relativistic mean field theory. The effect of the cluster-meson coupling constants on the dissolution density is discussed. Theoretical and experimental constraints are used to fix the cluster-meson couplings at T Almost-Equal-To 5 MeV.

  2. Relativistic mean-field models and nuclear matter constraints

    SciTech Connect

    Dutra, M.; Lourenco, O.; Carlson, B. V.; Delfino, A.; Menezes, D. P.; Avancini, S. S.; Stone, J. R.; Providencia, C.; Typel, S.

    2013-05-06

    This work presents a preliminary study of 147 relativistic mean-field (RMF) hadronic models used in the literature, regarding their behavior in the nuclear matter regime. We analyze here different kinds of such models, namely: (i) linear models, (ii) nonlinear {sigma}{sup 3}+{sigma}{sup 4} models, (iii) {sigma}{sup 3}+{sigma}{sup 4}+{omega}{sup 4} models, (iv) models containing mixing terms in the fields {sigma} and {omega}, (v) density dependent models, and (vi) point-coupling ones. In the finite range models, the attractive (repulsive) interaction is described in the Lagrangian density by the {sigma} ({omega}) field. The isospin dependence of the interaction is modeled by the {rho} meson field. We submit these sets of RMF models to eleven macroscopic (experimental and empirical) constraints, used in a recent study in which 240 Skyrme parametrizations were analyzed. Such constraints cover a wide range of properties related to symmetric nuclear matter (SNM), pure neutron matter (PNM), and both SNM and PNM.

  3. Relativistic mean field model based on realistic nuclear forces

    SciTech Connect

    Hirose, S.; Serra, M.; Ring, P.; Otsuka, T.; Akaishi, Y.

    2007-02-15

    In order to predict properties of asymmetric nuclear matter, we construct a relativistic mean field (RMF) model consisting of one-meson exchange (OME) terms and point coupling (PC) terms. In order to determine the density dependent parameters of this model, we use properties of isospin symmetric nuclear matter in combination with the information on nucleon-nucleon scattering data, which are given in the form of the density dependent G-matrix derived from Brueckner calculations based on the Tamagaki potential. We show that the medium- and long-range components of this G-matrix can be described reasonably well by our effective OME interaction. In order to take into account the short-range part of the nucleon-nucleon interaction, which cannot be described well in this manner, a point coupling term is added. Its analytical form is taken from a model based on chiral perturbation theory. It contains only one additional parameter, which does not depend on the density. It is, together with the parameters of the OME potentials adjusted to the equation of state of symmetric nuclear matter. We apply this model for the investigation of asymmetric nuclear matter and find that the results for the symmetry energy as well as for the equation of state of pure neutron matter are in good agreement with either experimental data or with presently adopted theoretical predictions. In order to test the model at higher density, we use its equation of state for an investigation of properties of neutron stars.

  4. Systematic study of bubble nuclei in relativistic mean field model

    NASA Astrophysics Data System (ADS)

    Shukla, A.; Åberg, S.; Bajpeyi, A.

    2016-01-01

    We have theoretically studied potential bubble nuclei (20,22O, 34,36Si, and 46Ar), which are experimentally accessible and have attracted several studies in the recent past. Relativistic mean field is employed in conjunction with the NL-SH parameter set. Our results show that among the possible candidates, 22Oand 34Si may be the most prominent candidates, showing significant depletion of density at the center, which could be verified experimentally in the near future with some of the experiments underway.

  5. Relativistic mean-field hadronic models under nuclear matter constraints

    NASA Astrophysics Data System (ADS)

    Dutra, M.; Lourenço, O.; Avancini, S. S.; Carlson, B. V.; Delfino, A.; Menezes, D. P.; Providência, C.; Typel, S.; Stone, J. R.

    2014-11-01

    Background: The microscopic composition and properties of infinite hadronic matter at a wide range of densities and temperatures have been subjects of intense investigation for decades. The equation of state (EoS) relating pressure, energy density, and temperature at a given particle number density is essential for modeling compact astrophysical objects such as neutron stars, core-collapse supernovae, and related phenomena, including the creation of chemical elements in the universe. The EoS depends not only on the particles present in the matter, but, more importantly, also on the forces acting among them. Because a realistic and quantitative description of infinite hadronic matter and nuclei from first principles in not available at present, a large variety of phenomenological models has been developed in the past several decades, but the scarcity of experimental and observational data does not allow a unique determination of the adjustable parameters. Purpose: It is essential for further development of the field to determine the most realistic parameter sets and to use them consistently. Recently, a set of constraints on properties of nuclear matter was formed and the performance of 240 nonrelativistic Skyrme parametrizations was assessed [M. Dutra et al., Phys. Rev. C 85, 035201 (2012), 10.1103/PhysRevC.85.035201] in describing nuclear matter up to about three times nuclear saturation density. In the present work we examine 263 relativistic-mean-field (RMF) models in a comparable approach. These models have been widely used because of several important aspects not always present in nonrelativistic models, such as intrinsic Lorentz covariance, automatic inclusion of spin, appropriate saturation mechanism for nuclear matter, causality, and, therefore, no problems related to superluminal speed of sound in medium. Method: Three different sets of constraints related to symmetric nuclear matter, pure neutron matter, symmetry energy, and its derivatives were used. The

  6. Relativistic mean field model for entrainment in general relativistic superfluid neutron stars

    NASA Astrophysics Data System (ADS)

    Comer, G. L.; Joynt, R.

    2003-07-01

    General relativistic superfluid neutron stars have a significantly more intricate dynamics than their ordinary fluid counterparts. Superfluidity allows different superfluid (and superconducting) species of particles to have independent fluid flows, a consequence of which is that the fluid equations of motion contain as many fluid element velocities as superfluid species. Whenever the particles of one superfluid interact with those of another, the momentum of each superfluid will be a linear combination of both superfluid velocities. This leads to the so-called entrainment effect whereby the motion of one superfluid will induce a momentum in the other superfluid. We have constructed a fully relativistic model for entrainment between superfluid neutrons and superconducting protons using a relativistic σ-ω mean field model for the nucleons and their interactions. In this context there are two notions of “relativistic”: relativistic motion of the individual nucleons with respect to a local region of the star (i.e. a fluid element containing, say, an Avogadro’s number of particles), and the motion of fluid elements with respect to the rest of the star. While it is the case that the fluid elements will typically maintain average speeds at a fraction of that of light, the supranuclear densities in the core of a neutron star can make the nucleons themselves have quite high average speeds within each fluid element. The formalism is applied to the problem of slowly rotating superfluid neutron star configurations, a distinguishing characteristic being that the neutrons can rotate at a rate different from that of the protons.

  7. Neutron-skin thickness of finite nuclei in relativistic mean-field models with chiral limits

    SciTech Connect

    Jiang Weizhou; Li Baoan; Chen Liewen

    2007-11-15

    We study several structure properties of finite nuclei using relativistic mean-field Lagrangians constructed according to the Brown-Rho scaling due to the chiral symmetry restoration at high densities. The models are consistent with current experimental constraints for the equations of state of symmetric matter at both normal and supranormal densities and of asymmetric matter at subsaturation densities. It is shown that these models can successfully describe the binding energies and charge radii of finite nuclei. Compared to calculations with usual relativistic mean-field models, these models give a reduced thickness of neutron skin in {sup 208}Pb between 0.17 fm and 0.21 fm. The reduction of the predicted neutron skin thickness is found to be due to not only the softening of the symmetry energy but also the scaling property of {rho} meson required by the partial restoration of chiral symmetry.

  8. Comparative study of nuclear masses in the relativistic mean-field model

    NASA Astrophysics Data System (ADS)

    Hua, XueMin; Heng, TaiHua; Niu, ZhongMing; Sun, BaoHua; Guo, JianYou

    2012-12-01

    With experimental masses updated from AME11, the predictive power of relativistic mean-field (RMF) mass model is carefully examined and compared with HFB-17, FRDM, WS*, and DZ28 mass models. In the relativistic mean-field model, the calculation with the PC-PK1 has improved significantly in describing masses compared to the TMA, especially for the neutron-deficient nuclei. The corresponding rms deviation with respect to the known masses falls to 1.4 MeV. Furthermore, it is found that the RMF mass model better describes the nuclei with large deformations. The rms deviation for nuclei with the absolute value of quadrupole deformation parameter greater than 0.25 falls to 0.93, crossing the 1 MeV accuracy threshold for the PC-PK1, which may indicate the new model is more suitable for those largely-deformed nuclei. In addition, the necessity of new high-precision experimental data to evaluate and develop the nuclear mass models is emphasized as well.

  9. Relativistic mean-field model with energy dependent self-energies

    SciTech Connect

    Antic, S.; Typel, S.

    2015-02-24

    Conventional relativistic mean-field theory is extended with the introduction of higher-order derivative couplings of nucleons with the meson fields. The Euler-Lagrange equations follow from the principle of stationary action. From invariance principles of the Lagrangian density the most general expressions for the conserved current and energy-momentum tensor are derived. The nucleon self-energies show the explicit dependence on the meson fields. They contain additional regulator functions which describe the energy dependence. The density dependence of meson-nucleon couplings causes the apperance of additional rearrangement contributions in the self-energies. The equation of state of infinite nuclear matter is obtained and the thermodynamical consistency of the model is demonstrated. This model is applied to the description of spherical, non-rotating stars in β-equilibrium. Stellar structure is calculated by solving the Tolman-Oppenheimer-Volkov (TOV) equations. The results for neutron stars are shown in terms of mass-radius relations.

  10. Nuclear reaction cross sections of exotic nuclei in the Glauber model for relativistic mean field densities

    SciTech Connect

    Patra, S. K.; Panda, R. N.; Arumugam, P.; Gupta, Raj K.

    2009-12-15

    We have calculated the total nuclear reaction cross sections of exotic nuclei in the framework of the Glauber model, using as inputs the standard relativistic mean field (RMF) densities and the densities obtained from the more recently developed effective-field-theory-motivated RMF (the E-RMF). Both light and heavy nuclei are taken as the representative targets, and the light neutron-rich nuclei as projectiles. We found the total nuclear reaction cross section to increase as a function of the mass number, for both the target and projectile nuclei. The differential nuclear elastic scattering cross sections are evaluated for some selected systems at various incident energies. We found a large dependence of the differential elastic scattering cross section on incident energy. Finally, we have applied the same formalism to calculate both the total nuclear reaction cross section and the differential nuclear elastic scattering cross section for the recently discussed superheavy nucleus with atomic number Z=122.

  11. Cluster radioactive decay within the preformed cluster model using relativistic mean-field theory densities

    SciTech Connect

    Singh, BirBikram; Patra, S. K.; Gupta, Raj K.

    2010-07-15

    We have studied the (ground-state) cluster radioactive decays within the preformed cluster model (PCM) of Gupta and collaborators [R. K. Gupta, in Proceedings of the 5th International Conference on Nuclear Reaction Mechanisms, Varenna, edited by E. Gadioli (Ricerca Scientifica ed Educazione Permanente, Milano, 1988), p. 416; S. S. Malik and R. K. Gupta, Phys. Rev. C 39, 1992 (1989)]. The relativistic mean-field (RMF) theory is used to obtain the nuclear matter densities for the double folding procedure used to construct the cluster-daughter potential with M3Y nucleon-nucleon interaction including exchange effects. Following the PCM approach, we have deduced empirically the preformation probability P{sub 0}{sup emp} from the experimental data on both the alpha- and exotic cluster-decays, specifically of parents in the trans-lead region having doubly magic {sup 208}Pb or its neighboring nuclei as daughters. Interestingly, the RMF-densities-based nuclear potential supports the concept of preformation for both the alpha and heavier clusters in radioactive nuclei. P{sub 0}{sup a}lpha{sup (emp)} for alpha decays is almost constant (approx10{sup -2}-10{sup -3}) for all the parent nuclei considered here, and P{sub 0}{sup c(emp)} for cluster decays of the same parents decrease with the size of clusters emitted from different parents. The results obtained for P{sub 0}{sup c(emp)} are reasonable and are within two to three orders of magnitude of the well-accepted phenomenological model of Blendowske-Walliser for light clusters.

  12. Relativistic mean-field theory

    NASA Astrophysics Data System (ADS)

    Meng, Jie; Ring, Peter; Zhao, Pengwei

    In this chapter, the covariant energy density functional is constructed with both the meson-exchange and the point-coupling pictures. Several widely used functionals with either nonlinear or density-dependent effective interactions are introduced. The applications of covariant density functional theory are demonstrated for infinite nuclear matter and finite nuclei with spherical symmetry, axially symmetric quadrupole deformation, and triaxial quadrupole shapes. Finally, a relativistic description of the nuclear landscape has been discussed, which is not only important for nuclear structure, but also important for nuclear astrophysics, where we are facing the problem of a reliable extrapolation to the very neutron-rich nuclei.

  13. Relativistic Mean Field description of exotic nuclei

    NASA Astrophysics Data System (ADS)

    Gambhir, Y. K.

    1994-03-01

    The Relativistic Mean Field (RMF) approach which essentially is an extension of the original σ — ω model of Walecka, has been applied to exotic nuclei as an illustration. We consider nuclei near Z = 34 in the very interesting 2p-1f region. The calculated binding energies, root mean square radii, deformations and other observables are very satisfactory and are in accordance with the experiment (where available) and also with the available empirical studies. Large deformations and shape co-existence are obtained for several cases.

  14. Vlasov formalism for extended relativistic mean field models: The crust-core transition and the stellar matter equation of state

    NASA Astrophysics Data System (ADS)

    Pais, Helena; Providência, Constança

    2016-07-01

    The Vlasov formalism is extended to relativistic mean field hadron models with nonlinear terms up to fourth order and applied to the calculation of the crust-core transition density. The effect of the nonlinear ω ρ and σ ρ coupling terms on the crust-core transition density and pressure and on the macroscopic properties of some families of hadronic stars is investigated. For that purpose, six families of relativistic mean field models are considered. Within each family, the members differ in the symmetry energy behavior. For all the models, the dynamical spinodals are calculated, and the crust-core transition density and pressure and the neutron star mass-radius relations are obtained. The effect on the star radius of the inclusion of a pasta calculation in the inner crust is discussed. The set of six models that best satisfy terrestrial and observational constraints predicts a radius of 13.6 ±0.3 km and a crust thickness of 1.36 ±0.06 km for a 1.4 M⊙ star.

  15. Relativistic mean-field models with scaled hadron masses and couplings: Hyperons and maximum neutron star mass

    NASA Astrophysics Data System (ADS)

    Maslov, K. A.; Kolomeitsev, E. E.; Voskresensky, D. N.

    2016-06-01

    An equation of state of cold nuclear matter with an arbitrary isotopic composition is studied within a relativistic mean-field approach with hadron masses and coupling constants depending self-consistently on the scalar mean-field. All hadron masses decrease universally with the scalar field growth, whereas meson-nucleon coupling constants can vary differently. More specifically we focus on two modifications of the KVOR model studied previously. One extension of the model (KVORcut) demonstrates that the equation of state stiffens if the increase of the scalar-field magnitude with the density is bounded from above at some value for baryon densities above the saturation nuclear density. This can be realized if the nucleon vector-meson coupling constant changes rapidly as a function of the scalar field slightly above the desired value. The other version of the model (MKVOR) utilizes a smaller value of the nucleon effective mass at the nuclear saturation density and a saturation of the scalar field in the isospin asymmetric matter induced by a strong variation of the nucleon isovector-meson coupling constant as function of the scalar field. A possibility of hyperonization of the matter in neutron star interiors is incorporated. Our equations of state fulfill majority of known empirical constraints including the pressure-density constraint from heavy-ion collisions, direct Urca constraint, gravitational-baryon mass constraint for the pulsar J0737-3039B, and the constraint on the maximum mass of the neutron stars.

  16. Stochastic kinetic mean field model

    NASA Astrophysics Data System (ADS)

    Erdélyi, Zoltán; Pasichnyy, Mykola; Bezpalchuk, Volodymyr; Tomán, János J.; Gajdics, Bence; Gusak, Andriy M.

    2016-07-01

    This paper introduces a new model for calculating the change in time of three-dimensional atomic configurations. The model is based on the kinetic mean field (KMF) approach, however we have transformed that model into a stochastic approach by introducing dynamic Langevin noise. The result is a stochastic kinetic mean field model (SKMF) which produces results similar to the lattice kinetic Monte Carlo (KMC). SKMF is, however, far more cost-effective and easier to implement the algorithm (open source program code is provided on

  17. The Brownian mean field model

    NASA Astrophysics Data System (ADS)

    Chavanis, Pierre-Henri

    2014-05-01

    We discuss the dynamics and thermodynamics of the Brownian mean field (BMF) model which is a system of N Brownian particles moving on a circle and interacting via a cosine potential. It can be viewed as the canonical version of the Hamiltonian mean field (HMF) model. The BMF model displays a second order phase transition from a homogeneous phase to an inhomogeneous phase below a critical temperature T c = 1 / 2. We first complete the description of this model in the mean field approximation valid for N → +∞. In the strong friction limit, the evolution of the density towards the mean field Boltzmann distribution is governed by the mean field Smoluchowski equation. For T < T c , this equation describes a process of self-organization from a non-magnetized (homogeneous) phase to a magnetized (inhomogeneous) phase. We obtain an analytical expression for the temporal evolution of the magnetization close to T c . Then, we take fluctuations (finite N effects) into account. The evolution of the density is governed by the stochastic Smoluchowski equation. From this equation, we derive a stochastic equation for the magnetization and study its properties both in the homogenous and inhomogeneous phase. We show that the fluctuations diverge at the critical point so that the mean field approximation ceases to be valid. Actually, the limits N → +∞ and T → T c do not commute. The validity of the mean field approximation requires N( T - T c ) → +∞ so that N must be larger and larger as T approaches T c . We show that the direction of the magnetization changes rapidly close to T c while its amplitude takes a long time to relax. We also indicate that, for systems with long-range interactions, the lifetime of metastable states scales as e N except close to a critical point. The BMF model shares many analogies with other systems of Brownian particles with long-range interactions such as self-gravitating Brownian particles, the Keller-Segel model describing the chemotaxis

  18. Relativistic mean field calculations in neutron-rich nuclei

    SciTech Connect

    Gangopadhyay, G.; Bhattacharya, Madhubrata; Roy, Subinit

    2014-08-14

    Relativistic mean field calculations have been employed to study neutron rich nuclei. The Lagrange's equations have been solved in the co-ordinate space. The effect of the continuum has been effectively taken into account through the method of resonant continuum. It is found that BCS approximation performs as well as a more involved Relativistic Continuum Hartree Bogoliubov approach. Calculations reveal the possibility of modification of magic numbers in neutron rich nuclei. Calculation for low energy proton scattering cross sections shows that the present approach reproduces the density in very light neutron rich nuclei.

  19. Finite- to zero-range relativistic mean-field interactions

    SciTech Connect

    Niksic, T.; Vretenar, D.; Lalazissis, G. A.; Ring, P.

    2008-03-15

    We study the relation between the finite-range (meson-exchange) and zero-range (point-coupling) representations of effective nuclear interactions in the relativistic mean-field framework. Starting from the phenomenological interaction DD-ME2 with density-dependent meson-nucleon couplings, we construct a family of point-coupling effective interactions for different values of the strength parameter of the isoscalar-scalar derivative term. In the meson-exchange picture this corresponds to different values of the {sigma}-meson mass. The parameters of the isoscalar-scalar and isovector-vector channels of the point-coupling interactions are adjusted to nuclear matter and ground-state properties of finite nuclei. By comparing results for infinite and semi-infinite nuclear matter, ground-state masses, charge radii, and collective excitations, we discuss constraints on the parameters of phenomenological point-coupling relativistic effective interaction.

  20. Spins and parities of the odd-A P isotopes within a relativistic mean-field model and elastic magnetic electron-scattering theory

    NASA Astrophysics Data System (ADS)

    Wang, Zaijun; Ren, Zhongzhou; Dong, Tiekuang; Xu, Chang

    2014-08-01

    The ground-state spins and parities of the odd-A phosphorus isotopes 25-47P are studied with the relativistic mean-field (RMF) model and relativistic elastic magnetic electron-scattering theory (REMES). Results of the RMF model with the NL-SH, TM2, and NL3 parameters show that the 2s1/2 and 1d3/2 proton level inversion may occur for the neutron-rich isotopes 37-47P, and, consequently, the possible spin-parity values of 37-47P may be 3/2+, which, except for P47, differs from those given by the NUBASE2012 nuclear data table by Audi et al. Calculations of the elastic magnetic electron scattering of 37-47P with the single valence proton in the 2s1/2 and 1d3/2 state show that the form factors have significant differences. The results imply that elastic magnetic electron scattering can be a possible way to study the 2s1/2 and 1d3/2 level inversion and the spin-parity values of 37-47P. The results can also provide new tests as to what extent the RMF model, along with its various parameter sets, is valid for describing the nuclear structures. In addition, the contributions of the upper and lower components of the Dirac four-spinors to the form factors and the isotopic shifts of the magnetic form factors are discussed.

  1. Effects of δ mesons in relativistic mean field theory

    NASA Astrophysics Data System (ADS)

    Singh, Shailesh K.; Biswal, S. K.; Bhuyan, M.; Patra, S. K.

    2014-04-01

    The effect of δ- and ω-ρ-meson cross couplings on asymmetry nuclear systems are analyzed in the framework of an effective field theory motivated relativistic mean field formalism. The calculations are done on top of the G2 parameter set, where these contributions are absent. To show the effect of δ meson on the nuclear system, we split the isospin coupling into two parts: (i) gρ due to ρ meson and (ii) gδ for δ meson. Thus, our investigation is based on varying the coupling strengths of the δ and ρ mesons to reproduce the binding energies of the nuclei Ca48 and Pb208. We calculate the root mean square radius, binding energy, single particle energy, density, and spin-orbit interaction potential for some selected nuclei and evaluate the Lsym and Esym coefficients for nuclear matter as function of δ- and ω-ρ-meson coupling strengths. As expected, the influence of these effects are negligible for the symmetric nuclear system, but substantial for the contribution with large isospin asymmetry.

  2. Hot and dense matter beyond relativistic mean field theory

    NASA Astrophysics Data System (ADS)

    Zhang, Xilin; Prakash, Madappa

    2016-05-01

    Properties of hot and dense matter are calculated in the framework of quantum hadrodynamics by including contributions from two-loop (TL) diagrams arising from the exchange of isoscalar and isovector mesons between nucleons. Our extension of mean field theory (MFT) employs the same five density-independent coupling strengths which are calibrated using the empirical properties at the equilibrium density of isospin-symmetric matter. Results of calculations from the MFT and TL approximations are compared for conditions of density, temperature, and proton fraction encountered in the study of core-collapse supernovae, young and old neutron stars, and mergers of compact binary stars. The TL results for the equation of state (EOS) of cold pure neutron matter at sub- and near-nuclear densities agree well with those of modern quantum Monte Carlo and effective field-theoretical approaches. Although the high-density EOS in the TL approximation for cold and β -equilibrated neutron-star matter is substantially softer than its MFT counterpart, it is able to support a 2 M⊙ neutron star required by recent precise determinations. In addition, radii of 1.4 M⊙ stars are smaller by ˜1 km than those obtained in MFT and lie in the range indicated by analysis of astronomical data. In contrast to MFT, the TL results also give a better account of the single-particle or optical potentials extracted from analyses of medium-energy proton-nucleus and heavy-ion experiments. In degenerate conditions, the thermal variables are well reproduced by results of Landau's Fermi-liquid theory in which density-dependent effective masses feature prominently. The ratio of the thermal components of pressure and energy density expressed as Γth=1 +(Pth/ɛth) , often used in astrophysical simulations, exhibits a stronger dependence on density than on proton fraction and temperature in both MFT and TL calculations. The prominent peak of Γth at supranuclear density found in MFT is, however, suppressed in

  3. Toward relativistic mean-field description of Nbar-nucleus reactions

    NASA Astrophysics Data System (ADS)

    Gaitanos, T.; Kaskulov, M.

    2015-08-01

    In this work we study the antinucleon-nucleus optical potential in the framework of the non-linear derivative (NLD) model with momentum dependent mean-fields. We apply the NLD model to interaction of antinucleons (Nbar) in nuclear matter and, in particular, to antiproton scattering on nuclei. In nuclear matter a strong suppression of the Nbar-optical potential at rest and at high kinetic energies is found and caused by the momentum dependence of relativistic mean-fields. The NLD results are consistent with known empirical Nbar-nucleus observations and agree well with antiproton-nucleus scattering data. This makes the NLD approach compatible with both, nucleon and antinucleon Dirac phenomenologies. Furthermore, in nuclear matter an effective mass splitting between nucleons and antinucleons is predicted.

  4. Microscopically constrained mean-field models from chiral nuclear thermodynamics

    NASA Astrophysics Data System (ADS)

    Rrapaj, Ermal; Roggero, Alessandro; Holt, Jeremy W.

    2016-06-01

    We explore the use of mean-field models to approximate microscopic nuclear equations of state derived from chiral effective field theory across the densities and temperatures relevant for simulating astrophysical phenomena such as core-collapse supernovae and binary neutron star mergers. We consider both relativistic mean-field theory with scalar and vector meson exchange as well as energy density functionals based on Skyrme phenomenology and compare to thermodynamic equations of state derived from chiral two- and three-nucleon forces in many-body perturbation theory. Quantum Monte Carlo simulations of symmetric nuclear matter and pure neutron matter are used to determine the density regimes in which perturbation theory with chiral nuclear forces is valid. Within the theoretical uncertainties associated with the many-body methods, we find that select mean-field models describe well microscopic nuclear thermodynamics. As an additional consistency requirement, we study as well the single-particle properties of nucleons in a hot/dense environment, which affect e.g., charged-current weak reactions in neutron-rich matter. The identified mean-field models can be used across a larger range of densities and temperatures in astrophysical simulations than more computationally expensive microscopic models.

  5. Systematic study of Bh isotopes in a relativistic mean field formalism

    NASA Astrophysics Data System (ADS)

    Mehta, M. S.; Raj, B. K.; Patra, S. K.; Gupta, Raj K.

    2002-10-01

    The binding energy, charge radius, and quadrupole deformation parameter for the isotopic chain of the superheavy element bohrium (107Bh), from proton to neutron drip line, are calculated by using an axially deformed relativistic mean field model. The potential energy surfaces for some of the selected nuclei are plotted and the various possible shapes are investigated. The rms radii, density distributions, and two-neutron separation energies are also evaluated and the single-particle energies for some illustrative cases are analyzed to see the magic structures. Furthermore, the α-decay rates are calculated and compared with the available experimental data for the recently observed new isotopes 266,267Bh.

  6. Invisible dynamo in mean-field models

    NASA Astrophysics Data System (ADS)

    Reshetnyak, M. Yu.

    2016-07-01

    The inverse problem in a spherical shell to find the two-dimensional spatial distributions of the α-effect and differential rotation in a mean-field dynamo model has been solved. The derived distributions lead to the generation of a magnetic field concentrated inside the convection zone. The magnetic field is shown to have no time to rise from the region of maximum generation located in the lower layers to the surface in the polarity reversal time due to magnetic diffusion. The ratio of the maximum magnetic energy in the convection zone to its value at the outer boundary reaches two orders of magnitude or more. This result is important in interpreting the observed stellar and planetary magnetic fields. The proposed method of solving the inverse nonlinear dynamo problem is easily adapted for a wide class of mathematical-physics problems.

  7. Nuclear matter properties in the relativistic mean-field theory at finite temperature with interaction between sigma-omega mesons

    SciTech Connect

    Costa, R. S.; Duarte, S. B.; Oliveira, J. C. T.; Chiapparini, M.

    2010-05-21

    We study the nuclear matter properties in the regime of high temperatures using a relativistic mean-field theory. Contrasting with the usual linear Walecka model, we include the sigma-omega meson coupling in order to investigate the role of this interaction in the nucleon effective mass behavior. Some numerical results are presented and discussed.

  8. Non-mean-field critical exponent in a mean-field model: dynamics versus statistical mechanics.

    PubMed

    Ogawa, Shun; Patelli, Aurelio; Yamaguchi, Yoshiyuki Y

    2014-03-01

    Mean-field theory tells us that the classical critical exponent of susceptibility is twice that of magnetization. However, linear response theory based on the Vlasov equation, which is naturally introduced by the mean-field nature, makes the former exponent half of the latter for families of quasistationary states having second order phase transitions in the Hamiltonian mean-field model and its variances, in the low-energy phase. We clarify that this strange exponent is due to the existence of Casimir invariants which trap the system in a quasistationary state for a time scale diverging with the system size. The theoretical prediction is numerically confirmed by N-body simulations for the equilibrium states and a family of quasistationary states. PMID:24730814

  9. Multiple chiral doublet candidate nucleus {sup 105}Rh in a relativistic mean-field approach

    SciTech Connect

    Li Jian; Zhang, S. Q.; Meng, J.

    2011-03-15

    Following the reports of two pairs of chiral doublet bands observed in {sup 105}Rh, the adiabatic and configuration-fixed constrained triaxial relativistic mean-field calculations are performed to investigate their triaxial deformations with the corresponding configuration and the possible multiple chiral doublet (M{chi}D) phenomenon. The existence of the M{chi}D phenomenon in {sup 105}Rh is highly expected.

  10. Antimagnetic rotation in 108,110In with tilted axis cranking relativistic mean-field approach

    NASA Astrophysics Data System (ADS)

    Sun, Wu-Ji; Xu, Hai-Dan; Li, Jian; Liu, Yong-Hao; Ma, Ke-Yan; Yang, Dong; Lu, Jing-Bing; Ma, Ying-Jun

    2016-08-01

    Based on tilted axis cranking relativistic mean-field theory within point-coupling interaction PC-PK1, the rotational structure and the characteristic features of antimagnetic rotation for ΔI = 2 bands in 108,110In are studied. Tilted axis cranking relativistic mean-field calculations reproduce the experimental energy spectrum well and are in agreement with the experimental I ∼ ω plot, although the calculated spin overestimates the experimental values. In addition, the two-shears-like mechanism in candidate antimagnetic rotation bands is clearly illustrated and the contributions from two-shears-like orbits, neutron (gd) orbits above Z = 50 shell and Z = 50, N = 50 core are investigated microscopically. The predicted B(E2), dynamic moment of inertia ℑ(2), deformation parameters β and γ, and ℑ(2)/B(E2) ratios in tilted axis cranking relativistic mean-field calculations are discussed and the characteristic features of antimagnetic rotation for the bands before and after alignment are shown. Supported by National Natural Science Foundation of China (11205068, 11205069, 11405072, 11475072, 11547308) and China Postdoctoral Science Foundation (2012M520667)

  11. Relativistic self-consistent mean-field description of Sm isotopes

    NASA Astrophysics Data System (ADS)

    Karim, Afaque; Ahmad, Shakeb

    2015-12-01

    The evolution of the shape from the spherical to the axially deformed shapes of the neutron-rich, even-even Sm-164144 transitional nuclei is investigated. The investigations are performed with explicit density-dependent meson-nucleon and point-coupling models within the framework of the covariant density functional theory. A nonlinear meson-nucleon coupling model represented by the NL3* parametrization of the relativistic mean-field Lagrangian has also been used. The bulk and the microscopic properties of these nuclei have been investigated to analyze the phase-transition region and the critical-point behavior. The microscopic and self-consistent quadrupole deformation-constrained calculations show a clear shape change for even-even Sm isotopes with N =82 -102 . The potential energy surfaces for 148Sm,150Sm , and 152Sm obtained using different interactions are found to be relatively flat, which may be the possible critical-point nuclei. By examining the single-particle spectra, it is found that these nuclei distribute more uniformly as compared to other isotopes. Investigations also support the proposed shell-closure properties of 162Sm. Overall good agreement is found within the different models used and between the calculated and experimental results wherever available.

  12. Application of relativistic mean field and effective field theory densities to scattering observables for Ca isotopes

    NASA Astrophysics Data System (ADS)

    Bhuyan, M.; Panda, R. N.; Routray, T. R.; Patra, S. K.

    2010-12-01

    In the framework of relativistic mean field (RMF) theory, we have calculated the density distribution of protons and neutrons for Ca40,42,44,48 with NL3 and G2 parameter sets. The microscopic proton-nucleus optical potentials for p+Ca40,42,44,48 systems are evaluated from the Dirac nucleon-nucleon scattering amplitude and the density of the target nucleus using relativistic-Love-Franey and McNeil-Ray-Wallace parametrizations. We have estimated the scattering observables, such as the elastic differential scattering cross section, analyzing power and the spin observables with the relativistic impulse approximation (RIA). The results have been compared with the experimental data for a few selective cases and we find that the use of density as well as the scattering matrix parametrizations are crucial for the theoretical prediction.

  13. Application of relativistic mean field and effective field theory densities to scattering observables for Ca isotopes

    SciTech Connect

    Bhuyan, M.; Panda, R. N.; Routray, T. R.; Patra, S. K.

    2010-12-15

    In the framework of relativistic mean field (RMF) theory, we have calculated the density distribution of protons and neutrons for {sup 40,42,44,48}Ca with NL3 and G2 parameter sets. The microscopic proton-nucleus optical potentials for p+{sup 40,42,44,48}Ca systems are evaluated from the Dirac nucleon-nucleon scattering amplitude and the density of the target nucleus using relativistic-Love-Franey and McNeil-Ray-Wallace parametrizations. We have estimated the scattering observables, such as the elastic differential scattering cross section, analyzing power and the spin observables with the relativistic impulse approximation (RIA). The results have been compared with the experimental data for a few selective cases and we find that the use of density as well as the scattering matrix parametrizations are crucial for the theoretical prediction.

  14. Extended Chiral ({sigma},{pi},{omega}) Mean-Field Model with Vacuum Fluctuation Corrections

    SciTech Connect

    Uechi, Schun T.; Uechi, Hiroshi

    2011-10-21

    Density-dependent relations among saturation properties of symmetric nuclear matter and properties of hadronic stars are discussed by applying the conserving chiral nonlinear ({sigma},{pi},{omega}) mean-field theory. The chiral nonlinear ({sigma},{pi},{omega}) mean-field theory is an extension of the conserving nonlinear (nonchiral){sigma}-{omega} mean-field theory, which is thermodynamically consistent, relativistic and Lorentz-covariant. In the extended chiral ({sigma},{pi},{omega}) mean-field model, all the masses of hadrons are produced by the spontaneous chiral symmetry breaking, which is different from conventional chiral partner models. By comparing both nonchiral and chiral mean-field approximations, the effects of the chiral symmetry breaking mechanism on the mass of {sigma}-meson, coefficients of nonlinear interactions and Fermi-liquid properties are investigated in nuclear matter and neutron stars.

  15. Asymmetric neutrino production in magnetized proto-neutron stars in fully relativistic mean-field approach

    SciTech Connect

    Maruyama, Tomoyuki; Kajino, Toshitaka; Hidaka, Jun; Takiwaki, Tomoya; Yasutake, Nobutoshi; Kuroda, Takami; Cheoun, Myung-Ki; Ryu, Chung-Yeol; Mathews, Grant J.

    2014-05-02

    We calculate the neutrino production cross-section in the proto-neutron-star matter under a strong magnetic field in the relativistic mean-field approach. We introduce a new parameter-set which can reproduce the 1.96 solar mass neutron star. We find that the production process increases emitted neutrinos along the direction parallel to the magnetic field and decrease those along its opposite direction. It means that resultant asymmetry due to the neutrino absorption and scattering process in the magnetic field becomes larger by the addition of the neutrino production process.

  16. Proton and neutron skins of light nuclei within the relativistic mean field theory

    NASA Astrophysics Data System (ADS)

    Geng, L. S.; Toki, H.; Ozawa, A.; Meng, J.

    2004-01-01

    The relativistic mean field (RMF) theory is applied to the analysis of ground-state properties of Ne, Na, Cl and Ar isotopes. In particular, we study the recently established proton skin in Ar isotopes and neutron skin in Na isotopes as a function of the difference between the proton and the neutron separation energy. We use the TMA effective interaction in the RMF Lagrangian, and describe pairing correlation by the density-independent delta-function interaction. We calculate single neutron and proton separation energies, quadrupole deformations, nuclear matter radii and differences between proton radii and neutron radii, and compare these results with the recent experimental data.

  17. Multidimensionally constrained relativistic mean-field study of triple-humped barriers in actinides

    NASA Astrophysics Data System (ADS)

    Zhao, Jie; Lu, Bing-Nan; Vretenar, Dario; Zhao, En-Guang; Zhou, Shan-Gui

    2015-01-01

    Background: Potential energy surfaces (PES's) of actinide nuclei are characterized by a two-humped barrier structure. At large deformations beyond the second barrier, the occurrence of a third barrier was predicted by macroscopic-microscopic model calculations in the 1970s, but contradictory results were later reported by a number of studies that used different methods. Purpose: Triple-humped barriers in actinide nuclei are investigated in the framework of covariant density functional theory (CDFT). Methods: Calculations are performed using the multidimensionally constrained relativistic mean field (MDC-RMF) model, with the nonlinear point-coupling functional PC-PK1 and the density-dependent meson exchange functional DD-ME2 in the particle-hole channel. Pairing correlations are treated in the BCS approximation with a separable pairing force of finite range. Results: Two-dimensional PES's of 226,228,230,232Th and 232,235,236,238U are mapped and the third minima on these surfaces are located. Then one-dimensional potential energy curves along the fission path are analyzed in detail and the energies of the second barrier, the third minimum, and the third barrier are determined. The functional DD-ME2 predicts the occurrence of a third barrier in all Th nuclei and 238U . The third minima in 230 ,232Th are very shallow, whereas those in 226 ,228Th and 238U are quite prominent. With the functional PC-PK1 a third barrier is found only in 226 ,228 ,230Th . Single-nucleon levels around the Fermi surface are analyzed in 226Th, and it is found that the formation of the third minimum is mainly due to the Z =90 proton energy gap at β20≈1.5 and β30≈0.7 . Conclusions: The possible occurrence of a third barrier on the PES's of actinide nuclei depends on the effective interaction used in multidimensional CDFT calculations. More pronounced minima are predicted by the DD-ME2 functional, as compared to the functional PC-PK1. The depth of the third well in Th isotopes decreases

  18. Strange hadronic stars in relativistic mean-field theory with the FSUGold parameter set

    SciTech Connect

    Wu Chen; Ren Zhongzhou

    2011-02-15

    Relativistic mean-field theory with parameter set FSUGold that includes the isoscalar-isovector cross interaction term is extended to study the properties of neutron star matter in {beta} equilibrium by including hyperons. The influence of the attractive and repulsive {Sigma} potential on the properties of neutron star matter and the maximum mass of neutron stars is examined. We also investigate the equations of state for pure neutron matter and for nonstrange hadronic matter for comparison. For a pure neutron star, the maximum mass is about 1.8M{sub sun}, while for a strange (nonstrange) hadronic star in {beta} equilibrium, the maximum mass is around 1.35M{sub sun} (1.7M{sub sun}).

  19. Green's function method for single-particle resonant states in relativistic mean field theory

    NASA Astrophysics Data System (ADS)

    Sun, T. T.; Zhang, S. Q.; Zhang, Y.; Hu, J. N.; Meng, J.

    2014-11-01

    Relativistic mean field theory is formulated with the Green's function method in coordinate space to investigate the single-particle bound states and resonant states on the same footing. Taking the density of states for free particles as a reference, the energies and widths of single-particle resonant states are extracted from the density of states without any ambiguity. As an example, the energies and widths for single-neutron resonant states in 120Sn are compared with those obtained by the scattering phase-shift method, the analytic continuation in the coupling constant approach, the real stabilization method, and the complex scaling method. Excellent agreements with these methods are found for the energies and widths of single-neutron resonant states.

  20. MEAN FIELD AND MONTE CARLO MODELING OF MULTIBLOCK COPOLYMERS

    SciTech Connect

    K. RASMUSSEN; ET AL

    2001-01-01

    The authors discuss and apply extensions needed to treat multiblock copolymers within the mean field theoretical framework for microphase separation in diblock copolymer metals, originally due to Leibler. The mean field calculations are complemented by lattice Monte Carlo realizations using the bond fluctuation model. They find that the microphase separation transition occurs at larger {sub {chi}}N as the number of blocks in increased beyond two (i.e., beyond diblock), and that the characteristic length scale of the emerging morphology decreases as the number of blocks increases. The latter prediction is in qualitative agreement with published experimental results due to Sontak and co-workers for model multiblock poly(styrene-isoprene) systems and recent results due to Hjelm and co-workers for a segmented poly(ester-urethane) relevant to Los Alamos interests. Additionally, the mean field predictions and bond fluctuation realizations yield consistent results.

  1. Incorporating spatial correlations into multispecies mean-field models

    NASA Astrophysics Data System (ADS)

    Markham, Deborah C.; Simpson, Matthew J.; Maini, Philip K.; Gaffney, Eamonn A.; Baker, Ruth E.

    2013-11-01

    In biology, we frequently observe different species existing within the same environment. For example, there are many cell types in a tumour, or different animal species may occupy a given habitat. In modeling interactions between such species, we often make use of the mean-field approximation, whereby spatial correlations between the locations of individuals are neglected. Whilst this approximation holds in certain situations, this is not always the case, and care must be taken to ensure the mean-field approximation is only used in appropriate settings. In circumstances where the mean-field approximation is unsuitable, we need to include information on the spatial distributions of individuals, which is not a simple task. In this paper, we provide a method that overcomes many of the failures of the mean-field approximation for an on-lattice volume-excluding birth-death-movement process with multiple species. We explicitly take into account spatial information on the distribution of individuals by including partial differential equation descriptions of lattice site occupancy correlations. We demonstrate how to derive these equations for the multispecies case and show results specific to a two-species problem. We compare averaged discrete results to both the mean-field approximation and our improved method, which incorporates spatial correlations. We note that the mean-field approximation fails dramatically in some cases, predicting very different behavior from that seen upon averaging multiple realizations of the discrete system. In contrast, our improved method provides excellent agreement with the averaged discrete behavior in all cases, thus providing a more reliable modeling framework. Furthermore, our method is tractable as the resulting partial differential equations can be solved efficiently using standard numerical techniques.

  2. Kaon Condensation and Lambda-Nucleon Loop in the Relativistic Mean-Field Approach

    SciTech Connect

    Tomoyuki Maruyama; Takumi Muto; Toshitaka Tatsumi; Kazuo Tsushima; Anthony W. Thomas

    2005-02-24

    The possibility of kaon condensation in high-density symmetric nuclear matter is investigated including both s- and p-wave kaon-baryon interactions within the relativistic mean-field (RMF) theory. Above a certain density, we have a collective K{sub s} state carrying the same quantum numbers as the antikaon. The appearance of the K{sub s} state is caused by the time component of the axial-vector interaction between kaons and baryons. It is shown that the system becomes unstable with respect to condensation of K-Kbar{sub s} pairs. We consider how the effective baryon masses affect the kaon self-energy coming from the time component of the axial-vector interaction. Also, the role of the spatial component of the axial-vector interaction on the possible existence of the collective kaonic states is discussed in connection with Lambda-mixing effects in the ground state of high-density matter. Implications of K-Kbar{sub s} condensation for high-energy heavy-ion collisions are briefly mentioned.

  3. Mean-field theory of a recurrent epidemiological model

    NASA Astrophysics Data System (ADS)

    Nagy, Viktor

    2009-06-01

    Our purpose is to provide a mean-field theory for the discrete time-step susceptible-infected-recovered-susceptible (SIRS) model on uncorrelated networks with arbitrary degree distributions. The effect of network structure, time delays, and infection rate on the stability of oscillating and fixed point solutions is examined through analysis of discrete time mean-field equations. Consideration of two scenarios for disease contagion demonstrates that the manner in which contagion is transmitted from an infected individual to a contacted susceptible individual is of primary importance. In particular, the manner of contagion transmission determines how the degree distribution affects model behavior. We find excellent agreement between our theoretical results and numerical simulations on networks with large average connectivity.

  4. Mean-field theory of a recurrent epidemiological model.

    PubMed

    Nagy, Viktor

    2009-06-01

    Our purpose is to provide a mean-field theory for the discrete time-step susceptible-infected-recovered-susceptible (SIRS) model on uncorrelated networks with arbitrary degree distributions. The effect of network structure, time delays, and infection rate on the stability of oscillating and fixed point solutions is examined through analysis of discrete time mean-field equations. Consideration of two scenarios for disease contagion demonstrates that the manner in which contagion is transmitted from an infected individual to a contacted susceptible individual is of primary importance. In particular, the manner of contagion transmission determines how the degree distribution affects model behavior. We find excellent agreement between our theoretical results and numerical simulations on networks with large average connectivity. PMID:19658562

  5. Disorder Chaos in the Spherical Mean-Field Model

    NASA Astrophysics Data System (ADS)

    Chen, Wei-Kuo; Hsieh, Hsi-Wei; Hwang, Chii-Ruey; Sheu, Yuan-Chung

    2015-07-01

    We study the problem of disorder chaos in the spherical mean-field model. It concerns the behavior of the overlap between two independently sampled spin configurations from two Gibbs measures with the same external parameters. The prediction states that if the disorders in the Hamiltonians are slightly decoupled, then the overlap will be concentrated near a constant value. Following Guerra's replica symmetry breaking scheme, we establish this at the levels of the free energy and the Gibbs measure.

  6. Mean-field approach in the multi-component gas of interacting particles applied to relativistic heavy-ion collisions

    NASA Astrophysics Data System (ADS)

    Anchishkin, D.; Vovchenko, V.

    2015-10-01

    A generalized mean-field approach for the thermodynamic description of relativistic single- and multi-component gas in the grand canonical ensemble is formulated. In the framework of the proposed approach, different phenomenological excluded-volume procedures are presented and compared to the existing ones. The mean-field approach is then used to effectively include hard-core repulsion in hadron-resonance gas model for description of chemical freeze-out in heavy-ion collisions. We calculate the collision energy dependence of several quantities for different values of hard-core hadron radius and for different excluded-volume procedures such as the van der Waals and Carnahan-Starling models. It is shown that a choice of the excluded-volume model becomes important for large particle densities. For large enough values of hadron radii (r≳ 0.9 fm) there can be a sizable difference between different excluded-volume procedures used to describe the chemical freeze-out in heavy-ion collisions. At the same time, for the smaller and more commonly used values of hard-core hadron radii (r≲ 0.5 fm), the precision of the van der Waals excluded-volume procedure is shown to be sufficient.

  7. Modeling asset price processes based on mean-field framework

    NASA Astrophysics Data System (ADS)

    Ieda, Masashi; Shiino, Masatoshi

    2011-12-01

    We propose a model of the dynamics of financial assets based on the mean-field framework. This framework allows us to construct a model which includes the interaction among the financial assets reflecting the market structure. Our study is on the cutting edge in the sense of a microscopic approach to modeling the financial market. To demonstrate the effectiveness of our model concretely, we provide a case study, which is the pricing problem of the European call option with short-time memory noise.

  8. A mechanical approach to mean field spin models

    NASA Astrophysics Data System (ADS)

    Genovese, Giuseppe; Barra, Adriano

    2009-05-01

    Inspired by the bridge pioneered by Guerra among statistical mechanics on lattice and analytical mechanics on 1+1 continuous Euclidean space time, we built a self-consistent method to solve for the thermodynamics of mean field models defined on lattice, whose order parameters self-average. We show the whole procedure by analyzing in full detail the simplest test case, namely, the Curie-Weiss model. Further, we report some applications also to models whose order parameters do not self-average by using the Sherrington-Kirkpatrick spin glass as a guide.

  9. Second relativistic mean field and virial equation of state for astrophysical simulations

    SciTech Connect

    Shen, G.; Horowitz, C. J.; O'Connor, E.

    2011-06-15

    We generate a second equation of state (EOS) of nuclear matter for a wide range of temperatures, densities, and proton fractions for use in supernovae, neutron star mergers, and black hole formation simulations. We employ full relativistic mean field (RMF) calculations for matter at intermediate density and high density, and the virial expansion of a nonideal gas for matter at low density. For this EOS we use the RMF effective interaction FSUGold, whereas our earlier EOS was based on the RMF effective interaction NL3. The FSUGold interaction has a lower pressure at high densities compared to the NL3 interaction. We calculate the resulting EOS at over 100 000 grid points in the temperature range T=0 to 80 MeV, the density range n{sub B}=10{sup -8} to 1.6 fm{sup -3}, and the proton fraction range Y{sub p}=0 to 0.56. We then interpolate these data points using a suitable scheme to generate a thermodynamically consistent equation of state table on a finer grid. We discuss differences between this EOS, our NL3-based EOS, and previous EOSs by Lattimer-Swesty and H. Shen et al. for the thermodynamic properties, composition, and neutron star structure. The original FSUGold interaction produces an EOS, which we call FSU1.7, that has a maximum neutron star mass of 1.7 solar masses. A modification in the high-density EOS is introduced to increase the maximum neutron star mass to 2.1 solar masses and results in a slightly different EOS that we call FSU2.1. The EOS tables for FSU1.7 and FSU2.1 are available for download.

  10. Mean field linear response within the elimination of the small component formalism to evaluate relativistic effects on magnetic properties

    NASA Astrophysics Data System (ADS)

    Roura, P. G.; Melo, J. I.; Ruiz de Azúa, M. C.; Giribet, C. G.

    2006-08-01

    The linear response within the elimination of the small component formalism is aimed at obtaining the leading order relativistic corrections to magnetic molecular properties in the context of the elimination of the small component approximation. In the present work we extend the method in order to include two-body effects in the form of a mean field one-body operator. To this end we consider the four-component Dirac-Hartree-Fock operator as the starting point in the evaluation of the second order relativistic expression of magnetic properties. The approach thus obtained is the fully consistent leading order approximation of the random phase approximation four-component formalism. The mean field effect on the relativistic corrections to both the diamagnetic and paramagnetic terms of magnetic properties taking into account both the Coulomb and Breit two-body interactions is considered.

  11. Tensor effects in shell evolution at Z, N=8, 20, and 28 using nonrelativistic and relativistic mean-field theory

    SciTech Connect

    Moreno-Torres, M.; Anguiano, M.; Grasso, M.; Van Giai, N.; Liang, H.; De Donno, V.

    2010-06-15

    Tensor effects in shell evolution are studied within the mean-field approach. Particular attention is paid to the analysis of the magic gaps in different regions of the nuclear chart, namely, Z,N=8, 20, and 28. Hartree-Fock calculations with Skyrme and Gogny interactions are performed where the tensor term has a zero and finite range, respectively. Results obtained with and without the tensor component are compared between them and with the experimental data, when available. To complete this analysis, the tensor effect is also investigated within the relativistic Hartree-Fock model, where the exchange of rho mesons and pions is taken into account. It turns out that the tensor effect in the evolution of the magic gaps can be more easily identified in the cases Z,N=8 and 20, whereas the interpretation of the effect is more complicated for Z or N= 28. Consequently, we indicate the regions defined by the magic numbers 8 and 20 as suitable for fitting the tensor parameters in a mean-field approach: We suggest to include explicitly the data associated to these gap evolutions in the fitting procedures. In general, with the parametrizations used in this work (which have not been fitted on these data), the mean-field results obtained with the tensor contribution do not reproduce the experimental trend, that is, the reduction of the gaps at 8 and 20 that is observed when going toward the drip lines. Since some of the considered nuclei have N=Z, a discussion will be devoted to the interpretation of the experimental data concerning these nuclei and to the Wigner-energy correction.

  12. Simulated Tempering and Swapping on Mean-Field Models

    NASA Astrophysics Data System (ADS)

    Bhatnagar, Nayantara; Randall, Dana

    2016-08-01

    Simulated and parallel tempering are families of Markov Chain Monte Carlo algorithms where a temperature parameter is varied during the simulation to overcome bottlenecks to convergence due to multimodality. In this work we introduce and analyze the convergence for a set of new tempering distributions which we call entropy dampening. For asymmetric exponential distributions and the mean field Ising model with an external field simulated tempering is known to converge slowly. We show that tempering with entropy dampening distributions mixes in polynomial time for these models. Examining slow mixing times of tempering more closely, we show that for the mean-field 3-state ferromagnetic Potts model, tempering converges slowly regardless of the temperature schedule chosen. On the other hand, tempering with entropy dampening distributions converges in polynomial time to stationarity. Finally we show that the slow mixing can be very expensive practically. In particular, the mixing time of simulated tempering is an exponential factor longer than the mixing time at the fixed temperature.

  13. Mean field bipartite spin models treated with mechanical techniques

    NASA Astrophysics Data System (ADS)

    Barra, Adriano; Galluzzi, Andrea; Guerra, Francesco; Pizzoferrato, Andrea; Tantari, Daniele

    2014-03-01

    Inspired by a continuously increasing interest in modeling and framing complex systems in a thermodynamic rationale, in this paper we continue our investigation in adapting well-known techniques (originally stemmed in fields of physics and mathematics far from the present) for solving for the free energy of mean field spin models in a statistical mechanics scenario. Focusing on the test cases of bipartite spin systems embedded with all the possible interactions (self and reciprocal), we show that both the fully interacting bipartite ferromagnet, as well as the spin glass counterpart, at least at the replica symmetric level, can be solved via the fundamental theorem of calculus, trough an analogy with the Hamilton-Jacobi theory and lastly with a mapping to a Fourier diffusion problem. All these technologies are shown symmetrically for ferromagnets and spin-glasses in full details and contribute as powerful tools in the investigation of complex systems.

  14. Driven-dissipative Ising model: Mean-field solution

    NASA Astrophysics Data System (ADS)

    Goldstein, G.; Aron, C.; Chamon, C.

    2015-11-01

    We study the fate of the Ising model and its universal properties when driven by a rapid periodic drive and weakly coupled to a bath at equilibrium. The far-from-equilibrium steady-state regime is accessed by means of a Floquet mean-field approach. We show that, depending on the details of the bath, the drive can strongly renormalize the critical temperature to higher temperatures, modify the critical exponents, or even change the nature of the phase transition from second to first order after the emergence of a tricritical point. Moreover, by judiciously selecting the frequency of the field and by engineering the spectrum of the bath, one can drive a ferromagnetic Hamiltonian to an antiferromagnetically ordered phase and vice versa.

  15. Classical mutual information in mean-field spin glass models

    NASA Astrophysics Data System (ADS)

    Alba, Vincenzo; Inglis, Stephen; Pollet, Lode

    2016-03-01

    We investigate the classical Rényi entropy Sn and the associated mutual information In in the Sherrington-Kirkpatrick (S-K) model, which is the paradigm model of mean-field spin glasses. Using classical Monte Carlo simulations and analytical tools we investigate the S-K model in the n -sheet booklet. This is achieved by gluing together n independent copies of the model, and it is the main ingredient for constructing the Rényi entanglement-related quantities. We find a glassy phase at low temperatures, whereas at high temperatures the model exhibits paramagnetic behavior, consistent with the regular S-K model. The temperature of the paramagnetic-glassy transition depends nontrivially on the geometry of the booklet. At high temperatures we provide the exact solution of the model by exploiting the replica symmetry. This is the permutation symmetry among the fictitious replicas that are used to perform disorder averages (via the replica trick). In the glassy phase the replica symmetry has to be broken. Using a generalization of the Parisi solution, we provide analytical results for Sn and In and for standard thermodynamic quantities. Both Sn and In exhibit a volume law in the whole phase diagram. We characterize the behavior of the corresponding densities, Sn/N and In/N , in the thermodynamic limit. Interestingly, at the critical point the mutual information does not exhibit any crossing for different system sizes, in contrast with local spin models.

  16. Hyperdeformation in the cranked relativistic mean field theory: The Z=40-58 region of the nuclear chart

    SciTech Connect

    Afanasjev, A. V.; Abusara, H.

    2008-07-15

    The systematic investigation of hyperdeformation (HD) at high spin in the Z=40-58 region of the nuclear chart was performed in the framework of the cranked relativistic mean-field theory. The properties of the moments of inertia of the HD bands, the role of the single-particle and necking degrees of freedom at HD, the spins at which the HD bands become yrast, the possibility to observe discrete HD bands, and so on are discussed in detail.

  17. Individual based and mean-field modeling of direct aggregation

    PubMed Central

    Burger, Martin; Haškovec, Jan; Wolfram, Marie-Therese

    2013-01-01

    We introduce two models of biological aggregation, based on randomly moving particles with individual stochasticity depending on the perceived average population density in their neighborhood. In the first-order model the location of each individual is subject to a density-dependent random walk, while in the second-order model the density-dependent random walk acts on the velocity variable, together with a density-dependent damping term. The main novelty of our models is that we do not assume any explicit aggregative force acting on the individuals; instead, aggregation is obtained exclusively by reducing the individual stochasticity in response to higher perceived density. We formally derive the corresponding mean-field limits, leading to nonlocal degenerate diffusions. Then, we carry out the mathematical analysis of the first-order model, in particular, we prove the existence of weak solutions and show that it allows for measure-valued steady states. We also perform linear stability analysis and identify conditions for pattern formation. Moreover, we discuss the role of the nonlocality for well-posedness of the first-order model. Finally, we present results of numerical simulations for both the first- and second-order model on the individual-based and continuum levels of description. PMID:24926113

  18. Multidimensionally-constrained relativistic mean-field study of spontaneous fission: Coupling between shape and pairing degrees of freedom

    NASA Astrophysics Data System (ADS)

    Zhao, Jie; Lu, Bing-Nan; Nikšić, Tamara; Vretenar, Dario; Zhou, Shan-Gui

    2016-04-01

    Background: Studies of fission dynamics, based on nuclear energy density functionals, have shown that the coupling between shape and pairing degrees of freedom has a pronounced effect on the nonperturbative collective inertia and, therefore, on dynamic (least-action) spontaneous fission paths and half-lives. Purpose: The aim is to analyze the effects of particle-number fluctuation degrees of freedom on symmetric and asymmetric spontaneous fission (SF) dynamics, and to compare the findings with the results of recent studies based on the self-consistent Hartree-Fock-Bogoliubov (HFB) method. Methods: Collective potentials and nonperturbative cranking collective inertia tensors are calculated using the multidimensionally-constrained relativistic-mean-field (MDC-RMF) model. Pairing correlations are treated in the BCS approximation using a separable pairing force of finite range. Pairing fluctuations are included as a collective variable using a constraint on particle-number dispersion. Fission paths are determined with the dynamic programming method by minimizing the action in multidimensional collective spaces. Results: The dynamics of spontaneous fission of 264Fm and 250Fm are explored. Fission paths, action integrals, and corresponding half-lives computed in the three-dimensional collective space of shape and pairing coordinates, using the relativistic functional DD-PC1 and a separable pairing force of finite range, are compared with results obtained without pairing fluctuations. Results for 264Fm are also discussed in relation with those recently obtained using the HFB model. Conclusions: The inclusion of pairing correlations in the space of collective coordinates favors axially symmetric shapes along the dynamic path of the fissioning system, amplifies pairing as the path traverses the fission barriers, significantly reduces the action integral, and shortens the

  19. First principles based mean field model for oxygen reduction reaction.

    PubMed

    Jinnouchi, Ryosuke; Kodama, Kensaku; Hatanaka, Tatsuya; Morimoto, Yu

    2011-12-21

    A first principles-based mean field model was developed for the oxygen reduction reaction (ORR) taking account of the coverage- and material-dependent reversible potentials of the elementary steps. This model was applied to the simulation of single crystal surfaces of Pt, Pt alloy and Pt core-shell catalysts under Ar and O(2) atmospheres. The results are consistent with those shown by past experimental and theoretical studies on surface coverages under Ar atmosphere, the shape of the current-voltage curve for the ORR on Pt(111) and the material-dependence of the ORR activity. This model suggests that the oxygen associative pathway including HO(2)(ads) formation is the main pathway on Pt(111), and that the rate determining step (RDS) is the removal step of O(ads) on Pt(111). This RDS is accelerated on several highly active Pt alloys and core-shell surfaces, and this acceleration decreases the reaction intermediate O(ads). The increase in the partial pressure of O(2)(g) increases the surface coverage with O(ads) and OH(ads), and this coverage increase reduces the apparent reaction order with respect to the partial pressure to less than unity. This model shows details on how the reaction pathway, RDS, surface coverages, Tafel slope, reaction order and material-dependent activity are interrelated. PMID:22064886

  20. Multiagent model and mean field theory of complex auction dynamics

    NASA Astrophysics Data System (ADS)

    Chen, Qinghua; Huang, Zi-Gang; Wang, Yougui; Lai, Ying-Cheng

    2015-09-01

    Recent years have witnessed a growing interest in analyzing a variety of socio-economic phenomena using methods from statistical and nonlinear physics. We study a class of complex systems arising from economics, the lowest unique bid auction (LUBA) systems, which is a recently emerged class of online auction game systems. Through analyzing large, empirical data sets of LUBA, we identify a general feature of the bid price distribution: an inverted J-shaped function with exponential decay in the large bid price region. To account for the distribution, we propose a multi-agent model in which each agent bids stochastically in the field of winner’s attractiveness, and develop a theoretical framework to obtain analytic solutions of the model based on mean field analysis. The theory produces bid-price distributions that are in excellent agreement with those from the real data. Our model and theory capture the essential features of human behaviors in the competitive environment as exemplified by LUBA, and may provide significant quantitative insights into complex socio-economic phenomena.

  1. Modeling distributed axonal delays in mean-field brain dynamics

    NASA Astrophysics Data System (ADS)

    Roberts, J. A.; Robinson, P. A.

    2008-11-01

    The range of conduction delays between connected neuronal populations is often modeled as a single discrete delay, assumed to be an effective value averaging over all fiber velocities. This paper shows the effects of distributed delays on signal propagation. A distribution acts as a linear filter, imposing an upper frequency cutoff that is inversely proportional to the delay width. Distributed thalamocortical and corticothalamic delays are incorporated into a physiologically based mean-field model of the cortex and thalamus to illustrate their effects on the electroencephalogram (EEG). The power spectrum is acutely sensitive to the width of the thalamocortical delay distribution, and more so than the corticothalamic distribution, because all input signals must travel along the thalamocortical pathway. This imposes a cutoff frequency above which the spectrum is overly damped. The positions of spectral peaks in the resting EEG depend primarily on the distribution mean, with only weak dependences on distribution width. Increasing distribution width increases the stability of fixed point solutions. A single discrete delay successfully approximates a distribution for frequencies below a cutoff that is inversely proportional to the delay width, provided that other model parameters are moderately adjusted. A pair of discrete delays together having the same mean, variance, and skewness as the distribution approximates the distribution over the same frequency range without needing parameter adjustment. Delay distributions with large fractional widths are well approximated by low-order differential equations.

  2. Metabifurcation analysis of a mean field model of the cortex

    NASA Astrophysics Data System (ADS)

    Frascoli, Federico; van Veen, Lennaert; Bojak, Ingo; Liley, David T. J.

    2011-05-01

    Mean field models (MFMs) of cortical tissue incorporate salient, average features of neural masses in order to model activity at the population level, thereby linking microscopic physiology to macroscopic observations, e.g., with the electroencephalogram (EEG). One of the common aspects of MFM descriptions is the presence of a high-dimensional parameter space capturing neurobiological attributes deemed relevant to the brain dynamics of interest. We study the physiological parameter space of a MFM of electrocortical activity and discover robust correlations between physiological attributes of the model cortex and its dynamical features. These correlations are revealed by the study of bifurcation plots, which show that the model responses to changes in inhibition belong to two archetypal categories or “families”. After investigating and characterizing them in depth, we discuss their essential differences in terms of four important aspects: power responses with respect to the modeled action of anesthetics, reaction to exogenous stimuli such as thalamic input, and distributions of model parameters and oscillatory repertoires when inhibition is enhanced. Furthermore, while the complexity of sustained periodic orbits differs significantly between families, we are able to show how metamorphoses between the families can be brought about by exogenous stimuli. We here unveil links between measurable physiological attributes of the brain and dynamical patterns that are not accessible by linear methods. They instead emerge when the nonlinear structure of parameter space is partitioned according to bifurcation responses. We call this general method “metabifurcation analysis”. The partitioning cannot be achieved by the investigation of only a small number of parameter sets and is instead the result of an automated bifurcation analysis of a representative sample of 73,454 physiologically admissible parameter sets. Our approach generalizes straightforwardly and is well

  3. Relativistic Energy Density Functionals: beyond mean-field description of exotic structures

    SciTech Connect

    Vretenar, D.; Niksic, T.; Ring, P.; Lalazissis, G. A.

    2009-01-28

    The framework of relativistic energy density functionals is extended to include correlations related to the restoration of broken symmetries and to fluctuations of collective variables. The generator coordinate method is used to perform configuration mixing of angular-momentum and particle-number projected relativistic wave functions. This approach enables a quantitative description of the evolution of shell-structure, deformation and shape coexistence phenomena in nuclei with soft potential energy surfaces, and singular properties of excitation spectra and transition rates at critical points of quantum shape phase transitions.

  4. Relativistic mean-field study of the properties of Z=117 nuclei and the decay chains of the {sup 293,294}117 isotopes

    SciTech Connect

    Bhuyan, M.; Patra, S. K.; Gupta, Raj K.

    2011-07-15

    We have calculated the binding energy, root-mean-square radius, and quadrupole deformation parameter for the recently synthesized superheavy element Z=117, using the axially deformed relativistic mean-field (RMF) model. The calculation is extended to various isotopes of the Z=117 element, starting from A=286 till A=310. We predict almost spherical structures in the ground state for almost all the isotopes. A shape transition appears at about A=292 from a prolate to an oblate shape structure of the Z=117 nucleus in our mean-field approach. The most stable isotope (largest binding energy per nucleon) is found to be the {sup 288}117 nucleus. Also, the Q{sub {alpha}} values and the half-life T{sub 1/2}{sup {alpha}} for the {alpha}-decay chains of {sup 293}117 and {sup 294}117 are calculated, supporting the magic numbers at N=172 and/or 184.

  5. Neural Population Dynamics Modeled by Mean-Field Graphs

    NASA Astrophysics Data System (ADS)

    Kozma, Robert; Puljic, Marko

    2011-09-01

    In this work we apply random graph theory approach to describe neural population dynamics. There are important advantages of using random graph theory approach in addition to ordinary and partial differential equations. The mathematical theory of large-scale random graphs provides an efficient tool to describe transitions between high- and low-dimensional spaces. Recent advances in studying neural correlates of higher cognition indicate the significance of sudden changes in space-time neurodynamics, which can be efficiently described as phase transitions in the neuropil medium. Phase transitions are rigorously defined mathematically on random graph sequences and they can be naturally generalized to a class of percolation processes called neuropercolation. In this work we employ mean-field graphs with given vertex degree distribution and edge strength distribution. We demonstrate the emergence of collective oscillations in the style of brains.

  6. Computer program for the relativistic mean field description of the ground state properties of even-even axially deformed nuclei

    NASA Astrophysics Data System (ADS)

    Ring, P.; Gambhir, Y. K.; Lalazissis, G. A.

    1997-09-01

    We present a Fortran program for the calculation of the ground state properties of axially deformed even-even nuclei in the framework of Relativistic Mean Field Theory (RMF). In this approach a set of coupled partial differentials has to be solved self-consistently: the Dirac equation for the nucleons moving in self-consistent fields and the Klein-Gordon equations for the meson fields and the electromagnetic field, whose sources are scalar and vector densities determined of the nucleons. For this purpose the Dirac spinors as well as the meson fields are expanded in terms of anisotropic oscillator wave functions in cylindrical coordinates. This requires a matrix diagonalization for the solution of the Dirac equations and the solution of an inhomogeneous matrix equation for the meson fields. For the determination of the Coulomb field the Greens function method is used.

  7. Single-proton resonant states and the isospin dependence investigated by Green’s function relativistic mean field theory

    NASA Astrophysics Data System (ADS)

    Sun, T. T.; Niu, Z. M.; Zhang, S. Q.

    2016-08-01

    The relativistic mean field theory formulated with Green’s function method (RMF-GF) is applied to investigate single-proton resonant states and isospin dependence. The calculated energies and widths for the single-proton resonant states in {}120{{Sn}} are in good agreement with previous investigations. The single-proton resonant states of the Sn isotopes and the N = 82 isotones are systematically studied and it is shown that the calculated energies and widths decrease monotonically with the increase of neutron number while increase monotonically with the increase of proton number. To further examine the evolutions of the single-proton resonant states, their dependence on the depth, radius and diffuseness of nuclear potential is investigated with the help of an analytic Woods-Saxon potential, and it is found that the increase of radius plays the most important role in the cross phenomenon appearing in the single-proton resonant states of the Sn isotopes.

  8. Origin of the neutron skin thickness of {sup 208}Pb in nuclear mean-field models

    SciTech Connect

    Centelles, M.; Roca-Maza, X.; Vinas, X.; Warda, M.

    2010-11-15

    We study whether the neutron skin thickness {Delta}r{sub np} of {sup 208}Pb originates from the bulk or from the surface of the nucleon density distributions, according to the mean-field models of nuclear structure, and find that it depends on the stiffness of the nuclear symmetry energy. The bulk contribution to {Delta}r{sub np} arises from an extended sharp radius of neutrons, whereas the surface contribution arises from different widths of the neutron and proton surfaces. Nuclear models where the symmetry energy is stiff, as typical of relativistic models, predict a bulk contribution in {Delta}r{sub np} of {sup 208}Pb about twice as large as the surface contribution. In contrast, models with a soft symmetry energy like common nonrelativistic models predict that {Delta}r{sub np} of {sup 208}Pb is divided similarly into bulk and surface parts. Indeed, if the symmetry energy is supersoft, the surface contribution becomes dominant. We note that the linear correlation of {Delta}r{sub np} of {sup 208}Pb with the density derivative of the nuclear symmetry energy arises from the bulk part of {Delta}r{sub np}. We also note that most models predict a mixed-type (between halo and skin) neutron distribution for {sup 208}Pb. Although the halo-type limit is actually found in the models with a supersoft symmetry energy, the skin-type limit is not supported by any mean-field model. Finally, we compute parity-violating electron scattering in the conditions of the {sup 208}Pb parity radius experiment (PREX) and obtain a pocket formula for the parity-violating asymmetry in terms of the parameters that characterize the shape of the {sup 208}Pb nucleon densities.

  9. Density-Dependent Properties of Hadronic Matter in the self-consistent Chiral ({sigma},{pi},{omega}) Mean-Field Model

    SciTech Connect

    Uechi, Schun T.; Uechi, Hiroshi

    2011-05-06

    Density-dependent relations among saturation properties of symmetric nuclear matter and properties of hadronic stars are discussed by applying the conserving chiral nonlinear ({sigma},{pi},{omega}) hadronic mean-field theory. The chiral nonlinear ({sigma},{pi},{omega}) mean-field theory is an extension of the conserving nonlinear (nonchiral) {sigma}-{omega} hadronic mean-field theory which is thermodynamically consistent, relativistic and is a Lorentz-covariant mean-field theory of hadrons. In the extended chiral ({sigma},{pi},{omega}) mean-field model, all the masses of hadrons are produced by the breaking of chiral symmetry, which is different from other conventional chiral partner models. By comparing both nonchiral and chiral mean-field approximations, the effects of the chiral symmetry breaking mechanism on the mass of {sigma}-meson, coefficients of nonlinear interactions and Fermi-liquid properties are investigated in nuclear matter and neutron stars.

  10. Corrected mean-field models for spatially dependent advection-diffusion-reaction phenomena

    NASA Astrophysics Data System (ADS)

    Simpson, Matthew J.; Baker, Ruth E.

    2011-05-01

    In the exclusion-process literature, mean-field models are often derived by assuming that the occupancy status of lattice sites is independent. Although this assumption is questionable, it is the foundation of many mean-field models. In this work we develop methods to relax the independence assumption for a range of discrete exclusion-process-based mechanisms motivated by applications from cell biology. Previous investigations that focused on relaxing the independence assumption have been limited to studying initially uniform populations and ignored any spatial variations. By ignoring spatial variations these previous studies were greatly simplified due to translational invariance of the lattice. These previous corrected mean-field models could not be applied to many important problems in cell biology such as invasion waves of cells that are characterized by moving fronts. Here we propose generalized methods that relax the independence assumption for spatially inhomogeneous problems, leading to corrected mean-field descriptions of a range of exclusion-process-based models that incorporate (i) unbiased motility, (ii) biased motility, and (iii) unbiased motility with agent birth and death processes. The corrected mean-field models derived here are applicable to spatially variable processes including invasion wave-type problems. We show that there can be large deviations between simulation data and traditional mean-field models based on invoking the independence assumption. Furthermore, we show that the corrected mean-field models give an improved match to the simulation data in all cases considered.

  11. Configuration Interaction as an Impurity Solver: Benchmark Dynamical Mean-Field Theory for the Hubbard Model

    NASA Astrophysics Data System (ADS)

    Go, Ara; Millis, Andrew J.

    2013-03-01

    The configuration interaction technique has been widely used in quantum chemistry to solve quantum many body systems with lower computational costs than exact diagonalization and was introduced by Dominika Zgid, Emanuel Gull, and Garnet Kin-Lic Chan [Phys. Rev. B 86, 165128 (2012)] as a solver for the impurity models of dynamical mean field theory. We extend their work, demonstrating for the one and two dimensional Hubbard model how the method reproduces the known results and allows convergence with bath size to be studied in cluster dynamical mean field theory. As an example of the power of the method, cluster dynamical mean field studies of the three band copper-oxygen model are presented. This work was supported by the CMCSN program of the US Department of Energy.

  12. {alpha}-decay properties of superheavy elements Z=113-125 in the relativistic mean-field theory with vector self-coupling of {omega} meson

    SciTech Connect

    Sharma, M.M.; Farhan, A.R.; Muenzenberg, G.

    2005-05-01

    We have investigated properties of {alpha}-decay chains of recently produced superheavy elements Z=115 and Z=113 using the new Lagrangian model NL-SV1 with inclusion of the vector self-coupling of the {omega} meson in the framework of relativistic mean-field theory. It is shown that the experimentally observed {alpha}-decay energies and half-lives are reproduced well by this Lagrangian model. Further calculations for the heavier elements with Z=117-125 show that these nuclei are superdeformed with a prolate shape in the ground state. A superdeformed shell closure at Z=118 lends an additional binding and an extra stability to nuclei in this region. Consequently, it is predicted that the corresponding Q{sub {alpha}} values provide {alpha}-decay half-lives for heavier superheavy nuclei within experimentally feasible conditions. The results are compared with those of macroscopic-microscopic approaches. A perspective of the difference in shell effects among various approaches is presented and its consequences for superheavy nuclei are discussed.

  13. a Mean-Field Version of the Ssb Model for X-Chromosome Inactivation

    NASA Astrophysics Data System (ADS)

    Gaeta, Giuseppe

    Nicodemi and Prisco recently proposed a model for X-chromosome inactivation in mammals, explaining this phenomenon in terms of a spontaneous symmetry-breaking mechanism [{\\it Phys. Rev. Lett.} 99 (2007), 108104]. Here we provide a mean-field version of their model.

  14. Magnetohydrodynamic Simulation-driven Kinematic Mean Field Model of the Solar Cycle

    NASA Astrophysics Data System (ADS)

    Simard, Corinne; Charbonneau, Paul; Bouchat, Amélie

    2013-05-01

    We construct a series of kinematic axisymmetric mean-field dynamo models operating in the αΩ, α2Ω and α2 regimes, all using the full α-tensor extracted from a global magnetohydrodynamical simulation of solar convection producing large-scale magnetic fields undergoing solar-like cyclic polarity reversals. We also include an internal differential rotation profile produced in a purely hydrodynamical parent simulation of solar convection, and a simple meridional flow profile described by a single cell per meridional quadrant. An α2Ω mean-field model, presumably closest to the mode of dynamo action characterizing the MHD simulation, produces a spatiotemporal evolution of magnetic fields that share some striking similarities with the zonally-averaged toroidal component extracted from the simulation. Comparison with α2 and αΩ mean-field models operating in the same parameter regimes indicates that much of the complexity observed in the spatiotemporal evolution of the large-scale magnetic field in the simulation can be traced to the turbulent electromotive force. Oscillating α2 solutions are readily produced, and show some similarities with the observed solar cycle, including a deep-seated toroidal component concentrated at low latitudes and migrating equatorward in the course of the solar cycle. Various numerical experiments performed using the mean-field models reveal that turbulent pumping plays an important role in setting the global characteristics of the magnetic cycles.

  15. Mean-field models for heterogeneous networks of two-dimensional integrate and fire neurons

    PubMed Central

    Nicola, Wilten; Campbell, Sue Ann

    2013-01-01

    We analytically derive mean-field models for all-to-all coupled networks of heterogeneous, adapting, two-dimensional integrate and fire neurons. The class of models we consider includes the Izhikevich, adaptive exponential and quartic integrate and fire models. The heterogeneity in the parameters leads to different moment closure assumptions that can be made in the derivation of the mean-field model from the population density equation for the large network. Three different moment closure assumptions lead to three different mean-field systems. These systems can be used for distinct purposes such as bifurcation analysis of the large networks, prediction of steady state firing rate distributions, parameter estimation for actual neurons and faster exploration of the parameter space. We use the mean-field systems to analyze adaptation induced bursting under realistic sources of heterogeneity in multiple parameters. Our analysis demonstrates that the presence of heterogeneity causes the Hopf bifurcation associated with the emergence of bursting to change from sub-critical to super-critical. This is confirmed with numerical simulations of the full network for biologically reasonable parameter values. This change decreases the plausibility of adaptation being the cause of bursting in hippocampal area CA3, an area with a sizable population of heavily coupled, strongly adapting neurons. PMID:24416013

  16. MAGNETOHYDRODYNAMIC SIMULATION-DRIVEN KINEMATIC MEAN FIELD MODEL OF THE SOLAR CYCLE

    SciTech Connect

    Simard, Corinne; Charbonneau, Paul; Bouchat, Amelie E-mail: paulchar@astro.umontreal.ca

    2013-05-01

    We construct a series of kinematic axisymmetric mean-field dynamo models operating in the {alpha}{Omega}, {alpha}{sup 2}{Omega} and {alpha}{sup 2} regimes, all using the full {alpha}-tensor extracted from a global magnetohydrodynamical simulation of solar convection producing large-scale magnetic fields undergoing solar-like cyclic polarity reversals. We also include an internal differential rotation profile produced in a purely hydrodynamical parent simulation of solar convection, and a simple meridional flow profile described by a single cell per meridional quadrant. An {alpha}{sup 2}{Omega} mean-field model, presumably closest to the mode of dynamo action characterizing the MHD simulation, produces a spatiotemporal evolution of magnetic fields that share some striking similarities with the zonally-averaged toroidal component extracted from the simulation. Comparison with {alpha}{sup 2} and {alpha}{Omega} mean-field models operating in the same parameter regimes indicates that much of the complexity observed in the spatiotemporal evolution of the large-scale magnetic field in the simulation can be traced to the turbulent electromotive force. Oscillating {alpha}{sup 2} solutions are readily produced, and show some similarities with the observed solar cycle, including a deep-seated toroidal component concentrated at low latitudes and migrating equatorward in the course of the solar cycle. Various numerical experiments performed using the mean-field models reveal that turbulent pumping plays an important role in setting the global characteristics of the magnetic cycles.

  17. Self-consistent chaos in a mean-field Hamiltonian model of fluids and plasmas

    NASA Astrophysics Data System (ADS)

    del-Castillo-Negrete, D.; Firpo, Marie-Christine

    2002-11-01

    We present a mean-field Hamiltonian model that describes the collective dynamics of marginally stable fluids and plasmas. In plasmas, the model describes the self-consistent evolution of electron holes and clumps in phase space. In fluids, the model describes the dynamics of vortices with negative and positive circulation in shear flows. The mean-field nature of the system makes it a tractable model to study the dynamics of large degrees-of-freedom, coupled Hamiltonian systems. Here we focus in the role of self-consistent chaos in the formation and destruction of phase space coherent structures. Numerical simulations in the finite N and in the Narrow kinetic limit (where N is the number of particles) show the existence of coherent, rotating dipole states. We approximate the dipole as two macroparticles, and show that the N = 2 limit has a family of rotating integrable solutions described by a one degree-of-freedom nontwist Hamiltonian. The coherence of the dipole is explained in terms of a parametric resonance between the rotation frequency of the macroparticles and the oscillation frequency of the self-consistent mean field. For a class of initial conditions, the mean field exhibits a self-consistent, elliptic-hyperbolic bifurcation that leads to the destruction of the dipole and violent mixing of the phase space.

  18. Mean-field cluster model for the critical behaviour of ferromagnets

    NASA Astrophysics Data System (ADS)

    Chamberlin, Ralph V.

    2000-11-01

    Two separate theories are often used to characterize the paramagnetic properties of ferromagnetic materials. At temperatures T well above the Curie temperature, TC (where the transition from paramagnetic to ferromagnetic behaviour occurs), classical mean-field theory yields the Curie-Weiss law for the magnetic susceptibility: χ( T) ~ 1/(T - Θ), where Θ is the Weiss constant. Close to TC, however, the standard mean-field approach breaks down so that better agreement with experimental data is provided by critical scaling theory: χ(T) ~ 1/(T - TC)γ , where γ is a scaling exponent. But there is no known model capable of predicting the measured values of γ nor its variation among different substances. Here I use a mean-field cluster model based on finite-size thermostatistics to extend the range of mean-field theory, thereby eliminating the need for a separate scaling regime. The mean-field approximation is justified by using a kinetic-energy term to maintain the microcanonical ensemble. The model reproduces the Curie-Weiss law at high temperatures, but the classical Weiss transition at TC = Θ is suppressed by finite-size effects. Instead, the fraction of clusters with a specific amount of order diverges at T C, yielding a transition that is mathematically similar to Bose-Einstein condensation. At all temperatures above T C, the model matches the measured magnetic susceptibilities of crystalline EuO, Gd, Co and Ni, thus providing a unified picture for both the critical-scaling and Curie-Weiss regimes.

  19. Oscillating epidemics in a dynamic network model: stochastic and mean-field analysis.

    PubMed

    Szabó-Solticzky, András; Berthouze, Luc; Kiss, Istvan Z; Simon, Péter L

    2016-04-01

    An adaptive network model using SIS epidemic propagation with link-type-dependent link activation and deletion is considered. Bifurcation analysis of the pairwise ODE approximation and the network-based stochastic simulation is carried out, showing that three typical behaviours may occur; namely, oscillations can be observed besides disease-free or endemic steady states. The oscillatory behaviour in the stochastic simulations is studied using Fourier analysis, as well as through analysing the exact master equations of the stochastic model. By going beyond simply comparing simulation results to mean-field models, our approach yields deeper insights into the observed phenomena and help better understand and map out the limitations of mean-field models. PMID:26063525

  20. Mean field theory of the linear sigma-model: chiral solitons

    SciTech Connect

    Kahana, S.; Ripka, G.

    1983-01-01

    The mean field theory of the chiral invariant sigma-model is outlined. bound states (solitons) of valence quarks are obtained self-consistently using a hedgehog shape for the pion field. A schematic model for the coupled fermion-boson fields is presented. Renormalization is worked out for the fermion one-loop corrections and numerical results presented for the purely scalar-field case. The interpretation of the baryon number of the perturbed vacuum is considered.

  1. Spin-orbit and orbit-orbit strengths for the radioactive neutron-rich doubly magic nucleus {sup 132}Sn in relativistic mean-field theory

    SciTech Connect

    Liang Haozhao; Zhao Pengwei; Li Lulu; Meng Jie

    2011-01-15

    Relativistic mean-field (RMF) theory is applied to investigate the properties of the radioactive neutron-rich doubly magic nucleus {sup 132}Sn and the corresponding isotopes and isotones. The two-neutron and two-proton separation energies are well reproduced by the RMF theory. In particular, the RMF results agree with the experimental single-particle spectrum in {sup 132}Sn as well as the Nilsson spin-orbit parameter C and orbit-orbit parameter D thus extracted, but remarkably differ from the traditional Nilsson parameters. Furthermore, the present results provide a guideline for the isospin dependence of the Nilsson parameters.

  2. Modeling of coherent ultrafast magneto-optical experiments: Light-induced molecular mean-field model

    SciTech Connect

    Hinschberger, Y.; Hervieux, P.-A.

    2015-12-28

    We present calculations which aim to describe coherent ultrafast magneto-optical effects observed in time-resolved pump-probe experiments. Our approach is based on a nonlinear semi-classical Drude-Voigt model and is used to interpret experiments performed on nickel ferromagnetic thin film. Within this framework, a phenomenological light-induced coherent molecular mean-field depending on the polarizations of the pump and probe pulses is proposed whose microscopic origin is related to a spin-orbit coupling involving the electron spins of the material sample and the electric field of the laser pulses. Theoretical predictions are compared to available experimental data. The model successfully reproduces the observed experimental trends and gives meaningful insight into the understanding of magneto-optical rotation behavior in the ultrafast regime. Theoretical predictions for further experimental studies are also proposed.

  3. Core-halo quasi-stationary states in the Hamiltonian mean-field model

    NASA Astrophysics Data System (ADS)

    Konishi, Eiji

    2016-04-01

    A characteristic feature of long-range interacting systems is that they become trapped in a non-equilibrium and long-lived quasi-stationary state (QSS) during the early stages of their development. We present a comprehensive review of recent studies of the core-halo structure of QSSs, in the Hamiltonian mean-field model (HMF), which is a mean-field model of mutually coupled ferromagnetic XY spins located at a point, obtained by starting from various unsteady rectangular water-bag type initial phase-space distributions. The main result exposed in this review is that the core-halo structure can be described by the superposition of two independent Lynden-Bell distributions. We discuss the completeness of collisionless relaxation of this double Lynden-Bell distribution by using both of Lynden-Bell entropy and double Lynden-Bell entropy for the systems at low energies per particle.

  4. Avalanche-size distributions in mean-field plastic yielding models.

    PubMed

    Jagla, E A

    2015-10-01

    We discuss the size distribution N(S) of avalanches occurring at the yielding transition of mean-field (i.e., Hebraud-Lequeux) models of amorphous solids. The size distribution follows a power law dependence of the form N(S)∼S(-τ). However (contrary to what is found in its depinning counterpart), the value of τ depends on details of the dynamic protocol used. For random triggering of avalanches we recover the τ=3/2 exponent typical of mean-field models, which, in particular, is valid for the depinning case. However, for the physically relevant case of external loading through a quasistatic increase of applied strain, a smaller exponent (close to 1) is obtained. This result is rationalized by mapping the problem to an effective random walk in the presence of a moving absorbing boundary. PMID:26565196

  5. Synchronization and Spin-Flop Transitions for a Mean-Field XY Model in Random Field

    NASA Astrophysics Data System (ADS)

    Collet, Francesca; Ruszel, Wioletta

    2016-08-01

    We characterize the phase space for the infinite volume limit of a ferromagnetic mean-field XY model in a random field pointing in one direction with two symmetric values. We determine the stationary solutions and detect possible phase transitions in the interaction strength for fixed random field intensity. We show that at low temperature magnetic ordering appears perpendicularly to the field. The latter situation corresponds to a spin-flop transition.

  6. Inverse magnetic catalysis in Nambu-Jona-Lasinio model beyond mean field

    NASA Astrophysics Data System (ADS)

    Mao, Shijun

    2016-07-01

    We study inverse magnetic catalysis in the Nambu-Jona-Lasinio model beyond mean field approximation. The feed-down from mesons to quarks is embedded in an effective coupling constant at finite temperature and magnetic field. While the magnetic catalysis is still the dominant effect at low temperature, the meson dressed quark mass drops down with increasing magnetic field at high temperature due to the dimension reduction of the Goldstone mode in the Pauli-Villars regularization scheme.

  7. An exactly solvable spherical mean-field plus extended monopole pairing model

    NASA Astrophysics Data System (ADS)

    Pan, Feng; Ding, Xiaoxue; Launey, Kristina D.; Li, Hui; Xu, Xinxin; Draayer, Jerry P.

    2016-03-01

    An extended pairing Hamiltonian that describes pairing interactions among monopole nucleon pairs up to an infinite order in a spherical mean field, such as the spherical shell model, is proposed based on the local E˜2 algebraic structure, which includes the extended pairing interaction within a deformed mean-field theory (Pan et al., 2004) [19] as a special case. The advantage of the model lies in the fact that numerical solutions of the model can be obtained more easily and with less computational time than the solutions to the standard pairing model. Thus, open-shell large-scale calculations within the model become feasible. As an example of the application, pairing contribution to the binding energy of 12-28O is estimated in the present model with neutron pairs allowed to occupy a no-core shell model space of 11 j-orbits up to the fifth major harmonic oscillator shell including excitations up to 14 ħω for 12O and up to 40 ħω for 28O. The results for 12O are also compared and found to be in agreement with those of ab initio calculations. It is shown that the pairing energy per particle in 12-28O ranges from 0.4 to 1.8 MeV/A with the strongest one observed for a small number of pairs.

  8. Dynamical Aspects of Mean Field Plane Rotators and the Kuramoto Model

    NASA Astrophysics Data System (ADS)

    Bertini, Lorenzo; Giacomin, Giambattista; Pakdaman, Khashayar

    2010-02-01

    The Kuramoto model has been introduced in order to describe synchronization phenomena observed in groups of cells, individuals, circuits, etc. We look at the Kuramoto model with white noise forces: in mathematical terms it is a set of N oscillators, each driven by an independent Brownian motion with a constant drift, that is each oscillator has its own frequency, which, in general, changes from one oscillator to another (these frequencies are usually taken to be random and they may be viewed as a quenched disorder). The interactions between oscillators are of long range type (mean field). We review some results on the Kuramoto model from a statistical mechanics standpoint: we give in particular necessary and sufficient conditions for reversibility and we point out a formal analogy, in the N→∞ limit, with local mean field models with conservative dynamics (an analogy that is exploited to identify in particular a Lyapunov functional in the reversible set-up). We then focus on the reversible Kuramoto model with sinusoidal interactions in the N→∞ limit and analyze the stability of the non-trivial stationary profiles arising when the interaction parameter K is larger than its critical value K c . We provide an analysis of the linear operator describing the time evolution in a neighborhood of the synchronized profile: we exhibit a Hilbert space in which this operator has a self-adjoint extension and we establish, as our main result, a spectral gap inequality for every K> K c .

  9. Mean-field approximation for a limit order driven market model.

    PubMed

    Slanina, F

    2001-11-01

    A mean-field variant of the model of limit order driven market introduced recently by Maslov is formulated and solved. The agents do not have any strategies and the memory of the system is kept within the order book. We show that the evolution of the order book is governed by a matrix multiplicative process. The resulting stationary distribution of step-to-step price changes is calculated. It exhibits a power-law tail with exponent 2. We obtain also the price autocorrelation function, which agrees qualitatively with the experimentally observed negative autocorrelation for short times. PMID:11736043

  10. Bose-Hubbard models in confining potentials: Inhomogeneous mean-field theory

    NASA Astrophysics Data System (ADS)

    Pai, Ramesh V.; Kurdestany, Jamshid Moradi; Sheshadri, K.; Pandit, Rahul

    2012-06-01

    We present an extensive study of Mott insulator (MI) and superfluid (SF) shells in Bose-Hubbard (BH) models for bosons in optical lattices with harmonic traps. For this we apply the inhomogeneous mean-field theory developed by Sheshadri [Phys. Rev. Lett.PRLTAO0031-900710.1103/PhysRevLett.75.4075 75, 4075 (1995)]. Our results for the BH model with one type of spinless bosons agree quantitatively with quantum Monte Carlo simulations. Our approach is numerically less intensive than such simulations, so we are able to perform calculations on experimentally realistic, large three-dimensional systems, explore a wide range of parameter values, and make direct contact with a variety of experimental measurements. We also extend our inhomogeneous mean-field theory to study BH models with harmonic traps and (a) two species of bosons or (b) spin-1 bosons. With two species of bosons, we obtain rich phase diagrams with a variety of SF and MI phases and associated shells when we include a quadratic confining potential. For the spin-1 BH model, we show, in a representative case, that the system can display alternating shells of polar SF and MI phases, and we make interesting predictions for experiments in such systems.

  11. Mean-field thalamocortical modeling of longitudinal EEG acquired during intensive meditation training.

    PubMed

    Saggar, Manish; Zanesco, Anthony P; King, Brandon G; Bridwell, David A; MacLean, Katherine A; Aichele, Stephen R; Jacobs, Tonya L; Wallace, B Alan; Saron, Clifford D; Miikkulainen, Risto

    2015-07-01

    Meditation training has been shown to enhance attention and improve emotion regulation. However, the brain processes associated with such training are poorly understood and a computational modeling framework is lacking. Modeling approaches that can realistically simulate neurophysiological data while conforming to basic anatomical and physiological constraints can provide a unique opportunity to generate concrete and testable hypotheses about the mechanisms supporting complex cognitive tasks such as meditation. Here we applied the mean-field computational modeling approach using the scalp-recorded electroencephalogram (EEG) collected at three assessment points from meditating participants during two separate 3-month-long shamatha meditation retreats. We modeled cortical, corticothalamic, and intrathalamic interactions to generate a simulation of EEG signals recorded across the scalp. We also present two novel extensions to the mean-field approach that allow for: (a) non-parametric analysis of changes in model parameter values across all channels and assessments; and (b) examination of variation in modeled thalamic reticular nucleus (TRN) connectivity over the retreat period. After successfully fitting whole-brain EEG data across three assessment points within each retreat, two model parameters were found to replicably change across both meditation retreats. First, after training, we observed an increased temporal delay between modeled cortical and thalamic cells. This increase provides a putative neural mechanism for a previously observed reduction in individual alpha frequency in these same participants. Second, we found decreased inhibitory connection strength between the TRN and secondary relay nuclei (SRN) of the modeled thalamus after training. This reduction in inhibitory strength was found to be associated with increased dynamical stability of the model. Altogether, this paper presents the first computational approach, taking core aspects of physiology and

  12. Economic dynamics with financial fragility and mean-field interaction: A model

    NASA Astrophysics Data System (ADS)

    Di Guilmi, C.; Gallegati, M.; Landini, S.

    2008-06-01

    Following Aoki’s statistical mechanics methodology [Masanao Aoki, New Approaches to Macroeconomic Modeling, Cambridge University Press, 1996; Masanao Aoki, Modeling Aggregate Behaviour and Fluctuations in Economics, Cambridge University Press, 2002; Masanao Aoki, and Hiroshi Yoshikawa, Reconstructing Macroeconomics, Cambridge University Press, 2006], we provide some insights into the well-known works of [Bruce Greenwald, Joseph Stiglitz, Macroeconomic models with equity and credit rationing, in: R. Hubbard (Ed.), Information, Capital Markets and Investment, Chicago University Press, Chicago, 1990; Bruce Greenwald, Joseph Stiglitz, Financial markets imperfections and business cycles, Quarterly journal of Economics (1993)]. Specifically, we reach analytically a closed form solution of their models overcoming the aggregation problem. The key idea is to represent the economy as an evolving complex system, composed by heterogeneous interacting agents, that can be partitioned into a space of macroscopic states. This meso level of aggregation permits to adopt mean-field interaction modeling and master equation techniques.

  13. Mean-Field Theory of the Symmetry Breaking Model for X Chromosome Inactivation

    NASA Astrophysics Data System (ADS)

    Scialdone, A.; Barbieri, M.; Pallotti, D.; Nicodemi, M.

    X Chromosome Inactivation (XCI) is the process in mammal femalecells whereby one of the X chromosomes is silenced to compensate dosage with respect to males. It is still mysterious how precisely one X chromosome is randomly chosen for inactivation. We discuss here a mean-field theory of the Symmetry Breaking (SB) model of XCI, a Statistical Mechanics model introduced to explain that process. The SB model poses that a single regulatory factor, an aggregate of molecules, is produced which acts to preserve from inactivation one of the X's. The model illustrates a physical mechanism, originating from a thermodynamic phase transition, for the self-assembling of such a single super-molecular aggregate which can spontaneously break the binding symmetry of equivalent targets. This results in a sharp, yet stochastic, regulatory mechanism of XCI. In particular, we focus here on how the model can predict the effects of genetic deletions.

  14. Octupole degree of freedom for the critical-point candidate nucleus {sup 152}Sm in a reflection-asymmetric relativistic mean-field approach

    SciTech Connect

    Zhang, W.; Li, Z. P.; Zhang, S. Q.; Meng, J.

    2010-03-15

    The potential energy surfaces of even-even {sup 146-156}Sm are investigated in the constrained reflection-asymmetric relativistic mean-field approach with parameter set PK1. It is shown that the critical-point candidate nucleus {sup 152}Sm marks the shape/phase transition not only from U(5) to SU(3) symmetry, but also from the octupole-deformed ground state in {sup 150}Sm to the quadrupole-deformed ground state in {sup 154}Sm. By including the octupole degree of freedom, an energy gap near the Fermi surface for single-particle levels in {sup 152}Sm with beta{sub 2}=0.14approx0.26 is found and the important role of the octupole deformation driving pair nu2f{sub 7/2} and nu1i{sub 13/2} is demonstrated.

  15. Structural Properties of the Disordered Spherical and Other Mean Field Spin Models

    NASA Astrophysics Data System (ADS)

    Sanctis, Luca De

    2007-03-01

    We extend the approach of Aizenman, Sims and Starr for the SK-type models to their spherical versions. Such an extension has already been performed for diluted spin glasses. The factorization property of the optimal structures found by Guerra for the SK model, which holds for diluted models as well, is verified also in the case of spherical systems, with the due modifications. Hence we show that there are some common structural features in various mean field spin models. These similarities seem to be quite paradigmatic, and we summarize the various techniques typically used to prove the structural analogies and to tackle the computation of the free energy per spin in the thermodynamic limit.

  16. Bayesian Learning in Sparse Graphical Factor Models via Variational Mean-Field Annealing

    PubMed Central

    Yoshida, Ryo; West, Mike

    2010-01-01

    We describe a class of sparse latent factor models, called graphical factor models (GFMs), and relevant sparse learning algorithms for posterior mode estimation. Linear, Gaussian GFMs have sparse, orthogonal factor loadings matrices, that, in addition to sparsity of the implied covariance matrices, also induce conditional independence structures via zeros in the implied precision matrices. We describe the models and their use for robust estimation of sparse latent factor structure and data/signal reconstruction. We develop computational algorithms for model exploration and posterior mode search, addressing the hard combinatorial optimization involved in the search over a huge space of potential sparse configurations. A mean-field variational technique coupled with annealing is developed to successively generate “artificial” posterior distributions that, at the limiting temperature in the annealing schedule, define required posterior modes in the GFM parameter space. Several detailed empirical studies and comparisons to related approaches are discussed, including analyses of handwritten digit image and cancer gene expression data. PMID:20890391

  17. Mean-field approximation for the Sznajd model in complex networks

    NASA Astrophysics Data System (ADS)

    Araújo, Maycon S.; Vannucchi, Fabio S.; Timpanaro, André M.; Prado, Carmen P. C.

    2015-02-01

    This paper studies the Sznajd model for opinion formation in a population connected through a general network. A master equation describing the time evolution of opinions is presented and solved in a mean-field approximation. Although quite simple, this approximation allows us to capture the most important features regarding the steady states of the model. When spontaneous opinion changes are included, a discontinuous transition from consensus to polarization can be found as the rate of spontaneous change is increased. In this case we show that a hybrid mean-field approach including interactions between second nearest neighbors is necessary to estimate correctly the critical point of the transition. The analytical prediction of the critical point is also compared with numerical simulations in a wide variety of networks, in particular Barabási-Albert networks, finding reasonable agreement despite the strong approximations involved. The same hybrid approach that made it possible to deal with second-order neighbors could just as well be adapted to treat other problems such as epidemic spreading or predator-prey systems.

  18. Modeling and computation of mean field equilibria in producers' game with emission permits trading

    NASA Astrophysics Data System (ADS)

    Zhang, Shuhua; Wang, Xinyu; Shanain, Aleksandr

    2016-08-01

    In this paper, we present a mean field game to model the production behaviors of a very large number of producers, whose carbon emissions are regulated by government. Especially, an emission permits trading scheme is considered in our model, in which each enterprise can trade its own permits flexibly. By means of the mean field equilibrium, we obtain a Hamilton-Jacobi-Bellman (HJB) equation coupled with a Kolmogorov equation, which are satisfied by the adjoint state and the density of producers (agents), respectively. Then, we propose a so-called fitted finite volume method to solve the HJB equation and the Kolmogorov equation. The efficiency and the usefulness of this method are illustrated by the numerical experiments. Under different conditions, the equilibrium states as well as the effects of the emission permits price are examined, which demonstrates that the emission permits trading scheme influences the producers' behaviors, that is, more populations would like to choose a lower rather than a higher emission level when the emission permits are expensive.

  19. Mean-field diffusion-limited aggregation: a "density" model for viscous fingering phenomena.

    PubMed

    Bogoyavlenskiy, V A

    2001-12-01

    We explore a universal "density" formalism to describe nonequilibrium growth processes, specifically, the immiscible viscous fingering in Hele-Shaw cells (usually referred to as the Saffman-Taylor problem). For that we develop an alternative approach to the viscous fingering phenomena, whose basic concepts have been recently published in a Rapid Communication [Phys. Rev. E 63, 045305(R) (2001)]. This approach uses the diffusion-limited aggregation (DLA) paradigm as a core: we introduce a mean-field DLA generalization in stochastic and deterministic formulations. The stochastic model, a quasicontinuum DLA, simulates Monte Carlo patterns, which demonstrate a striking resemblance to natural Hele-Shaw fingers and, for steady-state growth regimes, follow precisely the Saffman-Taylor analytical solutions in channel and sector configurations. The relevant deterministic theory, a complete set of differential equations for a time development of density fields, is derived from that stochastic model. As a principal conclusion, we prove an asymptotic equivalency of both the stochastic and deterministic mean-field DLA formulations to the classic Saffman-Taylor hydrodynamics in terms of an interface evolution. PMID:11736272

  20. Mean-field model for the growth and coarsening of stoichiometric precipitates at grain boundaries

    NASA Astrophysics Data System (ADS)

    Kozeschnik, E.; Svoboda, J.; Radis, R.; Fischer, F. D.

    2010-01-01

    In this paper, a model for growth and coarsening of precipitates at grain boundaries is developed. The concept takes into account that the evolution of grain boundary precipitates involves fast short-circuit diffusion along grain boundaries as well as slow bulk diffusion of atoms from the grain interior to the grain boundaries. The mathematical formalism is based on a mean-field approximation, utilizing the thermodynamic extremal principle. The model is applied to the precipitation of aluminum nitrides in microalloyed steel in austenite, where precipitation occurs predominately at the austenite grain boundaries. It is shown that the kinetics of precipitation predicted by the proposed model differs significantly from that calculated for randomly distributed precipitates with spherical diffusion fields. Good agreement of the numerical solution is found with experimental observations as well as theoretical treatment of precipitate coarsening.

  1. Ground state phase transition in the Nilsson mean-field plus standard pairing model

    NASA Astrophysics Data System (ADS)

    Guan, Xin; Xu, Haocheng; Zhang, Yu; Pan, Feng; Draayer, Jerry P.

    2016-08-01

    The ground state phase transition in Nd, Sm, and Gd isotopes is investigated by using the Nilsson mean-field plus standard pairing model based on the exact solutions obtained from the extended Heine-Stieltjes correspondence. The results of the model calculations successfully reproduce the critical phenomena observed experimentally in the odd-even mass differences, odd-even differences of two-neutron separation energy, and the α -decay and double β--decay energies of these isotopes. Since the odd-even effects are the most important signatures of pairing interactions in nuclei, the model calculations yield microscopic insight into the nature of the ground state phase transition manifested by the standard pairing interaction.

  2. Quantum correlated cluster mean-field theory applied to the transverse Ising model

    NASA Astrophysics Data System (ADS)

    Zimmer, F. M.; Schmidt, M.; Maziero, Jonas

    2016-06-01

    Mean-field theory (MFT) is one of the main available tools for analytical calculations entailed in investigations regarding many-body systems. Recently, there has been a surge of interest in ameliorating this kind of method, mainly with the aim of incorporating geometric and correlation properties of these systems. The correlated cluster MFT (CCMFT) is an improvement that succeeded quite well in doing that for classical spin systems. Nevertheless, even the CCMFT presents some deficiencies when applied to quantum systems. In this article, we address this issue by proposing the quantum CCMFT (QCCMFT), which, in contrast to its former approach, uses general quantum states in its self-consistent mean-field equations. We apply the introduced QCCMFT to the transverse Ising model in honeycomb, square, and simple cubic lattices and obtain fairly good results both for the Curie temperature of thermal phase transition and for the critical field of quantum phase transition. Actually, our results match those obtained via exact solutions, series expansions or Monte Carlo simulations.

  3. Effects of anisotropies in turbulent magnetic diffusion in mean-field solar dynamo models

    SciTech Connect

    Pipin, V. V.; Kosovichev, A. G.

    2014-04-10

    We study how anisotropies of turbulent diffusion affect the evolution of large-scale magnetic fields and the dynamo process on the Sun. The effect of anisotropy is calculated in a mean-field magnetohydrodynamics framework assuming that triple correlations provide relaxation to the turbulent electromotive force (so-called the 'minimal τ-approximation'). We examine two types of mean-field dynamo models: the well-known benchmark flux-transport model and a distributed-dynamo model with a subsurface rotational shear layer. For both models, we investigate effects of the double- and triple-cell meridional circulation, recently suggested by helioseismology and numerical simulations. To characterize the anisotropy effects, we introduce a parameter of anisotropy as a ratio of the radial and horizontal intensities of turbulent mixing. It is found that the anisotropy affects the distribution of magnetic fields inside the convection zone. The concentration of the magnetic flux near the bottom and top boundaries of the convection zone is greater when the anisotropy is stronger. It is shown that the critical dynamo number and the dynamo period approach to constant values for large values of the anisotropy parameter. The anisotropy reduces the overlap of toroidal magnetic fields generated in subsequent dynamo cycles, in the time-latitude 'butterfly' diagram. If we assume that sunspots are formed in the vicinity of the subsurface shear layer, then the distributed dynamo model with the anisotropic diffusivity satisfies the observational constraints from helioseismology and is consistent with the value of effective turbulent diffusion estimated from the dynamics of surface magnetic fields.

  4. Mean-field behavior of the negative-weight percolation model on random regular graphs.

    PubMed

    Melchert, Oliver; Hartmann, Alexander K; Mézard, Marc

    2011-10-01

    We investigate both analytically and numerically the ensemble of minimum-weight loops in the negative-weight percolation model on random graphs with fixed connectivity and bimodal weight distribution. This allows us to study the mean-field behavior of this model. The analytical study is based on a conjectured equivalence with the problem of self-avoiding walks in a random medium. The numerical study is based on a mapping to a standard minimum-weight matching problem for which fast algorithms exist. Both approaches yield results that are in agreement on the location of the phase transition, on the value of critical exponents, and on the absence of any sizable indications of a glass phase. By these results, the previously conjectured upper critical dimension of d(u)=6 is confirmed. PMID:22181086

  5. Dynamics and Self-consistent Chaos in a Mean Field Hamiltonian Model

    NASA Astrophysics Data System (ADS)

    del-Castillo-Negrete, Diego

    We study a mean field Hamiltonian model that describes the collective dynamics of marginally stable fluids and plasmas in the finite N and N-> infty kinetic limit (where N is the number of particles). The linear stability of equilibria in the kinetic model is studied as well as the initial value problem including Landau damping . Numerical simulations show the existence of coherent, rotating dipole states. We approximate the dipole as two macroparticles and show that the N=2 limit has a family of rotating integrable solutions that provide an accurate description of the dynamics. We discuss the role of self-consistent Hamiltonian chaos in the formation of coherent structures, and discuss a mechanism of "violent" mixing caused by a self-consistent elliptic-hyperbolic bifurcation in phase space.

  6. Angular momentum projection for a Nilsson mean-field plus pairing model

    NASA Astrophysics Data System (ADS)

    Wang, Yin; Pan, Feng; Launey, Kristina D.; Luo, Yan-An; Draayer, J. P.

    2016-06-01

    The angular momentum projection for the axially deformed Nilsson mean-field plus a modified standard pairing (MSP) or the nearest-level pairing (NLP) model is proposed. Both the exact projection, in which all intrinsic states are taken into consideration, and the approximate projection, in which only intrinsic states with K = 0 are taken in the projection, are considered. The analysis shows that the approximate projection with only K = 0 intrinsic states seems reasonable, of which the configuration subspace considered is greatly reduced. As simple examples for the model application, low-lying spectra and electromagnetic properties of 18O and 18Ne are described by using both the exact and approximate angular momentum projection of the MSP or the NLP, while those of 20Ne and 24Mg are described by using the approximate angular momentum projection of the MSP or NLP.

  7. Mean-field dynamic criticality and geometric transition in the Gaussian core model

    NASA Astrophysics Data System (ADS)

    Coslovich, Daniele; Ikeda, Atsushi; Miyazaki, Kunimasa

    2016-04-01

    We use molecular dynamics simulations to investigate dynamic heterogeneities and the potential energy landscape of the Gaussian core model (GCM). Despite the nearly Gaussian statistics of particles' displacements, the GCM exhibits giant dynamic heterogeneities close to the dynamic transition temperature. The divergence of the four-point susceptibility is quantitatively well described by the inhomogeneous version of the mode-coupling theory. Furthermore, the potential energy landscape of the GCM is characterized by large energy barriers, as expected from the lack of activated, hopping dynamics, and display features compatible with a geometric transition. These observations demonstrate that all major features of mean-field dynamic criticality can be observed in a physically sound, three-dimensional model.

  8. Mean-field dynamic criticality and geometric transition in the Gaussian core model.

    PubMed

    Coslovich, Daniele; Ikeda, Atsushi; Miyazaki, Kunimasa

    2016-04-01

    We use molecular dynamics simulations to investigate dynamic heterogeneities and the potential energy landscape of the Gaussian core model (GCM). Despite the nearly Gaussian statistics of particles' displacements, the GCM exhibits giant dynamic heterogeneities close to the dynamic transition temperature. The divergence of the four-point susceptibility is quantitatively well described by the inhomogeneous version of the mode-coupling theory. Furthermore, the potential energy landscape of the GCM is characterized by large energy barriers, as expected from the lack of activated, hopping dynamics, and display features compatible with a geometric transition. These observations demonstrate that all major features of mean-field dynamic criticality can be observed in a physically sound, three-dimensional model. PMID:27176347

  9. Thermodynamics of the one-dimensional parallel Kawasaki model: Exact solution and mean-field approximations

    NASA Astrophysics Data System (ADS)

    Pazzona, Federico G.; Demontis, Pierfranco; Suffritti, Giuseppe B.

    2014-08-01

    The adsorption isotherm for the recently proposed parallel Kawasaki (PK) lattice-gas model [Phys. Rev. E 88, 062144 (2013), 10.1103/PhysRevE.88.062144] is calculated exactly in one dimension. To do so, a third-order difference equation for the grand-canonical partition function is derived and solved analytically. In the present version of the PK model, the attraction and repulsion effects between two neighboring particles and between a particle and a neighboring empty site are ruled, respectively, by the dimensionless parameters ϕ and θ. We discuss the inflections induced in the isotherms by situations of high repulsion, the role played by finite lattice sizes in the emergence of substeps, and the adequacy of the two most widely used mean-field approximations in lattice gases, namely, the Bragg-Williams and the Bethe-Peierls approximations.

  10. Stability of inhomogeneous superstructures from renormalized mean-field theory of the t-J model

    NASA Astrophysics Data System (ADS)

    Poilblanc, Didier

    2005-08-01

    Using the t-J model (which can also include Coulomb repulsion) and the “plain vanilla” renormalized mean-field theory of Zhang, [Supercond. Sci. Technol. 1, 36 (1988)], stability of inhomogeneous 4a×4a superstructures, such as those observed in cuprates superconductors around 1/8 hole doping is investigated. We find a nonuniform 4a×4a bond order wave involving simultaneously small (˜10-2t) inhomogeneous staggered plaquette currents as well as a small charge-density modulation similar to pair density wave order. On the other hand, no supersolid phase involving a decoupling in the superconducting particle-particle channel is found.

  11. Improving Solar Cycle Prediction Using Variational Data Assimilation in a Mean-Field Dynamo Model

    NASA Astrophysics Data System (ADS)

    Fournier, A.; Hung, C. P.; Jouve, L.; Brun, S.

    2014-12-01

    We present our recent effort to implement modern variational data assimilation techniques into a mean field solar dynamo code. This work extends the work of (Jouve et al. 2011, ApJ) to take into account the correct spherical geometry and meridional circulation into so-called Babcock-Leigthon flux transport dynamo models. Based on twin-experiments, in which we observe our dynamo simulations, and on a well defined cost function using toroidal and poloidal field observations we are able to recover the main attributes of the dynamo solution used to test our data assimilation algorithm. By assimilating solar data (such as Wolf number or butterfly diagram) we are starting to deduce the profile and temporal variations of key ingredients of the solar dynamo. We find that the data sampling and the temporal window are key to get reliable results. We show how such a powerful technique can be used to improve our ability to predict the solar magnetic activity.

  12. Metastates in Mean-Field Models with Random External Fields Generated by Markov Chains

    NASA Astrophysics Data System (ADS)

    Formentin, M.; Külske, C.; Reichenbachs, A.

    2012-01-01

    We extend the construction by Külske and Iacobelli of metastates in finite-state mean-field models in independent disorder to situations where the local disorder terms are a sample of an external ergodic Markov chain in equilibrium. We show that for non-degenerate Markov chains, the structure of the theorems is analogous to the case of i.i.d. variables when the limiting weights in the metastate are expressed with the aid of a CLT for the occupation time measure of the chain. As a new phenomenon we also show in a Potts example that for a degenerate non-reversible chain this CLT approximation is not enough, and that the metastate can have less symmetry than the symmetry of the interaction and a Gaussian approximation of disorder fluctuations would suggest.

  13. Mean field study of the topological Haldane-Hubbard model of spin-1/2 fermions

    NASA Astrophysics Data System (ADS)

    Arun, V. S.; Sohal, R.; Hickey, C.; Paramekanti, A.

    2016-03-01

    Motivated by exploring the effect of interactions on Chern insulators, and by recent experiments realizing topological bands for ultracold atoms in synthetic gauge fields, we study the honeycomb lattice Haldane-Hubbard model of spin-1/2 fermions. Using an unrestricted mean field approach, we map out the instability of the topological band insulator towards magnetically ordered insulators which emerge with increasing Hubbard repulsion. In addition to the topological Néel phase, we recover various chiral noncoplanar magnetic orders previously identified within a strong-coupling approach. We compute the band structure of these ordered phases, showing that the triple-Q tetrahedral phase harbors topological Chern bands with large Chern numbers.

  14. Mean field theory for biology inspired duplication-divergence network model

    NASA Astrophysics Data System (ADS)

    Cai, Shuiming; Liu, Zengrong; Lee, H. C.

    2015-08-01

    The duplication-divergence network model is generally thought to incorporate key ingredients underlying the growth and evolution of protein-protein interaction networks. Properties of the model have been elucidated through numerous simulation studies. However, a comprehensive theoretical study of the model is lacking. Here, we derived analytic expressions for quantities describing key characteristics of the network—the average degree, the degree distribution, the clustering coefficient, and the neighbor connectivity—in the mean-field, large-N limit of an extended version of the model, duplication-divergence complemented with heterodimerization and addition. We carried out extensive simulations and verified excellent agreement between simulation and theory except for one partial case. All four quantities obeyed power-laws even at moderate network size ( N ˜104 ), except the degree distribution, which had an additional exponential factor observed to obey power-law. It is shown that our network model can lead to the emergence of scale-free property and hierarchical modularity simultaneously, reproducing the important topological properties of real protein-protein interaction networks.

  15. Mean field theory for biology inspired duplication-divergence network model.

    PubMed

    Cai, Shuiming; Liu, Zengrong; Lee, H C

    2015-08-01

    The duplication-divergence network model is generally thought to incorporate key ingredients underlying the growth and evolution of protein-protein interaction networks. Properties of the model have been elucidated through numerous simulation studies. However, a comprehensive theoretical study of the model is lacking. Here, we derived analytic expressions for quantities describing key characteristics of the network-the average degree, the degree distribution, the clustering coefficient, and the neighbor connectivity-in the mean-field, large-N limit of an extended version of the model, duplication-divergence complemented with heterodimerization and addition. We carried out extensive simulations and verified excellent agreement between simulation and theory except for one partial case. All four quantities obeyed power-laws even at moderate network size ( N∼10(4)), except the degree distribution, which had an additional exponential factor observed to obey power-law. It is shown that our network model can lead to the emergence of scale-free property and hierarchical modularity simultaneously, reproducing the important topological properties of real protein-protein interaction networks. PMID:26328557

  16. Algebraic Mean Field Theory

    NASA Astrophysics Data System (ADS)

    Dankova, T. S.; Rosensteel, G.

    1998-10-01

    Mean field theory has an unexpected group theoretic mathematical foundation. Instead of representation theory which applies to most group theoretic quantum models, Hartree-Fock and Hartree-Fock-Bogoliubov have been formulated in terms of coadjoint orbits for the groups U(n) and O(2n). The general theory of mean fields is formulated for an arbitrary Lie algebra L of fermion operators. The moment map provides the correspondence between the Hilbert space of microscopic wave functions and the dual space L^* of densities. The coadjoint orbits of the group in the dual space are phase spaces on which time-dependent mean field theory is equivalent to a classical Hamiltonian dynamical system. Indeed it forms a finite-dimensional Lax system. The mean field theories for the Elliott SU(3) and symplectic Sp(3,R) algebras are constructed explicitly in the coadjoint orbit framework.

  17. A stochastic mean field model for an excitatory and inhibitory synaptic drive cortical neuronal network.

    PubMed

    Hui, Qing; Haddad, Wassim M; Bailey, James M; Hayakawa, Tomohisa

    2014-04-01

    With the advances in biochemistry, molecular biology, and neurochemistry there has been impressive progress in understanding the molecular properties of anesthetic agents. However, there has been little focus on how the molecular properties of anesthetic agents lead to the observed macroscopic property that defines the anesthetic state, that is, lack of responsiveness to noxious stimuli. In this paper, we develop a mean field synaptic drive firing rate cortical neuronal model and demonstrate how the induction of general anesthesia can be explained using multistability; the property whereby the solutions of a dynamical system exhibit multiple attracting equilibria under asymptotically slowly changing inputs or system parameters. In particular, we demonstrate multistability in the mean when the system initial conditions or the system coefficients of the neuronal connectivity matrix are random variables. Uncertainty in the system coefficients is captured by representing system uncertain parameters by a multiplicative white noise model wherein stochastic integration is interpreted in the sense of Itô. Modeling a priori system parameter uncertainty using a multiplicative white noise model is motivated by means of the maximum entropy principle of Jaynes and statistical analysis. PMID:24807952

  18. Self-consistent mean-field model for palmitoyloleoylphosphatidylcholine-palmitoyl sphingomyelin-cholesterol lipid bilayers

    NASA Astrophysics Data System (ADS)

    Tumaneng, Paul W.; Pandit, Sagar A.; Zhao, Guijun; Scott, H. L.

    2011-03-01

    The connection between membrane inhomogeneity and the structural basis of lipid rafts has sparked interest in the lateral organization of model lipid bilayers of two and three components. In an effort to investigate anisotropic lipid distribution in mixed bilayers, a self-consistent mean-field theoretical model is applied to palmitoyloleoylphosphatidylcholine (POPC)-palmitoyl sphingomyelin (PSM)-cholesterol mixtures. The compositional dependence of lateral organization in these mixtures is mapped onto a ternary plot. The model utilizes molecular dynamics simulations to estimate interaction parameters and to construct chain conformation libraries. We find that at some concentration ratios the bilayers separate spatially into regions of higher and lower chain order coinciding with areas enriched with PSM and POPC, respectively. To examine the effect of the asymmetric chain structure of POPC on bilayer lateral inhomogeneity, we consider POPC-lipid interactions with and without angular dependence. Results are compared with experimental data and with results from a similar model for mixtures of dioleoylphosphatidylcholine, steroyl sphingomyelin, and cholesterol.

  19. Higgs-Yukawa model with higher dimension operators via extended mean field theory

    NASA Astrophysics Data System (ADS)

    Akerlund, Oscar; de Forcrand, Philippe

    2016-02-01

    Using extended mean field theory (EMFT) on the lattice, we study properties of the Higgs-Yukawa model as an approximation of the standard model Higgs sector, and the effect of higher dimension operators. We remark, as has been noted before, that the discussion of vacuum stability is completely modified in the presence of a ϕ6 term, and that the Higgs mass no longer appears fine tuned. We also study the finite temperature transition. Without higher dimension operators the transition is found to be second order (crossover with gauge fields) for the experimental value of the Higgs mass Mh=125 GeV . By taking a ϕ6 interaction in the Higgs potential as a proxy for a UV completion of the standard model, the transition becomes stronger and turns first order if the scale of new physics, i.e., the mass of the lightest mediator particle, is around 1.5 TeV. This implies that electroweak baryogenesis may be viable in models which introduce new particles around that scale.

  20. Spin-orbit interaction in relativistic nuclear structure models

    NASA Astrophysics Data System (ADS)

    Ebran, J.-P.; Mutschler, A.; Khan, E.; Vretenar, D.

    2016-08-01

    Relativistic self-consistent mean-field (SCMF) models naturally account for the coupling of the nucleon spin to its orbital motion, whereas nonrelativistic SCMF methods necessitate a phenomenological ansatz for the effective spin-orbit potential. Recent experimental studies aim to explore the isospin properties of the effective spin-orbit interaction in nuclei. SCMF models are very useful in the interpretation of the corresponding data; however, standard relativistic mean-field and nonrelativistic Hartree-Fock models use effective spin-orbit potentials with different isovector properties, mainly because exchange contributions are not treated explicitly in the former. The impact of exchange terms on the effective spin-orbit potential in relativistic mean-field models is analyzed, and it is shown that it leads to an isovector structure similar to the one used in standard nonrelativistic Hartree-Fock models. Data on the isospin dependence of spin-orbit splittings in spherical nuclei could be used to constrain the isovector-scalar channel of relativistic mean-field models. The reproduction of the empirical kink in the isotope shifts of even Pb nuclei by relativistic effective interactions points to the occurrence of pseudospin symmetry in the single-neutron spectra in these nuclei.

  1. Novel dynamics and thermodynamics of a new Hamiltonian mean field model

    NASA Astrophysics Data System (ADS)

    Curilef, Sergio; Atenas, Boris

    Statistical systems are idealized by the hypothesis that the particles do not interact among them, or the range of interactions is short enough, reaching very fast the statistical state that we know as equilibrium. However, systems with long-range interactions are common in nature because of they are observed from the atomic scale to the astronomical scale, exhibiting some anomalies as inequivalence of ensembles, negative heat capacity, ergodicity breaking, non equilibrium phase transitions, quasi-stationarity, anomalous diffusion, etc. We present in this contribution a new Hamiltonian mean field model whose potential is inspired in the dipole-dipole interactions. The equilibrium is analytically studied in the canonical ensemble and coincides with the one obtained from molecular dynamics simulations (microcanonical ensemble). We notice, this model presents a kind of inequivalence of ensembles in long-standing states before arriving at equilibrium. However, the novelty, compared to other models presented in recent literature, is that two quasi-stationary states appear in the behavior of this system. The first quasi-stationary state decays to a second one, which is different to the first, before going to the equilibrium.We characterize them by its dynamics and thermodynamics. We acknowledge partial financial support by Anillo ACT-1204, VRIDT-UCN105/2015. We appreciate the computational assistance of A. Pluchino.

  2. Critical behavior of a tumor growth model: directed percolation with a mean-field flavor.

    PubMed

    Lipowski, Adam; Ferreira, António Luis; Wendykier, Jacek

    2012-10-01

    We examine the critical behavior of a lattice model of tumor growth where supplied nutrients are correlated with the distribution of tumor cells. Our results support the previous report [Ferreira et al., Phys. Rev. E 85, 010901(R) (2012)], which suggested that the critical behavior of the model differs from the expected directed percolation (DP) universality class. Surprisingly, only some of the critical exponents (β, α, ν([perpendicular]), and z) take non-DP values while some others (β', ν(||), and spreading-dynamics exponents Θ, δ, z') remain very close to their DP counterparts. The obtained exponents satisfy the scaling relations β=αν(||), β'=δν(||), and the generalized hyperscaling relation Θ+α+δ=d/z, where the dynamical exponent z is, however, used instead of the spreading exponent z'. Both in d=1 and d=2 versions of our model, the exponent β most likely takes the mean-field value β=1, and we speculate that it might be due to the roulette-wheel selection, which is used to choose the site to supply a nutrient. PMID:23214560

  3. Mean-field theory for Bose-Hubbard model under a magnetic field

    SciTech Connect

    Oktel, M. Oe.; Tanatar, B.; Nita, M.

    2007-01-15

    We consider the superfluid-insulator transition for cold bosons under an effective magnetic field. We investigate how the applied magnetic field affects the Mott transition within mean-field theory and find that the critical hopping strength (t/U){sub c} increases with the applied field. The increase in the critical hopping follows the bandwidth of the Hofstadter butterfly at the given value of the magnetic field. We also calculate the magnetization and superfluid density within mean-field theory.

  4. General model of phospholipid bilayers in fluid phase within the single chain mean field theory

    SciTech Connect

    Guo, Yachong; Baulin, Vladimir A.; Pogodin, Sergey

    2014-05-07

    Coarse-grained model for saturated phospholipids: 1,2-didecanoyl-sn-glycero-3-phosphocholine (DCPC), 1,2-dilauroyl-sn-glycero-3-phosphocholine (DLPC), 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC), 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC), 1,2-distearoyl-sn-glycero-3-phosphocholine (DSPC) and unsaturated phospholipids: 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC), 1,2- dioleoyl-sn-glycero-3-phosphocholine (DOPC) is introduced within the single chain mean field theory. A single set of parameters adjusted for DMPC bilayers gives an adequate description of equilibrium and mechanical properties of a range of saturated lipid molecules that differ only in length of their hydrophobic tails and unsaturated (POPC, DOPC) phospholipids which have double bonds in the tails. A double bond is modeled with a fixed angle of 120°, while the rest of the parameters are kept the same as saturated lipids. The thickness of the bilayer and its hydrophobic core, the compressibility, and the equilibrium area per lipid correspond to experimentally measured values for each lipid, changing linearly with the length of the tail. The model for unsaturated phospholipids also fetches main thermodynamical properties of the bilayers. This model is used for an accurate estimation of the free energies of the compressed or stretched bilayers in stacks or multilayers and gives reasonable estimates for free energies. The proposed model may further be used for studies of mixtures of lipids, small molecule inclusions, interactions of bilayers with embedded proteins.

  5. Hamiltonian mean field model: Effect of network structure on synchronization dynamics.

    PubMed

    Virkar, Yogesh S; Restrepo, Juan G; Meiss, James D

    2015-11-01

    The Hamiltonian mean field model of coupled inertial Hamiltonian rotors is a prototype for conservative dynamics in systems with long-range interactions. We consider the case where the interactions between the rotors are governed by a network described by a weighted adjacency matrix. By studying the linear stability of the incoherent state, we find that the transition to synchrony begins when the coupling constant K is inversely proportional to the largest eigenvalue of the adjacency matrix. We derive a closed system of equations for a set of local order parameters to study the effect of network heterogeneity on the synchronization of the rotors. When K is just beyond the transition to synchronization, we find that the degree of synchronization is highly dependent on the network's heterogeneity, but that for large K the degree of synchronization is robust to changes in the degree distribution. Our results are illustrated with numerical simulations on Erdös-Renyi networks and networks with power-law degree distributions. PMID:26651739

  6. Correlation patterns of NIKKEI index constituents. Towards a mean-field model

    NASA Astrophysics Data System (ADS)

    Hayashi, Katsuhiko; Kaizoji, Taisei; Pichl, Lukáš

    2007-09-01

    An analysis of minute-tick data from the Japanese stock index market is reported for a three-year period of 2000/7/4-2003/6/30. Correlation patterns and principal component distributions were determined for 180 constituents of the NIKKEI 225 index, excluding the effects of after-hours trading and constituent revisions. The first principal component describes about 30% of the total variance in constituent log returns (subject to slow decrease with the size of the correlation window), suggesting that a small number of physical parameters may describe the internal dynamics of the index, allowing for an adiabatic representation of index dynamics, and a self-consistent mean-field model of its constituents. Finally, it is shown that the introduction of a time gap into minute-tick data significantly improves the correlations of the price-weighed index with its constituents, even when such gap inserts are strictly penalized. This phenomenon corresponds to a heterogenous response time of index constituents to the adiabatic collective motion and also demonstrates the inhomogeneous nature of equidistant time ticks in financial trading.

  7. Neutron stars in the Bogomol'nyi-Prasad-Sommerfield Skyrme model: Mean-field limit versus full field theory

    NASA Astrophysics Data System (ADS)

    Adam, C.; Naya, C.; Sanchez-Guillen, J.; Vazquez, R.; Wereszczynski, A.

    2015-08-01

    Using a solitonic model of nuclear matter, the Bogomol'nyi-Prasad-Sommerfield (BPS) Skyrme model, we compare neutron stars obtained in the full field theory, where gravitational backreaction is completely taken into account, with calculations in a mean-field approximation using the Tolman-Oppenheimer-Volkoff approach. In the latter case, a mean-field-theory equation of state is derived from the original BPS field theory. We show that in the full field theory, where the energy density is nonconstant even at equilibrium, there is no universal and coordinate-independent equation of state of nuclear matter, in contrast to the mean-field approximation. We also study how neutron star properties are modified by going beyond mean-field theory and find that the differences between mean-field theory and exact results can be considerable. Further, we compare both exact and mean-field results with some theoretical and phenomenological constraints on neutron star properties, demonstrating thus the relevance of our model even in its most simple version.

  8. Structural versus dynamical origins of mean-field behavior in a self-organized critical model of neuronal avalanches.

    PubMed

    Moosavi, S Amin; Montakhab, Afshin

    2015-11-01

    Critical dynamics of cortical neurons have been intensively studied over the past decade. Neuronal avalanches provide the main experimental as well as theoretical tools to consider criticality in such systems. Experimental studies show that critical neuronal avalanches show mean-field behavior. There are structural as well as recently proposed [Phys. Rev. E 89, 052139 (2014)] dynamical mechanisms that can lead to mean-field behavior. In this work we consider a simple model of neuronal dynamics based on threshold self-organized critical models with synaptic noise. We investigate the role of high-average connectivity, random long-range connections, as well as synaptic noise in achieving mean-field behavior. We employ finite-size scaling in order to extract critical exponents with good accuracy. We conclude that relevant structural mechanisms responsible for mean-field behavior cannot be justified in realistic models of the cortex. However, strong dynamical noise, which can have realistic justifications, always leads to mean-field behavior regardless of the underlying structure. Our work provides a different (dynamical) origin than the conventionally accepted (structural) mechanisms for mean-field behavior in neuronal avalanches. PMID:26651741

  9. Structural versus dynamical origins of mean-field behavior in a self-organized critical model of neuronal avalanches

    NASA Astrophysics Data System (ADS)

    Moosavi, S. Amin; Montakhab, Afshin

    2015-11-01

    Critical dynamics of cortical neurons have been intensively studied over the past decade. Neuronal avalanches provide the main experimental as well as theoretical tools to consider criticality in such systems. Experimental studies show that critical neuronal avalanches show mean-field behavior. There are structural as well as recently proposed [Phys. Rev. E 89, 052139 (2014), 10.1103/PhysRevE.89.052139] dynamical mechanisms that can lead to mean-field behavior. In this work we consider a simple model of neuronal dynamics based on threshold self-organized critical models with synaptic noise. We investigate the role of high-average connectivity, random long-range connections, as well as synaptic noise in achieving mean-field behavior. We employ finite-size scaling in order to extract critical exponents with good accuracy. We conclude that relevant structural mechanisms responsible for mean-field behavior cannot be justified in realistic models of the cortex. However, strong dynamical noise, which can have realistic justifications, always leads to mean-field behavior regardless of the underlying structure. Our work provides a different (dynamical) origin than the conventionally accepted (structural) mechanisms for mean-field behavior in neuronal avalanches.

  10. Modeling relativistic nuclear collisions.

    SciTech Connect

    Anderlik, C.; Magas, V.; Strottman, D.; Csernai, L. P.

    2001-01-01

    Modeling Ultra-Relativistic Heavy Ion Collisioiis at RHIC and LHC energies using a Multi Module Model is presented. The first Module is the Effective String Rope Model for the calculation of the initial stages of the reaction; the output of this module is used as the initial state for the subsequent one-fluid hydrodynainical calculation module. It is shown that such an initial state leads to the creation of the third flow component. The hydrodynamical evolution of the energy density distribution is presented for RHIC energies. The final module describing the Freeze Out; and Hadronization is also discussed.

  11. Dynamical instabilities in density-dependent hadronic relativistic models

    SciTech Connect

    Santos, A. M.; Brito, L.; Providencia, C.

    2008-04-15

    Unstable modes in asymmetric nuclear matter (ANM) at subsaturation densities are studied in the framework of relativistic mean-field density-dependent hadron models. The size of the instabilities that drive the system are calculated and a comparison with results obtained within the nonlinear Walecka model is presented. The distillation and antidistillation effects are discussed.

  12. Mean-Field Theory is Exact for the Random-Field Model with Long-Range Interactions

    NASA Astrophysics Data System (ADS)

    Tsuda, Junichi; Nishimori, Hidetoshi

    2014-07-01

    We study the classical spin model in random fields with long-range interactions and show the exactness of the mean-field theory under certain mild conditions. This is a generalization of the result of Mori for the non-random and spin-glass cases. To treat random fields, we evoke the self-averaging property of a function of random fields, without recourse to the replica method. The result is that the mean-field theory gives the exact expression of the canonical free energy for systems with power-decaying interactions if the power is smaller than or equal to the spatial dimension.

  13. Mean-field plus various types of pairing models and an exact boson mapping of the standard pairing model

    SciTech Connect

    Pan Feng; Wang Yin; Guan Xin; Jia Lu; Chen Xiangrong; Draayer, J. P.

    2011-06-28

    Exact solutions of Nilsson mean-field with various pairing interactions are reviewed. Some even-odd mass differences and moments of inertia of low-lying states for rare earth and actinide nuclei are calculated for the nearest-orbit pairing approximation as well as for the extended pairing model and compared to available experimental data. An exact boson mapping of the standard pairing Hamiltonian is also reported. Under the mapping, fermion pair operators are mapped exactly onto corresponding bosons. The image of the mapping is a Bose-Hubbard model with orbit-dependent hopping.

  14. Beyond-mean-field boson-fermion model for odd-mass nuclei

    NASA Astrophysics Data System (ADS)

    Nomura, K.; Nikšić, T.; Vretenar, D.

    2016-05-01

    A novel method for calculating spectroscopic properties of medium-mass and heavy atomic nuclei with an odd number of nucleons is introduced, based on the framework of nuclear energy density functional theory and the particle-core coupling scheme. The deformation energy surface of the even-even core, as well as the spherical single-particle energies and occupation probabilities of the odd particle(s), are obtained in a self-consistent mean-field calculation determined by the choice of the energy density functional and pairing interaction. This method uniquely determines the parameters of the Hamiltonian of the boson core, and only the strength of the particle-core coupling is specifically adjusted to selected data for a particular nucleus. The approach is illustrated in a systematic study of low-energy excitation spectra and transition rates of axially deformed odd-mass Eu isotopes.

  15. Analytic models of regularly branched polymer brushes using the self-consistent mean field theory

    NASA Astrophysics Data System (ADS)

    LeSher, Daniel

    2015-10-01

    Polymer brushes consist of multiple monomers connected together with one of the polymer chain's ends attached to a surface. Polymer brushes have shown great promise for a wide variety of applications including drug delivery dendrimer systems and as tunable brushes that can change their shape and physical properties in response to changes in their environment. Regularly branched polymer brushes which are structured as a function of their chemical indices are investigated here using the self-consistent mean field theory for electrically neutral polymers. The brushes were described using weighting functions, f(n), were n was the fewest number of monomers from a specified location to a free end. Brushes with weighting functions of the form f(n)=nb, f(n)=ebn, as well as f(n)=dan when d 2 and alpha > 2 were found to match the parabolic free chain end profile expected, while it was determined that polymer brushes described using f(n)=n b must be very small in order to remain in equilibrium. However, brushes described by f(n)=2G(N-n) N and f(n)2n were found to be unstable for real, positive values of the potential of the system.

  16. Mean Field Analysis of Stochastic Neural Network Models with Synaptic Depression

    NASA Astrophysics Data System (ADS)

    Yasuhiko Igarashi,; Masafumi Oizumi,; Masato Okada,

    2010-08-01

    We investigated the effects of synaptic depression on the macroscopic behavior of stochastic neural networks. Dynamical mean field equations were derived for such networks by taking the average of two stochastic variables: a firing-state variable and a synaptic variable. In these equations, the average product of thesevariables is decoupled as the product of their averages because the two stochastic variables are independent. We proved the independence of these two stochastic variables assuming that the synaptic weight Jij is of the order of 1/N with respect to the number of neurons N. Using these equations, we derived macroscopic steady-state equations for a network with uniform connections and for a ring attractor network with Mexican hat type connectivity and investigated the stability of the steady-state solutions. An oscillatory uniform state was observed in the network with uniform connections owing to a Hopf instability. For the ring network, high-frequency perturbations were shown not to affect system stability. Two mechanisms destabilize the inhomogeneous steady state, leading to two oscillatory states. A Turing instability leads to a rotating bump state, while a Hopf instability leads to an oscillatory bump state, which was previously unreported. Various oscillatory states take place in a network with synaptic depression depending on the strength of the interneuron connections.

  17. MEAN-FIELD MODELING OF AN α{sup 2} DYNAMO COUPLED WITH DIRECT NUMERICAL SIMULATIONS OF RIGIDLY ROTATING CONVECTION

    SciTech Connect

    Masada, Youhei; Sano, Takayoshi E-mail: sano@ile.osaka-u.ac.jp

    2014-10-10

    The mechanism of large-scale dynamos in rigidly rotating stratified convection is explored by direct numerical simulations (DNS) in Cartesian geometry. A mean-field dynamo model is also constructed using turbulent velocity profiles consistently extracted from the corresponding DNS results. By quantitative comparison between the DNS and our mean-field model, it is demonstrated that the oscillatory α{sup 2} dynamo wave, excited and sustained in the convection zone, is responsible for large-scale magnetic activities such as cyclic polarity reversal and spatiotemporal migration. The results provide strong evidence that a nonuniformity of the α-effect, which is a natural outcome of rotating stratified convection, can be an important prerequisite for large-scale stellar dynamos, even without the Ω-effect.

  18. Brownian regime of finite-N corrections to particle motion in the XY Hamiltonian mean field model

    NASA Astrophysics Data System (ADS)

    Ribeiro, Bruno V.; Amato, Marco A.; Elskens, Yves

    2016-08-01

    We study the dynamics of the N-particle system evolving in the XY Hamiltonian mean field (HMF) model for a repulsive potential, when no phase transition occurs. Starting from a homogeneous distribution, particles evolve in a mean field created by the interaction with all others. This interaction does not change the homogeneous state of the system, and particle motion is approximately ballistic with small corrections. For initial particle data approaching a waterbag, it is explicitly proved that corrections to the ballistic velocities are in the form of independent Brownian noises over a time scale diverging not slower than {N}2/5 as N\\to ∞ , which proves the propagation of molecular chaos. Molecular dynamics simulations of the XY-HMF model confirm our analytical findings.

  19. Relativistic dynamical collapse model

    NASA Astrophysics Data System (ADS)

    Pearle, Philip

    2015-05-01

    A model is discussed where all operators are constructed from a quantum scalar field whose energy spectrum takes on all real values. The Schrödinger picture wave function depends upon space and time coordinates for each particle, as well as an inexorably increasing evolution parameter s which labels a foliation of spacelike hypersurfaces. The model is constructed to be manifestly Lorentz invariant in the interaction picture. Free particle states and interactions are discussed in this framework. Then, the formalism of the continuous spontaneous localization (CSL) theory of dynamical collapse is applied. The collapse-generating operator is chosen to be the particle number space-time density. Unlike previous relativistically invariant models, the vacuum state is not excited. The collapse dynamics depends upon two parameters, a parameter Λ which represents the collapse rate/volume and a scale factor ℓ. A common example of collapse dynamics, involving a clump of matter in a superposition of two locations, is analyzed. The collapse rate is shown to be identical to that of nonrelativistic CSL when the GRW-CSL choice of ℓ=a =1 0-5 cm , is made, along with Λ =λ /a3 (GRW-CSL choice λ =1 0-16s-1). The collapse rate is also satisfactory with the choice ℓ as the size of the Universe, with Λ =λ /ℓa2. Because the collapse narrows wave functions in space and time, it increases a particle's momentum and energy, altering its mass. It is shown that, with ℓ=a , the change of mass of a nucleon is unacceptably large but, when ℓ is the size of the Universe, the change of mass over the age of the Universe is acceptably small.

  20. Chain architecture and micellization: A mean-field coarse-grained model for poly(ethylene oxide) alkyl ether surfactants

    SciTech Connect

    García Daza, Fabián A.; Mackie, Allan D.; Colville, Alexander J.

    2015-03-21

    Microscopic modeling of surfactant systems is expected to be an important tool to describe, understand, and take full advantage of the micellization process for different molecular architectures. Here, we implement a single chain mean field theory to study the relevant equilibrium properties such as the critical micelle concentration (CMC) and aggregation number for three sets of surfactants with different geometries maintaining constant the number of hydrophobic and hydrophilic monomers. The results demonstrate the direct effect of the block organization for the surfactants under study by means of an analysis of the excess energy and entropy which can be accurately determined from the mean-field scheme. Our analysis reveals that the CMC values are sensitive to branching in the hydrophilic head part of the surfactant and can be observed in the entropy-enthalpy balance, while aggregation numbers are also affected by splitting the hydrophobic tail of the surfactant and are manifested by slight changes in the packing entropy.

  1. Microstructure for ferroelastic transitions from strain pseudo-spin clock models in two and three dimensions: a mean field analysis

    SciTech Connect

    Lookman, Turab; Vasseur, Romain

    2009-01-01

    We obtain the microstructure of ferroelastic transitions in two and three dimensions from the solution of their corresponding discrete pseudo-spin models. In two dimensions we consider two transitions each from the high symmetry square and triangle symmetries: square-to-rectangle (SR), square-to-oblique (SO), triangle-to-centered rectangle (TR) and triangle-to-oblique (TO). In three dimensions we study the corresponding spin model for the cubic to tetragonal transition. The Landau free energies for these transitions result in N+ I states clock models (Z{sub N}) with long range interactions and we derive mean-field self-consistency equations for the clock model Hamiltonians. The microstructures from the mean-field solutions of the models are very similar to those obtained from the original continuum models or Monte Carlo simulations on the spin models (in the SR case), illustrating that these discrete models capture the salient physics. The models, in the presence of disorder, provide the basis for the study of the strain glass phase observed in martensitic alloys.

  2. A Mean-Field Photoreaction Model for the Pretilt Generation of a Liquid Crystal on Photopolymer Layers upon Ultraviolet Exposure

    NASA Astrophysics Data System (ADS)

    Na, Jun-Hee; Pae, Hyungwu; Kim, Jinyool; Yu, Chang-Jae; Lee, Sin-Doo

    2011-03-01

    We present a mean-field photoreaction model for the pretilt generation of a nematic liquid crystal (NLC) on the surfaces of photopolymers upon the exposure of ultraviolet (UV) light. The angular distribution function of photopolymer side chains, calculated in the photoreaction model, is used for determining the surface anchoring energy within the Rapini-Papoular approximation. The pretilt angle transition from the homeotropic alignment to the planar alignment of the NLC is demonstrated in two photopolymers with different alkyl chain lengths as a function of the UV exposure time. The main features of the experimental results agree well with theoretical predictions.

  3. About a solvable mean field model of a Gaussian spin glass

    NASA Astrophysics Data System (ADS)

    Barra, Adriano; Genovese, Giuseppe; Guerra, Francesco; Tantari, Daniele

    2014-04-01

    In a series of papers, we have studied a modified Hopfield model of a neural network, with learned words characterized by a Gaussian distribution. The model can be represented as a bipartite spin glass, with one party described by dichotomic Ising spins, and the other party by continuous spin variables, with an a priori Gaussian distribution. By application of standard interpolation methods, we have found it useful to compare the neural network model (bipartite) from one side, with two spin glass models, each monopartite, from the other side. Of these, the first is the usual Sherrington-Kirkpatrick model, the second is a spin glass model, with continuous spins and inbuilt highly nonlinear smooth cut-off interactions. This model is an invaluable laboratory for testing all techniques which have been useful in the study of spin glasses. The purpose of this paper is to give a synthetic description of the most peculiar aspects, by stressing the necessary novelties in the treatment. In particular, it will be shown that the control of the infinite volume limit, according to the well-known Guerra-Toninelli strategy, requires in addition one to consider the involvement of the cut-off interaction in the interpolation procedure. Moreover, the control of the ergodic region, the annealed case, cannot be directly achieved through the standard application of the Borel-Cantelli lemma, but requires previous modification of the interaction. This remark could find useful application in other cases. The replica symmetric expression for the free energy can be easily reached through a suitable version of the doubly stochastic interpolation technique. However, this model shares the unique property that the fully broken replica symmetry ansatz can be explicitly calculated. A very simple sum rule connects the general expression of the fully broken free energy trial function with the replica symmetric one. The definite sign of the error term shows that the replica solution is optimal. Then

  4. MEAN-FIELD SOLAR DYNAMO MODELS WITH A STRONG MERIDIONAL FLOW AT THE BOTTOM OF THE CONVECTION ZONE

    SciTech Connect

    Pipin, V. V.; Kosovichev, A. G.

    2011-09-01

    This paper presents a study of kinematic axisymmetric mean-field dynamo models for the case of meridional circulation with a deep-seated stagnation point and a strong return flow at the bottom of the convection zone. This kind of circulation follows from mean-field models of the angular momentum balance in the solar convection zone. The dynamo models include turbulent sources of the large-scale poloidal magnetic field production due to kinetic helicity and a combined effect due to the Coriolis force and large-scale electric current. In these models the toroidal magnetic field, which is responsible for sunspot production, is concentrated at the bottom of the convection zone and is transported to low-latitude regions by a meridional flow. The meridional component of the poloidal field is also concentrated at the bottom of the convection zone, while the radial component is concentrated in near-polar regions. We show that it is possible for this type of meridional circulation to construct kinematic dynamo models that resemble in some aspects the sunspot magnetic activity cycle. However, in the near-equatorial regions the phase relation between the toroidal and poloidal components disagrees with observations. We also show that the period of the magnetic cycle may not always monotonically decrease with the increase of the meridional flow speed. Thus, for further progress it is important to determine the structure of the meridional circulation, which is one of the critical properties, from helioseismology observations.

  5. Limit theorems in the imitative monomer-dimer mean-field model via Stein's method

    NASA Astrophysics Data System (ADS)

    Chen, Wei-Kuo

    2016-08-01

    We consider the imitative monomer-dimer model on the complete graph introduced in the work of Alberici et al. [J. Math. Phys. 55, 063301-1-063301-27 (2014)]. It was shown that this model is described by the monomer density and has a phase transition along certain coexistence curve, where the monomer and dimer phases coexist. More recently, it was understood [D. Alberici et al., Commun. Math. Phys. (published online, 2016)] that the monomer density exhibits the central limit theorem away from the coexistence curve and enjoys a non-normal limit theorem at criticality with normalized exponent 3/4. By reverting the model to a weighted Curie-Weiss model with hard core interaction, we establish the complete description of the fluctuation properties of the monomer density on the full parameter space via Stein's method of exchangeable pairs. Our approach recovers what were established in the work of Alberici et al. [Commun. Math. Phys. (published online, 2016)] and furthermore allows to obtain the conditional central limit theorems along the coexistence curve. In all these results, the Berry-Esseen inequalities for the Kolmogorov-Smirnov distance are given.

  6. Microstructure from ferroelastic transitions using strain pseudospin clock models in two and three dimensions: A local mean-field analysis

    NASA Astrophysics Data System (ADS)

    Vasseur, Romain; Lookman, Turab; Shenoy, Subodh R.

    2010-09-01

    We show how microstructure can arise in first-order ferroelastic structural transitions, in two and three spatial dimensions, through a local mean-field approximation of their pseudospin Hamiltonians, that include anisotropic elastic interactions. Such transitions have symmetry-selected physical strains as their NOP -component order parameters, with Landau free energies that have a single zero-strain “austenite” minimum at high temperatures, and spontaneous-strain “martensite” minima of NV structural variants at low temperatures. The total free energy also has gradient terms, and power-law anisotropic effective interactions, induced by “no-dislocation” St Venant compatibility constraints. In a reduced description, the strains at Landau minima induce temperature dependent, clocklike ZNV+1 Hamiltonians, with NOP -component strain-pseudospin vectors S⃗ pointing to NV+1 discrete values (including zero). We study elastic texturing in five such first-order structural transitions through a local mean-field approximation of their pseudospin Hamiltonians, that include the power-law interactions. As a prototype, we consider the two-variant square/rectangle transition, with a one-component pseudospin taking NV+1=3 values of S=0,±1 , as in a generalized Blume-Capel model. We then consider transitions with two-component (NOP=2) pseudospins: the equilateral to centered rectangle (NV=3) ; the square to oblique polygon (NV=4) ; the triangle to oblique (NV=6) transitions; and finally the three-dimensional (3D) cubic to tetragonal transition (NV=3) . The local mean-field solutions in two-dimensional and 3D yield oriented domain-wall patterns as from continuous-variable strain dynamics, showing the discrete-variable models capture the essential ferroelastic texturings. Other related Hamiltonians illustrate that structural transitions in materials science can be the source of interesting spin models in statistical mechanics.

  7. Dynamics and termination cost of spatially coupled mean-field models

    NASA Astrophysics Data System (ADS)

    Caltagirone, Francesco; Franz, Silvio; Morris, Richard G.; Zdeborová, Lenka

    2014-01-01

    This work is motivated by recent progress in information theory and signal processing where the so-called spatially coupled design of systems leads to considerably better performance. We address relevant open questions about spatially coupled systems through the study of a simple Ising model. In particular, we consider a chain of Curie-Weiss models that are coupled by interactions up to a certain range. Indeed, it is well known that the pure (uncoupled) Curie-Weiss model undergoes a first-order phase transition driven by the magnetic field, and furthermore in the spinodal region such systems are unable to reach equilibrium in subexponential time if initialized in the metastable state. In contrast, the spatially coupled system is instead able to reach the equilibrium even when initialized to the metastable state. The equilibrium phase propagates along the chain in the form of a traveling wave. Here we study the speed of the wave front and the so-called termination cost—i.e., the conditions necessary for the propagation to occur. We reach several interesting conclusions about optimization of the speed and the cost.

  8. Periodic solutions to a mean-field model for electrocortical activity

    NASA Astrophysics Data System (ADS)

    van Veen, L.; Green, K. R.

    2014-12-01

    We consider a continuum model of electrical signals in the human cortex, which takes the form of a system of semilinear, hyperbolic partial differential equations for the inhibitory and excitatory membrane potentials and the synaptic inputs. The coupling of these components is represented by sigmoidal and quadratic nonlinearities. We consider these equations on a square domain with periodic boundary conditions, in the vicinity of the primary transition from a stable equilibrium to time-periodic motion through an equivariant Hopf bifurcation. We compute part of a family of standing wave solutions, emanating from this point.

  9. Analytical approaches to modelling panspermia - beyond the mean-field paradigm

    NASA Astrophysics Data System (ADS)

    Lingam, Manasvi

    2016-01-01

    We model the process of panspermia by adopting two different approaches. The first method conceives it as a self-replication process, endowed with non-local creation and extinction. We show that some features suggestive of universal behaviour emerge, such as exponential decay or growth, and a power spectral density that displays a power-law behaviour in a particular regime. We also present a special case wherein the number density of the planets seeded through panspermia approaches a finite asymptotic distribution. The power spectral density for the independent and spontaneous emergence of life is investigated in conjunction with its counterpart for panspermia. The former exhibits attributes characteristic of a noise spectrum, including the resemblance to white noise in a certain regime. These features are absent in panspermia, suggesting that the power spectral density could be utilized as a future tool for differentiating between the two processes. Our second approach adopts the machinery of Markov processes and diffusion, and we show that the power spectral density exhibits a power-law tail in some domains, as earlier, suggesting that this behaviour may be fairly robust. We comment on a generalization of the diffusive model, and also indicate how the methods and results developed herein could be used to analyse other phenomena.

  10. Folding of small knotted proteins: Insights from a mean field coarse-grained model

    SciTech Connect

    Najafi, Saeed; Potestio, Raffaello

    2015-12-28

    A small but relevant number of proteins whose native structure is known features nontrivial topology, i.e., they are knotted. Understanding the process of folding from a swollen unknotted state to the biologically relevant native conformation is, for these proteins, particularly difficult, due to their rate-limiting topological entanglement. To shed some light into this conundrum, we introduced a structure-based coarse-grained model of the protein, where the information about the folded conformation is encoded in bonded angular interactions only, which do not favor the formation of native contacts. A stochastic search scheme in parameter space is employed to identify a set of interactions that maximizes the probability to attain the knotted state. The optimal knotting pathways of the two smallest knotted proteins, obtained through this approach, are consistent with the results derived by means of coarse-grained as well as full atomistic simulations.

  11. Mott transitions in a three-component Falicov-Kimball model: A slave boson mean-field study

    NASA Astrophysics Data System (ADS)

    Le, Duc-Anh; Tran, Minh-Tien

    2015-05-01

    Metal-insulator transitions in a three-component Falicov-Kimball model are investigated within the Kotliar-Ruckenstein slave boson mean-field approach. The model describes a mixture of two interacting fermion atom species loaded into an optical lattice at ultralow temperature. One species is two-component atoms, which can hop in the optical lattice, and the other is single-component atoms, which are localized. Different correlation-driven metal-insulator transitions are observed depending on the atom filling conditions and local interactions. These metal-insulator transitions are classified by the band renormalization factors and the double occupancies of the atom species. The filling conditions and the critical value of the local interactions for these metal-insulator transitions are also analytically established. The obtained results not only are in good agreement with the dynamical mean-field theory for the three-component Falicov-Kimball model but also clarify the nature and properties of the metal-insulator transitions in a simple physics picture.

  12. Mean-field dispersion-induced spatial synchrony, oscillation and amplitude death, and temporal stability in an ecological model.

    PubMed

    Banerjee, Tanmoy; Dutta, Partha Sharathi; Gupta, Anubhav

    2015-05-01

    One of the most important issues in spatial ecology is to understand how spatial synchrony and dispersal-induced stability interact. In the existing studies it is shown that dispersion among identical patches results in spatial synchrony; on the other hand, the combination of spatial heterogeneity and dispersion is necessary for dispersal-induced stability (or temporal stability). Population synchrony and temporal stability are thus often thought of as conflicting outcomes of dispersion. In contrast to the general belief, in this present study we show that mean-field dispersion is conducive to both spatial synchrony and dispersal-induced stability even in identical patches. This simultaneous occurrence of rather conflicting phenomena is governed by the suppression of oscillation states, namely amplitude death (AD) and oscillation death (OD). These states emerge through spatial synchrony of the oscillating patches in the strong-coupling strength. We present an interpretation of the mean-field diffusive coupling in the context of ecology and identify that, with increasing mean-field density, an open ecosystem transforms into a closed ecosystem. We report on the occurrence of OD in an ecological model and explain its significance. Using a detailed bifurcation analysis we show that, depending on the mortality rate and carrying capacity, the system shows either AD or both AD and OD. We also show that the results remain qualitatively the same for a network of oscillators. We identify a new transition scenario between the same type of oscillation suppression states whose geneses differ. In the parameter-mismatched case, we further report on the direct transition from OD to AD through a transcritical bifurcation. We believe that this study will lead to a proper interpretation of AD and OD in ecology, which may be important for the conservation and management of several communities in ecosystems. PMID:26066241

  13. Mean-field dispersion-induced spatial synchrony, oscillation and amplitude death, and temporal stability in an ecological model

    NASA Astrophysics Data System (ADS)

    Banerjee, Tanmoy; Dutta, Partha Sharathi; Gupta, Anubhav

    2015-05-01

    One of the most important issues in spatial ecology is to understand how spatial synchrony and dispersal-induced stability interact. In the existing studies it is shown that dispersion among identical patches results in spatial synchrony; on the other hand, the combination of spatial heterogeneity and dispersion is necessary for dispersal-induced stability (or temporal stability). Population synchrony and temporal stability are thus often thought of as conflicting outcomes of dispersion. In contrast to the general belief, in this present study we show that mean-field dispersion is conducive to both spatial synchrony and dispersal-induced stability even in identical patches. This simultaneous occurrence of rather conflicting phenomena is governed by the suppression of oscillation states, namely amplitude death (AD) and oscillation death (OD). These states emerge through spatial synchrony of the oscillating patches in the strong-coupling strength. We present an interpretation of the mean-field diffusive coupling in the context of ecology and identify that, with increasing mean-field density, an open ecosystem transforms into a closed ecosystem. We report on the occurrence of OD in an ecological model and explain its significance. Using a detailed bifurcation analysis we show that, depending on the mortality rate and carrying capacity, the system shows either AD or both AD and OD. We also show that the results remain qualitatively the same for a network of oscillators. We identify a new transition scenario between the same type of oscillation suppression states whose geneses differ. In the parameter-mismatched case, we further report on the direct transition from OD to AD through a transcritical bifurcation. We believe that this study will lead to a proper interpretation of AD and OD in ecology, which may be important for the conservation and management of several communities in ecosystems.

  14. Relativistic Continuum Shell Model

    NASA Astrophysics Data System (ADS)

    Grineviciute, Janina; Halderson, Dean

    2011-04-01

    The R-matrix formalism of Lane and Thomas has been extended to the relativistic case so that the many-coupled channels problem may be solved for systems in which binary breakup channels satisfy a relative Dirac equation. The formalism was previously applied to the relativistic impulse approximation RIA and now we applied it to Quantum Hadrodynamics QHD in the continuum Tamm-Dancoff approximation TDA with the classical meson fields replaced by one-meson exchange potentials. None of the published QHD parameters provide a decent fit to the 15 N + p elastic cross section. The deficiency is also evident in inability of the QHD parameters with the one meson exchange potentials to reproduce the QHD single particle energies. Results with alternate parameters sets are presented. A. M. Lane and R. G. Thomas, R-Matrix Theory of Nuclear Reactions, Reviews of Modern Physics, 30 (1958) 257

  15. Approximate analytical solution for nuclear matter in a mean-field Walecka model and Coester line behavior

    SciTech Connect

    Delfino, A.; Silva, J.B.; Malheiro, M.

    2006-03-15

    We study nuclear matter, at the mean-field approximation, by considering as equal the values of the scalar and the vector density in the Walecka model, which is a very reasonable approximation up to the nuclear matter saturation density. It turns out that the model has an analytical solution for the scalar and vector couplings as functions only of the nuclear matter density and binding energy. The nuclear matter properties are very close to the original version of the model. This solution allows us to show that the correlation between the binding energy and the saturation density is Coester line like. The liquid-gas phase transition is also studied and the critical and flash temperatures are again very similar to the original ones.

  16. Modeling MHD accretion-ejection: episodic ejections of jets triggered by a mean-field disk dynamo

    SciTech Connect

    Stepanovs, Deniss; Fendt, Christian; Sheikhnezami, Somayeh E-mail: fendt@mpia.de

    2014-11-20

    We present MHD simulations exploring the launching, acceleration, and collimation of jets and disk winds. The evolution of the disk structure is consistently taken into account. Extending our earlier studies, we now consider the self-generation of the magnetic field by an α{sup 2}Ω mean-field dynamo. The disk magnetization remains on a rather low level, which helps to evolve the simulations for T > 10, 000 dynamical time steps on a domain extending 1500 inner disk radii. We find the magnetic field of the inner disk to be similar to the commonly found open field structure, favoring magneto-centrifugal launching. The outer disk field is highly inclined and predominantly radial. Here, differential rotation induces a strong toroidal component, which plays a key role in outflow launching. These outflows from the outer disk are slower, denser, and less collimated. If the dynamo action is not quenched, magnetic flux is continuously generated, diffuses outward through the disk, and fills the entire disk. We have invented a toy model triggering a time-dependent mean-field dynamo. The duty cycles of this dynamo lead to episodic ejections on similar timescales. When the dynamo is suppressed as the magnetization falls below a critical value, the generation of the outflows and also accretion is inhibited. The general result is that we can steer episodic ejection and large-scale jet knots by a disk-intrinsic dynamo that is time-dependent and regenerates the jet-launching magnetic field.

  17. Parametric Smoothness and Self- Scaling of the Statistical Properties of a Minimal Climate Model: what beyond mean field theories?

    NASA Astrophysics Data System (ADS)

    Lucarini, V.; Speranza, A.; Vitolo, R.

    2009-04-01

    A quasi-geostrophic intermediate complexity model of the mid-latitude atmospheric circulation is considered, featuring simplified baroclinic conversion and barotropic convergence processes. The model undergoes baroclinic forcing towards a given latitudinal temperature profile controlled by the forced equator-to-pole temperature difference Te. When Te increases, a transition takes place from a stationary regime-Hadley equilibrium-to a periodic regime, and eventually to a chaotic regime where evolution takes place on a strange attractor. The attractor dimension, metric entropy, and bounding box volume in phase space have a smooth dependence on Te which results in power-law scaling properties. Power-law scalings are detected also for the statistical properties of global physical observables — the total energy of the system and the averaged zonal wind. The scaling laws, which constitute the main novel result of the present work, can be thought to result from the presence of a statistical process of baroclinic adjustment, which tends to decrease the equator-to-pole temperature difference and determines the properties of the attractor of the system. The self-similarity could be of great help in setting up a theory for the overall statistical properties of the general circulation of the atmosphere and in guiding-on a heuristic basis-both data analysis and realistic simulations, going beyond the unsatisfactory mean field theories and /brute force/ approaches. A leading example for this would be the possibility of estimating the sensitivity of the output of the system with respect to changes in the parameters. Ref: Valerio Lucarini, Antonio Speranza, Renato Vitolo, Parametric smoothness and self-scaling of the statistical properties of a minimal climate model: What beyond the mean field theories?, Physica D, 234 (2007), 105-123

  18. Incommensurate phase of a triangular frustrated Heisenberg model studied via Schwinger-boson mean-field theory

    NASA Astrophysics Data System (ADS)

    Li, Peng; Su, Haibin; Dong, Hui-Ning; Shen, Shun-Qing

    2009-08-01

    We study a triangular frustrated antiferromagnetic Heisenberg model with nearest-neighbor interactions J1 and third-nearest-neighbor interactions J3 by means of Schwinger-boson mean-field theory. By setting an antiferromagnetic J3 and varying J1 from positive to negative values, we disclose the low-temperature features of its interesting incommensurate phase. The gapless dispersion of quasiparticles leads to the intrinsic T2 law of specific heat. The magnetic susceptibility is linear in temperature. The local magnetization is significantly reduced by quantum fluctuations. We address possible relevance of these results to the low-temperature properties of NiGa2S4. From a careful analysis of the incommensurate spin wavevector, the interaction parameters are estimated as J1≈-3.8755 K and J3≈14.0628 K, in order to account for the experimental data.

  19. A Note on the Guerra and Talagrand Theorems for Mean Field Spin Glasses: The Simple Case of Spherical Models

    NASA Astrophysics Data System (ADS)

    Franz, Silvio; Tria, Francesca

    2006-01-01

    The aim of this paper is to discuss the main ideas of the Talagrand proof of the Parisi Ansatz for the free-energy of Mean Field Spin Glasses with a physicist's approach. We consider the case of the spherical p-spin model, which has the following advantages: (1) the Parisi Ansatz takes the simple "one step replica symmetry breaking form," (2) the replica free-energy as a function of the order parameters is simple enough to allow for numerical maximization with arbitrary precision. We present the essential ideas of the proof, we stress its connections with the theory of effective potentials for glassy systems, and we reduce the technically more difficult part of the Talagrand's analysis to an explicit evaluation of the solution of a variational problem.

  20. Continuum random-phase approximation for relativistic point coupling models

    SciTech Connect

    Daoutidis, J.; Ring, P.

    2009-08-15

    Continuum relativistic random-phase approximation (CRPA) is used to investigate collective excitation phenomena in several spherical nuclei along the periodic table. We start from relativistic mean-field calculations based on a covariant density functional with density-dependent zero-range forces. From the same functional an effective interaction is obtained as the second derivative with respect to the density. This interaction is used in relativistic CRPA calculations for the investigation of isoscalar monopole, isovector dipole, and isoscalar quadrupole resonances of spherical nuclei. In particular we study the low-lying E1 strength in the vicinity of the neutron evaporation threshold. The properties of the resonances, such as centroid energies and strengths distributions are compared with results of discrete RPA calculations for the same model as well as with experimental data.

  1. Flocking dynamics and mean-field limit in the Cucker-Smale-type model with topological interactions

    NASA Astrophysics Data System (ADS)

    Haskovec, Jan

    2013-10-01

    We introduce a Cucker-Smale-type model for flocking, where the strength of interaction between two agents depends on their relative separation (called “topological distance” in previous works), which is the number of intermediate individuals separating them. This makes the model scale-free and is motivated by recent extensive observations of starling flocks, suggesting that the interaction ruling animal collective behavior depends on topological rather than the metric distance. We study the conditions leading to asymptotic flocking in the topological model, defined as the convergence of the agents’ velocities to a common vector. The shift from metric to topological interactions requires development of new analytical methods, taking into account the graph-theoretical nature of the problem. Moreover, we provide a rigorous derivation of the mean-field limit of large populations, recovering kinetic and hydrodynamic descriptions. In particular, we introduce the novel concept of relative separation in continuum descriptions, which is applicable to a broad variety of models of collective behavior. As an example, we shortly discuss a topological modification of the attraction-repulsion model and illustrate with numerical simulations that the modified model produces interesting new pattern dynamics.

  2. Adapting Poisson-Boltzmann to the self-consistent mean field theory: Application to protein side-chain modeling

    NASA Astrophysics Data System (ADS)

    Koehl, Patrice; Orland, Henri; Delarue, Marc

    2011-08-01

    We present an extension of the self-consistent mean field theory for protein side-chain modeling in which solvation effects are included based on the Poisson-Boltzmann (PB) theory. In this approach, the protein is represented with multiple copies of its side chains. Each copy is assigned a weight that is refined iteratively based on the mean field energy generated by the rest of the protein, until self-consistency is reached. At each cycle, the variational free energy of the multi-copy system is computed; this free energy includes the internal energy of the protein that accounts for vdW and electrostatics interactions and a solvation free energy term that is computed using the PB equation. The method converges in only a few cycles and takes only minutes of central processing unit time on a commodity personal computer. The predicted conformation of each residue is then set to be its copy with the highest weight after convergence. We have tested this method on a database of hundred highly refined NMR structures to circumvent the problems of crystal packing inherent to x-ray structures. The use of the PB-derived solvation free energy significantly improves prediction accuracy for surface side chains. For example, the prediction accuracies for χ1 for surface cysteine, serine, and threonine residues improve from 68%, 35%, and 43% to 80%, 53%, and 57%, respectively. A comparison with other side-chain prediction algorithms demonstrates that our approach is consistently better in predicting the conformations of exposed side chains.

  3. Exact chiral spin liquids and mean-field perturbations of gamma matrix models on the ruby lattice

    NASA Astrophysics Data System (ADS)

    Whitsitt, Seth; Chua, Victor; Fiete, Gregory A.

    2012-11-01

    We theoretically studied an exactly solvable gamma matrix generalization of the Kitaev spin model on the ruby lattice, which is a honeycomb lattice with ‘expanded’ vertices and links. We find that this model displays an exceptionally rich phase diagram that includes (i) gapless phases with stable spin Fermi surfaces, (ii) gapless phases with low-energy Dirac cones and quadratic band touching points and (iii) gapped phases with finite Chern numbers possessing the values ±4,±3,±2 and ±1. The model is then generalized to include Ising-like interactions that break the exact solvability of the model in a controlled manner. When these terms are dominant, they lead to a trivial Ising ordered phase which is shown to be adiabatically connected to a large coupling limit of the exactly solvable phase. In the limit where these interactions are weak, we treat them within mean-field theory and present the resulting phase diagrams. We discuss the nature of the transitions between various phases. Our results show the richness of possible ground states in closely related magnetic systems.

  4. Large pseudocounts and L2-norm penalties are necessary for the mean-field inference of Ising and Potts models

    NASA Astrophysics Data System (ADS)

    Barton, J. P.; Cocco, S.; De Leonardis, E.; Monasson, R.

    2014-07-01

    The mean-field (MF) approximation offers a simple, fast way to infer direct interactions between elements in a network of correlated variables, a common, computationally challenging problem with practical applications in fields ranging from physics and biology to the social sciences. However, MF methods achieve their best performance with strong regularization, well beyond Bayesian expectations, an empirical fact that is poorly understood. In this work, we study the influence of pseudocount and L2-norm regularization schemes on the quality of inferred Ising or Potts interaction networks from correlation data within the MF approximation. We argue, based on the analysis of small systems, that the optimal value of the regularization strength remains finite even if the sampling noise tends to zero, in order to correct for systematic biases introduced by the MF approximation. Our claim is corroborated by extensive numerical studies of diverse model systems and by the analytical study of the m-component spin model for large but finite m. Additionally, we find that pseudocount regularization is robust against sampling noise and often outperforms L2-norm regularization, particularly when the underlying network of interactions is strongly heterogeneous. Much better performances are generally obtained for the Ising model than for the Potts model, for which only couplings incoming onto medium-frequency symbols are reliably inferred.

  5. Screening and nonlocal correlations in the extended Hubbard model from self-consistent combined GW and dynamical mean field theory

    NASA Astrophysics Data System (ADS)

    Ayral, Thomas; Biermann, Silke; Werner, Philipp

    2013-03-01

    We describe a recent implementation of the combined GW and dynamical mean field method (GW+DMFT) for the two-dimensional Hubbard model with onsite and nearest-neighbor repulsion. We clarify the relation of the GW+DMFT scheme to alternative approaches in the literature, and discuss the corresponding approximations to the free-energy functional of the model. Furthermore, we describe a numerically exact technique for the solution of the GW+DMFT equations, namely, the hybridization expansion continuous-time algorithm for impurity models with retarded interactions. We compute the low-temperature phase diagram of the half-filled extended Hubbard model, addressing the metal-insulator transition at small intersite interactions and the transition to a charge-ordered state for stronger intersite repulsions. GW+DMFT introduces a nontrivial momentum dependence into the many-body self-energy and polarization. We find that the charge fluctuations included in the present approach have a larger impact on the latter than on the former. Finally, within the GW+DMFT framework, as in extended DMFT, the intersite repulsion translates into a frequency dependence of the local effective interaction. We analyze this dependence and show how it affects the local spectral function.

  6. Mean-field model of the von Kármán sodium dynamo experiment using soft iron impellers

    NASA Astrophysics Data System (ADS)

    Nore, C.; Léorat, J.; Guermond, J.-L.; Giesecke, A.

    2015-01-01

    It has been observed that dynamo action occurs in the von-Kármán-Sodium (VKS) experiment only when the rotating disks and the blades are made of soft iron. The purpose of this paper is to numerically investigate the role of soft iron in the VKS dynamo scenario. This is done by using a mean-field model based on an axisymmetric mean flow, a localized permeability distribution, and a localized α effect modeling the action of the small velocity scales between the blades. The action of the rotating blades is modeled by an axisymmetric effective permeability field. Key properties of the flow giving to the numerical magnetic field a geometric structure similar to that observed experimentally are identified. Depending on the permeability of the disks and the effective permeability of the blades, the dynamo that is obtained is either oscillatory or stationary. Our numerical results confirm the leading role played by the ferromagnetic impellers. A scenario for the VKS dynamo is proposed.

  7. Mean-field model of the von Kármán sodium dynamo experiment using soft iron impellers.

    PubMed

    Nore, C; Léorat, J; Guermond, J-L; Giesecke, A

    2015-01-01

    It has been observed that dynamo action occurs in the von-Kármán-Sodium (VKS) experiment only when the rotating disks and the blades are made of soft iron. The purpose of this paper is to numerically investigate the role of soft iron in the VKS dynamo scenario. This is done by using a mean-field model based on an axisymmetric mean flow, a localized permeability distribution, and a localized α effect modeling the action of the small velocity scales between the blades. The action of the rotating blades is modeled by an axisymmetric effective permeability field. Key properties of the flow giving to the numerical magnetic field a geometric structure similar to that observed experimentally are identified. Depending on the permeability of the disks and the effective permeability of the blades, the dynamo that is obtained is either oscillatory or stationary. Our numerical results confirm the leading role played by the ferromagnetic impellers. A scenario for the VKS dynamo is proposed. PMID:25679709

  8. Primitive-path statistics of entangled polymers: mapping multi-chain simulations onto single-chain mean-field models

    NASA Astrophysics Data System (ADS)

    Steenbakkers, Rudi J. A.; Tzoumanekas, Christos; Li, Ying; Liu, Wing Kam; Kröger, Martin; Schieber, Jay D.

    2014-01-01

    We present a method to map the full equilibrium distribution of the primitive-path (PP) length, obtained from multi-chain simulations of polymer melts, onto a single-chain mean-field ‘target’ model. Most previous works used the Doi-Edwards tube model as a target. However, the average number of monomers per PP segment, obtained from multi-chain PP networks, has consistently shown a discrepancy of a factor of two with respect to tube-model estimates. Part of the problem is that the tube model neglects fluctuations in the lengths of PP segments, the number of entanglements per chain and the distribution of monomers among PP segments, while all these fluctuations are observed in multi-chain simulations. Here we use a recently proposed slip-link model, which includes fluctuations in all these variables as well as in the spatial positions of the entanglements. This turns out to be essential to obtain qualitative and quantitative agreement with the equilibrium PP-length distribution obtained from multi-chain simulations. By fitting this distribution, we are able to determine two of the three parameters of the model, which govern its equilibrium properties. This mapping is executed for four different linear polymers and for different molecular weights. The two parameters are found to depend on chemistry, but not on molecular weight. The model predicts a constant plateau modulus minus a correction inversely proportional to molecular weight. The value for well-entangled chains, with the parameters determined ab initio, lies in the range of experimental data for the materials investigated.

  9. Dispersal-induced synchrony, temporal stability, and clustering in a mean-field coupled Rosenzweig-MacArthur model

    NASA Astrophysics Data System (ADS)

    Arumugam, Ramesh; Dutta, Partha Sharathi; Banerjee, Tanmoy

    2015-10-01

    In spatial ecology, dispersal among a set of spatially separated habitats, named as metapopulation, preserves the diversity and persistence by interconnecting the local populations. Understanding the effects of several variants of dispersion in metapopulation dynamics and to identify the factors which promote population synchrony and population stability are important in ecology. In this paper, we consider the mean-field dispersion among the habitats in a network and study the collective dynamics of the spatially extended system. Using the Rosenzweig-MacArthur model for individual patches, we show that the population synchrony and temporal stability, which are believed to be of conflicting outcomes of dispersion, can be simultaneously achieved by oscillation quenching mechanisms. Particularly, we explore the more natural coupling configuration where the rates of dispersal of different habitats are disparate. We show that asymmetry in dispersal rate plays a crucial role in determining inhomogeneity in an otherwise homogeneous metapopulation. We further identify an unusual emergent state in the network, namely, a multi-branch clustered inhomogeneous steady state, which arises due to the intrinsic parameter mismatch among the patches. We believe that the present study will shed light on the cooperative behavior of spatially structured ecosystems.

  10. Dispersal-induced synchrony, temporal stability, and clustering in a mean-field coupled Rosenzweig-MacArthur model.

    PubMed

    Arumugam, Ramesh; Dutta, Partha Sharathi; Banerjee, Tanmoy

    2015-10-01

    In spatial ecology, dispersal among a set of spatially separated habitats, named as metapopulation, preserves the diversity and persistence by interconnecting the local populations. Understanding the effects of several variants of dispersion in metapopulation dynamics and to identify the factors which promote population synchrony and population stability are important in ecology. In this paper, we consider the mean-field dispersion among the habitats in a network and study the collective dynamics of the spatially extended system. Using the Rosenzweig-MacArthur model for individual patches, we show that the population synchrony and temporal stability, which are believed to be of conflicting outcomes of dispersion, can be simultaneously achieved by oscillation quenching mechanisms. Particularly, we explore the more natural coupling configuration where the rates of dispersal of different habitats are disparate. We show that asymmetry in dispersal rate plays a crucial role in determining inhomogeneity in an otherwise homogeneous metapopulation. We further identify an unusual emergent state in the network, namely, a multi-branch clustered inhomogeneous steady state, which arises due to the intrinsic parameter mismatch among the patches. We believe that the present study will shed light on the cooperative behavior of spatially structured ecosystems. PMID:26520087

  11. Majority-vote model on spatially embedded networks: Crossover from mean-field to Ising universality classes

    NASA Astrophysics Data System (ADS)

    Sampaio Filho, C. I. N.; dos Santos, T. B.; Moreira, A. A.; Moreira, F. G. B.; Andrade, J. S.

    2016-05-01

    We study through Monte Carlo simulations and finite-size scaling analysis the nonequilibrium phase transitions of the majority-vote model taking place on spatially embedded networks. These structures are built from an underlying regular lattice over which directed long-range connections are randomly added according to the probability Pi j˜r-α , where ri j is the Manhattan distance between nodes i and j , and the exponent α is a controlling parameter [J. M. Kleinberg, Nature (London) 406, 845 (2000), 10.1038/35022643]. Our results show that the collective behavior of this system exhibits a continuous order-disorder phase transition at a critical parameter, which is a decreasing function of the exponent α . Precisely, considering the scaling functions and the critical exponents calculated, we conclude that the system undergoes a crossover among distinct universality classes. For α ≤3 the critical behavior is described by mean-field exponents, while for α ≥4 it belongs to the Ising universality class. Finally, in the region where the crossover occurs, 3 <α <4 , the critical exponents are dependent on α .

  12. Improving solar 11yr magnetic cycle prediction by using variational data assimilation in a mean field dynamo model

    NASA Astrophysics Data System (ADS)

    Hung, Ching Pui; Jouve, Laurène; Brun, Allan-Sacha; Fournier, Alexandre; Talagrand, Olivier

    2015-04-01

    We present our recent effort to implement modern variational data assimilation techniques into a 2.5 D mean field solar dynamo code. This work extend the work of (Jouve et al. 2011, ApJ) to take into account the correct spherical geometry and meridional circulation into so-called Babccok-Leigthon flux transport dynamo models. Based on twin-experiments, in which we observe our dynamo simulations, and on a well defined cost function using toroidal and poloidal field observations we are able to recover the main attributes of the dynamo solution used to test our data assimilation algorithm. By assimilating solar data (such as Wolf number or butterfly diagram) we are starting to deduce the profile and temporal variations of key ingredients of the solar dynamo. We find that the data sampling and the temporal window are key to get reliable results. We show how such powerful technique can be used to improve our ability to predict the solar magnetic activity. This work is supported by Idex Sorbonne Paris Cite via the DAMSE project.

  13. Stochastic multi-scale models of competition within heterogeneous cellular populations: Simulation methods and mean-field analysis.

    PubMed

    Cruz, Roberto de la; Guerrero, Pilar; Spill, Fabian; Alarcón, Tomás

    2016-10-21

    We propose a modelling framework to analyse the stochastic behaviour of heterogeneous, multi-scale cellular populations. We illustrate our methodology with a particular example in which we study a population with an oxygen-regulated proliferation rate. Our formulation is based on an age-dependent stochastic process. Cells within the population are characterised by their age (i.e. time elapsed since they were born). The age-dependent (oxygen-regulated) birth rate is given by a stochastic model of oxygen-dependent cell cycle progression. Once the birth rate is determined, we formulate an age-dependent birth-and-death process, which dictates the time evolution of the cell population. The population is under a feedback loop which controls its steady state size (carrying capacity): cells consume oxygen which in turn fuels cell proliferation. We show that our stochastic model of cell cycle progression allows for heterogeneity within the cell population induced by stochastic effects. Such heterogeneous behaviour is reflected in variations in the proliferation rate. Within this set-up, we have established three main results. First, we have shown that the age to the G1/S transition, which essentially determines the birth rate, exhibits a remarkably simple scaling behaviour. Besides the fact that this simple behaviour emerges from a rather complex model, this allows for a huge simplification of our numerical methodology. A further result is the observation that heterogeneous populations undergo an internal process of quasi-neutral competition. Finally, we investigated the effects of cell-cycle-phase dependent therapies (such as radiation therapy) on heterogeneous populations. In particular, we have studied the case in which the population contains a quiescent sub-population. Our mean-field analysis and numerical simulations confirm that, if the survival fraction of the therapy is too high, rescue of the quiescent population occurs. This gives rise to emergence of resistance

  14. Hierarchical mean-field approach to the J1-J2 Heisenberg model on a square lattice

    NASA Astrophysics Data System (ADS)

    Isaev, Leonid; Ortiz, Gerardo; Dukelsky, Jorge

    2009-03-01

    We study the quantum phase diagram and excitation spectrum of the frustrated J1-J2 spin-1/2 Heisenberg Hamiltonian. A hierarchical mean-field approach, at the heart of which lies the idea of identifying relevant degrees of freedom, is developed. Thus, by performing educated, manifestly symmetry preserving mean-field approximations, we unveil fundamental properties of the system. We then compare various coverings of the square lattice with plaquettes, dimers and other degrees of freedom, and show that only the symmetric plaquette covering, which reproduces the original Bravais lattice, leads to the known phase diagram. The intermediate quantum paramagnetic phase is shown to be a (singlet) plaquette crystal, connected with the neighbouring N'eel phase by a continuous phase transition. We also introduce fluctuations around the hierarchical mean-field solutions, and demonstrate that in the paramagnetic phase the ground and first excited states are separated by a finite gap, which closes in the N'eel and columnar phases. Our results suggest that the quantum phase transition between N'eel and paramagnetic phases can be properly described within the Ginzburg-Landau-Wilson paradigm.

  15. Hierarchical mean-field approach to the J1-J2 Heisenberg model on a square lattice

    NASA Astrophysics Data System (ADS)

    Isaev, L.; Ortiz, G.; Dukelsky, J.

    2009-01-01

    We study the quantum phase diagram and excitation spectrum of the frustrated J1-J2 spin-1/2 Heisenberg Hamiltonian. A hierarchical mean-field approach, at the heart of which lies the idea of identifying relevant degrees of freedom, is developed. Thus, by performing educated, manifestly symmetry-preserving mean-field approximations, we unveil fundamental properties of the system. We then compare various coverings of the square lattice with plaquettes, dimers, and other degrees of freedom, and show that only the symmetric plaquette covering, which reproduces the original Bravais lattice, leads to the known phase diagram. The intermediate quantum paramagnetic phase is shown to be a (singlet) plaquette crystal, connected with the neighboring Néel phase by a continuous phase transition. We also introduce fluctuations around the hierarchical mean-field solutions, and demonstrate that in the paramagnetic phase the ground and first excited states are separated by a finite gap, which closes in the Néel and columnar phases. Our results suggest that the quantum phase transition between Néel and paramagnetic phases can be properly described within the Ginzburg-Landau-Wilson paradigm.

  16. Continuous Time Finite State Mean Field Games

    SciTech Connect

    Gomes, Diogo A.; Mohr, Joana Souza, Rafael Rigao

    2013-08-01

    In this paper we consider symmetric games where a large number of players can be in any one of d states. We derive a limiting mean field model and characterize its main properties. This mean field limit is a system of coupled ordinary differential equations with initial-terminal data. For this mean field problem we prove a trend to equilibrium theorem, that is convergence, in an appropriate limit, to stationary solutions. Then we study an N+1-player problem, which the mean field model attempts to approximate. Our main result is the convergence as N{yields}{infinity} of the mean field model and an estimate of the rate of convergence. We end the paper with some further examples for potential mean field games.

  17. Statistics of Dislocation Slip Avalanches in Nanosized Single Crystals Show Tuned Critical Behavior Predicted by a Simple Mean Field Model

    NASA Astrophysics Data System (ADS)

    Friedman, Nir; Jennings, Andrew T.; Tsekenis, Georgios; Kim, Ju-Young; Tao, Molei; Uhl, Jonathan T.; Greer, Julia R.; Dahmen, Karin A.

    2012-08-01

    We show that slowly sheared metallic nanocrystals deform via discrete strain bursts (slips), whose size distributions follow power laws with stress-dependent cutoffs. We show for the first time that plasticity reflects tuned criticality, by collapsing the stress-dependent slip-size distributions onto a predicted scaling function. Both power-law exponents and scaling function agree with mean-field theory predictions. Our study of 7 materials and 2 crystal structures, at various deformation rates, stresses, and crystal sizes down to 75 nm, attests to the universal characteristics of plasticity.

  18. Periodic mean-field solutions and the spectra of discrete bosonic fields: Trace formula for Bose-Hubbard models.

    PubMed

    Engl, Thomas; Urbina, Juan Diego; Richter, Klaus

    2015-12-01

    We consider the many-body spectra of interacting bosonic quantum fields on a lattice in the semiclassical limit of large particle number N. We show that the many-body density of states can be expressed as a coherent sum over oscillating long-wavelength contributions given by periodic, nonperturbative solutions of the, typically nonlinear, wave equation of the classical (mean-field) limit. To this end, we construct the semiclassical approximation for both the smooth and oscillatory parts of the many-body density of states in terms of a trace formula starting from the exact path integral form of the propagator between many-body quadrature states. We therefore avoid the use of a complexified classical limit characteristic of the coherent state representation. While quantum effects such as vacuum fluctuations and gauge invariance are exactly accounted for, our semiclassical approach captures quantum interference and therefore is valid well beyond the Ehrenfest time where naive quantum-classical correspondence breaks down. Remarkably, due to a special feature of harmonic systems with incommensurable frequencies, our formulas are generically valid also in the free-field case of noninteracting bosons. PMID:26764774

  19. Three-site Bose-Hubbard model subject to atom losses: Boson-pair dissipation channel and failure of the mean-field approach

    SciTech Connect

    Shchesnovich, V. S.; Mogilevtsev, D. S.

    2010-10-15

    We employ the perturbation series expansion for derivation of the reduced master equations for the three-site Bose-Hubbard model subject to strong atom losses from the central site. The model describes a condensate trapped in a triple-well potential subject to externally controlled removal of atoms. We find that the {pi}-phase state of the coherent superposition between the side wells decays via two dissipation channels, the single-boson channel (similar to the externally applied dissipation) and the boson-pair channel. The quantum derivation is compared to the classical adiabatic elimination within the mean-field approximation. We find that the boson-pair dissipation channel is not captured by the mean-field model, whereas the single-boson channel is described by it. Moreover, there is a matching condition between the zero-point energy bias of the side wells and the nonlinear interaction parameter which separates the regions where either the single-boson or the boson-pair dissipation channel dominate. Our results indicate that the M-site Bose-Hubbard models, for M>2, subject to atom losses may require an analysis which goes beyond the usual mean-field approximation for correct description of their dissipative features. This is an important result in view of the recent experimental works on the single-site addressability of condensates trapped in optical lattices.

  20. Information geometry of mean-field approximation.

    PubMed

    Tanaka, T

    2000-08-01

    I present a general theory of mean-field approximation based on information geometry and applicable not only to Boltzmann machines but also to wider classes of statistical models. Using perturbation expansion of the Kullback divergence (or Plefka expansion in statistical physics), a formulation of mean-field approximation of general orders is derived. It includes in a natural way the "naive" mean-field approximation and is consistent with the Thouless-Anderson-Palmer (TAP) approach and the linear response theorem in statistical physics. PMID:10953246

  1. Superconductivity in the single-band Hubbard model: Mean-field treatment of slave-boson pairing

    SciTech Connect

    Kopp, T.; Seco, F.J.; Schiller, S.; Woelfle, P.

    1988-12-01

    We investigate the possibility of superconductivity in the large-U limit of the Hubbard model near half filling, using a slave-boson representation. The local constraint is treated in a novel way, incorporating short-range fluctuations in addition to the global constraint. We find slave-boson pairing and a superconducting extended s-wave phase at hole concentrations delta above 4% and below 32%.

  2. Construction of traveling clusters in the Hamiltonian mean-field model by nonequilibrium statistical mechanics and Bernstein-Greene-Kruskal waves.

    PubMed

    Yamaguchi, Yoshiyuki Y

    2011-07-01

    Traveling clusters are ubiquitously observed in the Hamiltonian mean-field model for a wide class of initial states, which are not predicted to become spatially inhomogeneous states by nonequilibrium statistical mechanics and by nonlinear Landau damping. To predict such a cluster state from a given initial state, we combine nonequilibrium statistical mechanics and a construction method of Bernstein-Greene-Kruskal (BGK) waves with the aid of phenomenological assumptions. The phenomenological theory is partially successful, and the theoretically constructed cluster states are in good agreement with N-body simulations. Robustness of the theory is also discussed for unsuccessful initial states. PMID:21867277

  3. Equivalence between the semirelativistic limit of the Dirac-Maxwell equations and the Breit-Pauli model in the mean-field approximation

    NASA Astrophysics Data System (ADS)

    Hinschberger, Y.; Dixit, A.; Manfredi, G.; Hervieux, P.-A.

    2015-01-01

    We demonstrate the equivalence between (i) the semirelativistic limit (up to second order in the inverse of the speed of light) of the self-consistent Dirac-Maxwell equations and (ii) the Breit-Pauli equations in the mean-field (Hartree-like) approximation. We explain how the charge and current densities that act as sources in the Dirac-Maxwell equations are related to the microscopic two-electron interactions of the Breit-Pauli model (spin orbit, spin-other-orbit, and spin-spin). The key role played by the second-order corrections to the charge density is clarified.

  4. Statics and dynamics of the ten-state mean-field Potts glass model: a Monte Carlo study

    NASA Astrophysics Data System (ADS)

    Brangian, Claudio; Kob, Walter; Binder, Kurt

    2002-01-01

    We investigate by means of Monte Carlo simulations the fully connected p-state Potts model for different system sizes in order to see how the static and dynamic properties of a finite model compare with the, exactly known, behaviour of the system in the thermodynamic limit. Using p = 10 we are able to study the equilibrium dynamics for system sizes as large as N = 2560. We find that the static quantities, such as the energy, the entropy, the spin glass susceptibility as well as the distribution of the order parameter P(q) show very strong finite-size effects. From P(q) we calculate the fourth-order cumulant g4(N,T) and the Guerra parameter G(N,T) and show that these quantities cannot be used to locate the static transition temperature for the system sizes investigated. Also the spin-autocorrelation function C(t) shows strong finite-size effects in that it does not show a plateau even for temperatures around the dynamical critical temperature TD. We show that the dependence on N and T of the α-relaxation time can be understood by means of a dynamical finite-size scaling ansatz. C(t) does not obey the time-temperature superposition principle for temperatures around TD, but does so for significantly lower T. Finally we study the relaxation dynamics of the individual spins and show that their dependence on time depends strongly on the chosen spin, i.e. that the system is dynamically very heterogeneous, which explains the non-exponentiality of C(t).

  5. Scale-free memory model for multiagent reinforcement learning. Mean field approximation and rock-paper-scissors dynamics

    NASA Astrophysics Data System (ADS)

    Lubashevsky, I.; Kanemoto, S.

    2010-07-01

    A continuous time model for multiagent systems governed by reinforcement learning with scale-free memory is developed. The agents are assumed to act independently of one another in optimizing their choice of possible actions via trial-and-error search. To gain awareness about the action value the agents accumulate in their memory the rewards obtained from taking a specific action at each moment of time. The contribution of the rewards in the past to the agent current perception of action value is described by an integral operator with a power-law kernel. Finally a fractional differential equation governing the system dynamics is obtained. The agents are considered to interact with one another implicitly via the reward of one agent depending on the choice of the other agents. The pairwise interaction model is adopted to describe this effect. As a specific example of systems with non-transitive interactions, a two agent and three agent systems of the rock-paper-scissors type are analyzed in detail, including the stability analysis and numerical simulation. Scale-free memory is demonstrated to cause complex dynamics of the systems at hand. In particular, it is shown that there can be simultaneously two modes of the system instability undergoing subcritical and supercritical bifurcation, with the latter one exhibiting anomalous oscillations with the amplitude and period growing with time. Besides, the instability onset via this supercritical mode may be regarded as “altruism self-organization”. For the three agent system the instability dynamics is found to be rather irregular and can be composed of alternate fragments of oscillations different in their properties.

  6. Constraining relativistic models through heavy ion collisions

    SciTech Connect

    Menezes, D. P.; Providencia, C.; Chiapparini, M.; Bracco, M. E.; Delfino, A.; Malheiro, M.

    2007-12-15

    Relativistic models can be successfully applied to the description of compact star properties in nuclear astrophysics as well as to nuclear matter and finite nuclei properties, these studies taking place at low and moderate temperatures. Nevertheless, all results are model dependent, and so far it is unclear whether some of them should be discarded. Moreover, in the regime of hot hadronic matter, very few calculations exist using these relativistic models, in particular when applied to particle yields in heavy ion collisions. A very important investigation is the simulation of a supernova explosion that is based on the construction of an adequate equation of state that needs to be valid within very large ranges of temperatures (0 to 100 MeV at least) and densities (very low to ten times the nuclear saturation density at least). In the present work, we comment on the known constraints that can help the selection of adequate models in this wide regime and investigate the main differences that arise when the particle production during a Au+Au collision at the BNL Relativistic Heavy Ion Collider is calculated with different relativistic models. We conclude that most of the models investigated in the present work give a very good overall description of the data and make predictions for not yet measured particle ratios.

  7. Impurity-induced magnetic moments on the graphene-lattice Hubbard model: An inhomogeneous cluster dynamical mean-field theory study

    NASA Astrophysics Data System (ADS)

    Charlebois, M.; Sénéchal, D.; Gagnon, A.-M.; Tremblay, A.-M. S.

    2015-01-01

    Defect-induced magnetic moments are at the center of the research effort on spintronic applications of graphene. Here, we study the problem of a nonmagnetic impurity in graphene with a new theoretical method, inhomogeneous cluster dynamical mean-field theory (I-CDMFT), which takes into account interaction-induced short-range correlations while allowing long-range inhomogeneities. The system is described by a Hubbard model on the honeycomb lattice. The impurity is modeled by a local potential. For a large enough potential, interactions induce local antiferromagnetic correlations around the impurity and a net total spin 1/2 appears, in agreement with Lieb's theorem. Bound states caused by the impurity are visible in the local density of states (LDOS) and have their energies shifted by interactions in a spin-dependent way, leading to the antiferromagnetic correlations. Our results take into account dynamical correlations; nevertheless they qualitatively agree with previous mean-field and density functional theory (DFT) studies. Moreover, they provide a relation between impurity potential and on-site repulsion U that could in principle be used to determine experimentally the value of U .

  8. Anomalous diffusion in the evolution of soccer championship scores: Real data, mean-field analysis, and an agent-based model

    NASA Astrophysics Data System (ADS)

    da Silva, Roberto; Vainstein, Mendeli H.; Gonçalves, Sebastián; Paula, Felipe S. F.

    2013-08-01

    Statistics of soccer tournament scores based on the double round robin system of several countries are studied. Exploring the dynamics of team scoring during tournament seasons from recent years we find evidences of superdiffusion. A mean-field analysis results in a drift velocity equal to that of real data but in a different diffusion coefficient. Along with the analysis of real data we present the results of simulations of soccer tournaments obtained by an agent-based model which successfully describes the final scoring distribution [da Silva , Comput. Phys. Commun.CPHCBZ0010-465510.1016/j.cpc.2012.10.030 184, 661 (2013)]. Such model yields random walks of scores over time with the same anomalous diffusion as observed in real data.

  9. A stochastically forced time delay solar dynamo model: Self-consistent recovery from a maunder-like grand minimum necessitates a mean-field alpha effect

    SciTech Connect

    Hazra, Soumitra; Nandy, Dibyendu; Passos, Dário E-mail: dariopassos@ist.utl.pt

    2014-07-01

    Fluctuations in the Sun's magnetic activity, including episodes of grand minima such as the Maunder minimum have important consequences for space and planetary environments. However, the underlying dynamics of such extreme fluctuations remain ill-understood. Here, we use a novel mathematical model based on stochastically forced, non-linear delay differential equations to study solar cycle fluctuations in which time delays capture the physics of magnetic flux transport between spatially segregated dynamo source regions in the solar interior. Using this model, we explicitly demonstrate that the Babcock-Leighton poloidal field source based on dispersal of tilted bipolar sunspot flux, alone, cannot recover the sunspot cycle from a grand minimum. We find that an additional poloidal field source effective on weak fields—e.g., the mean-field α effect driven by helical turbulence—is necessary for self-consistent recovery of the sunspot cycle from grand minima episodes.

  10. Neutron star cooling: A challenge to the nuclear mean field

    SciTech Connect

    Hoang Sy Than; Nguyen Van Giai

    2009-12-15

    The two recent density-dependent versions of the finite-range M3Y interaction (CDM3Yn and M3Y-Pn) have been probed against the bulk properties of asymmetric nuclear matter (NM) in the nonrelativistic Hartree-Fock (HF) formalism. The same HF study has also been done with the famous Skyrme (SLy4) and Gogny (D1S and D1N) interactions that were well tested in the nuclear structure calculations. Our HF results are compared with those given by other many-body calculations like the Dirac-Brueckner Hartree-Fock approach or ab initio variational calculations using free nucleon-nucleon interaction and by both the nonrelativistic and relativistic mean-field studies using different model parameters. Although the two considered density-dependent versions of the M3Y interaction were proven to be quite realistic in the nuclear structure or reaction studies, they give two distinct behaviors of the NM symmetry energy at high densities, like the Asy-soft and Asy-stiff scenarios found earlier with other mean-field interactions. As a consequence, we obtain two different behaviors of the proton fraction in the {beta}-equilibrium that in turn can imply two drastically different mechanisms for the neutron star cooling. While some preference of the Asy-stiff scenario was found based on predictions of the latest microscopic many-body calculations or empirical NM pressure and isospin diffusion data deduced from heavy-ion collisions, a consistent mean-field description of nuclear structure database is more often given by some Asy-soft type interaction like the Gogny or M3Y-Pn ones. Such a dilemma poses an interesting challenge to the modern mean-field approaches.

  11. Spatial Correlation in the Three-band Copper Oxide Model: Dynamical Mean-field Study with Configuration Interaction Based Impurity Solver

    NASA Astrophysics Data System (ADS)

    Go, Ara; Millis, Andrew J.

    2014-03-01

    The three-band copper oxide model is studied using the single-site and four-site dynamical mean-field theory with configuration interaction based impurity solver. Comparison of the single and four site approximations shows that short ranged antiferromagnetic correlations are crucial to the physics. In the undoped case, they increase the gap size, shift the metal-insulator phase boundary and enhance the conductivity at the gap edge. The relation of antiferromagnetism and the pseudogap is discussed for the doped case. The new solver permits the inclusion of more bath orbitals which are crucial for accurate studies of spectral properties near the gap edge. This work was supported by the US Department of Energy under Grants No. DOE FG02-04ER46169 and DE-SC0006613.

  12. Spinodal instabilities and the distillation effect in relativistic hadronic models

    SciTech Connect

    Avancini, S. S.; Menezes, D. P.; Brito, L.; Provide circumflex ncia, C.; Chomaz, Ph.

    2006-08-15

    Liquid-gas phase transitions in asymmetric nuclear matter give rise to a distillation effect that corresponds to the formation of droplets of high-density symmetric matter in a background of a neutron gas possibly with a very small fraction of protons. In the present work we test the model dependence of this effect. We study the spinodal instabilities of asymmetric nuclear matter within six different mean-field relativistic models with both constant and density-dependent coupling parameters. We also consider the effects of introducing the {delta} meson and the nonlinear {omega}-{rho} coupling. It is shown that the distillation effect within density-dependent models is not so efficient and is comparable to results obtained for nonrelativistic models. Thermodynamical instabilities of nuclear matter neutralized by electrons as found in stellar matter are also investigated. The high Fermi energy of electrons completely erases the instability of density-dependent models. The other models still show a small region of instability but the distillation effect completely disappears because the electron presence freezes the proton fluctuations.

  13. A mean field approach to watershed hydrology

    NASA Astrophysics Data System (ADS)

    Bartlett, Mark; Porporato, Amilcare

    2016-04-01

    Mean field theory (also known as self-consistent field theory) is commonly used in statistical physics when modeling the space-time behavior of complex systems. The mean field theory approximates a complex multi-component system by considering a lumped (or average) effect for all individual components acting on a single component. Thus, the many body problem is reduced to a one body problem. For watershed hydrology, a mean field theory reduces the numerous point component effects to more tractable watershed averages, resulting in a consistent method for linking the average watershed fluxes to the local fluxes at each point. We apply this approach to the spatial distribution of soil moisture, and as a result, the numerous local interactions related to lateral fluxes of soil water are parameterized in terms of the average soil moisture. The mean field approach provides a basis for unifying and extending common event-based models (e.g. Soil Conservation Service curve number (SCS-CN) method) with more modern semi-distributed models (e.g. Variable Infiltration Capacity (VIC) model, the Probability Distributed (PDM) model, and TOPMODEL). We obtain simple equations for the fractions of the different source areas of runoff, the spatial variability of runoff, and the average runoff value (i.e., the so-called runoff curve). The resulting space time distribution of soil moisture offers a concise description of the variability of watershed fluxes.

  14. Quantum phase transition in the spherical mean-field plus quadrupole-quadrupole and pairing model in a single-j shell

    NASA Astrophysics Data System (ADS)

    Li, Bo; Pan, Feng; Draayer, J. P.

    2016-04-01

    The quantum-phase-transitional behavior of the spherical shell-model mean field plus the geometric quadrupole-quadrupole and standard pairing model within a single-j shell is analyzed in detail. Various quantities, such as low-lying energy levels, some typical energy ratios, the overlaps of the excited states with those of the corresponding limiting cases, B (E 2 ) values and electric quadrupole moments of some low-lying states and their ratios, as functions of the control parameter of the model in a j =15 /2 shell are calculated as an example, in which only a crossover occurs due to the Pauli exclusion. The results show that there are noticeable changes within the crossover region of the rotation-like to the pair-excitation (superconducting) phase transition, especially in the half-filling case. As an application to realistic nuclear systems, a chain of isotones 212Rn-213Fr-214Ra-215Ac is chosen to be described by the model with valence protons in the 1 h9 /2 shell. As far as the low-lying energy levels, the experimentally observed B (E 2 ) values, and the electric quadrupole moment within the yrast band are concerned, these nuclei seem fitted reasonably well. The results indicate that these nuclei are all within the rotation-like to the pair-excitation phase transition near the crossover point.

  15. Shell corrections, magic numbers, and mean field

    SciTech Connect

    Denisov, V. Yu.

    2007-02-15

    It is shown that the positions of deep local minima of shell corrections associated with magic numbers in the region of superheavy nuclei depend on the parameters of the central and spin-orbit mean-field potentials. The accuracy of nuclear-mass predictions made within various models for superheavy nuclei is analyzed.

  16. Systematics of light nuclei in a relativistic model

    SciTech Connect

    Price, C.E.

    1988-01-01

    The results of relativistic mean field calculations for non-spherical nuclei are presented and discussed. The need for non-linear scalar meson self-couplings in order to describe the properties of s-d shell nuclei is emphasized along with the importance of self-consistency in calculations of magnetic moments of odd-mass nuclei. 16 refs., 3 figs., 2 tabs.

  17. "Phase diagram" of a mean field game

    NASA Astrophysics Data System (ADS)

    Swiecicki, Igor; Gobron, Thierry; Ullmo, Denis

    2016-01-01

    Mean field games were introduced by J-M. Lasry and P-L. Lions in the mathematical community, and independently by M. Huang and co-workers in the engineering community, to deal with optimization problems when the number of agents becomes very large. In this article we study in detail a particular example called the "seminar problem" introduced by O. Guéant, J-M. Lasry, and P-L. Lions in 2010. This model contains the main ingredients of any mean field game but has the particular feature that all agents are coupled only through a simple random event (the seminar starting time) that they all contribute to form. In the mean field limit, this event becomes deterministic and its value can be fixed through a self consistent procedure. This allows for a rather thorough understanding of the solutions of the problem, through both exact results and a detailed analysis of various limiting regimes. For a sensible class of initial configurations, distinct behaviors can be associated to different domains in the parameter space. For this reason, the "seminar problem" appears to be an interesting toy model on which both intuition and technical approaches can be tested as a preliminary study toward more complex mean field game models.

  18. A general relativistic model for free-fall absolute gravimeters

    NASA Astrophysics Data System (ADS)

    Tan, Yu-Jie; Shao, Cheng-Gang; Li, Jia; Hu, Zhong-Kun

    2016-04-01

    Although the relativistic manifestations of gravitational fields in gravimetry were first studied 40 years ago, the relativistic effects combined with free-fall absolute gravimeters have rarely been considered. In light of this, we present a general relativistic model for free-fall absolute gravimeters in a local-Fermi coordinates system, where we focus on effects related to the measuring devices: relativistic transverse Doppler effects, gravitational redshift effects and Earth’s rotation effects. Based on this model, a general relativistic expression of the measured gravity acceleration is obtained.

  19. Finite nuclei in relativistic models with a light chiral scalar meson

    SciTech Connect

    Serot, B.D.; Furnstahl, R.J.

    1993-10-01

    Relativistic chiral models with a light scalar, meson appear to provide an economical marriage of successful relativistic mean-field theories and chiral symmetry. In these models, the scalar meson serves as both the chiral partner of the pion and the mediator of the intermediate-range nucleon-nucleon (NN) attraction. However, while some of these models can reproduce the empirical nuclear matter saturation point, they fail to reproduce observed properties of finite nuclei, such as spin-orbit splittings, shell structure, charge densities, and surface energetics. There deficiencies imply that this realization of chiral symmetry is incorrect. An alternative scenario for chiral hadronic models, which features a heavy chiral scalar and dynamical generation of the NN attraction, is discussed.

  20. Periodic ordering of clusters and stripes in a two-dimensional lattice model. I. Ground state, mean-field phase diagram and structure of the disordered phases

    SciTech Connect

    Pekalski, J.; Ciach, A.; Almarza, N. G.

    2014-03-21

    The short-range attraction and long-range repulsion between nanoparticles or macromolecules can lead to spontaneous pattern formation on solid surfaces, fluid interfaces, or membranes. In order to study the self-assembly in such systems we consider a triangular lattice model with nearest-neighbor attraction and third-neighbor repulsion. At the ground state of the model (T = 0) the lattice is empty for small values of the chemical potential μ, and fully occupied for large μ. For intermediate values of μ periodically distributed clusters, bubbles, or stripes appear if the repulsion is sufficiently strong. At the phase coexistences between the vacuum and the ordered cluster phases and between the cluster and the lamellar (stripe) phases the entropy per site does not vanish. As a consequence of this ground state degeneracy, disordered fluid phases consisting of clusters or stripes are stable, and the surface tension vanishes. For T > 0 we construct the phase diagram in the mean-field approximation and calculate the correlation function in the self-consistent Brazovskii-type field theory.

  1. Dynamic mean field theory for lattice gas models of fluids confined in porous materials: Higher order theory based on the Bethe-Peierls and path probability method approximations

    SciTech Connect

    Edison, John R.; Monson, Peter A.

    2014-07-14

    Recently we have developed a dynamic mean field theory (DMFT) for lattice gas models of fluids in porous materials [P. A. Monson, J. Chem. Phys. 128(8), 084701 (2008)]. The theory can be used to describe the relaxation processes in the approach to equilibrium or metastable states for fluids in pores and is especially useful for studying system exhibiting adsorption/desorption hysteresis. In this paper we discuss the extension of the theory to higher order by means of the path probability method (PPM) of Kikuchi and co-workers. We show that this leads to a treatment of the dynamics that is consistent with thermodynamics coming from the Bethe-Peierls or Quasi-Chemical approximation for the equilibrium or metastable equilibrium states of the lattice model. We compare the results from the PPM with those from DMFT and from dynamic Monte Carlo simulations. We find that the predictions from PPM are qualitatively similar to those from DMFT but give somewhat improved quantitative accuracy, in part due to the superior treatment of the underlying thermodynamics. This comes at the cost of greater computational expense associated with the larger number of equations that must be solved.

  2. A New Relativistic Jet Model of Blazars

    NASA Astrophysics Data System (ADS)

    Webb, James; Benitez, Erika; Howard, Emily

    1998-11-01

    The subclass of Active galaxies called Blazars encompass the most intrinsically luminous and rapidly variable sources known to astrophysicists. Attempts to model these sources has largely been frustrated due in part to observational difficulties, but also due to the lack of theoretical models capable of explaining the different characteristics of the observed sources. Leading candidate models all incorporate a massive, rotating black hole which is accreting galactic material, with some of this material being ejected out the ratational axis of the hole in the form of relativistically expanding jets. These jets are thought to emit energy via the synchrotron process across the entire spectrum from radio frequences all the way through the GEV (sometimes TEV) gamma-ray frequencies. Attempts to model these sources with single relativistic jets has proven difficult. We present a new model which features concentric interacting jets that do a much better job of explaining the types of Blazars we observe. We also discuss ways of testing this new model against multifreuqency observations.

  3. Magnetic fluctuations and dynamics in the vicinity of quantum spin liquids: Cluster dynamical mean-field study of the Kitaev model

    NASA Astrophysics Data System (ADS)

    Yoshitake, Junki; Nasu, Joji; Motome, Yukitoshi

    The quantum spin liquid, which does not show any long-range ordering down to the lowest temperature, has attracted broad interest as a new quantum state of matter. Since the ground state of the Kitaev model was shown to be a quantum spin liquid in two dimensions, there has been an explosion in both theoretical and experimental studies. Nevertheless, dynamical properties at finite temperatures remain a challenge, despite the relevance to analysis of recent experiments for Ir and Ru compounds. In this contribution, we address this problem by using the cluster dynamical mean-field approximation, which we newly develop on the basis of the Majorana fermion representation. Using the continuous-time quantum Monte Carlo method for the impurity solver, we calculate the magnetic susceptibility, dynamical spin structure factor, and relaxation time in the nuclear magnetic resonance. We find that these quantities show peculiar temperature dependences in the paramagnetic state when approaching the quantum spin liquid by decreasing temperature, which reflects the fractionalization of quantum spins. We will discuss the results while changing the anisotropy and sign (ferro/antiferro) of the exchange interactions, in comparison with experiments.

  4. Cluster dynamical mean field theory study of antiferromagnetic transition in the square-lattice Hubbard model: Optical conductivity and electronic structure

    NASA Astrophysics Data System (ADS)

    Sato, Toshihiro; Tsunetsugu, Hirokazu

    2016-08-01

    We numerically study optical conductivity σ (ω ) near the "antiferromagnetic" phase transition in the square-lattice Hubbard model at half filling. We use a cluster dynamical mean field theory and calculate conductivity including vertex corrections and, to this end, we have reformulated the vertex corrections in the antiferromagnetic phase. We find that the vertex corrections change various important details in temperature and ω dependencies of conductivity in the square lattice, and this contrasts sharply the case of the Mott transition in the frustrated triangular lattice. Generally, the vertex corrections enhance variations in the ω dependence, and sharpen the Drude peak and a high-ω incoherent peak in the paramagnetic phase. They also enhance the dip in σ (ω ) at ω =0 in the antiferromagnetic phase. Therefore, the dc conductivity is enhanced in the paramagnetic phase and suppressed in the antiferromagnetic phase, but this change occurs slightly below the transition temperature. We also find a temperature region above the transition temperature in which the dc conductivity shows an insulating behavior but σ (ω ) retains the Drude peak, and this region is stabilized by the vertex corrections. We also investigate which fluctuations are important in the vertex corrections and analyze momentum dependence of the vertex function in detail.

  5. Field theory for size- and charge-asymmetric primitive model of ionic systems: mean-field stability analysis and pretransitional effects.

    PubMed

    Ciach, A; Góźdź, W T; Stell, G

    2007-05-01

    The primitive model of ionic systems is investigated within a field-theoretic description for the whole range of diameter-, lambda , and charge, Z ratios of the two ionic species. Two order parameters (OP) are identified. The relation of the OP's to physically relevant quantities is nontrivial. Each OP is a linear combination of the charge density and the number-density waves. Instabilities of the disordered phase associated with the two OP's are determined in the mean-field approximation (MF). In MF a gas-liquid separation occurs for any Z and lambda is not equal to 1 . In addition, an instability with respect to various types of periodic ordering of the two kinds of ions is found. Depending on lambda and Z , one or the other transition is metastable in different thermodynamic states. For sufficiently large size disparity we find a sequence of fluid-crystal-fluid transitions for the increasing volume fraction of ions, in agreement with experimental observations. The instabilities found in MF represent weak ordering of the most probable instantaneous states, and are identified with structural loci associated with pretransitional effects. PMID:17677071

  6. Interacting Turing-Hopf Instabilities Drive Symmetry-Breaking Transitions in a Mean-Field Model of the Cortex: A Mechanism for the Slow Oscillation

    NASA Astrophysics Data System (ADS)

    Steyn-Ross, Moira L.; Steyn-Ross, D. A.; Sleigh, J. W.

    2013-04-01

    Electrical recordings of brain activity during the transition from wake to anesthetic coma show temporal and spectral alterations that are correlated with gross changes in the underlying brain state. Entry into anesthetic unconsciousness is signposted by the emergence of large, slow oscillations of electrical activity (≲1Hz) similar to the slow waves observed in natural sleep. Here we present a two-dimensional mean-field model of the cortex in which slow spatiotemporal oscillations arise spontaneously through a Turing (spatial) symmetry-breaking bifurcation that is modulated by a Hopf (temporal) instability. In our model, populations of neurons are densely interlinked by chemical synapses, and by interneuronal gap junctions represented as an inhibitory diffusive coupling. To demonstrate cortical behavior over a wide range of distinct brain states, we explore model dynamics in the vicinity of a general-anesthetic-induced transition from “wake” to “coma.” In this region, the system is poised at a codimension-2 point where competing Turing and Hopf instabilities coexist. We model anesthesia as a moderate reduction in inhibitory diffusion, paired with an increase in inhibitory postsynaptic response, producing a coma state that is characterized by emergent low-frequency oscillations whose dynamics is chaotic in time and space. The effect of long-range axonal white-matter connectivity is probed with the inclusion of a single idealized point-to-point connection. We find that the additional excitation from the long-range connection can provoke seizurelike bursts of cortical activity when inhibitory diffusion is weak, but has little impact on an active cortex. Our proposed dynamic mechanism for the origin of anesthetic slow waves complements—and contrasts with—conventional explanations that require cyclic modulation of ion-channel conductances. We postulate that a similar bifurcation mechanism might underpin the slow waves of natural sleep and comment on the

  7. The Mean-Field Flux Pinning Theory

    NASA Astrophysics Data System (ADS)

    Stejic, George

    We develop the Mean-Field Flux Pinning Theory, designed to model the flux line lattice (FLL) as it interacts with itself, the flux pinning centers and the geometry of the superconductor. Like other mean-field theories, the mean-field flux pinning theory does not attempt to model the FLL completely. Instead, it utilizes a simplified model for the FLL, termed the mean-field FLL, in which the FLL is modelled as a continuous vector field rather than as discrete fluxons as in other theories. By so doing, the interactions of the FLL are greatly simplified and more easily modelled. One application of the mean-field flux pinning theory is to predict J_{c} from microstructural data, which we use to determine the optimal Nb-Ti microstructures with (1) alpha -Ti pinning centers and (2) Nb pinning centers. The microstructure is modelled on a grid in which the local values of T_{c} and kappa reflect the spatial distribution of the pinning centers and the superconductor. Using this model, we solve the G-L equations and calculate the pinning potential defined as the vortex free energy as a function of position. We conclude that the ideal Nb-Ti microstructure with alpha-Ti pinning centers would require 40 volume percent of alpha -Ti and have 6nm thick pinning centers. In the Nb pinning center case, the ideal microstructure requires 50 volume percent of Nb and would have 6nm pinning centers. Another application for the mean-field flux pinning theory is to model the FLL as it interacts with the penetrating magnetic fields within lambda of the superconducting surface. Using this theory, we study the effects of sample geometry on the FLL and J _{c} for the thin film geometry. We find that the FLL becomes increasingly distorted as the film thickness is reduced and that J_{c } increases sharply for dimensions less that lambda. These predictions are experimentally evaluated in Nb-Ti thin films. Our results show that J_{c} values as high as 1/3 of J_{d} and a strong orientational

  8. Modelling of relativistic laser-plasma interactions

    NASA Astrophysics Data System (ADS)

    Berwick, Stuart James

    In order to characterise the propagation and stability of linearly polarised laser pulses of arbitrary intensity interacting with underdense plasma, a one-dimensional, fully relativistic, covariant electron fluid model is derived. As a first step, the model is Lorentz transformed into a frame moving with the group velocity of the laser pulse. A linear instability analysis is undertaken which generates an infinite hierarchy of homogeneous mode-coupling equations describing the decay of the laser pump via stimulated Raman forward scattering (SRFS), stimulated Raman back scattering (SRBS) and the relativistic modulational instability (RMI). SRFS and RMI are seen to merge into a hybrid instability at high intensities (1>1018Wcm-2) and a 6-wave analysis (rather than the conventional 3 or 4-wave) is required to accurately predict growth. Next, an Eulerian fluid code is developed in order to evolve the full non- linear equations. The method of characteristics is used to integrate the electromagnetic wave equation and a predictor-corrector algorithm is used to integrate the equations of continuity and momentum. After testing, this code is used to simulate the propagation and stability of ultra-short (<200fs), 'table-top' and cos2 modulated laser pulses of relativistic intensities in underdense plasma. Comparison is made to the predictions of the dispersion relation and growth rates obtained in each case are reconciled. The spatiotemporal behaviour is discussed with reference to the results of a 3-wave WKB model of the interaction. The importance of seeding mechanisms, pulse shape and relativity on the evolution of the instabilities is also discussed.

  9. Finite nuclei in relativistic models with a light chiral scalar meson

    SciTech Connect

    Furnstahl, R.J. ); Serot, B.D. )

    1993-05-01

    Relativistic chiral models with a light scalar meson appear to provide an economical marriage of successful relativistic mean-field theories and chiral symmetry. The scalar meson serves as both the chiral partner of the pion and the mediator of the intermediate-range nucleon-nucleon ([ital NN]) attraction. However, while some of these models can reproduce the empirical nuclear matter saturation point, they fail to reproduce observed properties of finite nuclei, such as spin-orbit splittings, shell structure, charge densities, and surface energetics. These deficiencies imply that this realization of chiral symmetry is incorrect. An alternative scenario, which features a heavy chiral scalar and dynamical generation of the [ital NN] attraction, is discussed.

  10. Relativistic interactions and realistic applications

    SciTech Connect

    Hoch, T.; Madland, D.; Manakos, P.; Mannel, T.; Nikolaus, B.A.; Strottman, D. |

    1992-12-31

    A four-fermion-coupling Lagrangian (relativistic Skyrme-type) interaction has been proposed for relativistic nuclear structure calculations. This interaction, which has the merit of simplicity, is from the outset tailored as an effective interaction for relativistic Hartree-Fock calculations. Various extensions of such a model are discussed and compared with Walecka`s meson-nucleon mean field approach. We also present results of the calculation of nuclear ground state properties with an extended (density dependent) version of the four fermion interaction in a relativistic Hartree-Fock approximation.

  11. Exact evolution of discrete relativistic cosmological models

    SciTech Connect

    Clifton, Timothy; Tavakol, Reza; Gregoris, Daniele; Rosquist, Kjell E-mail: danielegregoris@libero.it E-mail: r.tavakol@qmul.ac.uk

    2013-11-01

    We study the effects of inhomogeneities on the evolution of the Universe, by considering a range of cosmological models with discretized matter content. This is done using exact and fully relativistic methods that exploit the symmetries in and about submanifolds of spacetimes that themselves possess no continuous global symmetries. These methods allow us to follow the evolution of our models throughout their entire history, far beyond what has previously been possible. We find that while some space-like curves collapse to anisotropic singularities in finite time, others remain non-singular forever. The resulting picture is of a cosmological spacetime in which some behaviour remains close to Friedmann-like, while other behaviours deviate radically. In particular, we find that large-scale acceleration is possible without any violation of the energy conditions.

  12. Relativistic constituent quark model with infrared confinement

    SciTech Connect

    Branz, Tanja; Faessler, Amand; Gutsche, Thomas; Lyubovitskij, Valery E.; Ivanov, Mikhail A.; Koerner, Juergen G.

    2010-02-01

    We refine the relativistic constituent quark model developed in our previous papers to include the confinement of quarks. It is done, first, by introducing the scale integration in the space of {alpha} parameters, and, second, by cutting this scale integration on the upper limit which corresponds to an infrared cutoff. In this manner one removes all possible thresholds present in the initial quark diagram. The cutoff parameter is taken to be the same for all physical processes. We adjust other model parameters by fitting the calculated quantities of the basic physical processes to available experimental data. As an application, we calculate the electromagnetic form factors of the pion and the transition form factors of the {omega} and {eta} Dalitz decays.

  13. Mean-field sparse optimal control

    PubMed Central

    Fornasier, Massimo; Piccoli, Benedetto; Rossi, Francesco

    2014-01-01

    We introduce the rigorous limit process connecting finite dimensional sparse optimal control problems with ODE constraints, modelling parsimonious interventions on the dynamics of a moving population divided into leaders and followers, to an infinite dimensional optimal control problem with a constraint given by a system of ODE for the leaders coupled with a PDE of Vlasov-type, governing the dynamics of the probability distribution of the followers. In the classical mean-field theory, one studies the behaviour of a large number of small individuals freely interacting with each other, by simplifying the effect of all the other individuals on any given individual by a single averaged effect. In this paper, we address instead the situation where the leaders are actually influenced also by an external policy maker, and we propagate its effect for the number N of followers going to infinity. The technical derivation of the sparse mean-field optimal control is realized by the simultaneous development of the mean-field limit of the equations governing the followers dynamics together with the Γ-limit of the finite dimensional sparse optimal control problems. PMID:25288818

  14. Solutions to the relativistic precession model

    NASA Astrophysics Data System (ADS)

    Ingram, Adam; Motta, Sara

    2014-11-01

    The relativistic precession model (RPM) can be used to obtain a precise measurement of the mass and spin of a black hole when the appropriate set of quasi-periodic oscillations is detected in the power-density spectrum of an accreting black hole. However, in previous studies, the solution of the RPM equations could be obtained only through numerical methods at a price of an intensive computational effort. Here, we demonstrate that the RPM system of equations can be solved analytically, drastically reducing the computational load, now limited to the Monte Carlo simulation necessary to estimate the uncertainties. The analytical method not only provides an easy solution to the RPM system when three oscillations are detected, but in all the cases where the detection of two simultaneous oscillations is coupled with an independent mass measurement. We also present a computationally inexpensive method to place limits on the black hole mass and spin when only two oscillations are observed.

  15. One Model of Singular Relativistic Lagrangian

    NASA Astrophysics Data System (ADS)

    Laserra, E.; Strianese, M.; Pavlotsky, I. P.

    The so-called no-interaction theorem of D.G. Currie, T.F. Jordan, E.C. Sudarshan, H. Leutwyler, G. Marmo and N. Mukunda makes it possible to construct relativistic quasiclassical particle dynamics in the post-Galilean approximation only. It is obtained that in this approximation the Lagrangians are singular on some surfaces of the phase space. These peculiarities have different physical and mathematical nature from the ones studied by P.M.A. Dirac, where Hessians vanish in the whole phase-space. The dynamical properties are essentially peculiar on the studied singular surfaces. These properties, in our opinion, can be related with some physical phenomena. In the present paper the general situation is first described. As the whole problem is large and difficult we examined in detail just typical one- and two-dimensional models in the second part of the paper.

  16. Lagrangian approach to the semirelativistic electron dynamics in the mean-field approximation

    NASA Astrophysics Data System (ADS)

    Dixit, Anant; Hinschberger, Yannick; Zamanian, Jens; Manfredi, Giovanni; Hervieux, Paul-Antoine

    2013-09-01

    We derive a mean-field model that is based on a two-component Pauli-like equation and incorporates quantum, spin, and relativistic effects up to second order in 1/c. Using a Lagrangian approach, we obtain the self-consistent charge and current densities that act as sources in the Maxwell equations. A physical interpretation is provided for the second-order corrections to the sources. The Maxwell equations are also expanded to the same order. The resulting self-consistent model constitutes a suitable semirelativistic approximation to the full Dirac-Maxwell equations.

  17. Kinetic analysis of thermally relativistic flow with dissipation. II. Relativistic Boltzmann equation versus its kinetic models

    NASA Astrophysics Data System (ADS)

    Yano, Ryosuke; Matsumoto, Jun; Suzuki, Kojiro

    2011-06-01

    Thermally relativistic flow with dissipation was analyzed by solving the rarefied supersonic flow of thermally relativistic matter around a triangle prism by Yano and Suzuki [Phys. Rev. DPRVDAQ1550-7998 83, 023517 (2011)10.1103/PhysRevD.83.023517], where the Anderson-Witting (AW) model was used as a solver. In this paper, we solve the same problem, which was analyzed by Yano and Suzuki, using the relativistic Boltzmann equation (RBE). To solve the RBE, the conventional direct simulation Monte Carlo method for the nonrelativistic Boltzmann equation is extended to a new direct simulation Monte Carlo method for the RBE. Additionally, we solve the modified Marle (MM) model proposed by Yano-Suzuki-Kuroda for comparisons. The solution of the thermally relativistic shock layer around the triangle prism obtained using the relativistic Boltzmann equation is considered by focusing on profiles of macroscopic quantities, such as the density, velocity, temperature, heat flux and dynamic pressure along the stagnation streamline (SSL). Differences among profiles of the number density, velocity and temperature along the SSL obtained using the RBE, the AW and MM. models are described in the framework of the relativistic Navier-Stokes-Fourier law. Finally, distribution functions on the SSL obtained using the RBE are compared with those obtained using the AW and MM models. The distribution function inside the shock wave obtained using the RBE does not indicate a bimodal form, which is obtained using the AW and MM models, but a smooth deceleration of thermally relativistic matter inside a shock wave.

  18. Lattice Boltzmann model for resistive relativistic magnetohydrodynamics.

    PubMed

    Mohseni, F; Mendoza, M; Succi, S; Herrmann, H J

    2015-08-01

    In this paper, we develop a lattice Boltzmann model for relativistic magnetohydrodynamics (MHD). Even though the model is derived for resistive MHD, it is shown that it is numerically robust even in the high conductivity (ideal MHD) limit. In order to validate the numerical method, test simulations are carried out for both ideal and resistive limits, namely the propagation of Alfvén waves in the ideal MHD and the evolution of current sheets in the resistive regime, where very good agreement is observed comparing to the analytical results. Additionally, two-dimensional magnetic reconnection driven by Kelvin-Helmholtz instability is studied and the effects of different parameters on the reconnection rate are investigated. It is shown that the density ratio has a negligible effect on the magnetic reconnection rate, while an increase in shear velocity decreases the reconnection rate. Additionally, it is found that the reconnection rate is proportional to σ-1/2, σ being the conductivity, which is in agreement with the scaling law of the Sweet-Parker model. Finally, the numerical model is used to study the magnetic reconnection in a stellar flare. Three-dimensional simulation suggests that the reconnection between the background and flux rope magnetic lines in a stellar flare can take place as a result of a shear velocity in the photosphere. PMID:26382548

  19. Lattice Boltzmann model for resistive relativistic magnetohydrodynamics

    NASA Astrophysics Data System (ADS)

    Mohseni, F.; Mendoza, M.; Succi, S.; Herrmann, H. J.

    2015-08-01

    In this paper, we develop a lattice Boltzmann model for relativistic magnetohydrodynamics (MHD). Even though the model is derived for resistive MHD, it is shown that it is numerically robust even in the high conductivity (ideal MHD) limit. In order to validate the numerical method, test simulations are carried out for both ideal and resistive limits, namely the propagation of Alfvén waves in the ideal MHD and the evolution of current sheets in the resistive regime, where very good agreement is observed comparing to the analytical results. Additionally, two-dimensional magnetic reconnection driven by Kelvin-Helmholtz instability is studied and the effects of different parameters on the reconnection rate are investigated. It is shown that the density ratio has a negligible effect on the magnetic reconnection rate, while an increase in shear velocity decreases the reconnection rate. Additionally, it is found that the reconnection rate is proportional to σ-1 / 2, σ being the conductivity, which is in agreement with the scaling law of the Sweet-Parker model. Finally, the numerical model is used to study the magnetic reconnection in a stellar flare. Three-dimensional simulation suggests that the reconnection between the background and flux rope magnetic lines in a stellar flare can take place as a result of a shear velocity in the photosphere.

  20. Mean Field Theory for Nonequilibrium Network Reconstruction

    NASA Astrophysics Data System (ADS)

    Roudi, Yasser; Hertz, John

    2011-01-01

    There has been recent progress on inferring the structure of interactions in complex networks when they are in stationary states satisfying detailed balance, but little has been done for nonequilibrium systems. Here we introduce an approach to this problem, considering, as an example, the question of recovering the interactions in an asymmetrically coupled, synchronously updated Sherrington-Kirkpatrick model. We derive an exact iterative inversion algorithm and develop efficient approximations based on dynamical mean-field and Thouless-Anderson-Palmer equations that express the interactions in terms of equal-time and one-time-step-delayed correlation functions.

  1. Optical-model potential in a relativistic quantum field model

    NASA Astrophysics Data System (ADS)

    Jaminon, M.; Mahaux, C.; Rochus, P.

    1980-11-01

    The average nucleon-nucleus potential at low and medium energy is investigated in the framework of a relativistic quantum field model. Using the same input parameters as Brockmann in his recent study of nuclear ground states, we calculate the self-consistent relativistic Hartree potential at positive energy in the case of infinite nuclear matter and of 16O and 40Ca. This potential is the sum of a scalar operator and of the fourth component of a vector operator. We construct its Schrödinger-equivalent potential by eliminating the small component of the Dirac spinor. The central part of this Schrödinger-equivalent potential is in fair agreement with empirical values at low and intermediate energy. Particular attention is paid to the intermediate energy domain, in which the calculated potential is repulsive in the nuclear interior and attractive at the nuclear surface. This is in keeping with some empirical evidence and is similar to results found in the framework of the nonrelativistic Brueckner-Hartree-Fock approximation. The spin-orbit potential of the relativistic Hartree model is also in good agreement with empirical values. NUCLEAR REACTIONS Calculated average nuclear field of nuclear matter, 16O and 40Ca at positive energy from relativistic Hartree approximation.

  2. Hot and dense hadronic matter in an effective mean-field approach

    SciTech Connect

    Lavagno, A.

    2010-04-15

    We investigate the equation of state of hadronic matter at finite values of baryon density and temperature reachable in high-energy heavy-ion collisions. The analysis is performed by requiring the Gibbs conditions on the global conservation of baryon number, electric charge fraction, and zero net strangeness. We consider an effective relativistic mean-field model with the inclusion of DELTA isobars, hyperons, and the lightest pseudoscalar and vector meson degrees of freedom. In this context, we study the influence of the DELTA-isobar degrees of freedom in the hadronic equation of state and, in connection, the behavior of different particle-antiparticle ratios and strangeness production.

  3. Chemical potential beyond the quasiparticle mean field

    SciTech Connect

    Dinh Dang, N.; Hung, N. Quang

    2010-03-15

    The effects of quantal and thermal fluctuations beyond the BCS quasiparticle mean field on the chemical potential are studied within a model, which consists of N particles distributed amongst {Omega} doubly folded equidistant levels interacting via a pairing force with parameter G. The results obtained at zero and finite temperatures T within several approaches, which include the fluctuations beyond the BCS theory, are compared with the exact results. The chemical potential, defined as the Lagrangian multiplier to preserve the average number of particles, is compared with the corresponding quantity, which includes the effect from fluctuations of particle and quasiparticle numbers beyond the BCS quasiparticle mean field. The analysis of the results shows that the latter differs significantly from the former as functions of G and T. The chemical potential loses its physical meaning in the system with a fixed number of particles or after eliminating quantal fluctuations of particle (quasiparticle) numbers by means of particle number projection. The validity of the criterion for the signature of the transition to Bose-Einstein condensation, which occurs in infinite systems when the chemical potential hits the bottom of the energy spectrum, is reexamined for the finite multilevel model.

  4. Quark mean field approach with derivative coupling for nuclear matter

    SciTech Connect

    Kawabata, M.; Akiyama, S.; Futami, Y.; Nakasone, T.; Yukino, T.

    2008-05-15

    We propose the quark mean field model including derivative coupling between quarks and scalar mesons in nuclear matter. This model concisely interprets an increasing size of the nucleon as well as a modification of coupling constant in the nuclear environment.

  5. Global relativistic folding optical potential and the relativistic Green's function model

    NASA Astrophysics Data System (ADS)

    Ivanov, M. V.; Vignote, J. R.; Álvarez-Rodríguez, R.; Meucci, A.; Giusti, C.; Udías, J. M.

    2016-07-01

    Optical potentials provide critical input for calculations on a wide variety of nuclear reactions, in particular, for neutrino-nucleus reactions, which are of great interest in the light of the new neutrino oscillation experiments. We present the global relativistic folding optical potential (GRFOP) fits to elastic proton scattering data from 12C nucleus at energies between 20 and 1040 MeV. We estimate observables, such as the differential cross section, the analyzing power, and the spin rotation parameter, in elastic proton scattering within the relativistic impulse approximation. The new GRFOP potential is employed within the relativistic Green's function model for inclusive quasielastic electron scattering and for (anti)neutrino-nucleus scattering at MiniBooNE kinematics.

  6. Machine Learning for Dynamical Mean Field Theory

    NASA Astrophysics Data System (ADS)

    Arsenault, Louis-Francois; Lopez-Bezanilla, Alejandro; von Lilienfeld, O. Anatole; Littlewood, P. B.; Millis, Andy

    2014-03-01

    Machine Learning (ML), an approach that infers new results from accumulated knowledge, is in use for a variety of tasks ranging from face and voice recognition to internet searching and has recently been gaining increasing importance in chemistry and physics. In this talk, we investigate the possibility of using ML to solve the equations of dynamical mean field theory which otherwise requires the (numerically very expensive) solution of a quantum impurity model. Our ML scheme requires the relation between two functions: the hybridization function describing the bare (local) electronic structure of a material and the self-energy describing the many body physics. We discuss the parameterization of the two functions for the exact diagonalization solver and present examples, beginning with the Anderson Impurity model with a fixed bath density of states, demonstrating the advantages and the pitfalls of the method. DOE contract DE-AC02-06CH11357.

  7. Relativistic Landau models and generation of fuzzy spheres

    NASA Astrophysics Data System (ADS)

    Hasebe, Kazuki

    2016-07-01

    Noncommutative geometry naturally emerges in low energy physics of Landau models as a consequence of level projection. In this work, we proactively utilize the level projection as an effective tool to generate fuzzy geometry. The level projection is specifically applied to the relativistic Landau models. In the first half of the paper, a detail analysis of the relativistic Landau problems on a sphere is presented, where a concise expression of the Dirac-Landau operator eigenstates is obtained based on algebraic methods. We establish SU(2) “gauge” transformation between the relativistic Landau model and the Pauli-Schrödinger nonrelativistic quantum mechanics. After the SU(2) transformation, the Dirac operator and the angular momentum operators are found to satisfy the SO(3, 1) algebra. In the second half, the fuzzy geometries generated from the relativistic Landau levels are elucidated, where unique properties of the relativistic fuzzy geometries are clarified. We consider mass deformation of the relativistic Landau models and demonstrate its geometrical effects to fuzzy geometry. Super fuzzy geometry is also constructed from a supersymmetric quantum mechanics as the square of the Dirac-Landau operator. Finally, we apply the level projection method to real graphene system to generate valley fuzzy spheres.

  8. Fractional Relativistic Yamaleev Oscillator Model and Its Dynamical Behaviors

    NASA Astrophysics Data System (ADS)

    Luo, Shao-Kai; He, Jin-Man; Xu, Yan-Li; Zhang, Xiao-Tian

    2016-07-01

    In the paper we construct a new kind of fractional dynamical model, i.e. the fractional relativistic Yamaleev oscillator model, and explore its dynamical behaviors. We will find that the fractional relativistic Yamaleev oscillator model possesses Lie algebraic structure and satisfies generalized Poisson conservation law. We will also give the Poisson conserved quantities of the model. Further, the relation between conserved quantities and integral invariants of the model is studied and it is proved that, by using the Poisson conserved quantities, we can construct integral invariants of the model. Finally, the stability of the manifold of equilibrium states of the fractional relativistic Yamaleev oscillator model is studied. The paper provides a general method, i.e. fractional generalized Hamiltonian method, for constructing a family of fractional dynamical models of an actual dynamical system.

  9. Coulomb Glass: a Mean Field Study

    NASA Astrophysics Data System (ADS)

    Mandra, Salvatore; Palassini, Matteo

    2012-02-01

    We study the Coulomb glass model of disordered localized electrons with long-range Coulomb interaction, which describes systems such as disordered insulators, granular metals, amorphous semiconductors, or doped crystalline semiconductors. Long ago Efros and Shklovskii showed that the long-range repulsion induces a soft Coulomb gap in the single particle density of states at low temperatures. Recent works suggested that this gap is associated to a transition to a glass phase, similar to the Almeida-Thouless transition in spin glasses. In this work, we use a mean field approach to characterize several physical properties of the Coulomb glass. In particular, following a seminal work of Bray and Moore, we show that the Edward-Anderson parameter qEA and the spin glass susceptibility χSG are directly related to spectrum distribution of the Hessian matrix around free energy minima. Using this result, we show that no glass transition is associated to the gap formation.

  10. Stochastic mean-field polycrystal plasticity methods

    NASA Astrophysics Data System (ADS)

    Tonks, Michael R.

    To accommodate multiple length scales, mean-field polycrystal plasticity models treat each material point as an aggregate of N crystals. The crystal velocity gradients Lc are approximated and then used to evaluate the crystal stresses T c. The Tc are averaged to determine the material point stress T. Commonly, the Lc are approximated with the fully constrained model (FCM) based on the Taylor hypothesis which equates Lc to the macro-scale velocity gradient L. Herein, we present two stochastic models that relax the FCM constraint. Through various applications we show that these computationally efficient stochastic models provide realistic response predictions. We first investigate the texture evolution in a planar polycrystal with our stochastic Taylor model (STM), in which we define L c as a realization of a normal distribution with mean equal to L. Our STM predictions agree with crystal plasticity finite element method (CPFEM) predictions, demonstrating the development of a steady-state texture that is not predicted by the FCM. The computational cost of the STM is comparable to the FCM, i.e. substantially less than the CPFEM. We develop the STM for 3-D polycrystals based on CPFEM analysis results which show that Lc follows a normal distribution. In addition to the STM, we develop the stochastic no-constraints model (SNCM), which differs from the STM in the manner with which the Lc distribution means are determined. Calibration and validation of the models are performed using tantalum compression experiment data. Both models predict the compression textures more accurately than the FCM, and the SNCM predicts them more accurately than the STM. The STM is slightly more computationally expensive than the FCM, while the SNCM is three times more expensive. Finally, we incorporate the STM in a finite element simulation of the Taylor impact of two tantalum specimens. Our simulation predictions mimic the texture and deformation data measured from a powder metallurgy

  11. Isomeric state in {sup 53}Co: A mean field analysis

    SciTech Connect

    Patra, S. K.; Bhat, F. H.; Panda, R. N.; Arumugam, P.; Gupta, Raj K.

    2009-04-15

    We study the ground and the first excited intrinsic states of {sup 53}Co and its mirror nucleus {sup 53}Fe, within the frameworks of the relativistic and nonrelativistic mean field formalisms. The analysis of the single-particle energy spectra of these nuclei show a competition of spins 1/2{sup -} and 3/2{sup -} in a low-lying excited state, which agrees well with the recent experimental observation [D. Rudolph et al., Eur. Phys. J. A 36, 131 (2008)] of spin and parity J{sup {pi}}=3/2{sup -} for the isomeric configuration in {sup 53}Co.

  12. Relativistic models for Gaia at the (cross)check-point

    NASA Astrophysics Data System (ADS)

    Bertone, S.; Le Poncin-Lafitte, C.; Crosta, M.; Vecchiato, A.; Minazzoli, O.; Angonin, M.-C.

    2013-11-01

    Given the extreme accuracy of modern space astrometry, a precise relativistic modeling of observations is required. Concerning light propagation, most approaches rely on the solution of the null-geodesic equations. However, another approach based on the Time Transfer Functions (TTF) can be used to define an astrometric observation using an integral-based method derived from the Synge World Function. The availability of several models, formulated in different and independent ways, is indeed a security against the presence of systematic errors in the analysis of future experimental results. It is the case of the forthcoming Gaia mission. In this work, we review the modeling of a Gaia-like astrometric observation using the TTF and two other approaches, namely the Gaia RElativistic Model (GREM) and the Relativistic Astrometric MODel (RAMOD), and we provide explicit relations between their characteristic quantities.

  13. Mean Field Analysis of Quantum Annealing Correction

    NASA Astrophysics Data System (ADS)

    Matsuura, Shunji; Nishimori, Hidetoshi; Albash, Tameem; Lidar, Daniel A.

    2016-06-01

    Quantum annealing correction (QAC) is a method that combines encoding with energy penalties and decoding to suppress and correct errors that degrade the performance of quantum annealers in solving optimization problems. While QAC has been experimentally demonstrated to successfully error correct a range of optimization problems, a clear understanding of its operating mechanism has been lacking. Here we bridge this gap using tools from quantum statistical mechanics. We study analytically tractable models using a mean-field analysis, specifically the p -body ferromagnetic infinite-range transverse-field Ising model as well as the quantum Hopfield model. We demonstrate that for p =2 , where the phase transition is of second order, QAC pushes the transition to increasingly larger transverse field strengths. For p ≥3 , where the phase transition is of first order, QAC softens the closing of the gap for small energy penalty values and prevents its closure for sufficiently large energy penalty values. Thus QAC provides protection from excitations that occur near the quantum critical point. We find similar results for the Hopfield model, thus demonstrating that our conclusions hold in the presence of disorder.

  14. Mean Field Analysis of Quantum Annealing Correction.

    PubMed

    Matsuura, Shunji; Nishimori, Hidetoshi; Albash, Tameem; Lidar, Daniel A

    2016-06-01

    Quantum annealing correction (QAC) is a method that combines encoding with energy penalties and decoding to suppress and correct errors that degrade the performance of quantum annealers in solving optimization problems. While QAC has been experimentally demonstrated to successfully error correct a range of optimization problems, a clear understanding of its operating mechanism has been lacking. Here we bridge this gap using tools from quantum statistical mechanics. We study analytically tractable models using a mean-field analysis, specifically the p-body ferromagnetic infinite-range transverse-field Ising model as well as the quantum Hopfield model. We demonstrate that for p=2, where the phase transition is of second order, QAC pushes the transition to increasingly larger transverse field strengths. For p≥3, where the phase transition is of first order, QAC softens the closing of the gap for small energy penalty values and prevents its closure for sufficiently large energy penalty values. Thus QAC provides protection from excitations that occur near the quantum critical point. We find similar results for the Hopfield model, thus demonstrating that our conclusions hold in the presence of disorder. PMID:27314705

  15. Nonrelativistic approaches derived from point-coupling relativistic models

    SciTech Connect

    Lourenco, O.; Dutra, M.; Delfino, A.; Sa Martins, J. S.

    2010-03-15

    We construct nonrelativistic versions of relativistic nonlinear hadronic point-coupling models, based on new normalized spinor wave functions after small component reduction. These expansions give us energy density functionals that can be compared to their relativistic counterparts. We show that the agreement between the nonrelativistic limit approach and the Skyrme parametrizations becomes strongly dependent on the incompressibility of each model. We also show that the particular case A=B=0 (Walecka model) leads to the same energy density functional of the Skyrme parametrizations SV and ZR2, while the truncation scheme, up to order {rho}{sup 3}, leads to parametrizations for which {sigma}=1.

  16. Mean-Field Games for Marriage

    PubMed Central

    Bauso, Dario; Dia, Ben Mansour; Djehiche, Boualem; Tembine, Hamidou; Tempone, Raul

    2014-01-01

    This article examines mean-field games for marriage. The results support the argument that optimizing the long-term well-being through effort and social feeling state distribution (mean-field) will help to stabilize marriage. However, if the cost of effort is very high, the couple fluctuates in a bad feeling state or the marriage breaks down. We then examine the influence of society on a couple using mean-field sentimental games. We show that, in mean-field equilibrium, the optimal effort is always higher than the one-shot optimal effort. We illustrate numerically the influence of the couple’s network on their feeling states and their well-being. PMID:24804835

  17. Deterministic Mean-Field Ensemble Kalman Filtering

    DOE PAGESBeta

    Law, Kody J. H.; Tembine, Hamidou; Tempone, Raul

    2016-05-03

    The proof of convergence of the standard ensemble Kalman filter (EnKF) from Le Gland, Monbet, and Tran [Large sample asymptotics for the ensemble Kalman filter, in The Oxford Handbook of Nonlinear Filtering, Oxford University Press, Oxford, UK, 2011, pp. 598--631] is extended to non-Gaussian state-space models. In this paper, a density-based deterministic approximation of the mean-field limit EnKF (DMFEnKF) is proposed, consisting of a PDE solver and a quadrature rule. Given a certain minimal order of convergence κ between the two, this extends to the deterministic filter approximation, which is therefore asymptotically superior to standard EnKF for dimension d model, which arises in both deterministic and standard EnKF. Numerical results support and extend the theory.« less

  18. The relativistic feedback discharge model of terrestrial gamma ray flashes

    NASA Astrophysics Data System (ADS)

    Dwyer, Joseph R.

    2012-02-01

    As thunderclouds charge, the large-scale fields may approach the relativistic feedback threshold, above which the production of relativistic runaway electron avalanches becomes self-sustaining through the generation of backward propagating runaway positrons and backscattered X-rays. Positive intracloud (IC) lightning may force the large-scale electric fields inside thunderclouds above the relativistic feedback threshold, causing the number of runaway electrons, and the resulting X-ray and gamma ray emission, to grow exponentially, producing very large fluxes of energetic radiation. As the flux of runaway electrons increases, ionization eventually causes the electric field to discharge, bringing the field below the relativistic feedback threshold again and reducing the flux of runaway electrons. These processes are investigated with a new model that includes the production, propagation, diffusion, and avalanche multiplication of runaway electrons; the production and propagation of X-rays and gamma rays; and the production, propagation, and annihilation of runaway positrons. In this model, referred to as the relativistic feedback discharge model, the large-scale electric fields are calculated self-consistently from the charge motion of the drifting low-energy electrons and ions, produced from the ionization of air by the runaway electrons, including two- and three-body attachment and recombination. Simulation results show that when relativistic feedback is considered, bright gamma ray flashes are a natural consequence of upward +IC lightning propagating in large-scale thundercloud fields. Furthermore, these flashes have the same time structures, including both single and multiple pulses, intensities, angular distributions, current moments, and energy spectra as terrestrial gamma ray flashes, and produce large current moments that should be observable in radio waves.

  19. Local excitations in mean-field spin glasses

    NASA Astrophysics Data System (ADS)

    Krzakala, F.; Parisi, G.

    2004-06-01

    We address the question of geometrical as well as energetic properties of local excitations in mean-field Ising spin glasses. We study analytically the Random Energy Model and numerically a dilute mean-field model, first on tree-like graphs, equivalent to a replica-symmetric computation, and then directly on finite-connectivity random lattices. In the first model, characterized by a discontinuous replica symmetry breaking, we found that the energy of finite-volume excitation is infinite, whereas in the dilute mean-field model, described by a continuous replica symmetry breaking, it slowly decreases with sizes and saturates at a finite value, in contrast with what would be naively expected. The geometrical properties of these excitations are similar to those of lattice animals or branched polymers. We discuss the meaning of these results in terms of replica symmetry breaking and also possible relevance in finite-dimensional systems.

  20. Relativistic reflection: Review and recent developments in modeling

    NASA Astrophysics Data System (ADS)

    Dauser, T.; García, J.; Wilms, J.

    2016-05-01

    Measuring relativistic reflection is an important tool to study the innermost regions of the an accreting black hole system. In the following we present a brief review on the different aspects contributing to the relativistic reflection. The combined approach is for the first time incorporated in the new ``relxill'' model. The advantages of this more self-consistent approach are briefly summarized. A special focus is put on the new definition of the intrinsic reflection fraction in the lamp post geometry, which allows to draw conclusions about the primary source of radiation in these system. Additionally the influence of the high energy cutoff of the primary source on the reflection spectrum is motivated, revealing the remarkable capabilities of constraining E_cut by measuring relativistic reflection spectra from NuSTAR, preferably with lower energy coverage.

  1. Modeling terrestrial gamma ray flashes produced by relativistic feedback discharges

    NASA Astrophysics Data System (ADS)

    Liu, Ningyu; Dwyer, Joseph R.

    2013-05-01

    This paper reports a modeling study of terrestrial gamma ray flashes (TGFs) produced by relativistic feedback discharges. Terrestrial gamma ray flashes are intense energetic radiation originating from the Earth's atmosphere that has been observed by spacecraft. They are produced by bremsstrahlung interactions of energetic electrons, known as runaway electrons, with air atoms. An efficient physical mechanism for producing large fluxes of the runaway electrons to make the TGFs is the relativistic feedback discharge, where seed runaway electrons are generated by positrons and X-rays, products of the discharge itself. Once the relativistic feedback discharge becomes self-sustaining, an exponentially increasing number of relativistic electron avalanches propagate through the same high-field region inside the thundercloud until the electric field is partially discharged by the ionization created by the discharge. The modeling results indicate that the durations of the TGF pulses produced by the relativistic feedback discharge vary from tens of microseconds to several milliseconds, encompassing all durations of the TGFs observed so far. In addition, when a sufficiently large potential difference is available in thunderclouds, a self-propagating discharge known as the relativistic feedback streamer can be formed, which propagates like a conventional positive streamer. For the relativistic feedback streamer, the positive feedback mechanism of runaway electron production by the positrons and X-rays plays a similar role as the photoionization for the conventional positive streamer. The simulation results of the relativistic feedback streamer show that a sequence of TGF pulses with varying durations can be produced by the streamer. The relativistic streamer may initially propagate with a pulsed manner and turn into a continuous propagation mode at a later stage. Milliseconds long TGF pulses can be produced by the feedback streamer during its continuous propagation. However

  2. Universality in bipartite mean field spin glasses

    NASA Astrophysics Data System (ADS)

    Genovese, Giuseppe

    2012-12-01

    In this work, we give a proof of universality with respect to the choice of the statistical distribution of the quenched noise, for mean field bipartite spin glasses. We use mainly techniques of spin glasses theory, as Guerra's interpolation and the cavity approach.

  3. Instabilities in the Mean Field Limit

    NASA Astrophysics Data System (ADS)

    Han-Kwan, Daniel; Nguyen, Toan T.

    2016-03-01

    Consider a system of N particles interacting through Newton's second law with Coulomb interaction potential in one spatial dimension or a {C}^2 smooth potential in any dimension. We prove that in the mean field limit N → + ∞, the N particles system displays instabilities in times of order log N, for some configurations approximately distributed according to unstable homogeneous equilibria.

  4. Relativistic HD and MHD modelling for AGN jets

    NASA Astrophysics Data System (ADS)

    Keppens, R.; Porth, O.; Monceau-Baroux, R.; Walg, S.

    2013-12-01

    Relativistic hydro and magnetohydrodynamics (MHD) provide a continuum fluid description for plasma dynamics characterized by shock-dominated flows approaching the speed of light. Significant progress in its numerical modelling emerged in the last two decades; we highlight selected examples of modern grid-adaptive, massively parallel simulations realized by our open-source software MPI-AMRVAC (Keppens et al 2012 J. Comput. Phys. 231 718). Hydrodynamical models quantify how energy transfer from active galactic nuclei (AGN) jets to their surrounding interstellar/intergalactic medium (ISM/IGM) gets mediated through shocks and various fluid instability mechanisms (Monceau-Baroux et al 2012 Astron. Astrophys. 545 A62). With jet parameters representative for Fanaroff-Riley type-II jets with finite opening angles, we can quantify the ISM volumes affected by jet injection and distinguish the roles of mixing versus shock-heating in cocoon regions. This provides insight in energy feedback by AGN jets, usually incorporated parametrically in cosmological evolution scenarios. We discuss recent axisymmetric studies up to full 3D simulations for precessing relativistic jets, where synthetic radio maps can confront observations. While relativistic hydrodynamic models allow one to better constrain dynamical parameters like the Lorentz factor and density contrast between jets and their surroundings, the role of magnetic fields in AGN jet dynamics and propagation characteristics needs full relativistic MHD treatments. Then, we can demonstrate the collimating properties of an overal helical magnetic field backbone and study differences between poloidal versus toroidal field dominated scenarios (Keppens et al 2008 Astron. Astrophys. 486 663). Full 3D simulations allow one to consider the fate of non-axisymmetric perturbations on relativistic jet propagation from rotating magnetospheres (Porth 2013 Mon. Not. R. Astron. Soc. 429 2482). Self-stabilization mechanisms related to the detailed

  5. An Extended Magnetohydrodynamics Model for Relativistic Weakly Collisional Plasmas

    NASA Astrophysics Data System (ADS)

    Chandra, Mani; Gammie, Charles F.; Foucart, Francois; Quataert, Eliot

    2015-09-01

    Black holes that accrete far below the Eddington limit are believed to accrete through a geometrically thick, optically thin, rotationally supported plasma that we will refer to as a radiatively inefficient accretion flow (RIAF). RIAFs are typically collisionless in the sense that the Coulomb mean free path is large compared to {GM}/{c}2, and relativistically hot near the event horizon. In this paper we develop a phenomenological model for the plasma in RIAFs, motivated by the application to sources such as Sgr A* and M87. The model is derived using Israel–Stewart theory, which considers deviations up to second order from thermal equilibrium, but modified for a magnetized plasma. This leads to thermal conduction along magnetic field lines and a difference in pressure, parallel and perpendicular to the field lines (which is equivalent to anisotropic viscosity). In the non-relativistic limit, our model reduces to the widely used Braginskii theory of magnetized, weakly collisional plasmas. We compare our model to the existing literature on dissipative relativistic fluids, describe the linear theory of the plasma, and elucidate the physical meaning of the free parameters in the model. We also describe limits of the model when the conduction is saturated and when the viscosity implies a large pressure anisotropy. In future work, the formalism developed in this paper will be used in numerical models of RIAFs to assess the importance of non-ideal processes for the dynamics and radiative properties of slowly accreting black holes.

  6. Kinematic arguments against single relativistic shell models for GRBs

    NASA Technical Reports Server (NTRS)

    Fenimore, Edward E.; Ramirez, E.; Sumner, M. C.

    1997-01-01

    Two main types of models have been suggested to explain the long durations and multiple peaks of Gamma Ray Bursts (GRBs). In one, there is a very quick release of energy at a central site resulting in a single relativistic shell that produces peaks in the time history through its interactions with the ambient material. In the other, the central site sporadically releases energy over hundreds of seconds forming a peak with each burst of energy. The authors show that the average envelope of emission and the presence of gaps in GRBs are inconsistent with a single relativistic shell. They estimate that the maximum fraction of a single shell that can produce gamma-rays in a GRB with multiple peaks is 10(exp (minus)3), implying that single relativistic shells require 10(exp 3) times more energy than previously thought. They conclude that either the central site of a GRB must produce (approx)10(exp 51) erg/s(exp (minus)1) for hundreds of seconds, or the relativistic shell must have structure on a scales the order of (radical)(epsilon)(Gamma)(exp (minus)1), where (Gamma) is the bulk Lorentz factor ((approximately)10(exp 2) to 10(exp 3)) and (epsilon) is the efficiency.

  7. Relativistic modeling capabilities in PERSEUS extended MHD simulation code for HED plasmas

    SciTech Connect

    Hamlin, Nathaniel D.; Seyler, Charles E.

    2014-12-15

    We discuss the incorporation of relativistic modeling capabilities into the PERSEUS extended MHD simulation code for high-energy-density (HED) plasmas, and present the latest hybrid X-pinch simulation results. The use of fully relativistic equations enables the model to remain self-consistent in simulations of such relativistic phenomena as X-pinches and laser-plasma interactions. By suitable formulation of the relativistic generalized Ohm’s law as an evolution equation, we have reduced the recovery of primitive variables, a major technical challenge in relativistic codes, to a straightforward algebraic computation. Our code recovers expected results in the non-relativistic limit, and reveals new physics in the modeling of electron beam acceleration following an X-pinch. Through the use of a relaxation scheme, relativistic PERSEUS is able to handle nine orders of magnitude in density variation, making it the first fluid code, to our knowledge, that can simulate relativistic HED plasmas.

  8. On Mean Field Limits for Dynamical Systems

    NASA Astrophysics Data System (ADS)

    Boers, Niklas; Pickl, Peter

    2016-07-01

    We present a purely probabilistic proof of propagation of molecular chaos for N-particle systems in dimension 3 with interaction forces scaling like 1/\\vert q\\vert ^{3λ - 1} with λ smaller but close to one and cut-off at q = N^{-1/3}. The proof yields a Gronwall estimate for the maximal distance between exact microscopic and approximate mean-field dynamics. This can be used to show weak convergence of the one-particle marginals to solutions of the respective mean-field equation without cut-off in a quantitative way. Our results thus lead to a derivation of the Vlasov equation from the microscopic N-particle dynamics with force term arbitrarily close to the physically relevant Coulomb- and gravitational forces.

  9. Mean-Field Evolution of Fermionic Systems

    NASA Astrophysics Data System (ADS)

    Benedikter, Niels; Porta, Marcello; Schlein, Benjamin

    2014-11-01

    The mean field limit for systems of many fermions is naturally coupled with a semiclassical limit. This makes the analysis of the mean field regime much more involved, compared with bosonic systems. In this paper, we study the dynamics of initial data close to a Slater determinant, whose reduced one-particle density is an orthogonal projection ω N with the appropriate semiclassical structure. Assuming some regularity of the interaction potential, we show that the evolution of such an initial data remains close to a Slater determinant, with reduced one-particle density given by the solution of the Hartree-Fock equation with initial data ω N . Our result holds for all (semiclassical) times, and gives effective bounds on the rate of the convergence towards the Hartree-Fock dynamics.

  10. Derivation of manifestly covariant quantum models for spinning relativistic particles

    NASA Astrophysics Data System (ADS)

    Marnelius, Robert; Mårtensson, Ulf

    1990-05-01

    A method to construct manifestly covariant models for relativistic spinning particles is given. The models involve manifestly covariant internal variables leading to discrete state spaces apart from the coordinate and momentum variables. The precise form of the models are extracted from the Bargmann-Wigner conditions on the Pauli-Lubanski operator, together with some consistency conditions. Several simple models are derived and analysed, some of which are new. Also manifestly conformally invariant models for particles with arbitrary spins are derived using a condition of Bracken and Jessup.