A Multi-dimensional Code for Isothermal Magnetohydrodynamic Flows in Astrophysics

We present a multi-dimensional numerical code to solve isothermal magnetohydrodynamic (IMHD) equations for use in modeling astrophysical flows. First, we have built a one-dimensional code which is based on an explicit finite-difference method on an Eulerian grid, called the total variation diminishing (TVD) scheme. Recipes for building the one-dimensional IMHD code, including the normalized right and left eigenvectors of the IMHD Jacobian matrix, are presented. Then, we have extended the one-dimensional code to a multi-dimensional IMHD code through a Strang-type dimensional splitting. In the multi-dimensional code, an explicit cleaning step has been included to eliminate non-zero $\nabla\cdot B$ at every time step. To estimate the proformance of the code, one- and two-dimensional IMHD shock tube tests, and the decay test of a two-dimensional Alfv\'{e}n wave have been done. As an example of astrophysical applications, we have simulated the nonlinear evolution of the two-dimensional Parker instability under a uniform gravity.Comment: Accepted for publication in ApJ, using aaspp4.sty, 22 text pages with 10 figure

Analysis of Vaccine Effectiveness against COVID-19 and the Emergence of Delta and Other Variants of Concern in Utah

Since the emergence of SARS-CoV-2, vaccines have been heralded as the best way to curtail the COVID-19 pandemic. Clinical trials have shown SARS-CoV-2 vaccines to be highly efficacious against both disease and infection. However, those vaccines currently in use were primarily tested against early lineages. Data on vaccine effectiveness against variants of concern (VOCs) remains limited, including the Delta variant (B.1.617.2)

Antimatter spectra from a time-dependent modeling of supernova remnants

We calculate the energy spectra of cosmic rays (CR) and their secondaries produced in a supernova remnant (SNR), taking into account the time-dependence of the SNR shock. We model the trajectories of charged particles as a random walk with a prescribed diffusioncoefficient, accelerating the particles at each shock crossing. Secondary production by CRs colliding with gas is included as a Monte Carlo process. We find that SNRs produce less antimatter than suggested previously: The positron/electron ratio and the antiproton/proton ratio are a few percent and few $\times 10^{-5}$, respectively. Both ratios do not rise with energy.Comment: 4 pages, 4 eps figures; v2: results for time-dependent magnetic field adde

Co-Occurrence Patterns of Common and Rare Leaf-Litter Frogs, Epiphytic Ferns and Dung Beetles across a Gradient of Human Disturbance

Indicator taxa are commonly used to identify priority areas for conservation or to measure biological responses to environmental change. Despite their widespread use, there is no general consensus about the ability of indicator taxa to predict wider trends in biodiversity. Many studies have focused on large-scale patterns of species co-occurrence to identify areas of high biodiversity, threat or endemism, but there is much less information about patterns of species co-occurrence at local scales. In this study, we assess fine-scale co-occurrence patterns of three indicator taxa (epiphytic ferns, leaf litter frogs and dung beetles) across a remotely sensed gradient of human disturbance in the Ecuadorian Amazon. We measure the relative contribution of rare and common species to patterns of total richness in each taxon and determine the ability of common and rare species to act as surrogate measures of human disturbance and each other. We find that the species richness of indicator taxa changed across the human disturbance gradient but that the response differed among taxa, and between rare and common species. Although we find several patterns of co-occurrence, these patterns differed between common and rare species. Despite showing complex patterns of species co-occurrence, our results suggest that species or taxa can act as reliable indicators of each other but that this relationship must be established and not assumed

The mass and energy budget of Cassiopeia A

Further analysis of X-ray spectroscopy results recently obtained from the MOS CCD cameras on-board XMM-Newton provides a detailed description of the hot and cool X-ray emitting plasma in Cas A. Measurement of the Doppler broadening of the X-ray lines is consistent with the expected ion velocities, ~1500 km/s along the line of sight, in the post shock plasma. Assuming a constant total pressure throughout the remnant we estimate the total remnant mass as 10 Msun and the total thermal energy as 7E43 J. We derive the differential mass distribution as a function of ionisation age for both X-ray emitting components. This distribution is consistent with a hot component dominated by swept up mass heated by the primary shock and a cool component which are ablated clumpy ejecta material which were and are still being heated by interaction with the preheated swept up material. We calculate a balanced mass and energy budget for the supernova explosion giving 1E44 J in ejected mass; approximately 0.4 Msun of the ejecta were diffuse with an initial rms velocity of 15000 km/s while the remaining ~1.8 Msun were clumpy with an initial rms velocity of ~2400 km/s. Using the Doppler velocity measurements of the X-ray spectral lines we can project the mass into spherical coordinates about the remnant. This provides quantitative evidence for mass and energy beaming in the supernova explosion. The mass and energy occupy less than 4.5 sr (<40 % of the available solid angle) around the remnant and 64 % of the mass occurs in two jets within 45 degrees of a jet axis. We calculate a swept up mass of 7.9 Msun in the emitting plasma and estimate that the total mass lost from the progenitor prior to the explosion could be as high as ~20 Msun.Comment: 8 pages, 7 figures, submitted to Astronomy & Astrophysic

The evolution of mass loaded supernova remnants: II. Temperature dependent mass injection rates

