First-order quasilinear canonical representation of the characteristic formulation of the Einstein equations

We prescribe a choice of 18 variables in all that casts the equations of the fully nonlinear characteristic formulation of general relativity in first--order quasi-linear canonical form. At the analytical level, a formulation of this type allows us to make concrete statements about existence of solutions. In addition, it offers concrete advantages for numerical applications as it now becomes possible to incorporate advanced numerical techniques for first order systems, which had thus far not been applicable to the characteristic problem of the Einstein equations, as well as in providing a framework for a unified treatment of the vacuum and matter problems. This is of relevance to the accurate simulation of gravitational waves emitted in astrophysical scenarios such as stellar core collapse.Comment: revtex4, 7 pages, text and references added, typos corrected, to appear in Phys. Rev.

Gravitational waveforms with controlled accuracy

A partially first-order form of the characteristic formulation is introduced to control the accuracy in the computation of gravitational waveforms produced by highly distorted single black hole spacetimes. Our approach is to reduce the system of equations to first-order differential form on the angular derivatives, while retaining the proven radial and time integration schemes of the standard characteristic formulation. This results in significantly improved accuracy over the standard mixed-order approach in the extremely nonlinear post-merger regime of binary black hole collisions.Comment: Revised version, published in Phys. Rev. D, RevTeX, 16 pages, 4 figure

Energy spectrum of turbulent fluctuations in boundary driven reduced magnetohydrodynamics

The nonlinear dynamics of a bundle of magnetic flux ropes driven by stationary fluid motions at their endpoints is studied, by performing numerical simulations of the magnetohydrodynamic (MHD) equations. The development of MHD turbulence is shown, where the system reaches a state that is characterized by the ratio between the Alfven time (the time for incompressible MHD waves to travel along the field lines) and the convective time scale of the driving motions. This ratio of time scales determines the energy spectra and the relaxation toward different regimes ranging from weak to strong turbulence. A connection is made with phenomenological theories for the energy spectra in MHD turbulence.Comment: Published in Physics of Plasma

Flavor violating decays of the Higgs bosons in the THDM-III

We calculate the branching ratios for the decays of neutral Higgs bosons ($h^{0},H^{0},A^{0}$) into pairs of fermions, including flavor violating processes, in the context of the General Two Higgs Doublet Model III.Comment: 23 pages, 10 figures, 6 tables. Text clarifying equations and references added, typos correction

Pressure dependence of the melting mechanism at the limit of overheating in Lennard-Jones crystals

We study the pressure dependence of the melting mechanism of a surface free Lennard-Jones crystal by constant pressure Monte Carlo simulation. The difference between the overheating temperature($T_{OH}$) and the thermodynamical melting point($T_M$) increase for increasing pressure. When particles move into the repulsive part of the potential the properties at $T_{OH}$ change. There is a crossover pressure where the volume jump becomes pressure-independent. The overheating limit is pre-announced by thermal excitation of big clusters of defects. The temperature zone where the system is dominated by these big clusters of defects increases with increasing pressure. Beyond the crossover pressure we find that excitation of defects and clusters of them start at the same temperature scale related with $T_{OH}$.Comment: 6 pages, 5 figures. Accepted for publication in Physical Review

Montecarlo simulation of the role of defects as the melting mechanism

We study in this paper the melting transition of a crystal of fcc structure with the Lennard-Jones potential, by using isobaric-isothermal Monte Carlo simulations. Local and collective updates are sequentially used to optimize the convergence. We show the important role played by defects in the melting mechanism in favor of modern melting theories.Comment: 6 page, 10 figures included. Corrected version to appear in Phys. Rev.

Scalar field induced oscillations of neutron stars and gravitational collapse

We study the interaction of massless scalar fields with self-gravitating neutron stars by means of fully dynamic numerical simulations of the Einstein-Klein-Gordon perfect fluid system. Our investigation is restricted to spherical symmetry and the neutron stars are approximated by relativistic polytropes. Studying the nonlinear dynamics of isolated neutron stars is very effectively performed within the characteristic formulation of general relativity, in which the spacetime is foliated by a family of outgoing light cones. We are able to compactify the entire spacetime on a computational grid and simultaneously impose natural radiative boundary conditions and extract accurate radiative signals. We study the transfer of energy from the scalar field to the fluid star. We find, in particular, that depending on the compactness of the neutron star model, the scalar wave forces the neutron star either to oscillate in its radial modes of pulsation or to undergo gravitational collapse to a black hole on a dynamical timescale. The radiative signal, read off at future null infinity, shows quasi-normal oscillations before the setting of a late time power-law tail.Comment: 12 pages, 13 figures, submitted to Phys. Rev.

A Note On R-Parity Violation and Fermion Masses

We consider a class of supersymmetric SU(3)\times SU(2)\times U(1) multihiggs models in which R-parity is violated through bilinear Higgs-lepton interactions. The required, due to R-parity violation, higgs-lepton rotations introduce an alternative way to generate the phenomenologically desirable fermion mass matrix structures independently of the equality of Yukawas, possibly imposed by superstring or other unification.Comment: 8 pages, uses LaTeX2

Axisymmetric core collapse simulations using characteristic numerical relativity

We present results from axisymmetric stellar core collapse simulations in general relativity. Our hydrodynamics code has proved robust and accurate enough to allow for a detailed analysis of the global dynamics of the collapse. Contrary to traditional approaches based on the 3+1 formulation of the gravitational field equations, our framework uses a foliation based on a family of outgoing light cones, emanating from a regular center, and terminating at future null infinity. Such a coordinate system is well adapted to the study of interesting dynamical spacetimes in relativistic astrophysics such as stellar core collapse and neutron star formation. Perhaps most importantly this procedure allows for the unambiguous extraction of gravitational waves at future null infinity without any approximation, along with the commonly used quadrupole formalism for the gravitational wave extraction. Our results concerning the gravitational wave signals show noticeable disagreement when those are extracted by computing the Bondi news at future null infinity on the one hand and by using the quadrupole formula on the other hand. We have strong indication that for our setup the quadrupole formula on the null cone does not lead to physical gravitational wave signals. The Bondi gravitational wave signals extracted at infinity show typical oscillation frequencies of about 0.5 kHz.Comment: 17 pages, 18 figures, submitted to Phys. Rev.

The incorporation of matter into characteristic numerical relativity

A code that implements Einstein equations in the characteristic formulation in 3D has been developed and thoroughly tested for the vacuum case. Here, we describe how to incorporate matter, in the form of a perfect fluid, into the code. The extended code has been written and validated in a number of cases. It is stable and capable of contributing towards an understanding of a number of problems in black hole astrophysics.Comment: 15 pages + 4 (eps) figure