Search

The interaction of relativistic short-pulse lasers with matter produces fast electrons with over megaampere currents, which supposedly heats a solid target isochorically and forms a hot dense plasma. In a picosecond timescale, however, thermal diffusion from hot preformed plasma turns out to be the dominant process of isochoric heating. We describe a heating process, fast thermal diffusion, launched from the preformed plasma heated resistively by the fast electron current. We demonstrate the fast thermal diffusion in the keV range in a solid density plasma by a series of one-dimensional particle-in-cell simulations. A theoretical model of the fast thermal diffusion is developed and we derive the diffusion speed as a function of the laser amplitude and target density. Under continuous laser irradiation, the diffusion front propagates at a constant speed in uniform plasma. Our model can provide a guideline for fast isochoric heating using future kilojoule petawatt lasers.

International audience; Carbonate minerals, for example calcite and magnesite, exist on the planetary surfaces of the Earth, Mars, and Venus, and are subjected to hypervelocity collisions. The physical properties of planetary materials at extreme conditions are essential for understanding their dynamic behaviors at hypervelocity collisions and the mantle structure of rocky planets including Super-Earths. Here we report laboratory investigations of laser-shocked calcite at pressures of 200-960 GPa (impact velocities of 12-30 km/s and faster than escape velocity from the Earth) using decay shock techniques. Our measured temperatures above 200 GPa indicated a large difference from the previous theoretical models. The present shock Hugoniot data and temperature measurements, compared with the previous reports, indicate melting without decomposition at pressures of ~110 GPa to ~350 GPa and a bonded liquid up to 960 GPa from the calculated specific heat. Our temperature calculations of calcite at 1 atm adiabatically released from the Hugoniot points suggest that the released products vary depending on the shock pressures and affect the planetary atmosphere by the degassed species. The present results on calcite newly provide an important anchor for considering the theoretical EOS at the extreme conditions, where the model calculations show a significant diversity at present.

Laser experiments are becoming established as a new tool for astronomical research that complements observations and theoretical modeling. Localized strong magnetic fields have been observed at a shock front of supernova explosions. Experimental confirmation and identification of the physical mechanism for this observation are of great importance in understanding the evolution of the interstellar medium. However, it has been challenging to treat the interaction between hydrodynamic instabilities and an ambient magnetic field in the laboratory. Here, we developed an experimental platform to examine magnetized Richtmyer-Meshkov instability (RMI). The measured growth velocity was consistent with the linear theory, and the magnetic-field amplification was correlated with RMI growth. Our experiment validated the turbulent amplification of magnetic fields associated with the shock-induced interfacial instability in astrophysical conditions for the first time. Experimental elucidation of fundamental processes in magnetized plasmas is generally essential in various situations such as fusion plasmas and planetary sciences., 17 pages, 12 figures, 3 tables, accepted for publication in PRE

A counter-propagating laser-beam platform using a spherical plasma mirror was developed for the kilojoule-class petawatt LFEX laser. The temporal and spatial overlaps of the incoming and redirected beams were measured with an optical interferometer and an x-ray pinhole camera. The plasma mirror performance was evaluated by measuring fast electrons, ions, and neutrons generated in the counter-propagating laser interaction with a Cu-doped deuterated film on both sides. The reflectivity and peak intensity were estimated as ∼50% and ∼5 × 1018 W/cm2, respectively. The platform could enable studies of counter-streaming charged particles in high-energy-density plasmas for fundamental and inertial confinement fusion research.

Ablative Rayleigh-Taylor instability growth was investigated to elucidate the fundamental physics of thermal conduction suppression in a magnetic field. Experiments found that unstable modulation growth is faster in an external magnetic field. This result was reproduced by a magnetohydrodynamic simulation based on a Braginskii model of electron thermal transport. An external magnetic field reduces the electron thermal conduction across the magnetic field lines because the Larmor radius of the thermal electrons in the field is much shorter than the temperature scale length. Thermal conduction suppression leads to spatially nonuniform pressure and reduced thermal ablative stabilization, which in turn increases the growth of ablative Rayleigh-Taylor instability.

Nonlinear motion of vortex sheets with non-uniform current is investigated taking the magnetohydrodynamic Richtmyer–Meshkov instability (MHD RMI) and the magnetohydrodynamic Kelvin–Helmholtz instability (MHD KHI) as the examples. As the ratio of the magnetic force to the convective force increases, Alfven oscillations appear and the nonlinear growth of the interface as a vortex sheet is suppressed. We show that the turbulent energy possessing the interface flows into the magnetic energy, which causes the strong magnetic field amplification for both instabilities. We also discuss the difference of the temporal evolution of the interface between MHD RMI and MHD KHI.

Intense laser-plasma interactions are an essential tool for the laboratory study of ion acceleration at a collisionless shock. With two-dimensional particle-in-cell calculations of a multicomponent plasma we observe two electrostatic collisionless shocks at two distinct longitudinal positions when driven with a linearly polarized laser at normalized laser vector potential ${a}_{0}$ that exceeds 10. Moreover, these shocks, associated with protons and carbon ions, show a power-law dependence on ${a}_{0}$ and accelerate ions to different velocities in an expanding upstream with higher flux than in a single-component hydrogen or carbon plasma. This results from an electrostatic ion two-stream instability caused by differences in the charge-to-mass ratio of different ions. Particle acceleration in collisionless shocks in multicomponent plasma are ubiquitous in space and astrophysics, and these calculations identify the possibility for studying these complex processes in the laboratory.

