Collisional Simulations Research Articles

Effective Resistivity in Collisionless Magnetic Reconnection

An effective resistivity relevant to collisionless magnetic reconnection (MR) in plasma is presented. It is based on the argument that pitch angle scattering of electrons in the small electron diffusion region around the X line can lead to an effective, resistivity in collisionless plasma. The effective resistivity so obtained is in the form of a power law of the local plasma and magnetic field parameters. Its validity is confirmed by direct collisionless particle-in-cell (PIC) simulation. The result agrees very well with the resistivity (obtained from available data) of a large number of environments susceptible to MR: from the intergalactic and interstellar to solar and terrestrial to laboratory fusion plasmas. The scaling law can readily be incorporated into existing collisional magnetohydrodynamic simulation codes to investigate collisionless MR, as well as serve as a guide to ab initio theoretical investigations of the collisionless MR process.

Scaling of spoke rotation frequency within a Penning discharge

A rotating plasma spoke is shown to develop in two-dimensional full-sized kinetic simulations of a Penning discharge cross-section. Electron cross-field transport within the discharge is highly anomalous and correlates strongly with the spoke phase. Similarity between collisional and collisionless simulations demonstrates that ionization is not necessary for spoke formation. Parameter scans with discharge current Id, applied magnetic field strength B, and ion mass mi show that the spoke frequency scales with eErLn/mi, where Er is the radial electric field, Ln is the gradient length scale, and e is the fundamental charge. This scaling suggests that the spoke may develop as a non-linear phase of the collisionless Simon-Hoh instability.

Experiments on wake structures behind a microparticle in a magnetized plasma flow

The wake behind a spherical microparticle in a magnetized ion flow is studied experimentally by analyzing the arrangement of a pair of particles. It is shown that there are two stable particle arrangements at intermediate magnetic inductions, whereas only oblique (horizontal) particle configurations are found at the highest magnetic field. Self-consistent collisional molecular dynamics simulations of the particle system show that the underlying mechanism of these arrangements is the weakening of attractive wake forces by the increasing magnetic field. Plasma instabilities provide a trigger for the onset of the transition between the two different arrangements. Furthermore, the course of the transition is qualitatively explained by the charge variation of the downstream particle in the wake of the upstream particle. In addition, a thorough analysis of the sheath by means of particle-in-cell simulations in combination with particle resonance measurements yields consistent values of the particle mass and charge, as well as the levitating electric field and ion flow velocity.

ISR Spectra Simulations With Electron‐Ion Coulomb Collisions

AbstractIncoherent scatter radars (ISR) rely on Thomson scattering of very high frequency or ultrahigh frequency radio waves off electrons in the ionosphere and measure the backscattered power spectra in order to estimate altitude profiles of plasma density, electron temperature, ion temperature, and ion drift speed. These spectra result from the collective behavior of coupled ion and electron dynamics, and, for most cases, existing theories predict these well. However, when the radar points nearly perpendicular to the Earth's magnetic field, the motion of the plasma across the field lines becomes complex and Coulomb collisions between electrons and ions become important in interpreting ISR measurements. This paper presents the first fully kinetic, self‐consistent, particle‐in‐cell simulations of ISR spectra with electron‐ion Coulomb collisions. We implement a grid‐based Coulomb collision algorithm in the Electrostatic Parallel Particle‐in‐Cell simulator and obtain ISR spectra from simulations both with and without collisions. For radar directions greater than 5° away from perpendicular to the magnetic field, both sets of simulations match collisionless ISR theory well. For angles between 3° and 5°, the collisional simulation is well described by a simplified Brownian motion collision process. At angles less than 3° away from perpendicular the Brownian motion model fails, and the collisional simulation qualitatively agrees with previous single particle simulations. For radar directions exactly perpendicular to the magnetic field the simulated collisional spectra match those from the Brownian motion collision theory, in agreement with previous single particle simulations.

Clocking Femtosecond Collisional Dynamics via Resonant X-Ray Spectroscopy.

