Nucleosynthesis in the accretion disks of type II collapsars

We investigate nucleosynthesis inside the gamma-ray burst (GRB) accretion disks formed by the Type II collapsars. In these collapsars, the core collapse of massive stars first leads to the formation of a proto-neutron star. After that, an outward moving shock triggers a successful supernova. However, the supernova ejecta lacks momentum and within a few seconds the newly formed neutron star gets transformed to a stellar mass black hole via massive fallback. The hydrodynamics of such an accretion disk formed from the fallback material of the supernova ejecta has been studied extensively in the past. We use these well-established hydrodynamic models for our accretion disk in order to understand nucleosynthesis, which is mainly advection dominated in the outer regions. Neutrino cooling becomes important in the inner disk where the temperature and density are higher. The higher the accretion rate (M) over dot is, the higher the density and temperature are in the disks. We deal with accretion disks with relatively low accretion rates: 0.001 M-circle dot s(-1) less than or similar to (M) over dot less than or similar to 0.01 M-circle dot s(-1) and hence these disks are predominantly advection dominated. We use He-rich and Si-rich abundances as the initial condition of nucleosynthesis at the outer disk, and being equipped with the disk hydrodynamics and the nuclear network code, we study the abundance evolution as matter inflows and falls into the central object. We investigate the variation in the nucleosynthesis products in the disk with the change in the initial abundance at the outer disk and also with the change in the mass accretion rate. We report the synthesis of several unusual nuclei like P-31, K-39, Sc-43, Cl-35 and various isotopes of titanium, vanadium, chromium, manganese and copper. We also confirm that isotopes of iron, cobalt, nickel, argon, calcium, sulphur and silicon get synthesized in the disk, as shown by previous authors. Much of these heavy elements thus synthesized are ejected from the disk via outflows and hence they should leave their signature in observed data.

Nucleosynthesis in the outflows associated with accretion disks of type II collapsars

We investigate nucleosynthesis inside the outflows from gamma-ray burst (GRB) accretion disks formed by the Type II collapsars. In these collapsars, massive stars undergo core collapse to form a proto-neutron star initially, and a mild supernova (SN) explosion is driven. The SN ejecta lack momentum, and subsequently this newly formed neutron star gets transformed to a stellar mass black hole via massive fallback. The hydrodynamics and the nucleosynthesis in these accretion disks have been studied extensively in the past. Several heavy elements are synthesized in the disk, and much of these heavy elements are ejected from the disk via winds and outflows. We study nucleosynthesis in the outflows launched from these disks by using an adiabatic, spherically expanding outflow model, to understand which of these elements thus synthesized in the disk survive in the outflow. While studying this, we find that many new elements like isotopes of titanium, copper, zinc, etc., are present in the outflows. Ni-56 is abundantly synthesized in most of the cases in the outflow, which implies that the outflows from these disks in a majority of cases will lead to an observable SN explosion. It is mainly present when outflow is considered from the He-rich, Ni-56/Fe-54-rich zones of the disks. However, outflow from the Si-rich zone of the disk remains rich in silicon. Although emission lines of many of these heavy elements have been observed in the X-ray afterglows of several GRBs by Chandra, BeppoSAX, XMM-Newton, etc., Swift seems to have not yet detected these lines.

Concentration Fluctuations and Segmental Dynamics in Weakly Dynamic Asymmetric Blends in the Presence of Surface-Functionalized Multiwall Carbon Nanotubes

Surface-functionalized multiwall carbon nanotubes (MWCNTs) are incorporated in poly(methyl methacrylate)/styrene acrylonitrile (PMMA/SAN) blends and the pretransitional regime is monitored in situ by melt rheology and dielectric spectroscopy. As the blends exhibit weak dynamic asymmetry, the obvious transitions in the melt rheology due to thermal concentration fluctuations are weak. This is further supported by the weak temperature dependence of the correlation length ( approximate to 10-12 angstrom) in the vicinity of demixing. Hence, various rheological techniques in both the temperature and frequency domains are adopted to evaluate the demixing temperature. The spinodal decomposition temperature is manifested in an increase in the miscibility gap in the presence of MWCNTs. Furthermore, MWCNTs lead to a significant slowdown of the segmental dynamics in the blends. Thermally induced phase separation in the PMMA/SAN blends lead to selective localization of MWCNTs in the PMMA phase. This further manifests itself in a significant increase in the melt conductivity.

Revisiting some physics issues related to the new mass limit for magnetized white dwarfs

We clarify important physics issues related to the recently established new mass limit for magnetized white dwarfs which is significantly super-Chandrasekhar. The issues include, justification of high magnetic field and the corresponding formation of stable white dwarfs, contribution of the magnetic field to the total density and pressure, flux freezing, variation of magnetic field and related currents therein. We also attempt to address the observational connection of such highly magnetized white dwarfs.