Magnetohydrodynamical evolution of the Richtmyer-Meshkov instability (RMI) is investigated by two-dimensional MHD simulations. The RMI is suppressed by a strong magnetic field, whereas the RMI amplifies an ambient magnetic field by many orders of magnitude if the seed field is weak. We have found that the suppression and amplification processes can be evaluated continuously along with the amplitude of the Alfv\'en number $R_A$, which is defined as the ratio of the linear growth velocity of the RMI to the Alfv\'en speed at the interface. When the Alfv\'en number is less than unity, the Lorentz force acting on the fluid mitigates the unstable motion of the RMI significantly, and the interface oscillates stably in this limit. If $R_A \gtrsim 1$, on the other hand, the surface modulation increases due to the growth of the RMI. The maximum strength of the magnetic field is enhanced up to by a factor of $R_A$. This critical feature is universal and independent of the initial Mach number of the incident shock, the Atwood number, corrugation amplitude, and even the direction of the initial magnetic field., Comment: 11 pages, 8 figures, 1 table, accepted for publication in ApJ

We have developed a public code, Optab, that outputs Rosseland, Planck, and two-temperature Planck mean gas opacity tables for radiation hydrodynamics simulations in astrophysics. The code is developed for modern high-performance computing, being written in Fortran 90 and using Message Passing Interface and Hierarchical Data Format, Version 5. The purpose of this work is to provide a platform on which users can generate opacity tables for their own research purposes. Therefore, the code has been designed so that a user can easily modify, change, or add opacity sources in addition to those already implemented, which include bremsstrahlung, photoionization, Rayleigh scattering, line absorption, and collision-induced absorption. In this paper, we provide details of the opacity calculations in our code and present validation tests to evaluate the performance of our code., Astronomy and Astrophysics, in press

Calcite (CaCO3) as a planetary material is a source to the atmospheric carbon dioxide through degassing by high-velocity impact events. Revealing the behavior of calcite in the extreme pressure and temperature conditions is required to understand the impact-induced degassing phenomena. Here we report laboratory investigations of shock- compressed calcite beyond the impact velocity of 12 km/s (faster than escape velocity from the Earth). The present precise shock measurements elucidate the shape of the calcite Hugoniot curve continuously passing through the melting and metallization states up to a pressure of 1000 GPa (= 10-million atmospheres) or a corresponding impact velocity of 30 km/s, allowing us to predict the post-shock residual temperatures and the dominant carbon oxide species in the impact aftermath. These predictions suggest that CO emission is much more dominant than CO2 at the impact velocities of ∼10 km/s and above, affecting the planetary atmospheric chemistry, greenhouse processes, and environmental changes during planetary evolution.

Propagation features of circularly polarized (CP) electromagnetic waves in magnetized plasmas are determined by the plasma density and the magnetic field strength. This property can be applied to design a unique plasma photonic device for intense short-pulse lasers. We have demonstrated by numerical simulations that a thin plasma foil under an external magnetic field works as a polarizing plate to separate a linearly polarized laser into two CP waves traveling in the opposite direction. This plasma photonic device has an advantage for generating intense CP waves even with a relativistic amplitude. For various research purposes, intense CP lights are strongly required to create high energy density plasmas in the laboratory., 9 pages, 8 figures, accepted for publication in PRE

We have investigated the effects of a smooth transition layer at the contact discontinuity on the growth of the Richtmyer-Meshkov instability (RMI) by hydrodynamic numerical simulations and derived an empirical condition for the suppression of the instability. The transition layer has little influence on the RMI when the thickness $L$ is narrower than the wavelength of an interface modulation $\lambda$. However, if the transition layer becomes broader than $\lambda$, the perturbed velocity associated with the RMI is reduced considerably. The suppression condition is interpreted as the cases that the shock transit time through the transition layer is longer than the sound crossing time of the modulation wavelength. The fluctuation kinetic energy decreases as $L^{-p}$ with $p = 2.5$, which indicates that the growth velocity of the RMI decreases in proportion to $L^{-p/2}$ by the presence of the transition layer. This feature is found to be quite universal and appeared in a wide range of shock-interface interactions., Comment: 11 pages, 8 figures, 1 table, accepted for publication in PRE

Thermal fusion plasmas initiated by standing whistler waves are investigated numerically by two- and one-dimensional Particle-in-Cell simulations. When a standing whistler wave collapses due to the wave breaking of ion plasma waves, the energy of the electromagnetic waves transfers directly to the ion kinetic energy. Here we find that the ion heating by the standing whistler wave is operational even in multi-dimensional simulations of multi-ion species targets, such as deuterium-tritium (DT) ices and solid ammonia borane (H$_6$BN). The energy conversion efficiency to ions becomes as high as 15% of the injected laser energy, which depends significantly on the target thickness and laser pulse duration. The ion temperature could reach a few tens of keV or much higher if appropriate laser-plasma conditions are selected. DT fusion plasmas generated by this method must be useful as efficient neutron sources. Our numerical simulations suggest that the neutron generation efficiency exceeds 10$^9$ n/J per steradian, which is beyond the current achievements of the state-of-the-art laser experiments. The standing whistler wave heating would expand the experimental possibility for an alternative ignition design of magnetically confined laser fusion, and also for more difficult fusion reactions including the aneutronic proton-boron reaction., Comment: 16 pages, 11 figures, accepted for publication in PRE

