Particle-in-cell simulation study of a lower-hybrid shock

The expansion of a magnetized high-pressure plasma into a low-pressure ambient medium is examined with particle-in-cell simulations. The magnetic field points perpendicular to the plasma's expansion direction and binary collisions between particles are absent. The expanding plasma steepens into a quasi-electrostatic shock that is sustained by the lower-hybrid (LH) wave. The ambipolar electric field points in the expansion direction and it induces together with the background magnetic field a fast E cross B drift of electrons. The drifting electrons modify the background magnetic field, resulting in its pile-up by the LH shock. The magnetic pressure gradient force accelerates the ambient ions ahead of the LH shock, reducing the relative velocity between the ambient plasma and the LH shock to about the phase speed of the shocked LH wave, transforming the LH shock into a nonlinear LH wave. The oscillations of the electrostatic potential have a larger amplitude and wavelength in the magnetized plasma than in an unmagnetized one with otherwise identical conditions. The energy loss to the drifting electrons leads to a noticeable slowdown of the LH shock compared to that in an unmagnetized plasma.

Bare-tether sheath and current: comparison of asymptotic theory and kinetic simulations in stationary plasma

Analytical expressions for current to a cylindrical Langmuir probe at rest in unmagnetized plasma are compared with results from both steady-state Vlasov and particle-in-cell simulations. Probe bias potentials that are much greater than plasma temperature (assumed equal for ions and electrons), as of interest for bare conductive tethers, are considered. At a very high bias, both the electric potential and the attracted-species density exhibit complex radial profiles; in particular, the density exhibits a minimum well within the plasma sheath and a maximum closer to the probe. Excellent agreement is found between analytical and numerical results for values of the probe radiusR close to the maximum radius Rmax for orbital-motion-limited (OML) collection at a particular bias in the following number of profile features: the values and positions of density minimum and maximum, position of sheath boundary, and value of a radius characterizing the no-space-charge behavior of a potential near the high-bias probe. Good agreement between the theory and simulations is also found for parametric laws jointly covering the following three characteristic R ranges: sheath radius versus probe radius and bias for Rmax; density minimum versus probe bias for Rmax; and (weakly bias-dependent) current drop below the OML value versus the probe radius for R > Rmax.

Ion acceleration in underdense plasmas by ultra-short laser pulses

We report on the ion acceleration mechanisms that occur during the interaction of an intense and ultrashort laser pulse ( λ > μ I 2 1018 W cm−2 m2) with an underdense helium plasma produced from an ionized gas jet target. In this unexplored regime, where the laser pulse duration is comparable to the inverse of the electron plasma frequency ωpe, reproducible non-thermal ion bunches have been measured in the radial direction. The two He ion charge states present energy distributions with cutoff energies between 150 and 200 keV, and a striking energy gap around 50 keV appearing consistently for all the shots in a given density range. Fully electromagnetic particle-in-cell simulations explain the experimental behaviors. The acceleration results from a combination of target normal sheath acceleration and Coulomb explosion of a filament formed around the laser pulse propagation axis

Microscopic processes in global relativistic jets containing helical magnetic fields

In the study of relativistic jets one of the key open questions is their interaction with the environment on the microscopic level. Here, we study the initial evolution of both electron–proton (e−–p+) and electron–positron (e±) relativistic jets containing helical magnetic fields, focusing on their interaction with an ambient plasma. We have performed simulations of “global” jets containing helical magnetic fields in order to examine how helical magnetic fields affect kinetic instabilities such as the Weibel instability, the kinetic Kelvin-Helmholtz instability (kKHI) and the Mushroom instability (MI). In our initial simulation study these kinetic instabilities are suppressed and new types of instabilities can grow. In the e−–p+ jet simulation a recollimation-like instability occurs and jet electrons are strongly perturbed. In the e± jet simulation a recollimation-like instability occurs at early times followed by a kinetic instability and the general structure is similar to a simulation without helical magnetic field. Simulations using much larger systems are required in order to thoroughly follow the evolution of global jets containing helical magnetic fields.