We investigate the evolution of spherically symmetric supernova remnants in which mass loading takes place due to conductively driven evaporation of embedded clouds. Numerical simulations reveal significant differences between the evolution of conductively mass loaded and the ablatively mass loaded remnants studied in Paper I. A main difference is the way in which conductive mass loading is extinguished at fairly early times, once the interior temperature of the remnant falls below ~ 107 K. Thus, at late times remnants that ablatively mass load are dominated by loaded mass and thermal energy, while those that conductively mass load are dominated by swept-up mass and kinetic energy. Simple approximations to the remnant evolution, complementary to those in Paper I, are given

A Chandra Observation of Supernova Remnant G350.1-0.3 and Its Central Compact Object

We present a new Chandra observation of supernova remnant (SNR) G350.1-0.3. The high resolution X-ray data reveal previously unresolved filamentary structures and allow us to perform detailed spectroscopy in the diffuse regions of this SNR. Spectral analysis demonstrates that the region of brightest emission is dominated by hot, metal-rich ejecta while the ambient material along the perimeter of the ejecta region and throughout the remnant's western half is mostly low-temperature, shocked interstellar/circumstellar medium (ISM/CSM) with solar-type composition. The data reveal that the emission extends far to the west of the ejecta region and imply a lower limit of 6.6 pc on the diameter of the source (at a distance of 4.5 kpc). We show that G350.1-0.3 is likely in the free expansion (ejecta-dominated) stage and calculate an age of 600-1200 years. The derived relationship between the shock velocity and the electron/proton temperature ratio is found to be entirely consistent with that of other SNRs. We perform spectral fits on the X-ray source XMMU J172054.5-372652, a candidate central compact object (CCO), and find that its spectral properties fall within the typical range of other CCOs. We also present archival 24 um data of G350.1-0.3 taken with the Spitzer Space Telescope during the MIPSGAL galactic survey and find that the infrared and X-ray morphologies are well-correlated. These results help to explain this remnant's peculiar asymmetries and shed new light on its dynamics and evolution

Turbulent Mixing in the Interstellar Medium -- an application for Lagrangian Tracer Particles

We use 3-dimensional numerical simulations of self-gravitating compressible turbulent gas in combination with Lagrangian tracer particles to investigate the mixing process of molecular hydrogen (H2) in interstellar clouds. Tracer particles are used to represent shock-compressed dense gas, which is associated with H2. We deposit tracer particles in regions of density contrast in excess of ten times the mean density. Following their trajectories and using probability distribution functions, we find an upper limit for the mixing timescale of H2, which is of order 0.3 Myr. This is significantly smaller than the lifetime of molecular clouds, which demonstrates the importance of the turbulent mixing of H2 as a preliminary stage to star formation.Comment: 10 pages, 5 figures, conference proceedings "Turbulent Mixing and Beyond 2007

Ionization front-driven turbulence in the clumpy interstellar medium

We present 3D radiation-gasdynamical simulations of an ionization front running into a dense clump. In our setup, a B0 star irradiates an overdensity which is at a distance of 10 pc and modelled as a supercritical 100 M_sol Bonnor-Ebert sphere. The radiation from the star heats up the gas and creates a shock front that expands into the interstellar medium. The shock compresses the clump material while the ionizing radiation heats it up. The outcome of this "cloud-crushing" process is a fully turbulent gas in the wake of the clump. In the end, the clump entirely dissolves. We propose that this mechanism is very efficient in creating short-living supersonic turbulence in the vicinity of massive stars

The Nucleosynthetic Imprint of 15-40 Solar Mass Primordial Supernovae on Metal-Poor Stars

The inclusion of rotationally-induced mixing in stellar evolution can alter the structure and composition of presupernova stars. We survey the effects of progenitor rotation on nucleosynthetic yields in Population III and II supernovae using the new adaptive mesh refinement (AMR) code CASTRO. We examine spherical explosions in 15, 25 and 40 solar mass stars at Z = 0 and 10^-4 solar metallicity with three explosion energies and two rotation rates. Rotation in the Z = 0 models resulted in primary nitrogen production and a stronger hydrogen burning shell which led all models to die as red supergiants. On the other hand, the Z=10^-4 solar metallicity models that included rotation ended their lives as compact blue stars. Because of their extended structure, the hydrodynamics favors more mixing and less fallback in the metal free stars than the Z = 10^-4 models. As expected, higher energy explosions produce more enrichment and less fallback than do lower energy explosions, and less massive stars produce more enrichment and leave behind smaller remnants than do more massive stars. We compare our nucleosynthetic yields to the chemical abundances in the three most iron-poor stars yet found and reproduce the abundance pattern of one, HE 0557-4840, with a zero metallicity 15 solar mass, 2.4 x 10^51 erg supernova. A Salpeter IMF averaged integration of our yields for Z=0 models with explosion energies of 2.4x10^51 ergs or less is in good agreement with the abundances observed in larger samples of extremely metal-poor stars, provided 15 solar mass stars are included. Since the abundance patterns of extremely metal-poor stars likely arise from a representative sample of progenitors, our yields suggest that low-mass supernovae contributed the bulk of the metals to the early universe.Comment: 16 pages, 11 figures; submitted to Ap