International audience; We describe a platform developed on the LULI2000 laser facility to investigate the evolution of Rayleigh-Taylor instability (RTI) in scaled conditions relevant to young Supernovae Remnant (SNR) up to 200 years. An RT unstable interface is imaged with a short-pulse laser-driven (PICO2000) x-ray source, providing an unprecedented simultaneous high spatial (24 µm) and temporal (10 ps) resolution. This experiment provides relevant data to compare with astrophysical codes, as observational data on the development of RTI at the early stage of the SNR expansion are missing. A comparison is also performed with FLASH radiative magneto-hydrodynamic simulations.

International audience

A preliminary numerical experiment is conducted for laboratory experiments on the generation of magnetized collisionless shocks with high-power lasers by using one-dimensional particle-in-cell simulation. The present study deals with the interaction between a moving aluminum plasma and a nitrogen plasma at rest. In the numerical experiment, the nitrogen plasma is unmagnetized or magnetized by a weak external magnetic field. Since the previous study suggested the generation of a spontaneous magnetic field in the piston (aluminum) plasma due to the Biermann battery, the effect of the magnetic field is of interest. Sharp jumps of the electron density and magnetic field are observed around the interface between the two plasmas as long as one of the two plasmas is magnetized, which indicates the formation of tangential electron-magneto-hydro-dynamic discontinuity. When the aluminum plasma is magnetized, strong compression of both the density and the magnetic field takes place in the pure aluminum plasma during the gyration of nitrogen ions in the aluminum plasma region. The formation of a shock downstream is obtained from the shock jump condition. The results suggest that the spontaneous magnetic field in the piston (aluminum) plasma plays an essential role in the formation of a perpendicular collisionless shock., ファイル公開：2020-03-05

Efficient energy transfer from electromagnetic waves to ions has been demanded to control laboratory plasmas for various applications and could be useful to understand the nature of space and astrophysical plasmas. However, there exists a severe unsolved problem that most of the wave energy is converted quickly to electrons, but not to ions. Here, an energy conversion process to ions in overdense plasmas associated with whistler waves is investigated by numerical simulations and theoretical model. Whistler waves propagating along a magnetic field in space and laboratories often form the standing waves by the collision of counter-propagating waves or through the reflection. We find that ions in the standing whistler waves acquire a large amount of energy directly from the waves in a short timescale comparable to the wave oscillation period. Thermalized ion temperature increases in proportion to the square of the wave amplitude and becomes much higher than the electron temperature in a wide range of wave-plasma conditions. This efficient ion-heating mechanism applies to various plasma phenomena in space physics and fusion energy sciences., Comment: 11 pages, 9 figures, 1 table, accepted for publication in PRE

Fast isochoric laser heating is a scheme to heat a matter with relativistic-intensity ($>$ 10$^{18}$ W/cm$^2$) laser pulse or X-ray free electron laser pulse. The fast isochoric laser heating has been studied for creating efficiently ultra-high-energy-density (UHED) state. We demonstrate an fast isochoric heating of an imploded dense plasma using a multi-picosecond kJ-class petawatt laser with an assistance of externally applied kilo-tesla magnetic fields for guiding fast electrons to the dense plasma.The UHED state with 2.2 Peta-Pascal is achieved experimentally with 4.6 kJ of total laser energy that is one order of magnitude lower than the energy used in the conventional implosion scheme. A two-dimensional particle-in-cell simulation reveals that diffusive heating from a laser-plasma interaction zone to the dense plasma plays an essential role to the efficient creation of the UHED state., 8 pages, 4 figures, 1 table

Water, methane, and ammonia are commonly considered to be the key components of the interiors of Uranus and Neptune. Modelling the planets’ internal structure, evolution, and dynamo heavily relies on the properties of the complex mixtures with uncertain exact composition in their deep interiors. Therefore, characterising icy mixtures with varying composition at planetary conditions of several hundred gigapascal and a few thousand Kelvin is crucial to improve our understanding of the ice giants. In this work, pure water, a water-ethanol mixture, and a water-ethanol-ammonia “synthetic planetary mixture” (SPM) have been compressed through laser-driven decaying shocks along their principal Hugoniot curves up to 270, 280, and 260 GPa, respectively. Measured temperatures spanned from 4000 to 25000 K, just above the coldest predicted adiabatic Uranus and Neptune profiles (3000–4000 K) but more similar to those predicted by more recent models including a thermal boundary layer (7000–14000 K). The experiments were performed at the GEKKO XII and LULI2000 laser facilities using standard optical diagnostics (Doppler velocimetry and optical pyrometry) to measure the thermodynamic state and the shock-front reflectivity at two different wavelengths. The results show that water and the mixtures undergo a similar compression path under single shock loading in agreement with Density Functional Theory Molecular Dynamics (DFT-MD) calculations using the Linear Mixing Approximation (LMA). On the contrary, their shock-front reflectivities behave differently by what concerns both the onset pressures and the saturation values, with possible impact on planetary dynamos.

We developed an experimental platform for studying magnetic reconnection in an external magnetic field with simultaneous measurements of plasma imaging, flow velocity, and magnetic-field variation. Here, we investigate the stagnation and acceleration in counter-streaming plasmas generated by high-power laser beams. A plasma flow perpendicular to the initial flow directions is measured with laser Thomson scattering. The flow is, interestingly, accelerated toward the high-density region, which is opposite to the direction of the acceleration by pressure gradients. This acceleration is possibly interpreted by the interaction of two magnetic field loops initially generated by Biermann battery effect, resulting in a magnetic reconnection forming a single field loop and additional acceleration by a magnetic tension force., Comment: 6 pages, 4 figures, Physics of Plasmas, in press