Electron-ion collisional dynamics is of fundamental importance in determining plasma transport properties, nonequilibrium plasma evolution, and electron damage in diffraction imaging applications using bright x-ray free-electron lasers (FELs). Here we describe the first experimental measurements of ultrafast electron impact collisional ionization dynamics using resonant core-hole spectroscopy in a solid-density magnesium plasma, created and diagnosed with the Linac Coherent Light Source x-ray FEL. By resonantly pumping the 1s→2p transition in highly charged ions within an optically thin plasma, we have measured how off-resonance charge states are populated via collisional processes on femtosecond time scales. We present a collisional cross section model that matches our results and demonstrates how the cross sections are enhanced by dense-plasma effects including continuum lowering. Nonlocal thermodynamic equilibrium collisional radiative simulations show excellent agreement with the experimental results and provide new insight on collisional ionization and three-body-recombination processes in the dense-plasma regime.

Debris disc constraints on planetesimal formation

Two basic routes for planetesimal formation have been proposed over the last few decades. One is a classical "slow-growth" scenario. Another one is particle concentration models, in which small pebbles are concentrated locally and then collapse gravitationally to form planetesimals. Both types of models make certain predictions for the size spectrum and internal structure of newly-born planetesimals. We use these predictions as input to simulate collisional evolution of debris discs left after the gas dispersal. The debris disc emission as a function of a system's age computed in these simulations is compared with several Spitzer and Herschel debris disc surveys around A-type stars. We confirm that the observed brightness evolution for the majority of discs can be reproduced by classical models. Further, we find that it is equally consistent with the size distribution of planetesimals predicted by particle concentration models - provided the objects are loosely bound "pebble piles" as these models also predict. Regardless of the assumed planetesimal formation mechanism, explaining the brightest debris discs in the samples uncovers a "disc mass problem." To reproduce such discs by collisional simulations, a total mass of planetesimals of up to ~1000 Earth masses is required, which exceeds the total mass of solids available in the protoplanetary progenitors of debris discs. This may indicate that stirring was delayed in some of the bright discs, that giant impacts occurred recently in some of them, that some systems may be younger than previously thought, or that non-collisional processes contribute significantly to the dust production.

Characterization of Protonated Model Disaccharides from Tandem Mass Spectrometry and Chemical Dynamics Simulations.

The fragmentation mechanisms of prototypical disaccharides have been studied herein by coupling tandem mass spectrometry (MS) with collisional chemical dynamics simulations. These calculations were performed by explicitly considering the collisions between the protonated sugar and the neutral target gas, which led to an ensemble of trajectories for each system, from which it was possible to obtain reaction products and mechanisms without pre-imposing them. The β-aminoethyl and aminopropyl derivatives of cellobiose, maltose, and gentiobiose were studied to observe differences in both the stereochemistry and the location of the glycosidic linkage. Chemical dynamics simulations of MS/MS and MS/MS/MS were used to suggest some primary and secondary fragmentation mechanisms for some experimentally observed product ions. These simulations provided some new insights into the fundamentals of the unimolecular dissociation of protonated sugars under collisional induced dissociation conditions.

Effect of Radiation on Chromospheric Magnetic Reconnection: Reactive and Collisional Multi-fluid Simulations

Abstract We study magnetic reconnection under chromospheric conditions in five different ionization levels from 0.5% to 50% using a self-consistent two-fluid (ions + neutrals) model that accounts for compressibility, collisional effects, chemical inequilibrium, and anisotropic heat conduction. Results with and without radiation are compared, using two models for the radiative losses: an optically thin radiation loss function, and an approximation of the radiative losses of a plasma with photospheric abundances. The results without radiation show that reconnection occurs faster for the weakly ionized cases as a result of the effect of ambipolar diffusion and fast recombination. The tearing mode instability appears earlier in the low ionized cases and grows rapidly. We find that radiative losses have a stronger effect than was found in previous results as the cooling changes the plasma pressure and the concentration of ions inside the current sheet. This affects the ambipolar diffusion and the chemical equilibrium, resulting in thin current sheets and enhanced reconnection. The results quantify this complex nonlinear interaction by showing that a strong cooling produces faster reconnections than have been found in models without radiation. The results accounting for radiation show timescales and outflows comparable to spicules and chromospheric jets.

Kinetic simulations of gas breakdown in the dense plasma focus

The first fully kinetic, collisional, and electromagnetic simulations of the breakdown phase of a MA-scale dense plasma focus are described and shown to agree with measured electrical characteristics, including breakdown time. In the model, avalanche ionization is driven by cathode electron emission, and this results in incomplete gas breakdown along the insulator. This reinforces the importance of the conditioning process that creates a metallic layer on the insulator surface. The simulations, nonetheless, help explain the relationship between the gas pressure, the insulator length, and the coaxial gap width. Previously, researchers noted three breakdown patterns related to pressure. Simulation and analytical results show that at low pressures, long ionization path lengths lead to volumetric breakdown, while high pressures lead to breakdown across the relatively small coaxial electrode gap. In an intermediate pressure regime, ionization path lengths are comparable to the insulator length which promotes ideal breakdown along the insulator surface.

