An algorithm for computing the 2D structure of fast rotating stars

Stars may be understood as self-gravitating masses of a compressible fluid whose radiative cooling is compensated by nuclear reactions or gravitational contraction. The understanding of their time evolution requires the use of detailed models that account for a complex microphysics including that of opacities, equation of state and nuclear reactions. The present stellar models are essentially one-dimensional, namely spherically symmetric. However, the interpretation of recent data like the surface abundances of elements or the distribution of internal rotation have reached the limits of validity of one-dimensional models because of their very simplified representation of large-scale fluid flows. In this article, we describe the ESTER code, which is the first code able to compute in a consistent way a two-dimensional model of a fast rotating star including its large-scale flows. Compared to classical 1D stellar evolution codes, many numerical innovations have been introduced to deal with this complex problem. First, the spectral discretization based on spherical harmonics and Chebyshev polynomials is used to represent the 2D axisymmetric fields. A nonlinear mapping maps the spheroidal star and allows a smooth spectral representation of the fields. The properties of Picard and Newton iterations for solving the nonlinear partial differential equations of the problem are discussed. It turns out that the Picard scheme is efficient on the computation of the simple polytropic stars, but Newton algorithm is unsurpassed when stellar models include complex microphysics. Finally, we discuss the numerical efficiency of our solver of Newton iterations. This linear solver combines the iterative Conjugate Gradient Squared algorithm together with an LU-factorization serving as a preconditionner of the Jacobian matrix.Comment: 40 pages, 12 figures, accepted in J. Comput. Physic

Physical processes leading to surface inhomogeneities: the case of rotation

In this lecture I discuss the bulk surface heterogeneity of rotating stars, namely gravity darkening. I especially detail the derivation of the omega-model of Espinosa Lara & Rieutord (2011), which gives the gravity darkening in early-type stars. I also discuss the problem of deriving gravity darkening in stars owning a convective envelope and in those that are members of a binary system.Comment: 23 pages, 11 figure, Lecture given to the school on the cartography of the Sun and the stars (May 2014 in Besan\c{c}on), to appear in LNP, Neiner and Rozelot edts V2: typos correcte

Neutrino Masses in the Lee-Wick Standard Model

Recently, an extension of the standard model based on ideas of Lee and Wick has been discussed. This theory is free of quadratic divergences and hence has a Higgs mass that is stable against radiative corrections. Here, we address the question of whether or not it is possible to couple very heavy particles, with masses much greater than the weak scale, to the Lee-Wick standard model degrees of freedom and still preserve the stability of the weak scale. We show that in the LW-standard model the familiar see-saw mechanism for generating neutrino masses preserves the solution to the hierarchy puzzle provided by the higher derivative terms. The very heavy right handed neutrinos do not destabilize the Higgs mass. We give an example of new heavy degrees of freedom that would destabilize the hierarchy, and discuss a general mechanism for coupling other heavy degrees of freedom to the Higgs doublet while preserving the hierarchy.Comment: 7 pages, 1 figur

The baryogenesis window in the MSSM

Thermal two-loop QCD corrections associated with light stops have a dramatic effect on the strength of the MSSM electroweak phase transition, making it more strongly first order as required for the viability of electroweak baryogenesis. We perform a perturbative analysis of the transition strength in this model, including these important contributions, extending previous work to arbitrary values of the pseudoscalar Higgs boson mass, m_A. We find a strong enough transition in a region with 2 120 GeV, a light Higgs boson with nearly standard couplings, and mass below 85 GeV within the reach of LEP II, and one stop not much heavier than the top quark. In addition, we give a qualitative discussion of the parameter space dependence of the transition strength and comment on the possibility that the transition turns to a crossover for sufficiently large Higgs masses.Comment: 33 pages, latex2e, 5 figures, epsfig.sty. Final version to appear in Nuclear Physics

Moving embedded lattice solitons

It was recently proved that isolated unstable "embedded lattice solitons" (ELS) may exist in discrete systems. The discovery of these ELS gives rise to relevant questions such as the following: are there continuous families of ELS?, can ELS be stable?, is it possible for ELS to move along the lattice?, how do ELS interact?. The present work addresses these questions by showing that a novel differential-difference equation (a discrete version of a complex mKdV equation) has a two-parameter continuous family of exact ELS. The numerical tests reveal that these solitons are stable and robust enough to withstand collisions. The model may apply to the description of a Bose-Einstein condensate with dipole-dipole interactions between the atoms, trapped in a deep optical-lattice potential.Comment: 13 pages, 11 figure

Black holes and Higgs stability

We study the effect of primordial black holes on the classical rate of nucleation of AdS regions within the standard electroweak vacuum. We find that the energy barrier for transitions to the new vacuum, which characterizes the exponential suppression of the nucleation rate, can be reduced significantly in the black-hole background. A precise analysis is required in order to determine whether the the existence of primordial black holes is compatible with the form of the Higgs potential at high temperature or density in the Standard Model or its extensions.Comment: 27 pages, 10 figures, conclusions expanded, to appear in JCA

Economic analysis of potato, corn and wheat response to nitrogen and phosphorus application in the highlands of Ecuador

Cover title.Includes bibliographical references (page 39)

Condensation in an Economic Model with Brand Competition

We present a linear agent based model on brand competition. Each agent belongs to one of the two brands and interacts with its nearest neighbors. In the process the agent can decide to change to the other brand if the move is beneficial. The numerical simulations show that the systems always condenses into a state when all agents belong to a single brand. We study the condensation times for different parameters of the model and the influence of different mechanisms to avoid condensation, like anti monopoly rules and brand fidelity.Comment: Accepted in: International Journal of Modern Physics