A theoretical model is proposed to describe fully nonlinear dynamics of interfaces in two-dimensional MHD flows based on an idea of non-uniform current-vortex sheet. Application of vortex sheet model to MHD flows has a crucial difficulty because of non-conservative nature of magnetic tension. However, it is shown that when a magnetic field is initially parallel to an interface, the concept of vortex sheet can be extended to MHD flows (current-vortex sheet). Two-dimensional MHD flows are then described only by a one-dimensional Lagrange parameter on the sheet. It is also shown that bulk magnetic field and velocity can be calculated from their values on the sheet. The model is tested by MHD Richtmyer–Meshkov instability with sinusoidal vortex sheet strength. Two-dimensional ideal MHD simulations show that the nonlinear dynamics of a shocked interface with density stratification agrees fairly well with that for its corresponding potential flow. Numerical solutions of the model reproduce properly the results of the ideal MHD simulations, such as the roll-up of spike, exponential growth of magnetic field, and its saturation and oscillation. Nonlinear evolution of the interface is found to be determined by the Alfven and Atwood numbers. Some of their dependence on the sheet dynamics and magnetic field amplification are discussed. It is shown by the model that the magnetic field amplification occurs locally associated with the nonlinear dynamics of the current-vortex sheet. We expect that our model can be applicable to a wide variety of MHD shear flows.

International audience; The first laboratory evidence of a radiative shock (RS) decelerating during its free expansion phase in an optically thick medium is presented. A shock is generated in a multilayer solid target under the irradiation of a high-power laser at the GEKKO XII laser facility. The rear surface of the target is connected to a gas cell filled with Xe. Upon breakout, an RS, characterized by low Boltzmann number Bo MUCH LESS-THAN 1 and Mihalas number R 10, is generated. Experimental results reveal that radiative losses through the radiative precursor cause the shock to lose energy and decelerate. A model is developed that describes the shock propagation as a function of time. The model is in agreement with both numerical simulations and experimental results. These results have tremendous consequences for astrophysical systems, such as SN 1987A, where radiative deceleration may play a role in the formation of the observed hotspots in the circumstellar ring.

Magnetic reconnections play essential roles in space, astrophysical, and laboratory plasmas, where the anti-parallel magnetic field components re-connect and the magnetic energy is converted to the plasma energy as Alfvénic out flows. Although the electron dynamics is considered to be essential, it is highly challenging to observe electron scale reconnections. Here we show the experimental results on an electron scale reconnection driven by the electron dynamics in laser-produced plasmas. We apply a weak-external magnetic field in the direction perpendicular to the plasma propagation, where the magnetic field is directly coupled with only the electrons but not for the ions. Since the kinetic pressure of plasma is much larger than the magnetic pressure, the magnetic field is distorted and locally anti-parallel. We observe plasma collimations, cusp and plasmoid like features with optical diagnostics. The plasmoid propagates at the electron Alfvén velocity, indicating a reconnection driven by the electron dynamics., Magnetic reconnection is the process of releasing energy by magnetized and space plasma. Here the authors report experimental observation of magnetic reconnection in laser-produced plasma and the role of electron scaling on reconnection.

The physics of compressible turbulence in high energy density (HED) plasmas is an unchartered experimental area. Simulations of compressible and radiative flows relevant for astrophysics rely mainly on subscale parameters. Therefore, we plan to perform turbulent hydrodynamics experiments in HED plasmas (TurboHEDP) in order to improve our understanding of such important phenomena for interest in both communities: laser plasma physics and astrophysics. We will focus on the physics of supernovae remnants which are complex structures subject to fluid instabilities such as the Rayleigh–Taylor and Kelvin–Helmholtz instabilities. The advent of megajoule laser facilities, like the National Ignition Facility and the Laser Megajoule, creates novel opportunities in laboratory astrophysics, as it provides unique platforms to study turbulent mixing flows in HED plasmas. Indeed, the physics requires accelerating targets over larger distances and longer time periods than previously achieved. In a preparatory phase, scaling from experiments at lower laser energies is used to guarantee the performance of future MJ experiments. This subscale experiments allow us to develop experimental skills and numerical tools in this new field of research, and are stepping stones to achieve our objectives on larger laser facilities. We review first in this paper recent advances in high energy density experiments devoted to laboratory astrophysics. Then we describe the necessary steps forward to commission an experimental platform devoted to turbulent hydrodynamics on a megajoule laser facility. Recent novel experimental results acquired on LULI2000, as well as supporting radiative hydrodynamics simulations, are presented. Together with the development of LiF detectors as transformative X-ray diagnostics, these preliminary results are promising on the way to achieve micrometric spatial resolution in turbulent HED physics experiments in the near future.

The interaction of dense plasmas with an intense laser under a strong external magnetic field has been investigated. When the cyclotron frequency for the ambient magnetic field is higher than the laser frequency, the laser's electromagnetic field is converted to the whistler mode that propagates along the field line. Because of the nature of the whistler wave, the laser light penetrates into dense plasmas with no cutoff density, and produces superthermal electrons through cyclotron resonance. It is found that the cyclotron resonance absorption occurs effectively under the broadened conditions, or a wider range of the external field, which is caused by the presence of relativistic electrons accelerated by the laser field. The upper limit of the ambient field for the resonance increases in proportion to the square root of the relativistic laser intensity. The propagation of a large-amplitude whistler wave could raise the possibility for plasma heating and particle acceleration deep inside dense plasmas., Comment: 8 pages, 8 figures, accepted for publication in PRE

