Riemann solvers in relativistic astrophysics

AbstractOur contribution reviews High Resolution Shock Capturing methods (HRSC) in the field of relativistic hydrodynamics with special emphasis on Riemann solvers. HRSC techniques achieve highly accurate numerical approximations (formally second order or better) in smooth regions of the flow, and capture the motion of unresolved steep gradients without creating spurious oscillations. One objective of our contribution is to show how these techniques have been extended to relativistic hydrodynamics, making it possible to explore some challenging astrophysical scenarios. We will review recent literature concerning the main properties of different special relativistic Riemann solvers, and discuss several 1D and 2D test problems which are commonly used to evaluate the performance of numerical methods in relativistic hydrodynamics. Finally, we will illustrate the use of HRSC methods in several applications in special and general relativistic hydrodynamics

Yukawa Unification, b --> s gamma and Bino-Stau Coannihilation

The minimal supersymmetric standard model with universal boundary conditions and "asymptotic" Yukawa unification is considered. The full one-loop effective potential for radiative electroweak symmetry breaking as well as the one-loop corrections to the charged Higgs boson, b-quark and tau lepton masses are included. The CP-even Higgs boson masses are corrected to two-loops. The relic abundance of the lightest supersymmetric particle (bino) is calculated by including its coannihilations with the next-to-lightest supersymmetric particle (lightest stau) consistently with Yukawa unification. The branching ratio of b --> s gamma is evaluated by incorporating all the applicable next-to-leading order QCD corrections. The bino-stau coannihilations reduce the bino relic abundance below the upper bound from cold dark matter considerations in a sizable fraction of the parameter space allowed by b --> s gamma for mu>0. Thus, the mu>0 case, which also predicts an acceptable b-quark mass, is perfectly compatible with data.Comment: 16 pages including 3 figures, Revtex, major revisions are mad

Noise-induced flow in quasigeostrophic turbulence with bottom friction

Randomly-forced fluid flow in the presence of scale-unselective dissipation develops mean currents following topographic contours. Known mechanisms based on the scale-selective action of damping processes are not at work in this situation. Coarse-graining reveals that the phenomenon is a kind of noise-rectification mechanism, in which lack of detailed balance and the symmetry-breaking provided by topography play an important role.Comment: 8 pages Revtex, no figures. Related material at http://www.imedea.uib.es

Maude: specification and programming in rewriting logic

AbstractMaude is a high-level language and a high-performance system supporting executable specification and declarative programming in rewriting logic. Since rewriting logic contains equational logic, Maude also supports equational specification and programming in its sublanguage of functional modules and theories. The underlying equational logic chosen for Maude is membership equational logic, that has sorts, subsorts, operator overloading, and partiality definable by membership and equality conditions. Rewriting logic is reflective, in the sense of being able to express its own metalevel at the object level. Reflection is systematically exploited in Maude endowing the language with powerful metaprogramming capabilities, including both user-definable module operations and declarative strategies to guide the deduction process. This paper explains and illustrates with examples the main concepts of Maude's language design, including its underlying logic, functional, system and object-oriented modules, as well as parameterized modules, theories, and views. We also explain how Maude supports reflection, metaprogramming and internal strategies. The paper outlines the principles underlying the Maude system implementation, including its semicompilation techniques. We conclude with some remarks about applications, work on a formal environment for Maude, and a mobile language extension of Maude

Cold Dark Matter and b --> s gamma in the Horava-Witten Theory

The minimal supersymmetric standard model with complete, partial or no Yukawa unification and radiative electroweak breaking with boundary conditions from the Horava-Witten theory is considered. The parameters are restricted by constraining the lightest sparticle relic abundance by cold dark matter considerations and requiring the b-quark mass after supersymmetric corrections and the branching ratio of b --> s gamma to be compatible with data. Complete Yukawa unification can be excluded. Also, t-b Yukawa unification is strongly disfavored since it requires almost degenerate lightest and next-to-lightest sparticle masses. However, the b-tau or no Yukawa unification cases avoid this degeneracy. The latter with mu<0 is the most natural case. The lightest sparticle mass, in this case, can be as low as about 77 GeV.Comment: 16 pages including 4 figures, Revtex, version to appear in Phys. Lett.