Towards optical polarization control of laser-driven proton acceleration in foils undergoing relativistic transparency

Control of the collective response of plasma particles to intense laser light is intrinsic to relativistic optics, the development of compact laser-driven particle and radiation sources, as well as investigations of some laboratory astrophysics phenomena. We recently demonstrated that a relativistic plasma aperture produced in an ultra-thin foil at the focus of intense laser radiation can induce diffraction, enabling polarization-based control of the collective motion of plasma electrons. Here we show that under these conditions the electron dynamics are mapped into the beam of protons accelerated via strong charge-separation-induced electrostatic fields. It is demonstrated experimentally and numerically via 3D particle-in-cell simulations that the degree of ellipticity of the laser polarization strongly influences the spatial-intensity distribution of the beam of multi-MeV protons. The influence on both sheath-accelerated and radiation pressure-accelerated protons is investigated. This approach opens up a potential new route to control laser-driven ion sources.

Numerical and Analytical Methods for Laser-Plasma Acceleration Physics

Theories and numerical modeling are fundamental tools for understanding, optimizing and designing present and future laser-plasma accelerators (LPAs). Laser evolution and plasma wave excitation in a LPA driven by a weakly relativistically intense, short-pulse laser propagating in a preformed parabolic plasma channel, is studied analytically in 3D including the effects of pulse steepening and energy depletion. At higher laser intensities, the process of electron self-injection in the nonlinear bubble wake regime is studied by means of fully self-consistent Particle-in-Cell simulations. Considering a non-evolving laser driver propagating with a prescribed velocity, the geometrical properties of the non-evolving bubble wake are studied. For a range of parameters of interest for laser plasma acceleration, The dependence of the threshold for self-injection in the non-evolving wake on laser intensity and wake velocity is characterized. Due to the nonlinear and complex nature of the Physics involved, computationally challenging numerical simulations are required to model laser-plasma accelerators operating at relativistic laser intensities. The numerical and computational optimizations, that combined in the codes INF&RNO and INF&RNO/quasi-static give the possibility to accurately model multi-GeV laser wakefield acceleration stages with present supercomputing architectures, are discussed. The PIC code jasmine, capable of efficiently running laser-plasma simulations on Graphics Processing Units (GPUs) clusters, is presented. GPUs deliver exceptional performance to PIC codes, but the core algorithms had to be redesigned for satisfying the constraints imposed by the intrinsic parallelism of the architecture. The simulation campaigns, run with the code jasmine for modeling the recent LPA experiments with the INFN-FLAME and CNR-ILIL laser systems, are also presented.

Ultra intense laser/plasma interaction at normal incidence: Relativistic mirrors effects, high harmonics generation and absorption.

An analytical study of the relativistic interaction of a linearly-polarized laser-field of w frequency with highly overdense plasma is presented. Very intense high harmonics are generated produced by relativistic mirrors effects due to the relativistic electron plasma oscillation. Also, in agreement with 1D Particle-In-Cell Simulations (PICS), the model self-consistently explains the transition between the sheath inverse bremsstrahlung (SIB) absorption regime and the J×B heating (responsible for the 2w electron bunches), as well as the mean electron energy.

Efficient on the fly calculation of time correlation functions in computer simulations

Time correlation functions yield profound information about the dynamics of a physical system and hence are frequently calculated in computer simulations. For systems whose dynamics span a wide range of time, currently used methods require significant computer time and memory. In this paper, we discuss the multiple-tau correlator method for the efficient calculation of accurate time correlation functions on the fly during computer simulations. The multiple-tau correlator is efficacious in terms of computational requirements and can be tuned to the desired level of accuracy. Further, we derive estimates for the error arising from the use of the multiple-tau correlator and extend it for use in the calculation of mean-square particle displacements and dynamic structure factors. The method described here, in hardware implementation, is routinely used in light scattering experiments but has not yet found widespread use in computer simulations.

