Observation of Collisionless Shocks in Laser-Plasma Experiments

The propagation in a rarefied plasma (n(e)less than or similar to 10(15) cm(-3)) of collisionless shock waves and ion-acoustic solitons, excited following the interaction of a long (tau(L)similar to 470 ps) and intense (I similar to 10(15) W cm(-2)) laser pulse with solid targets, has been investigated via proton probing techniques. The shocks' structures and related electric field distributions were reconstructed with high spatial and temporal resolution. The experimental results were interpreted within the framework of the nonlinear wave description based on the Korteweg-de Vries-Burgers equation.

Evolution of elastic x-ray scattering in laser-shocked warm dense lithium

We have studied the dynamics of warm dense Li with near-elastic x-ray scattering. Li foils were heated and compressed using shock waves driven by 4-ns-long laser pulses. Separate 1-ns-long laser pulses were used to generate a bright source of 2.96 keV Cl Ly-alpha photons for x-ray scattering, and the spectrum of scattered photons was recorded at a scattering angle of 120 degrees using a highly oriented pyrolytic graphite crystal operated in the von Hamos geometry. A variable delay between the heater and backlighter laser beams measured the scattering time evolution. Comparison with radiation-hydrodynamics simulations shows that the plasma is highly coupled during the first several nanoseconds, then relaxes to a moderate coupling state at later times. Near-elastic scattering amplitudes have been successfully simulated using the screened one-component plasma model. Our main finding is that the near-elastic scattering amplitudes are quite sensitive to the mean ionization state Z and by extension to the choice of ionization model in the radiation-hydrodynamics simulations used to predict plasma properties within the shocked Li.

Oscillatory Behavior of Supersonic Impinging Jet Flows

Oscillatory flows of a choked underexpanded supersonic impinging jet issuing from a convergent nozzle have been computed using the axisymmetric unsteady Navier–Stokes system. This paper focuses on the oscillatory flow features associated with the variation of the nozzle-to-plate distance and nozzle pressure ratio. Frequencies of the surface pressure oscillation and flow structural changes from computational results have been analyzed. Staging behaviour of the oscillation frequency has been observed for both cases of nozzle-to-plate distance variation and pressure ratio variation. However, the staging behavior for each case exhibits different features. These two distinct staging behaviors of the oscillation frequency are found to correlate well if the frequency and the distance are normalized by the length of the shock cell. It is further found that the staging behaviour is strongly correlated with the change of the pressure wave pattern in the jet shear layer, but not with the shock cell structure.

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.

A two-armed pattern in flickering maps of the nova-like variable UU aquarii

We report the analysis of a uniform sample of 31 light curves of the nova-like variable UU Aqr with eclipse-mapping techniques. The data were combined to derive eclipse maps of the average steady-light component, the long-term brightness changes, and the low- and high-frequency flickering components. The long-term variability responsible for the ""low-brightness`` and ""high-brightness`` states is explained in terms of the response of a viscous disk to changes of 20%-50% in the mass transfer rate from the donor star. Low- and high-frequency flickering maps are dominated by emission from two asymmetric arcs reminiscent of those seen in the outbursting dwarf nova IP Peg, and they are similarly interpreted as manifestations of a tidally induced spiral shock wave in the outer regions of a large accretion disk. The asymmetric arcs are also seen in the map of the steady light aside from the broad brightness distribution of a roughly steady-state disk. The arcs account for 25% of the steady-light flux and are a long-lasting feature in the accretion disk of UU Aqr. We infer an opening angle of 10 degrees +/- 3 degrees for the spiral arcs. The results suggest that the flickering in UU Aqr is caused by turbulence generated after the collision of disk gas with the density-enhanced spiral wave in the accretion disk.

Revisiting 2D numerical models for the 19th century outbursts of eta Carinae