Collisions between Sintered Icy Aggregates

Abstract Collisions between sintered icy dust aggregates are numerically simulated. If the temperature of an icy aggregate is sufficiently high, sintering promotes molecular transport and a neck between adjacent grains grows. This growth changes the mechanical responses of the neck. We included this effect in a simulation code, and conducted collisional simulations. For porous aggregates, the critical velocity for growth, below which the mass of an aggregate increases, decreased from 50 m s−1 for the non-sintered case to 20 m s−1. For compacted aggregates, the main collisional outcome is bouncing. These results come from the fact that the strength of the neck is increased by sintering. The numerical results suggest that the collisional growth of icy grain aggregates is strongly affected by sintering.

HD 106906: A Case Study for External Perturbations of a Debris Disk

Abstract Models of debris disk morphology are often focused on the effects of a planet orbiting interior to or within the disk. Nonetheless, an exterior planetary-mass perturber can also excite eccentricities in a debris disk, via Laplace–Lagrange secular perturbations in the coplanar case or Kozai–Lidov perturbations for mutually inclined companions and disks. HD 106906 is an ideal example of such a a system, as it harbors a confirmed exterior companion at a projected separation of 650 au outside a resolved, asymmetric disk. We use collisional and dynamical simulations to investigate the interactions between the disk and the companion, and to use the disk’s observed morphology to place constraints on the companion’s orbit. We conclude that the disk’s observed morphology is consistent with perturbations from the observed exterior companion. Generalizing this result, we suggest that exterior perturbers, as well as interior planets, should be considered when investigating the cause of observed asymmetries in a debris disk.

Features of collisionless turbulence in the intracluster medium from simulated Faraday rotation maps – II. The effects of instabilities feedback

Statistical analysis of Faraday Rotation Measure (RM) maps of the intracluster medium (ICM) of galaxy clusters provides a unique tool to evaluate some spatial features of the magnetic fields there. Its combination with numerical simulations of magnetohydrodynamic (MHD) turbulence allows the diagnosis of the ICM turbulence. Being the ICM plasma weakly collisional, the thermal velocity distribution of the particles naturally develops anisotropies as a consequence of the large scale motions and the conservation of the magnetic moment of the charged particles. A previous study (Paper I) analyzed the impact of large scale thermal anisotropy on the statistics of RM maps synthesized from simulations of turbulence; these simulations employed a collisionless MHD model which considered a tensor pressure with uniform anisotropy. In the present work, we extend that analysis to a collisionless MHD model in which the thermal anisotropy develops according to the conservation of the magnetic moment of the thermal particles. We also consider the effect of anisotropy relaxation caused by the micro-scale mirror and firehose instabilities. We show that if the relaxation rate is fast enough to keep the anisotropy limited by the threshold values of the instabilities, the dispersion and power spectrum of the RM maps are indistinguishable from those obtained from collisional MHD. Otherwise, there is a reduction in the dispersion and steepening of the power spectrum of the RM maps (compared to the collisional case). Considering the first scenario, the use of collisional MHD simulations for modeling the RM statistics in the ICM becomes better justified.

Origin and evolution of two-component debris discs and an application to the q1Eridani system

Many debris discs reveal a two-component structure, with an outer Kuiper-belt analogue and a warm inner component whose origin is still a matter of debate. One possibility is that warm emission stems from an "asteroid belt" closer in to the star. We consider a scenario in which a set of giant planets is formed in an initially extended planetesimal disc. These planets carve a broad gap around their orbits, splitting up the disc into the outer and the inner belts. After the gas dispersal, both belts undergo collisional evolution in a steady-state regime. This scenario is explored with detailed collisional simulations involving realistic physics to describe a long-term collisional depletion of the two-component disc. We find that the inner disc may be able to retain larger amounts of material at older ages than thought before on the basis of simplified analytic models. We show that the proposed scenario is consistent with a suite of thermal emission and scattered light observational data for a bright two-temperature debris disc around a nearby solar-type star q$^1$ Eridani. This implies a Solar System-like architecture of the system, with an outer massive "Kuiper belt", an inner "asteroid belt", and a few Neptune- to Jupiter-mass planets in between.