Uniqueness in MHD in divergence form: right nullvectors and well-posedness

Magnetohydrodynamics in divergence form describes a hyperbolic system of covariant and constraint-free equations. It comprises a linear combination of an algebraic constraint and Faraday's equations. Here, we study the problem of well-posedness, and identify a preferred linear combination in this divergence formulation. The limit of weak magnetic fields shows the slow magnetosonic and Alfven waves to bifurcate from the contact discontinuity (entropy waves), while the fast magnetosonic wave is a regular perturbation of the hydrodynamical sound speed. These results are further reported as a starting point for characteristic based shock capturing schemes for simulations with ultra-relativistic shocks in magnetized relativistic fluids.Comment: To appear in J Math Phy

Reformulation of Smoothed Particle Hydrodynamics with Riemann Solver

Smoothed Particle Hydrodynamics is reformulated in terms of the convolution of the original hydrodynamics equations, and the new evolution equations for the particles are derived. The same evolution equation of motion is also derived using a new action principle. The force acting on each particle is determined by solving the Riemann problem. The use of the Riemann Solver strengthens the method, making it accurate for the study of phenomena with strong shocks. The prescription for the variable smoothing length is shown. These techniques are implemented in strict conservation form. The results of a few test problems are also shown.Comment: 30 pages, 13 figure

Investigating Rare Events by Transition Interface Sampling

We briefly review simulation schemes for the investigation of rare transitions and we resume the recently introduced Transition Interface Sampling, a method in which the computation of rate constants is recast into the computation of fluxes through interfaces dividing the reactant and product state.Comment: 12 pages, 1 figure, contributed paper to the proceedings of NEXT 2003, Second Sardinian International Conference on News and Expectations in Thermostatistics, 21-28 Sep 2003, Cagliari (Italy

AMRA: An Adaptive Mesh Refinement Hydrodynamic Code for Astrophysics

Implementation details and test cases of a newly developed hydrodynamic code, AMRA, are presented. The numerical scheme exploits the adaptive mesh refinement technique coupled to modern high-resolution schemes which are suitable for relativistic and non-relativistic flows. Various physical processes are incorporated using the operator splitting approach, and include self-gravity, nuclear burning, physical viscosity, implicit and explicit schemes for conductive transport, simplified photoionization, and radiative losses from an optically thin plasma. Several aspects related to the accuracy and stability of the scheme are discussed in the context of hydrodynamic and astrophysical flows.Comment: 41 pages, 21 figures (9 low-resolution), LaTeX, requires elsart.cls, submitted to Comp. Phys. Comm.; additional documentation and high-resolution figures available from http://www.camk.edu.pl/~tomek/AMRA/index.htm

Unconventional low-energy SUSY from warped geometry

Supersymmetric models with a warped fifth spatial dimension can solve the hierarchy problem, avoiding some shortcomings of non-supersymmetric constructions, and predict a plethora of new phenomena at typical scales Lambda not far from the electroweak scale (Lambda ~ a few TeV). In this paper we derive the low-energy effective theories of these models, valid at energies below Lambda. We find that, in general, such effective theories can deviate significantly from the Minimal Supersymmetric Standard Model (MSSM) or other popular extensions of it, like the NMSSM: they have non-minimal Kaehler potentials (even in the Mp -> \infty limit), and the radion is coupled to the visible fields, both in the superpotential and the Kaehler potential, in a non-trivial (and quite model-independent) fashion. The corresponding phenomenology is pretty unconventional, in particular the electroweak breaking occurs in a non-radiative way, with tan beta \simeq 1 as a quite robust prediction, while the mass of the lightest Higgs boson can be as high as ~ 700 GeV.Comment: 53 pages, 2 ps figure