Titanium hidden in dust Iyudin, A F; Müller, E; Obergaulinger, M
Monthly notices of the Royal Astronomical Society,
05/2019, Letnik:
485, Številka:
3
Journal Article
Recenzirano
Odprti dostop
Abstract
Cassiopeia A, one of the most intriguing Galactic supernova remnants, has been a target of many observational efforts including most recent observations by Atacama Large ...Millimeter/submillimeter Array (ALMA), Hubble, Herschel, Spitzer, NuSTAR, INTEGRAL, and other observatories. We use recent gamma-ray lines observations of the radioactive products of Cas A supernova explosive nucleosynthesis and spectral energy densities derived for Cas A at infrared wavelengths to speculate about the possibility of radioactive ^{44}_{}\mathrm{Ti}$ being locked into large dust grains. This suggestion is also supported by the possible observation of a pre-supernova outburst about 80 yr before the actual Cas A supernova explosion in 1671 ad by Italian astronomer G. D. Cassini. The plausibility of such a scenario is discussed also with reference to recent supernovae, and to the contribution of core-collapse supernovae to the overall dust production in the Galaxy.
Context. Possible effects of magnetic fields in core collapse supernovae rely on an efficient amplification of the weak pre-collapse fields. It has been suggested that the magneto-rotational ...instability (MRI) leads to a rapid growth for these weak seed fields. Although plenty of MRI studies exist for accretion disks, the application of their results to core collapse supernovae is inhibited as the physics of supernova cores is substantially different from that of accretion discs. Aims. We address the problem of growth and saturation of the MRI in core collapse supernovae by studying its evolution by means of semi-global simulations, which combine elements of global and local simulations by taking the presence of global background gradients into account and using a local computational grid. We investigate, in particular, the termination of the growth of the MRI and the properties of the turbulence in the saturated state. Methods. We analyze the dispersion relation of the MRI to identify different regimes of the instability. This analysis is complemented by semi-global ideal MHD simulations, where we consider core matter in a local computational box (size ∼$1\,$km) rotating at sub-Keplerian velocity and where we allow for the presence of a radial entropy gradient, but neglect neutrino radiation. Results. We identify six regimes of the MRI depending on the ratio of the entropy and angular velocity gradient. Our numerical models confirm the instability criteria and growth rates for all regimes relevant to core-collapse supernovae. The MRI grows exponentially on time scales of milliseconds, the flow and magnetic field geometries being dominated by channel flows. We find MHD turbulence and efficient transport of angular momentum. The MRI growth ceases once the channels are disrupted by resistive instabilities (stemming from to the finite conductivity of the numerical code), and MHD turbulence sets in. From an analysis of the growth rates of the resistive instabilities, we deduce scaling laws for the termination amplitude of the MRI, which agree well with our numerical models. We determine the dependence of the development of large-scale coherent flow structures in the saturated state on the aspect ratio of the simulation boxes. Conclusions. The MRI can grow rapidly under the conditions considered here, i.e., a rapidly rotating core in hydrostatic equilibrium, possibly endowed with a nonvanishing entropy gradient, leading to fields exceeding $10^{15}~\mathrm{G}$. More investigations are required to cover the parameter space more comprehensively and to include more physical effects.
Context. Global magnetohydrodynamic simulations show the growth of Kelvin-Helmholtz instabilities at the contact surface of two merging neutron stars. That region has been identified as the site of ...efficient amplification of magnetic fields. However, these global simulations, due to numerical limitations, were unable to determine the saturation level of the field strength, and thus the possible back-reaction of the magnetic field onto the flow. Aims. We investigate the amplification of initially weak magnetic fields in Kelvin-Helmholtz unstable shear flows, and the back-reaction of the field onto the flow. Methods. We use a high-resolution finite-volume ideal MHD code to perform 2D and 3D local simulations of hydromagnetic shear flows, both for idealized systems and simplified models of merger flows. Results. In 2D, the magnetic field is amplified on time scales of less than 0.01 ms until it reaches locally equipartition with the kinetic energy. Subsequently, it saturates due to resistive instabilities that disrupt the Kelvin-Helmholtz unstable vortex and decelerate the shear flow on a secular time scale. We determine scaling laws of the field amplification with the initial field strength and the grid resolution. In 3D, the hydromagnetic mechanism seen in 2D may be dominated by purely hydrodynamic instabilities leading to less filed amplification. We find maximum magnetic fields ~1016 G locally, and rms maxima within the box ~1015 G. However, due to the fast decay of the shear flow such strong fields exist only for a short period (<0.1 ms). In the saturated state of most models, the magnetic field is mainly oriented parallel to the shear flow for rather strong initial fields, while weaker initial fields tend to lead to a more balanced distribution of the field energy among the components. In all models the flow shows small-scale features. The magnetic field is at most in energetic equipartition with the decaying shear flow. Conclusions. The magnetic field may be amplified efficiently to very high field strengths, the maximum field energy reaching values of the order of the kinetic energy associated with the velocity components transverse to the interface between the two neutron stars. However, the dynamic impact of the field onto the flow is limited to the shear layer, and it may not be adequate to produce outflows, because the time during which the magnetic field stays close to its maximum value is short compared to the time scale for launching an outflow (i.e., a few milliseconds).
The Kelvin-Helmholtz instability serves as a simple, well-defined setup for assessing the accuracy of different numerical methods for solving the equations of hydrodynamics. We use it to extend our ...previous analysis of the convergence and the numerical dissipation in models of the propagation of waves and in the tearing-mode instability in magnetohydrodynamic models. To this end, we perform two-dimensional simulations with and without explicit physical viscosity at different resolutions. A comparison of the growth of the modes excited by our initial perturbations allows us to estimate the effective numerical viscosity of two spatial reconstruction schemes (fifth-order monotonicity preserving and second-order piecewise linear schemes).
Measuring the eye's mechanical properties in vivo and with minimally invasive techniques can be the key for individualized solutions to a number of eye pathologies. The development of such techniques ...largely relies on a computational modelling of the eyeball and, it optimally requires the synergic interplay between experimentation and numerical simulation. In Astrophysics and Geophysics the remote measurement of structural properties of the systems of their realm is performed on the basis of (helio-)seismic techniques. As a biomechanical system, the eyeball possesses normal vibrational modes encompassing rich information about its structure and mechanical properties. However, the integral analysis of the eyeball vibrational modes has not been performed yet. Here we develop a new finite difference method to compute both the spheroidal and, specially, the toroidal eigenfrequencies of the human eye. Using this numerical model, we show that the vibrational eigenfrequencies of the human eye fall in the interval 100 Hz-10 MHz. We find that compressible vibrational modes may release a trace on high frequency changes of the intraocular pressure, while incompressible normal modes could be registered analyzing the scattering pattern that the motions of the vitreous humour leave on the retina. Existing contact lenses with embebed devices operating at high sampling frequency could be used to register the microfluctuations of the eyeball shape we obtain. We advance that an inverse problem to obtain the mechanical properties of a given eye (e.g., Young's modulus, Poisson ratio) measuring its normal frequencies is doable. These measurements can be done using non-invasive techniques, opening very interesting perspectives to estimate the mechanical properties of eyes in vivo. Future research might relate various ocular pathologies with anomalies in measured vibrational frequencies of the eye.
Aims.We have performed a comprehensive parameter study of the collapse of rotating, strongly magnetized stellar cores in axisymmetry to determine their gravitational wave signature based on the ...Einstein quadrupole formula. Methods.We use a Newtonian explicit magnetohydrodynamic Eulerian code based on the relaxing-TVD method for the solution of the ideal MHD equations, and apply the constraint-transport method to guarantee a divergence-free evolution of the magnetic field. We neglect effects due to neutrino transport and employ a simplified equation of state. The initial models are polytropes in rotational equilibrium with a prescribed degree of differential rotation and rotational energy. The initial magnetic fields are purely poloidal the field strength ranging from $10^{10}~{\rm G}$ to $10^{13}~{\rm G}$. The evolution of the core is followed until a few ten milliseconds past core bounce. Results.The initial magnetic fields are amplified mainly by the differential rotation of the core giving rise to a strong toroidal field component with an energy comparable to the rotational energy. The poloidal field component grows by compression during collapse, but does not change significantly after core bounce. In large parts of the simulated cores the growth time of the magneto-rotational instability (MRI) is of the order of a few milliseconds. The saturation field strengths that can be reached both via a pure Ω dynamo or the MRI are of the order of $10^{15}~{\rm G}$ at the surface of the core. Sheet-like circulation flows which produce a strong poloidal field component transporting angular momentum outwards develop due to MRI, provided the initial field is not too weak. Weak initial magnetic fields ($\la$$10^{11}~{\rm G}$) have no significant effect on the dynamics of the core and the gravitational wave signal. Strong initial fields ($\ga$$10^{12}~{\rm G}$) cause considerable angular momentum transport whereby rotational energy is extracted from the collapsed core which loses centrifugal support and enters a phase of secular contraction. The gravitational wave amplitude at bounce changes by up to a few ten percent compared to the corresponding non-magnetic model. If the angular momentum losses are large, the post-bounce model. If the angular momentum losses are large the post-bounce equilibrium state of the core changes from a centrifugally to a pressure supported one. This transition imprints in the gravitational wave signal a reduction of the amplitude of the large-scale oscillations characteristic of cores bouncing due to centrifugal forces. In some models the quasi-periodic large-scale oscillations are replaced by higher frequency irregular oscillations. This pattern defines a new signal type which we call a type IV gravitational wave signal. Collimated bipolar outflows give rise to a unique feature that may allow their detection by means of gravitational wave astronomy: a large positive quadrupole wave amplitude of similar size as that of the bounce signal.
We have performed 2.5D and 3D simulations of conical jets driven by the rotation of an ordered, large-scale magnetic field in a stratified atmosphere. The simulations cover about three orders of ...magnitude in distance to capture the centrifugal acceleration as well as the evolution past the Alfvénsurface. We find that the jets develop kink instabilities, the characteristics of which depend on the velocity profile imposed at the base of the flow. The instabilities are especially pronounced with a rigid rotation profile, which induces a shearless magnetic field. The jet's expansion appears to be limiting the growth of Alfvénmode instabilities.
We study the evolution of the field on the surface of proto-neutron stars in the immediate aftermath of stellar core collapse by analyzing the results of self-consistent, axisymmetric simulations of ...the cores of rapidly rotating high-mass stars. To this end, we compare the field topology and the angular spectra of the poloidal and toroidal field components over a time of about one seconds for cores. Both components are characterized by a complex geometry with high power at intermediate angular scales. The structure is mostly the result of the accretion of magnetic flux embedded in the matter falling through the turbulent post-shock layer onto the PNS. Our results may help to guide further studies of the long-term magneto-thermal evolution of proto-neutron stars. We find that the accretion of stellar progenitor layers endowed with low or null magnetization bury the magnetic field on the PNS surface very effectively.
We continue our investigations of the magnetorotational collapse of stellar cores by discussing simulations performed with a modified Newtonian gravitational potential that mimics general ...relativistic effects. The approximate TOV gravitational potential used in our simulations captures several basic features of fully relativistic simulations quite well. In particular, it is able to correctly reproduce the behavior of models that show a qualitative change both of the dynamics and the gravitational wave signal when switching from Newtonian to fully relativistic simulations. For models where the dynamics and gravitational wave signals are already captured qualitatively correctly by a Newtonian potential, the results of the Newtonian and the approximate TOV models differ quantitatively. The collapse proceeds to higher densities with the approximate TOV potential, allowing for a more efficient amplification of the magnetic field by differential rotation. The strength of the saturation fields (${\sim} 10^{15} ~ \mathrm{G}$ at the surface of the inner core) is a factor of two to three higher than in Newtonian gravity. Due to the more efficient field amplification, the influence of magnetic fields is considerably more pronounced than in the Newtonian case for some of the models. As in the Newtonian case, sufficiently strong magnetic fields slow down the core's rotation and trigger a secular contraction phase to higher densities. More clearly than in Newtonian models, the collapsed cores of these models exhibit two different kinds of shock generation. Due to magnetic braking, a first shock wave created during the initial centrifugal bounce at subnuclear densities does not suffice for ejecting any mass, and the temporarily stabilized core continues to collapse to supranuclear densities. Another stronger shock wave is generated during the second bounce as the core exceeds nuclear matter density. The gravitational wave signal of these models does not fit into the standard classification. Therefore, in the first paper of this series we introduced a new type of gravitational wave signal, which we call type IV or “magnetic type”. This signal type is more frequent for the approximate relativistic potential than for the Newtonian one. Most of our weak-field models are marginally detectable with the current LIGO interferometer for a source located at a distance of 10 kpc. Strongly magnetized models emit a substantial fraction of their GW power at very low frequencies. A flat spectrum between 10 Hz and $\la $100 kHz denotes the generation of a jet-like hydromagnetic outflow.