2D simulations of the double-detonation model for thermonuclear transients from low-mass carbon-oxygen white dwarfs

Thermonuclear explosions may arise in binary star systems in which a carbon-oxygen (CO) white dwarf (WD) accretes helium-rich material from a companion star. If the accretion rate allows a sufficiently large mass of helium to accumulate prior to ignition of nuclear burning, the helium surface layer may detonate, giving rise to an astrophysical transient. Detonation of the accreted helium layer generates shock waves that propagate into the underlying CO WD. This might directly ignite a detonation of the CO WD at its surface (an edge-lit secondary detonation) or compress the core of the WD sufficiently to trigger a CO detonation near the centre. If either of these ignition mechanisms works, the two detonations (helium and CO) can then release sufficient energy to completely unbind the WD. These 'double-detonation' scenarios for thermonuclear explosion of WDs have previously been investigated as a potential channel for the production of Type Ia supernovae from WDs of ~ 1 M . Here we extend our 2D studies of the double-detonation model to significantly less massive CO WDs, the explosion of which could produce fainter, more rapidly evolving transients. We investigate the feasibility of triggering a secondary core detonation by shock convergence in low-mass CO WDs and the observable consequences of such a detonation. Our results suggest that core detonation is probable, even for the lowest CO core masses that are likely to be realized in nature. To quantify the observable signatures of core detonation, we compute spectra and light curves for models in which either an edge-lit or compression-triggered CO detonation is assumed to occur. We compare these to synthetic observables for models in which no CO detonation was allowed to occur. If significant shock compression of the CO WD occurs prior to detonation, explosion of the CO WD can produce a sufficiently large mass of radioactive iron-group nuclei to significantly affect the light curves. In particular, this can lead to relatively slow post-maximum decline. If the secondary detonation is edge-lit, however, the CO WD explosion primarily yields intermediate-mass elements that affect the observables more subtly. In this case, near-infrared observations and detailed spectroscopic analysis would be needed to determine whether a core detonation occurred. We comment on the implications of our results for understanding peculiar astrophysical transients including SN 2002bj, SN 2010X and SN 2005E. © 2012 The Authors Monthly Notices of the Royal Astronomical Society © 2012 RAS.

Cosmic-ray acceleration and escape from supernova remnants

Galactic cosmic-ray (CR) acceleration to the knee in the spectrum at a few PeV is only possible if the magnetic field ahead of a supernova remnant (SNR) shock is strongly amplified by CRs escaping the SNR. A model formulated in terms of the electric charge carried by escaping CRs predicts the maximum CR energy and the energy spectrum of CRs released into the surrounding medium. We find that historical SNRs such as Cas A, Tycho and Kepler may be expanding too slowly to accelerate CRs to the knee at the present time.

Cosmic ray acceleration at oblique shocks

We show that the diffusion approximation breaks down for particle acceleration at oblique shocks with velocities typical of young supernova remnants. Higher order anisotropies flatten the spectral index at quasi-parallel shocks and steepen the spectral index at quasi-perpendicular shocks. We compare the theory with observed spectral indices.

TWO-DIMENSIONAL PARTICLE-IN-CELL SIMULATIONS OF THE NONRESONANT, COSMIC-RAY-DRIVEN INSTABILITY IN SUPERNOVA REMNANT SHOCKS

In supernova remnants, the nonlinear amplification of magnetic fields upstream of collisionless shocks is essential for the acceleration of cosmic rays to the energy of the "knee" at 10(15.5) eV. A nonresonant instability driven by the cosmic ray current is thought to be responsible for this effect. We perform two-dimensional, particle-in-cell simulations of this instability. We observe an initial growth of circularly polarized nonpropagating magnetic waves as predicted in linear theory. It is demonstrated that in some cases the magnetic energy density in the growing waves can grow to at least 10 times its initial value. We find no evidence of competing modes, nor of significant modification by thermal effects. At late times, we observe saturation of the instability in the simulation, but the mechanism responsible is an artifact of the periodic boundary conditions and has no counterpart in the supernova-shock scenario.

Tracking magnetic bright point motions through the solar atmosphere

High-cadence, multiwavelength observations and simulations are employed for the analysis of solar photospheric magnetic bright points (MBPs) in the quiet Sun. The observations were obtained with the Rapid Oscillations in the Solar Atmosphere (ROSA) imager and the Interferometric Bidimensional Spectrometer at the Dunn Solar Telescope. Our analysis reveals that photospheric MBPs have an average transverse velocity of approximately 1 km s-1, whereas their chromospheric counterparts have a slightly higher average velocity of 1.4 km s-1. Additionally, chromospheric MBPs were found to be around 63 per cent larger than the equivalent photospheric MBPs. These velocity values were compared with the output of numerical simulations generated using the muram code. The simulated results were similar, but slightly elevated, when compared to the observed data. An average velocity of 1.3 km s-1 was found in the simulated G-band images and an average of 1.8 km s-1 seen in the velocity domain at a height of 500 km above the continuum formation layer. Delays in the change of velocities were also analysed. Average delays of ˜4 s between layers of the simulated data set were established and values of ˜29 s observed between G-band and Ca ii K ROSA observations. The delays in the simulations are likely to be the result of oblique granular shock waves, whereas those found in the observations are possibly the result of a semi-rigid flux tube.

