Transport of intense ion beams in plasmas: collimation and energy-loss reduction

We compare the transport properties of a well-characterized hydrogen plasma for low and high current ion beams. The energy-loss of low current beams can be well understood, within the framework of current stopping power models. However, for high current proton beams, significant energy-loss reduction and collimation is observed in the experiment. We have developed a new particle-in-cell code, which includes both collective electromagnetic effects and collisional interactions. Our simulations indicate that resistive magnetic fields, induced by the transport of an intense proton beam, act to collimate the proton beam and simultaneously deplete the local plasma density along the beam path. This in turn causes the energy-loss reduction detected in the experiment

Anomalous stopping of laser-accelerated intense proton beam in dense ionized matter

Ultrahigh-intensity lasers (10$^{18}$-10$^{22}$W/cm$^{2}$) have opened up new perspectives in many fields of research and application [1-5]. By irradiating a thin foil, an ultrahigh accelerating field (10$^{12}$ V/m) can be formed and multi-MeV ions with unprecedentedly high intensity (10$^{10}$A/cm$^2$) in short time scale ($\sim$ps) are produced [6-14]. Such beams provide new options in radiography [15], high-yield neutron sources [16], high-energy-density-matter generation [17], and ion fast ignition [18,19]. An accurate understanding of the nonlinear behavior of beam transport in matter is crucial for all these applications. We report here the first experimental evidence of anomalous stopping of a laser-generated high-current proton beam in well-characterized dense ionized matter. The observed stopping power is one order of magnitude higher than single-particle slowing-down theory predictions. We attribute this phenomenon to collective effects where the intense beam drives an decelerating electric field approaching 1GV/m in the dense ionized matter. This finding will have considerable impact on the future path to inertial fusion energy.Comment: 8 pages, 4 figure

Energy loss enhancement of very intense proton beams in dense matter due to the beam-density effect

Thoroughly understanding the transport and energy loss of intense ion beams in dense matter is essential for high-energy-density physics and inertial confinement fusion. Here, we report a stopping power experiment with a high-intensity laser-driven proton beam in cold, dense matter. The measured energy loss is one order of magnitude higher than the expectation of individual particle stopping models. We attribute this finding to the proximity of beam ions to each other, which is usually insignificant for relatively-low-current beams from classical accelerators. The ionization of the cold target by the intense ion beam is important for the stopping power calculation and has been considered using proper ionization cross section data. Final theoretical values agree well with the experimental results. Additionally, we extend the stopping power calculation for intense ion beams to plasma scenario based on Ohm's law. Both the proximity- and the Ohmic effect can enhance the energy loss of intense beams in dense matter, which are also summarized as the beam-density effect. This finding is useful for the stopping power estimation of intense beams and significant to fast ignition fusion driven by intense ion beams

Target density effects on charge tansfer of laser-accelerated carbon ions in dense plasma

We report on charge state measurements of laser-accelerated carbon ions in the energy range of several MeV penetrating a dense partially ionized plasma. The plasma was generated by irradiation of a foam target with laser-induced hohlraum radiation in the soft X-ray regime. We used the tri-cellulose acetate (C$_{9}$H$_{16}$O$_{8}$) foam of 2 mg/cm$^{-3}$ density, and $1$-mm interaction length as target material. This kind of plasma is advantageous for high-precision measurements, due to good uniformity and long lifetime compared to the ion pulse length and the interaction duration. The plasma parameters were diagnosed to be T$_{e}$=17 eV and n$_{e}$=4 $\times$ 10$^{20}$ cm$^{-3}$. The average charge states passing through the plasma were observed to be higher than those predicted by the commonly-used semiempirical formula. Through solving the rate equations, we attribute the enhancement to the target density effects which will increase the ionization rates on one hand and reduce the electron capture rates on the other hand. In previsous measurement with partially ionized plasma from gas discharge and z-pinch to laser direct irradiation, no target density effects were ever demonstrated. For the first time, we were able to experimentally prove that target density effects start to play a significant role in plasma near the critical density of Nd-Glass laser radiation. The finding is important for heavy ion beam driven high energy density physics and fast ignitions.Comment: 7 pages, 4 figures, 35 conference

Binding Mechanism of Inhibitors to BRD4 and BRD9 Decoded by Multiple Independent Molecular Dynamics Simulations and Deep Learning