International audience; Investigation of the iron phase diagram under high pressure and temperature is crucial for the determination of the composition of the cores of rocky planets and for better understanding the generation of planetary magnetic fields. Here we present X-ray diffraction results from laser-driven shock-compressed single-crystal and polycrystalline iron, indicating the presence of solid hexagonal close-packed iron up to pressure of at least 170 GPa along the principal Hugoniot, corresponding to a temperature of 4,150 K. This is confirmed by the agreement between the pressure obtained from the measurement of the iron volume in the sample and the inferred shock strength from velocimetry deductions. Results presented in this study are of the first importance regarding pure Fe phase diagram probed under dynamic compression and can be applied to study conditions that are relevant to Earth and super-Earth cores.

Shocked forsterite above 250 GPa indicates incongruent crystallization of MgO, its phase transition, and remelting., Forsterite (Mg2SiO4) is one of the major planetary materials, and its behavior under extreme conditions is important to understand the interior structure of large planets, such as super-Earths, and large-scale planetary impact events. Previous shock compression measurements of forsterite indicate that it may melt below 200 GPa, but these measurements did not go beyond 200 GPa. We report the shock response of forsterite above ~250 GPa, obtained using the laser shock wave technique. We simultaneously measured the Hugoniot and temperature of shocked forsterite and interpreted the results to suggest the following: (i) incongruent crystallization of MgO at 271 to 285 GPa, (ii) phase transition of MgO at 285 to 344 GPa, and (iii) remelting above ~470 to 500 GPa. These exothermic and endothermic reactions are seen to occur under extreme conditions of pressure and temperature. They indicate complex structural and chemical changes in the system MgO-SiO2 at extreme pressures and temperatures and will affect the way we understand the interior processes of large rocky planets as well as material transformation by impacts in the formation of planetary systems.

We report the first successful simulation of spontaneous formation of surface magnetic structures from a large-scale dynamo by strongly-stratified thermal convection in Cartesian geometry. The large-scale dynamo observed in our strongly-stratified model has physical properties similar to those in earlier weakly-stratified convective dynamo simulations, indicating that the $\alpha^2$-type mechanism is responsible for it. Additionally to the large-scale dynamo, we find that large-scale structures of the vertical magnetic field are spontaneously formed in the convection zone surface only for the case of strongly-stratified atmosphere. The organization of the vertical magnetic field proceeds in the upper convection zone within tens of convective turn-over time and band-like bipolar structures are recurrently-appeared in the dynamo-saturated stage. We examine possibilities of several candidates as the origin of the surface magnetic structure formation, and then suggest the existence of an as-yet-unknown mechanism for the self-organization of the large-scale magnetic structure, which should be inherent in the strongly-stratified convective atmosphere., Comment: Accepted for Publication in ApJ Letters (7 pages, 5 figure)

We investigate spherical collisionless shocks in the presence of an external magnetic field. Spherical collisionless shocks are common resultant of interactions between a expanding plasma and a surrounding plasma, such as the solar wind, stellar winds, and supernova remnants. Anisotropies often observed in shock propagations and their emissions, and it is widely believed a magnetic field plays a major role. Since the local observations of magnetic fields in astrophysical plasmas are not accessible, laboratory experiments provide unique capability to investigate such phenomena. We model the spherical shocks in the universe by irradiating a solid spherical target surrounded by a plasma in the presence of a magnetic field. We present preliminary results obtained by shadowgraphy.

One of the important and interesting problems in astrophysics and plasma physics is collimation of plasma jets. The collimation mechanism, which causes a plasma flow to propagate a long distance, has not been understood in detail. We have been investigating a model experiment to simulate astrophysical plasma jets with an external magnetic field [Nishio et al., EPJ. Web of Conferences 59, 15005 (2013)]. The experiment was performed by using Gekko XII HIPER laser system at Institute of Laser Engineering, Osaka University. We shot CH plane targets (3 mm × 3 mm × 10 μm) and observed rear-side plasma flows. A collimated plasma flow or plasma jet was generated by separating focal spots of laser beams. In this report, we measured plasma jet structure without an external magnetic field with shadowgraphy, and simultaneously measured the local parameters of the plasma jet, i.e., electron density, electron and ion temperatures, charge state, and drift velocity, with collective Thomson scattering.

A model experiment of magnetic field amplification (MFA) via the Richtmyer-Meshkov instability (RMI) in supernova remnants (SNRs) was performed using a high-power laser. In order to account for very-fast acceleration of cosmic rays observed in SNRs, it is considered that the magnetic field has to be amplified by orders of magnitude from its background level. A possible mechanism for the MFA in SNRs is stretching and mixing of the magnetic field via the RMI when shock waves pass through dense molecular clouds in interstellar media. In order to model the astrophysical phenomenon in laboratories, there are three necessary factors for the RMI to be operative: a shock wave, an external magnetic field, and density inhomogeneity. By irradiating a double-foil target with several laser beams with focal spot displacement under influence of an external magnetic field, shock waves were excited and passed through the density inhomogeneity. Radiative hydrodynamic simulations show that the RMI evolves as the density inhomogeneity is shocked, resulting in higher MFA.