Dust-acoustic shocks in strongly coupled dusty plasmas

Electrostatic dust-acoustic shock waves are investigated in a viscous, complex plasma consisting of dust particles, electrons, and ions. The system is modelled using the generalized hydrodynamic equations, with strong coupling between the dust particles being accounted for by employing the effective electrostatic temperature approach. Using a reductive perturbation method, it is demonstrated that this model predicts the existence of weakly nonlinear dust-acoustic shock waves, arising as solutions to Burgers's equation, in which the nonlinear forces are balanced by dissipative forces, in this case, associated with viscosity. The evolution and stability of dust-acoustic shocks is investigated via a series of numerical simulations, which confirms our analytical predictions on the shock characteristics.

Developed turbulence and nonlinear amplification of magnetic fields in laboratory and astrophysical plasmas

The visible matter in the universe is turbulent and magnetized. Turbulence in galaxy clusters is produced by mergers and by jets of the central galaxies and believed responsible for the amplification of magnetic fields. We report on experiments looking at the collision of two laser-produced plasma clouds, mimicking, in the laboratory, a cluster merger event. By measuring the spectrum of the density fluctuations, we infer developed, Kolmogorov-like turbulence. From spectral line broadening, we estimate a level of turbulence consistent with turbulent heating balancing radiative cooling, as it likely does in galaxy clusters. We show that the magnetic field is amplified by turbulent motions, reaching a nonlinear regime that is a precursor to turbulent dynamo. Thus, our experiment provides a promising platform for understanding the structure of turbulence and the amplification of magnetic fields in the universe.

Model reduction for nonlinear aerodynamics and aeroelasticity using a Discrete Empirical Interpolation Method

A novel surrogate model is proposed in lieu of Computational Fluid Dynamics (CFD) solvers, for fast nonlinear aerodynamic and aeroelastic modeling. A nonlinear function is identified on selected interpolation points by
a discrete empirical interpolation method (DEIM). The flow field is then reconstructed using a least square approximation of the flow modes extracted
by proper orthogonal decomposition (POD). The aeroelastic reduce order
model (ROM) is completed by introducing a nonlinear mapping function
between displacements and the DEIM points. The proposed model is investigated to predict the aerodynamic forces due to forced motions using
a N ACA 0012 airfoil undergoing a prescribed pitching oscillation. To investigate aeroelastic problems at transonic conditions, a pitch/plunge airfoil
and a cropped delta wing aeroelastic models are built using linear structural models. The presence of shock-waves triggers the appearance of limit
cycle oscillations (LCO), which the model is able to predict. For all cases
tested, the new ROM shows the ability to replicate the nonlinear aerodynamic forces, structural displacements and reconstruct the complete flow
field with sufficient accuracy at a fraction of the cost of full order CFD
model.

Nonlinear theory of solitary waves associated with longitudinal particle motion in lattices - Application to longitudinal grain oscillations in a dust crystal

The nonlinear aspects of longitudinal motion of interacting point masses in a lattice are revisited, with emphasis on the paradigm of charged dust grains in a dusty plasma (DP) crystal. Different types of localized excitations, predicted by nonlinear wave theories, are reviewed and conditions for their occurrence (and characteristics) in DP crystals are discussed. Making use of a general formulation, allowing for an arbitrary (e.g. the Debye electrostatic or else) analytic potential form phi(r) and arbitrarily long site-to-site range of interactions, it is shown that dust-crystals support nonlinear kink-shaped localized excitations propagating at velocities above the characteristic DP lattice sound speed v(0). Both compressive and rarefactive kink-type excitations are predicted, depending on the physical parameter values, which represent pulse- (shock-)like coherent structures for the dust grain relative displacement. Furthermore, the existence of breather-type localized oscillations, envelope-modulated wavepackets and shocks is established. The relation to previous results on atomic chains as well as to experimental results on strongly-coupled dust layers in gas discharge plasmas is discussed.

Electron-scale electrostatic solitary waves and shocks: the role of superthermal electrons