CONVERGENCE ANALYSIS OF THE LOWEST ORDER WEAKLY PENALIZED ADAPTIVE DISCONTINUOUS GALERKIN METHODS

In this article, we prove convergence of the weakly penalized adaptive discontinuous Galerkin methods. Unlike other works, we derive the contraction property for various discontinuous Galerkin methods only assuming the stabilizing parameters are large enough to stabilize the method. A central idea in the analysis is to construct an auxiliary solution from the discontinuous Galerkin solution by a simple post processing. Based on the auxiliary solution, we define the adaptive algorithm which guides to the convergence of adaptive discontinuous Galerkin methods.

Revised density of magnetized nuclear matter at the neutron drip line

We study the onset of the neutron drip in high-density matter in the presence of a magnetic field. It has been found that, for systems having only protons and electrons, in the presence of a magnetic field greater than or similar to 10(15) G, neutronization occurs at a density that is at least an order of magnitude higher compared to that in a nonmagnetic system. In a system with heavier ions, the effect of the magnetic field, however, starts arising at a much higher field, greater than or similar to 10(17) G. These results may have important implications for high-magnetic-field neutron stars and white dwarfs and, in general, in nuclear astrophysics when the system is embedded within a strong magnetic field.

Maximum mass of stable magnetized highly super-Chandrasekhar white dwarfs: stable solutions with varying magnetic fields

We address the issue of stability of recently proposed significantly super-Chandrasekhar white dwarfs. We present stable solutions of magnetostatic equilibrium models for super-Chandrasekhar white dwarfs pertaining to various magnetic field profiles. This has been obtained by self-consistently including the effects of the magnetic pressure gradient and total magnetic density in a general relativistic framework. We estimate that the maximum stable mass of magnetized white dwarfs could be more than 3 solar mass. This is very useful to explain peculiar, overluminous type Ia supernovae which do not conform to the traditional Chandrasekhar mass-limit.