Materials at high pressures and temperatures are of great current interest for warm dense matter physics, planetary sciences, and inertial fusion energy research. Shock-compression equation-of-state data and optical reflectivities of the fluid dense oxide, Gd3Ga5O12 (GGG), were measured at extremely high pressures up to 2.6 TPa (26 Mbar) generated by high-power laser irradiation and magnetically-driven hypervelocity impacts. Above 0.75 TPa, the GGG Hugoniot data approach/reach a universal linear line of fluid metals, and the optical reflectivity most likely reaches a constant value indicating that GGG undergoes a crossover from fluid semiconductor to poor metal with minimum metallic conductivity (MMC). These results suggest that most fluid compounds, e.g., strong planetary oxides, reach a common state on the universal Hugoniot of fluid metals (UHFM) with MMC at sufficiently extreme pressures and temperatures. The systematic behaviors of warm dense fluid would be useful benchmarks for developing theoretical equation-of-state and transport models in the warm dense matter regime in determining computational predictions.

Large regions of protoplanetary discs are believed to be too weakly ionised to support magnetorotational instabilities, because abundant tiny dust grains soak up free electrons and reduce the conductivity of the gas. At the outer edge of this “dead zone”, the ionisation fraction increases gradually and the resistivity drops until the magnetorotational instability can develop turbulence. We identify a new viscous instability which operates in the semi-turbulent transition region between “dead” and “alive” zones. The strength of the saturated turbulence depends strongly on the local resistivity in this transition region. A slight increase (decrease) in dust density leads to a slight increase (decrease) in resistivity and a slight decrease (increase) in turbulent viscosity. Such spatial variation in the turbulence strength causes a mass pile-up where the turbulence is weak, leading to a run-away process where turbulence is weakened and mass continues to pile up. The final result is the appearance of high-amplitude pressure bumps and deep pressure valleys. Here we present a local linear stability analysis of weakly ionised accretion discs and identify the linear instability responsible for the pressure bumps. A paper in preparation concerns numerical results which confirm and expand the existence of the linear instability.

We apply ionization balance and MHD calculations to investigate whether magnetic activity moderated by recombination on dust can account for the mass accretion rates and the mid-infrared spectra and variability of protostellar disks. The MHD calculations use the stratified shearing-box approach and include grain settling and the feedback from the changing dust abundance on the resistivity of the gas. The two-decade spread in accretion rates among T Tauri stars is too large to result solely from variety in the grain size and stellar X-ray luminosity, but can be produced by varying these together with the disk magnetic flux. The diversity in the silicate bands can come from the coupling of grain settling to the distribution of the magneto-rotational turbulence, through three effects: (1) Recombination on grains yields a magnetically inactive dead zone extending above two scale heights, while turbulence in the magnetically active disk atmosphere overshoots the dead zone boundary by only about one scale height. (2) Grains deep in the dead zone oscillate vertically in waves driven by the turbulent layer above, but on average settle at the laminar rates, so the interior of the dead zone is a particle sink and the disk atmosphere becomes dust-depleted. (3) With sufficient depletion, the dead zone is thinner and mixing dredges grains off the midplane. The MHD results also show that the magnetic activity intermittently lifts clouds of dust into the atmosphere. The photosphere height changes by up to one-third over a few orbits, while the extinction along lines of sight grazing the disk surface varies by factors of two over times down to 0.1 orbit. We suggest that the changing shadows cast by the dust clouds on the outer disk are a cause of the daily to monthly mid-infrared variability in some young stars. (Abridged.), Comment: ApJ in press

Hydrogen at high pressure in the fluid state is of great interest for understanding interiors of gas giant planets. We newly obtained Hugoniot data for liquid hydrogen up to 55 GPa under laser-driven shock loading using impedance matching to a quartz standard. The shocked temperature was determined simultaneously from the brightness temperature. The compression and temperature are almost consistent with theoretical models. High reflectivity of hydrogen was observed at 40 GPa, which suggests the fluid becomes conducting.

The origins of extrasolar gas giant planets have been discussed, based on our understanding of the gas giant planets in the solar system, Jupiter and Saturn. However, how Jupiter and Saturn formed is still uncertain because of the uncertainty in their interiors, especially the core mass (Mc). The uncertainty in Mc is partly due to those in observational data such as gravitational moments (J2n), equatorial radius (Req) and 1-bar temperatures (T1bar). New frontiers mission to Jupiter by NASA (JUNO) launched in 2011 is expected to reduce the observational errors. However, it is not necessarily clear yet which observational uncertainty dominates and how accurate observation is needed to constrain Mc enough to know the origin of Jupiter. Thus, modeling the interior of Jupiter, we evaluate each effect on Mc and required precision. We have found that the observational error of 5% in T1bar yields an error of several M⊕ in Mc. We have also found that the values of J6 of our successful models are confined in a narrow range compared to its observational error. This implies that comparison between the values of J6 of our successful models and the J6 value obtained from JUNO mission helps us to know whether the present theoretical model is valid.

