Simulation of the interaction of high-energy C60 cluster ions with amorphous targets

Detailed simulations of the interaction of energetic C-60 beams with amorphous targets are presented here. The spatial evolution of the cluster components is calculated accounting for multiple scattering and Coulomb explosion by means of Monte Carlo and molecular dynamics, respectively. The charge states of the individual cluster components (atoms, atomic ions, fragment cluster ions) as a function of penetration depth are also calculated in tandem with the above calculations by means of the Monte Carlo method. The relative importance of scattering versus Coulomb repulsion is studied as a function of the C-60 cluster energy. The effect of the neighboring cluster constituents on the average charge state of the cluster atoms is calculated as a function of the depth of penetration for a C-60 cluster of 40 MeV. The calculation accounts for the increase in ionization energy of the atom due to the other constituents. Relative track radii are calculated as a function of penetration depth and good agreement with the experimental results is obtained for the interaction of a 30 MeV carbon cluster with silicon. Track splitting observed well into the target as measured by Dunlop in yttrium iron garnet is obtained in the simulations described here for the case of amorphous carbon, provided the Coulomb repulsion is screened by the four valence electrons. Collective energy deposition enhancement is calculated for the 720 MeV cluster. Here the cluster constituents are nearly fully ionized, thereby minimizing the ambiguity related to the value of the ionic charge in the calculation

Seismic modeling using the frozen Gaussian approximation

We adopt the frozen Gaussian approximation (FGA) for modeling seismic waves. The method belongs to the category of ray-based beam methods. It decomposes seismic wavefield into a set of Gaussian functions and propagates these Gaussian functions along appropriate ray paths. As opposed to the classic Gaussian-beam method, FGA keeps the Gaussians frozen (at a fixed width) during the propagation process and adjusts their amplitudes to produce an accurate approximation after summation. We perform the initial decomposition of seismic data using a fast version of the Fourier-Bros-Iagolnitzer (FBI) transform and propagate the frozen Gaussian beams numerically using ray tracing. A test using a smoothed Marmousi model confirms the validity of FGA for accurate modeling of seismic wavefields.Comment: 5 pages, 8 figure