Classifying general nonlinear force laws in cell-based models via the continuum limit

though discrete cell-based frameworks are now commonly used to simulate a whole range of biological phenomena, it is typically not obvious how the numerous different types of model are related to one another, nor which one is most appropriate in a given context. Here we demonstrate how individual cell movement on the discrete scale modeled using nonlinear force laws can be described by nonlinear diffusion coefficients on the continuum scale. A general relationship between nonlinear force laws and their respective diffusion coefficients is derived in one spatial dimension and, subsequently, a range of particular examples is considered. For each case excellent agreement is observed between numerical solutions of the discrete and corresponding continuum models. Three case studies are considered in which we demonstrate how the derived nonlinear diffusion coefficients can be used to (a) relate different discrete models of cell behavior; (b) derive discrete, intercell force laws from previously posed diffusion coefficients, and (c) describe aggregative behavior in discrete simulations.

A study on the inclusion of forest canopy morphology data in numerical simulations for the purpose of wind resource assessment

A series of numerical simulations of the flow over a forest stand have been conducted using two different turbulence closure models along with various levels of canopy morphology data. Simulations have been validated against Stereoscopic Particle Image Velocimetry measurements from a wind tunnel study using one hundred architectural model trees, the porosities of which have been assessed using a photographic technique. It has been found that an accurate assessment of the porosity of the canopy, and specifically the variability with height, improves simulation quality regardless of the turbulence closure model used or the level of canopy geometry included. The observed flow field and recovery of the wake is in line with characteristic canopy flows published in the literature and it was found that the shear stress transport turbulence model was best able to capture this detail numerically.

Characterisation Of The Membrane Transport Of Pilocarpine In Cell Suspension Cultures Of Pilocarpus Microphyllus.

Pilocarpine is an alkaloid obtained from the leaves of Pilocarpus genus, with important pharmaceutical applications. Previous reports have investigated the production of pilocarpine by Pilocarpus microphyllus cell cultures and tried to establish the alkaloid biosynthetic route. However, the site of pilocarpine accumulation inside of the cell and its exchange to the medium culture is still unknown. Therefore, the aim of this study was to determine the intracellular accumulation of pilocarpine and characterise its transport across membranes in cell suspension cultures of P. microphyllus. Histochemical analysis and toxicity assays indicated that pilocarpine is most likely stored in the vacuoles probably to avoid cell toxicity. Assays with exogenous pilocarpine supplementation to the culture medium showed that the alkaloid is promptly uptaken but it is rapidly metabolised. Treatment with specific ABC protein transporter inhibitors and substances that disturb the activity of secondary active transporters suppressed pilocarpine uptake and release suggesting that both proteins may participate in the traffic of pilocarpine to inside and outside of the cells. As bafilomicin A1, a specific V-type ATPase inhibitor, had little effect and NH4Cl (induces membrane proton gradient dissipation) had moderate effect, while cyclosporin A and nifedipine (ABC proteins inhibitors) strongly inhibited the transport of pilocarpine, it is believed that ABC proteins play a major role in the alkaloid transport across membranes but it is not the exclusive one. Kinetic studies supported these results.

Vimentin expression and the influence of Matrigel in cell lines of head and neck squamous cell carcinoma