Neutrino radiation takes a major role in the momentum, heat, and lepton transports in proto-neutron stars (PNSs). These diffusive processes affect the growth of magnetorotational instability (MRI) in PNSs. We perform a local linear analysis for the axisymmetric and nonaxisymmetric MRI including the effects of neutrino transports and ohmic dissipation. We find that the MRI can grow even in the multi-diffusive situations that are realized in neutrino loaded PNSs. When the toroidal magnetic component dominates over the poloidal one, nonaxisymmetric MRI modes grow much faster than axisymmetric modes. These results suggest the importance of the nonaxisymmetric MRI in PNSs. Thus the understandings of three-dimensional nonlinear evolutions of the MRI are necessary to reveal the explosion mechanism of core-collapse supernovae., Comment: Accepted for publication in ApJ, 24 pages,6 figures

We analyze the ionization state of the magnetohydrodynamically turbulent protoplanetary disks and propose a new mechanism of sustaining ionization. First, we show that in the quasi-steady state of turbulence driven by magnetorotational instability in a typical protoplanetary disk with dust grains, the amount of energy dissipation should be sufficient for providing the ionization energy that is required for activating magnetorotational instability. Second, we show that in the disk with dust grains the energetic electrons that compose electric currents in weakly ionized gas can provide collisional ionization, depending on the actual saturation state of magnetorotational turbulence. On the other hand, we show that in the protoplanetary disks with the reduced effect of dust grains, the turbulent motion can homogenize the ionization degree, leading to the activation of magnetorotational instability even in the absence of other ionization processes. The results in this Letter indicate that most of the regions in protoplanetary disks remain magnetically active, and we thus require a change in the theoretical modeling of planet formation., Comment: 11 pages, 2 figures. Accepted for publication in The Astrophysical Journal Letters

The saturation level of the magnetorotational instability (MRI) is investigated using three-dimensional MHD simulations. The shearing box approximation is adopted and the vertical component of gravity is ignored, so that the evolution of the MRI is followed in a small local part of the disk. We focus on the dependence of the saturation level of the stress on the gas pressure, which is a key assumption in the standard alpha disk model. From our numerical experiments it is found that there is a weak power-law relation between the saturation level of the Maxwell stress and the gas pressure in the nonlinear regime; the higher the gas pressure, the larger the stress. Although the power-law index depends slightly on the initial field geometry, the relationship between stress and gas pressure is independent of the initial field strength, and is unaffected by Ohmic dissipation if the magnetic Reynolds number is at least 10. The relationship is the same in adiabatic calculations, where pressure increases over time, and nearly-isothermal calculations, where pressure varies little with time. Our numerical results are qualitatively consistent with an idea that the saturation level of the MRI is determined by a balance between the growth of the MRI and the dissipation of the field through reconnection. The quantitative interpretation of the pressure-stress relation, however, may require advances in the theoretical understanding of non-steady magnetic reconnection., 45 pages, 5 tables, 17 figures, accepted for publication in ApJ

The evolution of the magnetorotational instability (MRI) during the transition from outburst to quiescence in a dwarf nova disk is investigated using three-dimensional MHD simulations. The shearing box approximation is adopted for the analysis so that the efficiency of angular momentum transport is studied in a small local patch of the disk: this is usually referred as to a one-zone model. To take account of the low ionization fraction of the disk, the induction equation includes both ohmic dissipation and the Hall effect. We induce a transition from outburst to quiescence by an instantaneous decrease of the temperature. The evolution of the MRI during the transition is found to be very sensitive to the temperature of the quiescent disk. As long as the temperature is higher than a critical value of about 2000 K, MHD turbulence and angular momentum transport is sustained by the MRI. However, MHD turbulence dies away within an orbital time if the temperature falls below this critical value. In this case, the stress drops off by more than 2 orders of magnitude and is dominated by the Reynolds stress associated with the remnant motions from the outburst. The critical temperature depends slightly on the distance from the central star and the local density of the disk.

In this paper, we present the recent results obtained regarding highly radiative shocks (RSs) generated in a low-density gas filled cell on the GEKKO XII laser facility. The RS was generated by using an ablator-pusher two-layer target (CH/Sn) and a propagation medium (Xe). High velocity RSs have been generated (100–140 km/s), while limiting as much as possible the preheating produced by the corona emission. Both self-emission and visible probe diagnostics highlighted a strong emission in the shock and an electron density in the downstream gas. The RS characteristics that depend on the initial conditions are described here as well as its precursor interaction with an aluminium foil used as an obstacle. The obtained results are discussed which show a strong extension of the radiative precursor (1 mm) leading to an expansion velocity of the obstacle up to ≈30 km/s compatible to a 20 eV temperature.

Convective dynamo simulations are performed in local Cartesian geometry. We report the first successful simulation of a large-scale oscillatory dynamo in rigidly rotating convection without stably stratified layers. A key requirement for exciting the large-scale dynamo is a sufficiently long integration time comparable to the ohmic diffusion time. By comparing two models with and without stably stratified layers, their effect on the large-scale dynamo is also studied. The spatiotemporal evolution of the large-scale magnetic field is similar in both models. However, it is intriguing that the magnetic cycle is much shorter in the model without the stable layer than with the stable layer. This suggests that the stable layer impedes the cyclic variations of the large-scale magnetic field., Comment: 5 pages, 4 figures, accepted to PASJ [Hinode Special Issue]

