Computer simulations of cosmic-ray diffusion near supernova remnant shock waves

A plasma simulation model was used to study the resonant interactions between streaming cosmic-ray ions and a self-consistent spectrum of Alfven waves, such as might exist in the interstellar medium upstream of a supernova remnant shock wave. The computational model is a hybrid one, in which the background interstellar medium is an MHD fluid and the cosmic-rays are discrete kinetic particles. The particle sources for the electromagnetic fields are obtained by averaging over the fast cyclotron motions. When the perturbed magnetic field is larger than 10 percent of the background field, the macro- and microphysics are no longer correctly predicted by quasi-linear theory. The particles are trapped by the waves and show sharp jumps in their pitch-angles relative to the background magnetic field, and the effective ninety-degree scattering time for diffusion parallel to the background magnetic field is reduced to between 5 and 30 cyclotron periods. Simulation results suggest that Type 1 supernova remnants may be the principal sites of cosmic ray acceleration

Acceleration in perpendicular relativistic shocks for plasmas consisting of leptons and hadrons

We investigate the acceleration of light particles in perpendicular shocks for plasmas consisting of a mixture of leptonic and hadronic particles. Starting from the full set of conservation equations for the mixed plasma constituents, we generalize the magneto-hydrodynamical jump conditions for a multi-component plasma, including information about the specific adiabatic constants for the different species. The impact of deviations from the standard model of an ideal gas is compared in theory and particle-in-cell simulations, showing that the standard-MHD model is a good approximation. The simulations of shocks in electron-positron-ion plasmas are for the first time multi-dimensional, transverse effects are small in this configuration and 1D simulations are a good representation if the initial magnetization is chosen high. 1D runs with a mass ratio of 1836 are performed, which identify the Larmor frequency \omega_{ci} as the dominant frequency that determines the shock physics in mixed component plasmas. The maximum energy in the non-thermal tail of the particle spectra evolves in time according to a power-law proportional to t^\alpha with \alpha in the range 1/3 < \alpha < 1, depending on the initial parameters. A connection is made with transport theoretical models by Drury (1983) and Gargate & Spitkovsky (2011), which predict an acceleration time proportional to \gamma and the theory for small wavelength scattering by Kirk & Reville (2010), which predicts a behavior rather as proportional to \gamma^2. Furthermore, we compare different magnetic field orientations with B_0 inside and out of the plane, observing qualitatively different particle spectra than in pure electron-ion shocks

Transverse quasilinear relaxation in inhomogeneous magnetic field

Transverse quasilinear relaxation of the cyclotron-Cherenkov instability in the inhomogeneous magnetic field of pulsar magnetospheres is considered. We find quasilinear states in which the kinetic cyclotron-Cherenkov instability of a beam propagating through strongly magnetized pair plasma is saturated by the force arising in the inhomogeneous field due to the conservation of the adiabatic invariant. The resulting wave intensities generally have nonpower law frequency dependence, but in a broad frequency range can be well approximated by the power law with the spectral index -2. The emergent spectra and fluxes are consistent with the one observed from pulsars.Comment: 14 Pages, 4 Figure

Understanding and Affecting Student Reasoning About Sound Waves

Student learning of sound waves can be helped through the creation of group-learning classroom materials whose development and design rely on explicit investigations into student understanding. We describe reasoning in terms of sets of resources, i.e. grouped building blocks of thinking that are commonly used in many different settings. Students in our university physics classes often used sets of resources that were different from the ones we wish them to use. By designing curriculum materials that ask students to think about the physics from a different view, we bring about improvement in student understanding of sound waves. Our curriculum modifications are specific to our own classes, but our description of student learning is more generally useful for teachers. We describe how students can use multiple sets of resources in their thinking, and raise questions that should be considered by both instructors and researchers.Comment: 23 pages, 4 figures, 3 tables, 28 references, 7 notes. Accepted for publication in the International Journal of Science Educatio

A Tale of Two Current Sheets

I outline a new model of particle acceleration in the current sheet separating the closed from the open field lines in the force-free model of pulsar magnetospheres, based on reconnection at the light cylinder and "auroral" acceleration occurring in the return current channel that connects the light cylinder to the neutron star surface. I discuss recent studies of Pulsar Wind Nebulae, which find that pair outflow rates in excess of those predicted by existing theories of pair creation occur, and use those results to point out that dissipation of the magnetic field in a pulsar's wind upstream of the termination shock is restored to life as a viable model for the solution of the "$\sigma$" problem as a consequence of the lower wind 4-velocity implied by the larger mass loading.Comment: 17 pages, 6 figures, Invited Review, Proceedings of the "ICREA Workshop on The High-Energy Emission from Pulsars and their Systems", Sant Cugat, Spain, April 12-16, 201

Theory of high-energy emission from the pulsar/Be-star system PSR 1259−-63 I: radiation mechanisms and interaction geometry

