Numerical experiments of adjusted BSSN systems for controlling constraint violations

We present our numerical comparisons between the BSSN formulation widely used in numerical relativity today and its adjusted versions using constraints. We performed three testbeds: gauge-wave, linear wave, and Gowdy-wave tests, proposed by the Mexico workshop on the formulation problem of the Einstein equations. We tried three kinds of adjustments, which were previously proposed from the analysis of the constraint propagation equations, and investigated how they improve the accuracy and stability of evolutions. We observed that the signature of the proposed Lagrange multipliers are always right and the adjustments improve the convergence and stability of the simulations. When the original BSSN system already shows satisfactory good evolutions (e.g., linear wave test), the adjusted versions also coincide with those evolutions; while in some cases (e.g., gauge-wave or Gowdy-wave tests) the simulations using the adjusted systems last 10 times as long as those using the original BSSN equations. Our demonstrations imply a potential to construct a robust evolution system against constraint violations even in highly dynamical situations.Comment: to be published in PR

Luminosity Dependent Evolution of Lyman Break Galaxies from redshift 5 to 3

In this contribution we briefly describe our recent results on the properties of Lyman break galaxies at z~5 obtained from deep and wide blank field surveys using Subaru telescope, and through the comparison with samples at lower redshift ranges we discuss the evolution of star-forming galaxies in the early universe.Comment: 2 pages, 1 figure, for the proceedings of the IAU Symposium 235, Galaxies Across the Hubble Time, J. Palous & F. Combes, ed

Non-axisymmetric instabilities of neutron star with toroidal magnetic fields

The aim of this paper is to clarify the stabilities of neutron stars with strong toroidal magnetic fields against non-axisymmetric perturbation. The motivation comes from the fact that super magnetized neutron stars of $\sim 10^{15}$G, magnetars, and magnetized proto-neutron stars born after the magnetically-driven supernovae are likely to have such strong toroidal magnetic fields. Long-term, three-dimensional general relativistic magneto-hydrodynamic simulations are performed, preparing isentropic neutron stars with toroidal magnetic fields in equilibrium as initial conditions. To explore the effects of rotations on the stability, simulations are done for both non-rotating and rigidly rotating models. We find the emergence of the Parker and/or Tayler instabilities in both the non-rotating and rotating models. For both non-rotating and rotating models, the Parker instability is the primary instability as predicted by the local linear perturbation analysis. The interchange instability also appears in the rotating models. It is found that rapid rotation is not enough to suppress the Parker instability, and this finding does not agree with the perturbation analysis. The reason for this is that rigidly and rapidly rotating stars are marginally stable, and hence, in the presence of stellar pulsations by which the rotational profile is deformed, unstable regions with negative gradient of angular momentum profile is developed. After the onset of the instabilities, a turbulence is excited. Contrary to the axisymmetric case, the magnetic fields never reach an equilibrium state after the development of the turbulence. This conclusion suggests that three-dimensional simulation is indispensable for exploring the formation of magnetars or prominence activities of magnetars such as giant flares.Comment: 19 pages, 11 figures, to be published in A&

Gravitational waves and neutrino emission from the merger of binary neutron stars

Numerical simulations for the merger of binary neutron stars are performed in full general relativity incorporating a finite-temperature (Shen's) equation of state (EOS) and neutrino cooling for the first time. It is found that for this stiff EOS, a hypermassive neutron star (HMNS) with a long lifetime ($\gg 10$ ms) is the outcome for the total mass \alt 3.0M_{\odot}. It is shown that the typical total neutrino luminosity of the HMNS is $\sim 3$--$8\times 10^{53}$ ergs/s and the effective amplitude of gravitational waves from the HMNS is 4--$6 \times 10^{-22}$ at $f=2.1$--2.5 kHz for a source distance of 100 Mpc. We also present the neutrino luminosity curve when a black hole is formed for the first time.Comment: 4 pages, 4 figures (Fig.2 is in low resolution), Accepted for publication in PR

The Influence of Thermal Pressure on Equilibrium Models of Hypermassive Neutron Star Merger Remnants

The merger of two neutron stars leaves behind a rapidly spinning hypermassive object whose survival is believed to depend on the maximum mass supported by the nuclear equation of state, angular momentum redistribution by (magneto-)rotational instabilities, and spindown by gravitational waves. The high temperatures (~5-40 MeV) prevailing in the merger remnant may provide thermal pressure support that could increase its maximum mass and, thus, its life on a neutrino-cooling timescale. We investigate the role of thermal pressure support in hypermassive merger remnants by computing sequences of spherically-symmetric and axisymmetric uniformly and differentially rotating equilibrium solutions to the general-relativistic stellar structure equations. Using a set of finite-temperature nuclear equations of state, we find that hot maximum-mass critically spinning configurations generally do not support larger baryonic masses than their cold counterparts. However, subcritically spinning configurations with mean density of less than a few times nuclear saturation density yield a significantly thermally enhanced mass. Even without decreasing the maximum mass, cooling and other forms of energy loss can drive the remnant to an unstable state. We infer secular instability by identifying approximate energy turning points in equilibrium sequences of constant baryonic mass parametrized by maximum density. Energy loss carries the remnant along the direction of decreasing gravitational mass and higher density until instability triggers collapse. Since configurations with more thermal pressure support are less compact and thus begin their evolution at a lower maximum density, they remain stable for longer periods after merger.Comment: 20 pages, 12 figures. Accepted for publication in Ap

Evaluation of Trapped Magnetic Field Properties in Superconducting MgB2_2 Bulk Magnets of Various Shapes by Finite Element Method

The trapped magnetic field properties of superconducting MgB2 bulk magnets with various shapes such as a triangular, a quadrangular, a hexangular bulk were calculated by the Finite Elements Method (FEM). The effect for the combination of several numbers of bulks was also investigated for several kinds of shapes to obtain large area of bulk surface in spite of one large bulk. In this calculation, the simple magnetization process replaced by the field-cool magnetization was used to obtain the equivalent distribution of the magnetic field, and the thermal equation in FEM was omitted. The trapped magnetic field for the triangular bulk by FEM was compared with the experimental result. It was found that the calculated results agreed well with the experimental result. The maximum trapped magnetic field was obtained in the cylindrical shape among several kinds of shapes. The trapped magnetic field was increased by the combination of multi-bulks. It was confirmed that the trapped magnetic field of the multi-bulks was larger than that of the single bulk. The trapped magnetic field increases with increasing the number of the bulks.28th International Symposium on Superconductivity(ISS 2015), November 16-18, 2015, Tokyo, Japa

Flux pinning characteristics and irreversibility line in high temperature superconductors

The flux pinning properties in high temperature superconductors are strongly influenced by thermally activated flux motion. The scaling relation of the pinning force density and the irreversibility line in various high temperature superconductors are numerically analyzed in terms of the flux creep model. The effect of two factors, i.e., the flux pinning strength and the dimensionality of the material, on these properties are investigated. It is speculated that the irreversibility line in Bi-2212 superconductors is one order of magnitude smaller than that in Y-123, even if the flux pinning strength in Bi-2212 is improved up to the level of Y-123. It is concluded that these two factors are equally important in determination of the flux pinning characteristics at high temperatures