We investigate the stability of uniformly magnetized accretion disks, including the effect of ohmic dissipation. The growth of axisymmetric local and global modes is examined using linear perturbation theory. A simple local analysis shows that the dissipation process generally suppresses the growth of magnetorotational instability when the magnetic Reynolds number is less than unity. The characteristic length scale of unstable modes becomes longer than that of the ideal MHD case, and its unstable growth rate is inversely proportional to the magnetic diffusivity. We perform a global linear analysis in which the vertical structure of the disk is considered. The growth rate of the magnetorotational instability is obtained by solving eigen equations numerically. We find that the conditions for existing unstable global modes are ?c1 and ?cRmc1, where ?c and Rmc are the plasma beta value and the magnetic Reynolds number defined at the midplane of the disk. The global stability criteria are approximately given by whether the minimum unstable wavelength expected by the local analysis would be shorter than the scale height of the disk or not. We also find unstable modes whose eigenfunctions of the perturbed velocities have amplitude localized near the surface layer of the disk. These unstable modes indicate layered accretion in the nonlinear regime. We apply the results of linear analysis to protoplanetary disks. For the case of the minimum-mass solar nebula, the magnetorotational instability occurs at the region farther out than 15 AU. This result suggests nonsteady accretion onto a central star in protoplanetary disks.

The nonlinear saturation of the magnetorotational instability due to ohmic dissipation is investigated with two-dimensional magnetohydrodynamic simulations of accretion disks. We adopt a local disk model in which the physical quantities are assumed to be spatially uniform except for the shear velocity in the azimuthal direction. Starting with a weak vertical field, the magnetorotational instability is saturated in the nonlinear regime when the magnetic Reynolds number Rm is less than unity. In the saturated turbulent state, the growth of the magnetic field for the instability and the damping by the ohmic dissipation are almost balanced. The efficiency of angular momentum transport α is of the order of 10-2 to 10-3, and α is found to be inversely proportional to the initial β value. When Rm1, channel solutions appear and the magnetic energy continues growing even with the dissipation process. We can roughly explain these features with the help of the results obtained by the linear analysis. Our conclusion is that the magnetic dissipation is one of the most important processes that determine the saturation level of the magnetorotational instability.

The critical strength of a magnetic field required for the suppression of the Richtmyer-Meshkov instability (RMI) is investigated numerically by using a two-dimensional single-mode analysis. For the cases of MHD parallel shocks, the RMI can be stabilized as a result of the extraction of vorticity from the interface. A useful formula describing a critical condition for MHD RMI has been introduced, and which is successfully confirmed by the direct numerical simulations. The critical field strength is found to be largely depending on the Mach number of the incident shock. If the shock is strong enough, even low-$\beta$ plasmas can be subject to the growth of the RMI., Comment: 5 pages, 3 figures, accepted for publication in PRL

We examine the conditions under which the disks of gas and dust orbiting young gas giant planets are sufficiently conducting to experience turbulence driven by the magneto-rotational instability. By modeling the ionization and conductivity in the disk around proto-Jupiter, we find that turbulence is possible if the X-rays emitted near the Sun reach the planet's vicinity and either (1) the gas surface densities are in the range of the minimum-mass models constructed by augmenting Jupiter's satellites to Solar composition, while dust is depleted from the disk atmosphere, or (2) the surface densities are much less, and in the range of gas-starved models fed with material from the Solar nebula, but not so low that ambipolar diffusion decouples the neutral gas from the plasma. The results lend support to both minimum-mass and gas-starved models of the protojovian disk: (1) The dusty minimum-mass models have internal conductivities low enough to prevent angular momentum transfer by magnetic forces, as required for the material to remain in place while the satellites form. (2) The gas-starved models have magnetically-active surface layers and a decoupled interior "dead zone". Similar active layers in the Solar nebula yield accretion stresses in the range assumed in constructing the circumjovian gas-starved models. Our results also point to aspects of both classes of models that can be further developed. Non-turbulent minimum-mass models will lose dust from their atmospheres by settling, enabling gas to accrete through a thin surface layer. For the gas-starved models it is crucial to learn whether enough stellar X-ray and ultraviolet photons reach the circumjovian disk. Additionally the stress-to-pressure ratio ought to increase with distance from the planet, likely leading to episodic accretion outbursts., 35 pages, 8 figures. Chemical network description expanded in sec. 4.1. Turbulent mixing discussion reorganized in secs. 6.3-6.4. Accepted to ApJ

We report our recent efforts on the experimental investigations related to the origins of cosmic rays. The origins of cosmic rays are long standing open issues in astrophysics. The galactic and extragalactic cosmic rays are considered to be accelerated in non-relativistic and relativistic collisionless shocks in the universe, respectively. However, the acceleration and transport processes of the cosmic rays are not well understood, and how the collisionless shocks are created is still under investigation. Recent high-power and high-intensity laser technologies allow us to simulate astrophysical phenomena in laboratories. We present our experimental results of collisionless shock formations in laser-produced plasmas. © 2012 IOP Publishing Ltd.

We report experimental results on Kelvin–Helmholtz (KH) instability and resultant vortices in laser-produced plasmas. By irradiating a double plane target with a laser beam, asymmetric counterstreaming plasmas are created. The interaction of the plasmas with different velocities and densities results in the formation of asymmetric shocks, where the shear flow exists along the contact surface and the KH instability is excited. We observe the spatial and temporal evolution of plasmas and shocks with time-resolved diagnostics over several shots. Our results clearly show the evolution of transverse fluctuations, wavelike structures, and circular features, which are interpreted as the KH instability and resultant vortices. The relevant numerical simulations demonstrate the time evolution of KH vortices and show qualitative agreement with experimental results. Shocks, and thus the contact surfaces, are ubiquitous in the universe; our experimental results show general consequences where two plasmas interact.