We study the physical processes of the PSR B1259-63 system containing a 47 ms pulsar orbiting around a Be star in a highly eccentric orbit. Motivated by the results of a multiwavelength campaign during the January 1994 periastron passage of PSR B1259-63, we discuss several issues regarding the mechanism of high-energy emission. Unpulsed power law emission from the this system was detected near periastron in the energy range 1-200 keV. We find that the observed high energy emission from the PSR B1259-63 system is not compatible with accretion or propeller-powered emission. Shock-powered high-energy emission produced by the pulsar/outflow interaction is consistent with all high energy observations. By studying the evolution of the pulsar cavity we constrain the magnitude and geometry of the mass outflow outflow of the Be star. The pulsar/outflow interaction is most likely mediated by a collisionless shock at the internal boundary of the pulsar cavity. The system shows all the characteristics of a {\it binary plerion} being {\it diffuse} and {\it compact} near apastron and periastron, respectively. The PSR B1259-63 cavity is subject to different radiative regimes depending on whether synchrotron or inverse Compton (IC) cooling dominates the radiation of electron/positron pairs advected away from the inner boundary of the pulsar cavity. The highly non-thermal nature of the observed X-ray/gamma-ray emission near periastron establishes the existence of an efficient particle acceleration mechanism within a timescale shown to be less than $\sim 10^2-10^3$ s. A synchrotron/IC model of emission of e\pm-pairs accelerated at the inner shock front of the pulsar cavity and adiabatically expanding in the MHD flow provides an excellent explanation of the observed time variableX-ray flux and spectrum from the PSRComment: 68 pages, accepted for publication in the Astrophys. J. on Aug. 26, 199

LISA, binary stars, and the mass of the graviton

We extend and improve earlier estimates of the ability of the proposed LISA (Laser Interferometer Space Antenna) gravitational wave detector to place upper bounds on the graviton mass, m_g, by comparing the arrival times of gravitational and electromagnetic signals from binary star systems. We show that the best possible limit on m_g obtainable this way is ~ 50 times better than the current limit set by Solar System measurements. Among currently known, well-understood binaries, 4U1820-30 is the best for this purpose; LISA observations of 4U1820-30 should yield a limit ~ 3-4 times better than the present Solar System bound. AM CVn-type binaries offer the prospect of improving the limit by a factor of 10, if such systems can be better understood by the time of the LISA mission. We briefly discuss the likelihood that radio and optical searches during the next decade will yield binaries that more closely approach the best possible case.Comment: ReVTeX 4, 6 pages, 1 figure, submitted to Phys Rev

Arrival Processes for Vessels in a Port Simulation

Ports provide jetty facilities for ships to load and unload their cargo. Jetty capacity is costly and therefore limited, causing delays for arriving ships. However, ship delays are also costly, so terminal operators attempt to min imize their number and duration. Here, simulation has proved to be a very suitable tool. However, in port simulation models, the impact of the arrival process of ships on the model outcomes tends to be underestimated. This report considers three arrival processes: stock-controlled, equidistant, and uncontrolled. We assess how their deployment in a port simulation model, based on data from a real case study, affects the efficiency of the loading and unloading process, making a case for careful modeling of arrival processes in port simulations. Uncontrolled, which is an assumed arrival process property in many client-oriented simulations, actually performs worst in terms of both ship delays and required storage capacity. Stock-controlled arrivals perform best with regard to large vessel delays and storage capacity. Additional control of the arrival process through the application of a priority scheme in processing ships further impacts efficiency in all three cases.Ports provide jetty facilities for ships to load and unload their cargo. Jetty capacity is costly and therefore limited, causing delays for arriving ships. However, ship delays are also costly, so terminal operators attempt to min imize their number and duration. Here, simulation has proved to be a very suitable tool. However, in port simulation models, the impact of the arrival process of ships on the model outcomes tends to be underestimated. This report considers three arrival processes: stock-controlled, equidistant, and uncontrolled. We assess how their deployment in a port simulation model, based on data from a real case study, affects the efficiency of the loading and unloading process, making a case for careful modeling of arrival processes in port simulations. Uncontrolled, which is an assumed arrival process property in many client-oriented simulations, actually performs worst in terms of both ship delays and required storage capacity. Stock-controlled arrivals perform best with regard to large vessel delays and storage capacity. Additional control of the arrival process through the application of a priority scheme in processing ships further impacts efficiency in all three cases

Statistics of Neutron Stars at the Stage of Supersonic Propeller

We analyze the statistical distribution of neutron stars at the stage of a supersonic propeller. An important point of our analysis is allowance for the evolution of the angle of inclination of the magnetic axis to the spin axis of the neutron star for the boundary of the transition to the supersonic propeller stage for two models: the model with hindered particle escape from the stellar surface and the model with free particle escape. As a result, we have shown that a consistent allowance for the evolution of the inclination angle in the region of extinct radio pulsars for the two models leads to an increase in the total number of neutron stars at the supersonic propeller stage. This increase stems from he fact that when allowing for the evolution of the inclination angle $\chi$ for neutron stars in the region of extinct radio pulsars and, hence, for the boundary of the transition to the propeller stage, this transition is possible at shorter spin periods (P~5-10 s) than assumed in the standard model.Comment: 15 pages, 6 figures; scale corrected for figures 3-