The dust grain size–stellar luminosity trend in debris discs

The cross section of material in debris discs is thought to be dominated by the smallest grains that can still stay in bound orbits despite the repelling action of stellar radiation pressure. Thus the minimum (and typical) grain size $s_\text{min}$ is expected to be close to the radiation pressure blowout size $s_\text{blow}$. Yet a recent analysis of a sample of Herschel-resolved debris discs showed the ratio $s_\text{min}/s_\text{blow}$ to systematically decrease with the stellar luminosity from about ten for solar-type stars to nearly unity in the discs around the most luminous A-type stars. Here we explore this trend in more detail, checking how significant it is and seeking to find possible explanations. We show that the trend is robust to variation of the composition and porosity of dust particles. For any assumed grain properties and stellar parameters, we suggest a recipe of how to estimate the "true" radius of a spatially unresolved debris disc, based solely on its spectral energy distribution. The results of our collisional simulations are qualitatively consistent with the trend, although additional effects may also be at work. In particular, the lack of grains with small $s_\text{min}/s_\text{blow}$ for lower luminosity stars might be caused by the grain surface energy constraint that should limit the size of the smallest collisional fragments. Also, a better agreement between the data and the collisional simulations is achieved when assuming debris discs of more luminous stars to have higher dynamical excitation than those of less luminous primaries. This would imply that protoplanetary discs of more massive young stars are more efficient in forming big planetesimals or planets that act as stirrers in the debris discs at the subsequent evolutionary stage.

FOMALHAUT b AS A DUST CLOUD: FREQUENT COLLISIONS WITHIN THE FOMALHAUT DISK

The planet candidate Fomalhaut b is bright in optical light but undetected in longer wavelengths, requiring a large, reflective dust cloud. The most recent observations find an extremely eccentric orbit (e ~ 0.8), indicating that Fomalhaut b cannot be the planet that is constraining the system's eccentric debris ring. An irregular satellite swarm around a super-Earth has been proposed, however, explaining the well-constrained debris ring requires an additional planet on an orbit that crosses that of the putative super-Earth. This paper expands upon a second theory: Fomalhaut b is a transient dust cloud produced by a catastrophic collision between planetesimals in the disk. We perform collisional probability simulations of the Fomalhaut debris disk based on the structure of our Kuiper belt, finding that the catastrophic disruption rate of d ~ 100 km bodies in the high-eccentricity scattering component is several per decade. This model paints a picture of the Fomalhaut system as having recently (within ~10-100 Myr) experienced a dynamical instability within its planetary system, which scattered a massive number of planetesimals onto large, high eccentricity orbits similar to that of Fom b. If Fomalhaut b is indeed a dust cloud produced by such a collision, we should soon see another appear, while Fomalhaut b will expand until it is either resolved or becomes too faint to be seen.

Transition in electron physics of magnetic reconnection in weakly collisional plasma

Using particle-in-cell (PIC) simulations with a Monte Carlo treatment of the Coulomb collision operator, we study the transition in electron dynamics of magnetic reconnection for various levels of collisionality. The weakly collisional cases considered all fall into the so-called Hall or kinetic regime. Nevertheless, collisions may still alter the electron kinetic physics characteristic of collisionless reconnection, where adiabatic trapping energizes the electrons and leads to strong anisotropy of the electron velocity distribution and pressure. This anisotropy can support extended current sheets, associated with secondary island formation and turbulent flux rope interactions in three dimensional systems. The collisional simulations demonstrate how weak collisions may modify or eliminate these electron structures in the kinetic regimes. While the reconnection rate is not sensitive to the collisionality in the range studied, we find that increasing collisionality reduces the level of electron energization near the reconnection site. Finally, the results provide guidance for new laboratory reconnection experiments that will access the weakly collisional regimes.

TRI3DYN – Collisional computer simulation of the dynamic evolution of 3-dimensional nanostructures under ion irradiation