Vimentin is a cytoeskeletal intermediate filament protein commonly observed in mesenchymal cells; however, it can also be found in malignant epithelial cells. It is demonstrated in several carcinomas, such as those of the cervix, breast and bladder, in which it is widely used as a marker of the epithelial to mesenchymal transition that takes place during embryogenesis and metastasis. Vimentin is associated with tumors that show a high degree of invasiveness, being detected in invasion front cells. Its expression seems to be influenced by the tumor microenvironment. The aim of this study was to evaluate vimentin expression in head and neck squamous cell carcinoma (HNSCC) cell lines, and to investigate the contribution of the microenvironment to its expression. HNSCC cell lines (HN6, HN30 and HN31) and an immortalized nontumorigenic cell line (HaCaT) were submitted to a three-dimensional assay with Matrigel. Cytoplasmatic staining of the HN6 cell line cultured without Matrigel and of the HN30 and HN31 cell lines cultured with Matrigel was demonstrated through immunohistochemistry. Western Blotting revealed a significant decrease in vimentin expression for the HN6 cell line and a significant increase for the HN30 and HN31 cell lines cultured with Matrigel. The results suggest that vimentin can be expressed in HNSCC cells and its presence is influenced by the microenvironment of a tumor.

Correlation between MMPs and their inhibitors in breast cancer tumor tissue specimens and in cell lines with different metastatic potential

Background: The metastatic disease rather than the primary tumor itself is responsible for death in most solid tumors, including breast cancer. The role of matrix metalloproteinases ( MMPs), tissue inhibitors of MMPs (TIMPs) and Reversion-inducing cysteine-rich protein with Kazal motifs ( RECK) in the metastatic process has previously been established. However, in all published studies only a limited number of MMPs/MMP inhibitors was analyzed in a limited number of cell lines. Here, we propose a more comprehensive approach by analyzing the expression levels of several MMPs (MMP-2, MMP-9 and MMP-14) and MMP inhibitors (TIMP-1, TIMP-2 and RECK) in different models ( five human breast cancer cell lines, 72 primary breast tumors and 30 adjacent normal tissues). Methods: We analyzed the expression levels of MMP-2, MMP-9 and MMP-14 and their inhibitors (TIMP-1, TIMP-2 and RECK) by quantitative RT-PCR (qRT-PCR) in five human breast cancer cell lines presenting increased invasiveness and metastatic potential, 72 primary breast tumors and 30 adjacent normal tissues. Moreover, the role of cell-extracellular matrix elements interactions in the regulation of expression and activity of MMPs and their inhibitors was analyzed by culturing these cell lines on plastic or on artificial ECM (Matrigel). Results: The results demonstrated that MMPs mRNA expression levels displayed a positive and statistically significant correlation with the transcriptional expression levels of their inhibitors both in the cell line models and in the tumor tissue samples. Furthermore, the expression of all MMP inhibitors was modulated by cell-Matrigel contact only in highly invasive and metastatic cell lines. The enzyme/inhibitor balance at the transcriptional level significantly favors the enzyme which is more evident in tumor than in adjacent non-tumor tissue samples. Conclusion: Our results suggest that the expression of MMPs and their inhibitors, at least at the transcriptional level, might be regulated by common factors and signaling pathways. Therefore, the multi-factorial analysis of these molecules could provide new and independent prognostic information contributing to the determination of more adequate therapy strategies for each patient.`

Chaotic communication based on the particle-in-a-box electronic circuit

A secure communication system based on the error-feedback synchronization of the electronic model of the particle-in-a-box system is proposed. This circuit allows a robust and simple electronic emulation of the mechanical behavior of the collisions of a particle inside a box, exhibiting rich chaotic behavior. The required nonlinearity to emulate the box walls is implemented in a simple way when compared with other analog electronic chaotic circuits. A master/slave synchronization of two circuits exhibiting a rich chaotic behavior demonstrates the potentiality of this system to secure communication. In this system, binary data stream information modulates the bifurcation parameter of the particle-in-a-box electronic circuit in the transmitter. In the receiver circuit, this parameter is estimated using Pecora-Carroll synchronization and error-feedback synchronization. The performance of the demodulation process is verified through the eye pattern technique applied on the recovered bit stream. During the demodulation process, the error-feedback synchronization presented better performance compared with the Pecora-Carroll synchronization. The application of the particle-in-a-box electronic circuit in a secure communication system is demonstrated.