We present here new results of two-dimensional hydrodynamical simulations of the eruptive events of the 1840s (the great) and the 1890s (the minor) eruptions suffered by the massive star eta Carinae (Car). The two bipolar nebulae commonly known as the Homunculus and the little Homunculus (LH) were formed from the interaction of these eruptive events with the underlying stellar wind. We assume here an interacting, non-spherical multiple-phase wind scenario to explain the shape and the kinematics of both Homunculi, but adopt a more realistic parametrization of the phases of the wind. During the 1890s eruptive event, the outflow speed decreased for a short period of time. This fact suggests that the LH is formed when the eruption ends, from the impact of the post-outburst eta Car wind (that follows the 1890s event) with the eruptive flow (rather than by the collision of the eruptive flow with the pre-outburst wind, as claimed in previous models; Gonzalez et al.). Our simulations reproduce quite well the shape and the observed expansion speed of the large Homunculus. The LH (which is embedded within the large Homunculus) becomes Rayleigh-Taylor unstable and develop filamentary structures that resemble the spatial features observed in the polar caps. In addition, we find that the interior cavity between the two Homunculi is partially filled by material that is expelled during the decades following the great eruption. This result may be connected with the observed double-shell structure in the polar lobes of the eta Car nebula. Finally, as in previous work, we find the formation of tenuous, equatorial, high-speed features that seem to be related to the observed equatorial skirt of eta Car.

Local star formation triggered by supernova shocks in magnetized diffuse neutral clouds

In this work, considering the impact of a supernova remnant (SNR) with a neutral magnetized cloud we derived analytically a set of conditions that are favourable for driving gravitational instability in the cloud and thus star formation. Using these conditions, we have built diagrams of the SNR radius, R(SNR), versus the initial cloud density, n(c), that constrain a domain in the parameter space where star formation is allowed. This work is an extension to previous study performed without considering magnetic fields (Melioli et al. 2006, hereafter Paper I). The diagrams are also tested with fully three-dimensional MHD radiative cooling simulations involving a SNR and a self-gravitating cloud and we find that the numerical analysis is consistent with the results predicted by the diagrams. While the inclusion of a homogeneous magnetic field approximately perpendicular to the impact velocity of the SNR with an intensity similar to 1 mu G within the cloud results only a small shrinking of the star formation zone in the diagram relative to that without magnetic field, a larger magnetic field (similar to 10 mu G) causes a significant shrinking, as expected. Though derived from simple analytical considerations these diagrams provide a useful tool for identifying sites where star formation could be triggered by the impact of a supernova blast wave. Applications of them to a few regions of our own Galaxy (e.g. the large CO shell in the direction of Cassiopeia, and the Edge Cloud 2 in the direction of the Scorpious constellation) have revealed that star formation in those sites could have been triggered by shock waves from SNRs for specific values of the initial neutral cloud density and the SNR radius. Finally, we have evaluated the effective star formation efficiency for this sort of interaction and found that it is generally smaller than the observed values in our own Galaxy (SFE similar to 0.01-0.3). This result is consistent with previous work in the literature and also suggests that the mechanism presently investigated, though very powerful to drive structure formation, supersonic turbulence and eventually, local star formation, does not seem to be sufficient to drive global star formation in normal star-forming galaxies, not even when the magnetic field in the neutral clouds is neglected.

SPH simulations of clumps formation by dissipative collision of molecular clouds - I. Non magnetic case

Computer experiments of interstellar cloud collisions were performed with a new smoothed-particle-hydrodynamics (SPH) code. The SPH quantities were calculated by using spatially adaptive smoothing lengths and the SPH fluid equations of motion were solved by means of a hierarchical multiple time-scale leapfrog. Such a combination of methods allows the code to deal with a large range of hydrodynamic quantities. A careful treatment of gas cooling by H, H(2), CO and H II, as well as a heating mechanism by cosmic rays and by H(2) production on grains surface, were also included in the code. The gas model reproduces approximately the typical environment of dark molecular clouds. The experiments were performed by impinging two dynamically identical spherical clouds onto each other with a relative velocity of 10 km s(-1) but with a different impact parameter for each case. Each object has an initial density profile obeying an r(-1)-law with a cutoff radius of 10 pc and with an initial temperature of 20 K. As a main result, cloud-cloud collision triggers fragmentation but in expense of a large amount of energy dissipated, which occurred in the head-on case only. Off-center collision did not allow remnants to fragment along the considered time (similar to 6 Myr). However, it dissipated a considerable amount of orbital energy. Structures as small as 0.1 pc, with densities of similar to 10(4) cm(-3), were observed in the more energetic collision.

A Novel Plasma Technique for Surface Treatment: The Plasma Expander

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)