The propagation of electron-acoustic solitary waves and shock structures is investigated in a plasma characterized by a superthermal electron population. A three-component plasma model configuration is employed, consisting of inertial (“cold”) electrons, inertialess ? (kappa) distributed superthermal (“hot”) electrons and stationary ions. A multiscale method is employed, leading to a Korteweg-de Vries (KdV) equation for the electrostatic potential (in the absence of dissipation). Taking into account dissipation, a hybrid Korteweg-de Vries-Burgers (KdVB) equation is derived. Exact negative-potential pulse- and kink-shaped solutions (shocks) are obtained. The relative strength among dispersion, nonlinearity and damping coefficients is discussed. Excitations formed in superthermal plasma (finite ?) are narrower and steeper, compared to the Maxwellian case (infinite ?). A series of numerical simulations confirms that energy initially stored in a solitary pulse which propagates in a stable manner for large ? (Maxwellian plasma) may break down to smaller structures or/and to random oscillations, when it encounters a small-? (nonthermal) region. On the other hand, shock structures used as initial conditions for numerical simulations were shown to be robust, essentially responding to changed in the environment by a simple profile change (in width).

Universal behaviour of shock precursors in the presence of efficient cosmic ray acceleration

The self-consistent interaction between energetic particles and self-generated hydromagnetic waves in a cosmic ray pressure dominated plasma is considered. Using a three-dimensional hybrid magnetohydrodynamics (MHD)-kinetic code, which utilizes a spherical harmonic expansion of the Vlasov-Fokker-Planck equation, high-resolution simulations of the magnetic field growth including feedback on the cosmic rays are carried out. It is found that for shocks with high cosmic ray acceleration efficiency, the magnetic fields become highly disorganized, resulting in near isotropic diffusion, independent of the initial orientation of the ambient magnetic field. The possibility of sub-Bohm diffusion is demonstrated for parallel shocks, while the diffusion coefficient approaches the Bohm limit from below for oblique shocks. This universal behaviour suggests that Bohm diffusion in the root-mean-squared field inferred from observation may provide a realistic estimate for the maximum energy acceleration time-scale in young supernova remnants. Although disordered, the magnetic field is not self-similar suggesting a non-uniform energy-dependent behaviour of the energetic particle transport in the precursor. Possible indirect radiative signatures of cosmic ray driven magnetic field amplification are discussed.

Diffusive shock acceleration at laser-driven shocks: studying cosmic-ray accelerators in the laboratory

The non-thermal particle spectra responsible for the emission from many astrophysical systems are thought to originate from shocks via a first order Fermi process otherwise known as diffusive shock acceleration. The same mechanism is also widely believed to be responsible for the production of high energy cosmic rays. With the growing interest in collisionless shock physics in laser produced plasmas, the possibility of reproducing and detecting shock acceleration in controlled laboratory experiments should be considered. The various experimental constraints that must be satisfied are reviewed. It is demonstrated that several currently operating laser facilities may fulfil the necessary criteria to confirm the occurrence of diffusive shock acceleration of electrons at laser produced shocks. Successful reproduction of Fermi acceleration in the laboratory could open a range of possibilities, providing insight into the complex plasma processes that occur near astrophysical sources of cosmic rays.

Modification of the formation of high-Mach number electrostatic shock-like structures by the ion acoustic instability

The formation of unmagnetized electrostatic shock-like structures with a high Mach number is examined with one- and two-dimensional particle-in-cell (PIC) simulations. The structures are generated through the collision of two identical plasma clouds, which consist of equally hot electrons and ions with a mass ratio of 250. The Mach number of the collision speed with respect to the initial ion acoustic speed of the plasma is set to 4.6. This high Mach number delays the formation of such structures by tens of inverse ion plasma frequencies. A pair of stable shock-like structures is observed after this time in the 1D simulation, which gradually evolve into electrostatic shocks. The ion acoustic instability, which can develop in the 2D simulation but not in the 1D one, competes with the nonlinear process that gives rise to these structures. The oblique ion acoustic waves fragment their electric field. The transition layer, across which the bulk of the ions change their speed, widens and their speed change is reduced. Double layer-shock hybrid structures develop.

The dynamics of neutrino-driven supernova explosions after shock revival in 2D and 3D

We study the growth of the explosion energy after shock revival in neutrino-driven explosions in two and three dimensions (2D/3D) using multi-group neutrino hydrodynamics simulations of an 11.2 M⊙ star. The 3D model shows a faster and steadier growth of the explosion energy and already shows signs of subsiding accretion after one second. By contrast, the growth of the explosion energy in 2D is unsteady, and accretion lasts for several seconds as confirmed by additional long-time simulations of stars of similar masses. Appreciable explosion energies can still be reached, albeit at the expense of rather high neutron star masses. In 2D, the binding energy at the gain radius is larger because the strong excitation of downward-propagating g modes removes energy from the freshly accreted material in the downflows. Consequently, the mass outflow rate is considerably lower in 2D than in 3D. This is only partially compensated by additional heating by outward-propagating acoustic waves in 2D. Moreover, the mass outflow rate in 2D is reduced because much of the neutrino energy deposition occurs in downflows or bubbles confined by secondary shocks without driving outflows. Episodic constriction of outflows and vertical mixing of colder shocked material and hot, neutrino-heated ejecta due to Rayleigh–Taylor instability further hamper the growth of the explosion energy in 2D. Further simulations will be necessary to determine whether these effects are generic over a wider range of supernova progenitors.