Reply to ``Comment on `Strongly magnetized cold degenerate electron gas: Mass-radius relation of the magnetized white dwarf'''

We show that the upper bound for the central magnetic field of a super-Chandrasekhar white dwarf calculated by Nityananda and Konar Phys. Rev. D 89, 103017 (2014)] and in the concerned comment, by the same authors, against our work U. Das and B. Mukhopadhyay, Phys. Rev. D 86, 042001 (2012)] is erroneous. This in turn strengthens the argument in favor of the stability of the recently proposed magnetized super-Chandrasekhar white dwarfs. We also point out several other numerical errors in their work. Overall we conclude that the arguments put forth by Nityananda and Konar are misleading.

Cross-correlation-aided transport in stochastically driven accretion flows

The origin of linear instability resulting in rotating sheared accretion flows has remained a controversial subject for a long time. While some explanations of such non-normal transient growth of disturbances in the Rayleigh stable limit were available for magnetized accretion flows, similar instabilities in the absence of magnetic perturbations remained unexplained. This dichotomy was resolved in two recent publications by Chattopadhyay and co-workers Mukhopadhyay and Chattopadhyay, J. Phys. A 46, 035501 (2013); Nath et al., Phys. Rev. E 88, 013010 (2013)] where it was shown that such instabilities, especially for nonmagnetized accretion flows, were introduced through interaction of the inherent stochastic noise in the system (even a ``cold'' accretion flow at 3000Kis too ``hot'' in the statistical parlance and is capable of inducing strong thermal modes) with the underlying Taylor-Couette flow profiles. Both studies, however, excluded the additional energy influx (or efflux) that could result from nonzero cross correlation of a noise perturbing the velocity flow, say, with the noise that is driving the vorticity flow (or equivalently the magnetic field and magnetic vorticity flow dynamics). Through the introduction of such a time symmetry violating effect, in this article we show that nonzero noise cross correlations essentially renormalize the strength of temporal correlations. Apart from an overall boost in the energy rate (both for spatial and temporal correlations, and hence in the ensemble averaged energy spectra), this results in mutual competition in growth rates of affected variables often resulting in suppression of oscillating Alfven waves at small times while leading to faster saturations at relatively longer time scales. The effects are seen to be more pronounced with magnetic field fluxes where the noise cross correlation magnifies the strength of the field concerned. Another remarkable feature noted specifically for the autocorrelation functions is the removal of energy degeneracy in the temporal profiles of fast growing non-normal modes leading to faster saturation with minimum oscillations. These results, including those presented in the previous two publications, now convincingly explain subcritical transition to turbulence in the linear limit for all possible situations that could now serve as the benchmark for nonlinear stability studies in Keplerian accretion disks.

GRMHD formulation of highly super-Chandrasekhar magnetized white dwarfs: stable configurations of non-spherical white dwarfs

The topic of magnetized super-Chandrasekhar white dwarfs is in the limelight, particularly in the last few years, since our proposal of their existence. By full-scale general relativistic magnetohydrodynamic (GRMHD) numerical analysis, we confirm in this work the existence of stable, highly magnetized, significantly super-Chandrasekhar white dwarfs with mass more than 3 solar mass. While a poloidal field geometry renders the white dwarfs oblate, a toroidal field makes them prolate retaining an overall quasi-spherical shape, as speculated in our earlier work. These white dwarfs are expected to serve as the progenitors of over-luminous type Ia supernovae.

Origin of nonlinearity and plausible turbulence by hydromagnetic transient growth in accretion disks: Faster growth rate than magnetorotational instability

We investigate the evolution of hydromagnetic perturbations in a small section of accretion disks. It is known that molecular viscosity is negligible in accretion disks. Hence, it has been argued that a mechanism, known as magnetorotational instability (MRI), is responsible for transporting matter in the presence of a weak magnetic field. However, there are some shortcomings, which question the effectiveness of MRI. Now the question arises, whether other hydromagnetic effects, e.g., transient growth (TG), can play an important role in bringing nonlinearity into the system, even at weak magnetic fields. In addition, it should be determined whether MRI or TG is primarily responsible for revealing nonlinearity in order to make the flow turbulent. Our results prove explicitly that the flows with a high Reynolds number (Re), which is the case for realistic astrophysical accretion disks, exhibit nonlinearity via TG of perturbation modes faster than that by modes producing MRI. For a fixed wave vector, MRI dominates over transient effects only at low Re, lower than the value expected to be in astrophysical accretion disks, and low magnetic fields. This calls into serious question the (overall) persuasiveness of MRI in astrophysical accretion disks.

GRMHD formulation of highly super-Chandrasekhar rotating magnetized white dwarfs: stable configurations of non-spherical white dwarfs

Here we extend the exploration of significantly super-Chandrasekhar magnetized white dwarfs by numerically computing axisymmetric stationary equilibria of differentially rotating magnetized polytropic compact stars in general relativity (GR), within the ideal magnetohydrodynamic regime. We use a general relativistic magnetohydrodynamic (GRMHD) framework that describes rotating and magnetized axisymmetric white dwarfs, choosing appropriate rotation laws and magnetic field profiles (toroidal and poloidal). The numerical procedure for finding solutions in this framework uses the 3 + 1 formalism of numerical relativity, implemented in the open source XNS code. We construct equilibrium sequences by varying different physical quantities in turn, and highlight the plausible existence of super-Chandrasekhar white dwarfs, with masses in the range of 2-3 solar mass, with central (deep interior) magnetic fields of the order of 10(14) G and differential rotation with surface time periods of about 1-10 s. We note that such white dwarfs are candidates for the progenitors of peculiar, overluminous Type Ia supernovae, to which observational evidence ascribes mass in the range 2.1-2.8 solar mass. We also present some interesting results related to the structure of such white dwarfs, especially the existence of polar hollows in special cases.

Weakly nonlinear acoustic wave propagation in a nonlinear orthotropic circular cylindrical waveguide

Nonlinear acoustic wave propagation is considered in an infinite orthotropic thin circular cylindrical waveguide. The modes are non-planar having small but finite amplitude. The fluid is assumed to be ideal and inviscid with no mean flow. The cylindrical waveguide is modeled using the Donnell's nonlinear theory for thin cylindrical shells. The approximate solutions for the acoustic velocity potential are found using the method of multiple scales (MMS) in space and time. The calculations are presented up to the third order of the small parameter. It is found that at some frequencies the amplitude modulation is governed by the Nonlinear Schrodinger Equation (NLSE). The first objective is to study the nonlinear term in the NLSE, as the sign of the nonlinear term determines the stability of the amplitude modulation. On the other hand, at other specific frequencies, interactions occur between the primary wave and its higher harmonics. Here, the objective is to identify the frequencies of the higher harmonic interactions. Lastly, the linear terms in the NLSE obtained using the MMS calculations are validated. All three objectives are met using an asymptotic analysis of the dispersion equation. (C) 2015 Acoustical Society of America.

Weakly corrected numerical solutions to stochastically driven nonlinear dynamical systems

A method to weakly correct the solutions of stochastically driven nonlinear dynamical systems, herein numerically approximated through the Eule-Maruyama (EM) time-marching map, is proposed. An essential feature of the method is a change of measures that aims at rendering the EM-approximated solution measurable with respect to the filtration generated by an appropriately defined error process. Using Ito's formula and adopting a Monte Carlo (MC) setup, it is shown that the correction term may be additively applied to the realizations of the numerically integrated trajectories. Numerical evidence, presently gathered via applications of the proposed method to a few nonlinear mechanical oscillators and a semi-discrete form of a 1-D Burger's equation, lends credence to the remarkably improved numerical accuracy of the corrected solutions even with relatively large time step sizes. (C) 2015 Elsevier Inc. All rights reserved.