A computer program is described which allows fully three-dimensional dynamic collisional simulations of ion irradiation effects in particular in nanosystems. The program is based on the well-known TRIDYN code which describes dynamic modifications along one dimension. In a computational volume composed of fixed and equisized voxels, local atomic density changes due to ion implantation or recoil relocation are relaxed to constant predefined atomic volumes. This is accomplished by interaction with neighboring voxels and by material transport towards and away from the surface. Arbitrary initial system shapes and elemental distributions can be defined as well as a wide range of irradiation conditions including non-uniform beam profiles. The paper addresses details of the involved algorithms and discusses potential artifacts with respect to surface sputtering, broadening of compositional profiles, and surface contours. It demonstrates the capabilities of the program by examples of focused ion beam erosion, self-organized surface pattering, and ion implantation and sputter-shaping of nanostructures.

MAGNETIC FIELD AMPLIFICATION AND EVOLUTION IN TURBULENT COLLISIONLESS MAGNETOHYDRODYNAMICS: AN APPLICATION TO THE INTRACLUSTER MEDIUM

The amplification and maintenance of the observed magnetic fields in the ICM are usually attributed to the turbulent dynamo action. This is generally derived employing a collisional MHD model. However, in the ICM the ion mean free path between collisions is of the order of the dynamical scales, thus requiring a collisionless MHD description. Unlike collisional MHD simulations, our study uses an anisotropic plasma pressure with respect to the direction of the local magnetic field, which brings the plasma within a parameter space where collisionless instabilities should take place. Within the adopted model these instabilities are contained at bay through the relaxation term of the pressure anisotropy which simulates the feedback of the mirror and firehose instabilities. This relaxation acts to get the plasma distribution function consistent with the empirical studies of collisionless plasmas. Our 3D numerical simulations of forced transonic turbulence motivated by modeling of the turbulent ICM are performed for different initial values of the magnetic field intensity, and different relaxation rates of the pressure anisotropy. We found that in the high beta plasma regime corresponding to the ICM conditions, a fast anisotropy relaxation rate gives results which are similar to the collisional-MHD model as far as the statistical properties of the turbulence are concerned. Also, the amplification of seed magnetic fields due to the turbulent dynamo action is similar to the collisional-MHD model. Our simulations that do not employ the anisotropy relaxation deviate significantly from the collisional-MHD results, and show more power at the small-scale fluctuations of both density and velocity representing the results of the instabilities. For these simulations the large scale fluctuations in the magnetic field are mostly suppressed and the turbulent dynamo fails in amplifying seed magnetic fields.

HERSCHEL's “COLD DEBRIS DISKS”: BACKGROUND GALAXIES OR QUIESCENT RIMS OF PLANETARY SYSTEMS?

(abridged) Infrared excesses associated with debris disk host stars detected so far peak at wavelengths around ~100{\mu}m or shorter. However, six out of 31 excess sources in the Herschel OTKP DUNES have been seen to show significant - and in some cases extended - excess emission at 160{\mu}m, which is larger than the 100{\mu}m excess. This excess emission has been suggested to stem from debris disks colder than those known previously. Using several methods, we re-consider whether some or even all of the candidates may be associated with unrelated galactic or extragalactic emission and conclude that it is highly unlikely that none of the candidates represents a true circumstellar disk. For true disks, both the dust temperatures inferred from the SEDs and the disk radii estimated from the images suggest that the dust is nearly as cold as a blackbody. This requires the grains to be larger than ~100{\mu}m, regardless of their material composition. To explain the dearth of small grains, we explore several conceivable scenarios: transport-dominated disks, disks of low dynamical excitation, and disks of unstirred primordial macroscopic grains. Our qualitative analysis and collisional simulations rule out the first two of these scenarios, but show the feasibility of the third one. We show that such disks can survive for gigayears, largely preserving the primordial size distribution. They should be composed of macroscopic solids larger than millimeters, but smaller than kilometers in size. Thus planetesimal formation, at least in the outer regions of the systems, has stopped before "cometary" or "asteroidal" sizes were reached.

Irradiation-induced tailoring of the magnetism of CoFeB/MgO ultrathin films

We study perpendicularly magnetized Ta/CoFeB/MgO films and investigate whether their irradiation with light ions can improve their properties by inducing a different crystallization dynamics. We report the magnetization, anisotropy, g-factor, and damping dependence upon irradiation fluence and discuss their evolutions with collisional mixing simulations and its expected consequence on magnetic properties. We show that after a short irradiation at 100 °C, the anisotropy increases close to the value obtained by conventional high temperature annealing. Higher irradiation-induced increase of anisotropy can be obtained but with a detrimental effect on the damping that can be understood from spin-orbit contributions.