Bromodomain 4 and 9 (BRD4 and BRD9) have been regarded as important targets of drug designs in regard to the treatment of multiple diseases. In our current study, molecular dynamics (MD) simulations, deep learning (DL) and binding free energy calculations are integrated to probe the binding modes of three inhibitors (H1B, JQ1 and TVU) to BRD4 and BRD9. The MD trajectory-based DL successfully identify significant functional function domains, such as BC-loop and ZA-loop. The information from the post-processing analysis of MD simulations indicates that inhibitor binding highly influences the structural flexibility and dynamic behavior of BRD4 and BRD9. The results of the MM-GBSA calculations not only suggest that the binding ability of H1B, JQ1 and TVU to BRD9 are stronger than to BRD4, but they also verify that van der Walls interactions are the primary forces responsible for inhibitor binding. The hot spots of BRD4 and BRD9 revealed by residue-based free energy estimation provide target sites of drug design in regard to BRD4 and BRD9. This work is anticipated to provide useful theoretical aids for the development of selective inhibitors over BRD family members

Identification Mechanism of BACE1 on Inhibitors Probed by Using Multiple Separate Molecular Dynamics Simulations and Comparative Calculations of Binding Free Energies

β-amyloid cleaving enzyme 1 (BACE1) is regarded as an important target of drug design toward the treatment of Alzheimer’s disease (AD). In this study, three separate molecular dynamics (MD) simulations and calculations of binding free energies were carried out to comparatively determine the identification mechanism of BACE1 for three inhibitors, 60W, 954 and 60X. The analyses of MD trajectories indicated that the presence of three inhibitors influences the structural stability, flexibility and internal dynamics of BACE1. Binding free energies calculated by using solvated interaction energy (SIE) and molecular mechanics generalized Born surface area (MM-GBSA) methods reveal that the hydrophobic interactions provide decisive forces for inhibitor–BACE1 binding. The calculations of residue-based free energy decomposition suggest that the sidechains of residues L91, D93, S96, V130, Q134, W137, F169 and I179 play key roles in inhibitor–BACE1 binding, which provides a direction for future drug design toward the treatment of AD

Binding Mechanism of Inhibitors to Heat Shock Protein 90 Investigated by Multiple Independent Molecular Dynamics Simulations and Prediction of Binding Free Energy

The heat shock protein (HSP90) has been an import target of drug design in the treatment of human disease. An exploration of the conformational changes in HSP90 can provide useful information for the development of efficient inhibitors targeting HSP90. In this work, multiple independent all-atom molecular dynamics (AAMD) simulations followed by calculations of the molecular mechanics generalized Born surface area (MM-GBSA) were performed to explore the binding mechanism of three inhibitors (W8Y, W8V, and W8S) to HSP90. The dynamics analyses verified that the presence of inhibitors impacts the structural flexibility, correlated movements, and dynamics behavior of HSP90. The results of the MM-GBSA calculations suggest that the selection of GB models and empirical parameters has important influences on the predicted results and verify that van der Waals interactions are the main forces that determine inhibitor–HSP90 binding. The contributions of separate residues to the inhibitor–HSP90 binding process indicate that hydrogen-bonding interactions (HBIs) and hydrophobic interactions play important roles in HSP90–inhibitor identifications. Moreover, residues L34, N37, D40, A41, D79, I82, G83, M84, F124, and T171 are recognized as hot spots of inhibitor–HSP90 binding and provide significant target sites of for the design of drugs related to HSP90. This study aims to contribute to the development of efficient inhibitors that target HSP90 by providing an energy-based and theoretical foundation

DataSheet1_Exploring the deactivation mechanism of human β2 adrenergic receptor by accelerated molecular dynamic simulations.PDF

The β2 adrenergic receptor (β2AR), one of important members of the G protein coupled receptors (GPCRs), has been suggested as an important target for cardiac and asthma drugs. Two replicas of Gaussian accelerated molecular dynamics (GaMD) simulations are performed to explore the deactivation mechanism of the active β2AR bound by three different substrates, including the agonist (P0G), antagonist (JTZ) and inverse agonist (JRZ). The simulation results indicate that the Gs protein is needed to stabilize the active state of the β2AR. Without the Gs protein, the receptor could transit from the active state toward the inactive state. During the transition process, helix TM6 moves toward TM3 and TM5 in geometric space and TM5 shrinks upwards. The intermediate state is captured during the transition process of the active β2AR toward the inactive one, moreover the changes in hydrophobic interaction networks between helixes TM3, TM5, and TM6 and the formation of a salt bridge between residues Arg3.50 and Glu6.30 drive the transition process. We expect that this finding can provide energetic basis and molecular mechanism for further understanding the function and target roles of the β2AR.</p