Source author record

Shiu-Hang Lee

Shiu-Hang Lee appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

19works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

19 published item(s)

preprint2022arXiv

A Multi-Wavelength View on the Rapidly-Evolving Supernova 2018ivc: An Analog of SN IIb 1993J but Powered Primarily by Circumstellar Interaction

SN 2018ivc is an unusual type II supernova (SN II). It is a variant of SNe IIL, which might represent a transitional case between SNe IIP with a massive H-rich envelope, and IIb with only a small amount of the H-rich envelope. However, SN 2018ivc shows an optical light curve evolution more complicated than canonical SNe IIL. In this paper, we present the results of prompt follow-up observations of SN 2018ivc with the Atacama Large Millimeter/submillimeter Array (ALMA). Its synchrotron emission is similar to that of SN IIb 1993J, suggesting that it is intrinsically an SN IIb-like explosion of a He star with a modest (~0.5 - 1 Msun) extended H-rich envelope. Its radio, optical, and X-ray light curves are explained primarily by the interaction between the SN ejecta and the circumstellar material (CSM); we thus suggest that it is a rare example (and the first involving the `canonical' SN IIb ejecta) for which the multi-wavelength emission is powered mainly by the SN-CSM interaction. The inner CSM density, reflecting the progenitor activity in the final decade, is comparable to that of SN IIb 2013cu that showed a flash spectral feature. The outer CSM density, and therefore the mass-loss rate in the final ~200 years, is larger than that of SN 1993J by a factor of ~5. We suggest that SN 2018ivc represents a missing link between SNe IIP and IIb/Ib/Ic in the binary evolution scenario.

preprint2022arXiv

Long-term evolution of a supernova remnant hosting a double neutron star binary

An ultra-stripped supernova (USSN) is a type of core-collapse SN explosion proposed to be a candidate formation site of a double neutron star (DNS) binary. We investigate the dynamical evolution of an ultra-stripped supernova remnant (USSNR), which should host a DNS at its center. By accounting for the mass-loss history of the progenitor binary using a model developed by a previous study, we construct the large-scale structure of the {circumstellar medium (CSM)} up to a radius $\sim 100\,{\rm pc}$, and simulate the explosion and subsequent evolution of a USSN surrounded by such a CSM environment. We find that the CSM encompasses an extended region characterized by a hot plasma with a temperature $\sim 10^8\,$K located around the termination shock of the wind from the progenitor binary ($\sim 10\,$pc), and the USSNR blastwave is drastically weakened while penetrating through this hot plasma. Radio continuum emission from a young USSNR is sufficiently bright to be detectable if it inhabits our Galaxy but faint compared to the observed Galactic SNRs, and thereafter declines in luminosity through adiabatic cooling. Within our parameter space, USSNRs typically exhibit a low radio luminosity and surface brightness compared to the known Galactic SNRs. Due to the small event rate of USSNe and their relatively short observable lifespan, we calculate that USSNRs account for only $\sim 0.1$-$1$ % of the total SNR population. This is consistent with the fact that no SNR hosting a DNS binary has been discovered in the Milky Way so far.

preprint2022arXiv

Long-term evolution of non-thermal emission from Type Ia and core-collapse supernova remnants in a diversified circumstellar medium

The contribution of galactic supernova remnants (SNRs) to the origin of cosmic rays (CRs) is an important open question in modern astrophysics. Broadband non-thermal emission is a useful proxy for probing the energy budget and production history of CRs in SNRs. We conduct hydrodynamic simulations to model the long-term SNR evolution from explosion all the way to the radiative phase (or $3\times10^5$ yrs at maximum), and compute the time evolution of the broadband non-thermal spectrum to explore its potential applications on constraining the surrounding environments as well as the natures and mass-loss histories of the SNR progenitors. A parametric survey is performed on the ambient environments separated into two main groups, namely a homogeneous medium with a uniform gas density and one with the presence of a circumstellar structure created by the stellar wind of a massive red-supergiant (RSG) progenitor star. Our results reveal a highly diverse evolution history of the non-thermal emission closely correlated to the environmental characteristics of a SNR. Up to the radiative phase, the roles of CR re-acceleration and ion-neutral wave damping on the spectral evolution are investigated. Finally, we make an assessment of the future prospect of SNR observations by the next-generation hard X-ray space observatory FORCE and predict what we can learn from their comparison with our evolution models.

preprint2021arXiv

Rapid Deceleration of Blast Waves Witnessed in Tycho's Supernova Remnant

In spite of their importance as standard candles in cosmology and as major major sites of nucleosynthesis in the Universe, what kinds of progenitor systems lead to type Ia supernovae (SN) remains a subject of considerable debate in the literature. This is true even for the case of Tycho's SN exploded in 1572 although it has been deeply studied both observationally and theoretically. Analyzing X-ray data of Tycho's supernova remnant (SNR) obtained with Chandra in 2003, 2007, 2009, and 2015, we discover that the expansion before 2007 was substantially faster than radio measurements reported in the past decades and then rapidly decelerated during the last ~ 15 years. The result is well explained if the shock waves recently hit a wall of dense gas surrounding the SNR. Such a gas structure is in fact expected in the so-called single-degenerate scenario, in which the progenitor is a binary system consisting of a white dwarf and a stellar companion, whereas it is not generally predicted by a competing scenario, the double-degenerate scenario, which has a binary of two white dwarfs as the progenitor. Our result thus favors the former scenario. This work also demonstrates a novel technique to probe gas environments surrounding SNRs and thus disentangle the two progenitor scenarios for Type Ia SNe.

preprint2021arXiv

Resurrection of Non-thermal Emissions from Type Ib/c Supernova Remnants

Supernova remnants (SNRs) are important objects in investigating the links among supernova (SN) explosion mechanism(s), progenitor stars, and cosmic-ray acceleration. Non-thermal emission from SNRs is an effective and promising tool for probing their surrounding circumstellar media (CSM) and, in turn, the stellar evolution and mass-loss mechanism(s) of massive stars. In this work, we calculate the time evolution of broadband non-thermal emissions from Type Ib/c SNRs whose CSM structures are derived from the mass-loss history of their progenitors. Our results predict that Type Ib/c SNRs make a transition of brightness in radio and $γ$-ray bands from an undetectable dark for a certain period to a re-brightening phase. This transition originates from their inhomogeneous CSM structures in which the SNRs are embedded within a low-density wind cavity surrounded by a high-density wind shell and the ambient interstellar medium (ISM). The "resurrection" in non-thermal luminosity happens at an age of ~1,000 yrs old for a Wolf-Rayet star progenitor evolved within a typical ISM density. Combining with the results of Type II SNR evolution recently reported by Yasuda et al. (2021), this result sheds light on a comprehensive understanding of non-thermal emissions from SNRs with different SN progenitor types and ages, which is made possible for the first time by the incorporation of realistic mass-loss histories of the progenitors.

preprint2020arXiv

Evidence of a Type Ia Progenitor for Supernova Remnant 3C 397

The explosive origin of the young supernova remnant (SNR) 3C 397 (G41.1-0.3) is debated. Its elongated morphology and proximity to a molecular cloud are suggestive of a core-collapse (CC) SN origin, yet recent X-ray studies of heavy metals show chemical yields and line centroid energies consistent with a Type Ia SN. In this paper, we analyze the full X-ray spectrum from 0.7-10 keV of 3C 397 observed with Suzaku and compare the line centroid energies, fluxes, and elemental abundances of intermediate-mass and heavy metals (Mg to Ni) to Type Ia and CC hydrodynamical model predictions. Based on the results, we conclude that 3C 397 likely arises from an energetic Type Ia explosion in a high-density ambient medium, and we show that the progenitor was a near Chandrasekhar mass white dwarf.

preprint2019arXiv

Radio Emission from Supernovae in the Very Early Phase: Implications for the Dynamical Mass Loss of Massive Stars

Recent high-cadence transient surveys and rapid follow-up observations indicate that some massive stars may dynamically lose their own mass within decades before supernovae (SNe). Such a mass-loss forms `confined' circumstellar medium (CSM); a high density material distributed only within a small radius ($\lesssim 10^{15}$ cm with the mass-loss rate of 0.01 $\sim 10^{-4} M_\odot$ yr$^{-1}$). While the SN shock should trigger particle acceleration and magnetic field amplification in the `confined' CSM, synchrotron emission may be masked in centimeter wavelengths due to free-free absorption; the millimeter range can however be a potential new window. We investigate the time evolution of synchrotron radiation from the system of a red super giant surrounded by the `confined' CSM, relevant to typical type II-P SNe. We show that synchrotron millimeter emission is generally detectable, and the signal can be used as a sensitive tracer of the nature of the `confined' CSM; it traces different CSM density parameter space than in the optical. Furthermore, our simulations show that the `confined' CSM efficiently produces secondary electrons and positrons through proton inelastic collisions, which can become main contributors to the synchrotron emission in several ten days since the SN. We predict that the synchrotron emission is detectable by ALMA, and suggest that it will provide a robust evidence of the existence of the `confined' CSM.

preprint2016arXiv

The Quiescent Intracluster Medium in the Core of the Perseus Cluster

Clusters of galaxies are the most massive gravitationally-bound objects in the Universe and are still forming. They are thus important probes of cosmological parameters and a host of astrophysical processes. Knowledge of the dynamics of the pervasive hot gas, which dominates in mass over stars in a cluster, is a crucial missing ingredient. It can enable new insights into mechanical energy injection by the central supermassive black hole and the use of hydrostatic equilibrium for the determination of cluster masses. X-rays from the core of the Perseus cluster are emitted by the 50 million K diffuse hot plasma filling its gravitational potential well. The Active Galactic Nucleus of the central galaxy NGC1275 is pumping jetted energy into the surrounding intracluster medium, creating buoyant bubbles filled with relativistic plasma. These likely induce motions in the intracluster medium and heat the inner gas preventing runaway radiative cooling; a process known as Active Galactic Nucleus Feedback. Here we report on Hitomi X-ray observations of the Perseus cluster core, which reveal a remarkably quiescent atmosphere where the gas has a line-of-sight velocity dispersion of 164+/-10 km/s in a region 30-60 kpc from the central nucleus. A gradient in the line-of-sight velocity of 150+/-70 km/s is found across the 60 kpc image of the cluster core. Turbulent pressure support in the gas is 4% or less of the thermodynamic pressure, with large scale shear at most doubling that estimate. We infer that total cluster masses determined from hydrostatic equilibrium in the central regions need little correction for turbulent pressure.

preprint2015arXiv

Are Models for Core-Collapse Supernova Progenitors Consistent with the Properties of Supernova Remnants?

The recent discovery that the Fe-K line luminosities and energy centroids observed in nearby SNRs are a strong discriminant of both progenitor type and circumstellar environment has implications for our understanding of supernova progenitor evolution. Using models for the chemical composition of core-collapse supernova ejecta, we model the dynamics and thermal X-ray emission from shocked ejecta and circumstellar material, modeled as an $r^{-2}$ wind, to ages of 3000 years. We compare the X-ray spectra expected from these models to observations made with the Suzaku satellite. We also model the dynamics and X-ray emission from Type Ia progenitor models. We find a clear distinction in Fe-K line energy centroid between core-collapse and Type Ia models. The core-collapse supernova models predict higher Fe-K line centroid energies than the Type Ia models, in agreement with observations. We argue that the higher line centroids are a consequence of the increased densities found in the circumstellar environment created by the expansion of the slow-moving wind from the massive progenitors.

preprint2015arXiv

Electron and Ion Acceleration in Relativistic Shocks with Applications to GRB Afterglows

We have modeled the simultaneous first-order Fermi shock acceleration of protons, electrons, and helium nuclei by relativistic shocks. By parameterizing the particle diffusion, our steady-state Monte Carlo simulation allows us to follow particles from particle injection at nonthermal thermal energies to above PeV energies, including the nonlinear smoothing of the shock structure due to cosmic-ray (CR) backpressure. We observe the mass-to-charge (A/Z) enhancement effect believed to occur in efficient Fermi acceleration in non-relativistic shocks and we parameterize the transfer of ion energy to electrons seen in particle-in-cell (PIC) simulations. For a given set of environmental and model parameters, the Monte Carlo simulation determines the absolute normalization of the particle distributions and the resulting synchrotron, inverse-Compton, and pion-decay emission in a largely self-consistent manner. The simulation is flexible and can be readily used with a wide range of parameters typical of gamma-ray burst (GRB) afterglows. We describe some preliminary results for photon emission from shocks of different Lorentz factors and outline how the Monte Carlo simulation can be generalized and coupled to hydrodynamic simulations of GRB blast waves. We assume Bohm diffusion for simplicity but emphasize that the nonlinear effects we describe stem mainly from an extended shock precursor where higher energy particles diffuse further upstream. Quantitative differences will occur with different diffusion models, particularly for the maximum CR energy and photon emission, but these nonlinear effects should be qualitatively similar as long as the scattering mean free path is an increasing function of momentum.

preprint2015arXiv

Evidence for Thermal X-Ray Line Emission from the Synchrotron-Dominated Supernova Remnant RX J1713.7-3946

We report the first detection of thermal X-ray line emission from the supernova remnant (SNR) RX J1713.7-3946, the prototype of the small class of synchrotron dominated SNRs. A softness-ratio map generated using XMM-Newton data shows that faint interior regions are softer than bright shell regions. Using Suzaku and deep XMM-Newton observations, we have extracted X-ray spectra from the softest area, finding clear line features at 1 keV and 1.35 keV. These lines can be best explained as Ne Ly-alpha and Mg He-alpha from a thermal emission component. Since the abundance ratios of metals to Fe are much higher than solar values in the thermal component, we attribute the thermal emission to reverse-shocked SN ejecta. The measured Mg/Ne, Si/Ne, and Fe/Ne ratios of 2.0-2.6, 1.5-2.0, and <0.05 solar suggest that the progenitor star of RX J1713.7-3946 was a relatively low-mass star (<~20 M_sun), consistent with a previous inference based on the effect of stellar winds of the progenitor star on the surrounding medium. Since the mean blastwave speed of ~6000 km/s (the radius of 9.6 pc divided by the age of 1600 yr) is relatively fast compared with other core-collapse SNRs, we propose that RX J1713.7-3946 is a result of a Type Ib/c supernova whose progenitor was a member of an interacting binary. While our analysis provides strong evidence for X-ray line emission, our interpretation of its nature as thermal emission from SN ejecta requires further confirmation especially through future precision spectroscopic measurements using ASTRO-H.

preprint2015arXiv

Modeling Bright Gamma-ray and Radio Emission at Fast Cloud Shocks

Recent observations by the Large Area Telescope (LAT) onboard the Fermi satellite have revealed bright gamma-ray emission from middle-aged supernova remnants (SNRs) inside our Galaxy. These remnants, which also possess bright non-thermal radio shells, are often found to be interacting directly with surrounding gas clouds. We explore the non-thermal emission mechanism at these dynamically evolved SNRs by constructing a hydrodynamical model. Two scenarios of particle acceleration, either a re-acceleration of Galactic cosmic rays (CRs) or an efficient nonlinear diffusive shock acceleration (NLDSA) of particles injected from downstream, are considered. Using parameters inferred from observations, our models are contrasted with the observed spectra of SNR W44. For the re-acceleration case, we predict a significant enhancement of radio and GeV emission as the SNR undergoes a transition into the radiative phase. If sufficiently strong magnetic turbulence is present in the molecular cloud, the re-acceleration scenario can explain the observed broadband spectral properties. The NLDSA scenario also succeeds in explaining the $γ$-ray spectrum but fails to reproduce the radio spectral index. Efficient NLDSA also results in a significant post-shock non-thermal pressure that limits the compression during cooling and prevents the formation of a prominent dense shell. Some other interesting differences between the two models in hydrodynamical behavior and resulting spectral features are illustrated.

preprint2015arXiv

Suzaku spectra of a Type II Supernova Remnant Kes 79

This paper reports results of a Suzaku observation of the supernova remnant (SNR) Kes 79 (G33.6+0.1). The X-ray spectrum is best fitted by a two-temperature model: a non-equilibrium ionization (NEI) plasma and a collisional ionization equilibrium (CIE) plasma. The NEI plasma is spatially confined within the inner radio shell with kT~0.8 keV, while the CIE plasma is found in more spatially extended regions associated with the outer radio shell with kT~0.2 keV and solar abundance. Therefore, the NEI plasma is attributable to the SN ejecta and the CIE plasma is forward shocked interstellar medium. In the NEI plasma, we discovered K-shell line of Al, Ar and Ca for the first time. The abundance pattern and estimated mass of the ejecta are consistent with the core-collapse supernova explosion of a ~30-40 solar mass progenitor star. An Fe line with center energy of ~6.4 keV is also found in the southeast (SE) portion of the SNR, a close peripheral region around dense molecular clouds. One possibility is that the line is associated with the ejecta. However, the centroid energy of ~6.4 keV and the spatial distribution of enhancement near the SE peripheral do not favor this scenario. Since the ~6.4 keV emitting region coincides to the molecular clouds, we propose another possibility that the Fe line is due to K-shell ionization of neutral Fe by the interaction of locally accelerated protons (LECRp) with the surrounding molecular cloud. Both these possibilities, heated ejecta or LECRp origin, are discussed based on the observational facts.

preprint2014arXiv

Reverse and Forward Shock X-ray Emission in an Evolutionary Model of Supernova Remnants undergoing Efficient Diffusive Shock Acceleration

We present new models for the forward and reverse shock thermal X-ray emission from core-collapse and Type Ia supernova remnants (SNRs) which include the efficient production of cosmic rays via non-linear diffusive shock acceleration (DSA). Our CR-hydro-NEI code takes into account non-equilibrium ionization (NEI), hydrodynamic effects of efficient CR production on the SNR evolution, and collisional temperature equilibration among heavy ions and electrons in both the shocked supernova (SN) ejecta and the shocked circumstellar material. While X-ray emission is emphasized here, our code self-consistently determines both thermal and non-thermal broadband emission from radio to TeV energies. We include Doppler broadening of the spectral lines by thermal motions of the ions and by the remnant expansion. We study, in general terms, the roles which the ambient environment, progenitor models, temperature equilibration, and processes related to DSA have on the thermal and non-thermal spectra. The study of X-ray line emission from young SNRs is a powerful tool for determining specific SN elemental contributions, and for providing critical information that helps to understand the type and energetics of the explosion, the composition of the ambient medium in which the SN exploded, and the ionization and dynamics of the hot plasma in the shocked SN ejecta and interstellar medium. With the approaching launch of the next-generation X-ray satellite Astro-H, observations of spectral lines with unprecedented high resolution will become a reality. Our self-consistent calculations of the X-ray spectra from various progenitors will help interpret future observations of SNRs.

preprint2014arXiv

Spectral and Polarization Properties of Photospheric Emission From Stratified Jets

We explore the spectral and polarization properties of photospheric emissions from stratified jets in which multiple components, separated by a sharp velocity shear regions, are distributed in lateral direction. Propagation of thermal photons injected at high optical depth region are calculated until they escape from the photosphere. It is found that presence of the lateral structure within the jet leads to non-thermal feature of the spectra and significant polarization signal in the resulting emission. The deviation from thermal spectra as well as the polarization degree tends to be enhanced as the velocity gradient in the shear region increases. In particular, we show that emissions from multi-component jet can reproduce the typical observed spectra of gamma-ray bursts (GRBs) irrespective to the position of the observer when a velocity shear region is closely spaced in various lateral ($θ$) positions. The degree of polarization associated in the emission is significant (> few%) at wide range of observer angles and can be higher than 30%.

preprint2013arXiv

A CR-hydro-NEI Model of Multi-wavelength Emission from the Vela Jr. Supernova Remnant (SNR RX J0852.0-4622)

Based largely on energy budget considerations and the observed cosmic-ray (CR) ionic composition, supernova remnant (SNR) blast waves are the most likely sources of CR ions with energies at least up to the "knee" near 3 PeV. Shocks in young shell-type TeV-bright SNRs are surely producing TeV particles, but the emission could be dominated by ions producing neutral pion-decay emission or electrons producing inverse-Compton gamma-rays. Unambiguously identifying the GeV-TeV emission process in a particular SNR will not only help pin down the origin of CRs, it will add significantly to our understanding of the diffusive shock acceleration (DSA) mechanism and improve our understanding of supernovae and the impact SNRs have on the circumstellar medium. In this study, we investigate the Vela Jr. SNR, an example of TeV-bright non-thermal SNRs. We perform hydrodynamic simulations coupled with non-linear DSA and non-equilibrium ionization near the forward shock (FS) to confront currently available multi-wavelength data. We find, with an analysis similar to that used earlier for SNR RX J1713.7-3946, that self-consistently modeling the thermal X-ray line emission with the non-thermal continuum in our one-dimensional model strongly constrains the fitting parameters, and this leads convincingly to a leptonic origin for the GeV-TeV emission for Vela Jr. This conclusion is further supported by applying additional constraints from observation, including the radial brightness profiles of the SNR shell in TeV gamma-rays, and the spatial variation of the X-ray synchrotron spectral index. We will discuss implications of our models on future observations by the next-generation telescopes.

preprint2013arXiv

Matter mixing in aspherical core-collapse supernovae: a search for possible conditions for conveying $^{56}$Ni into high velocity regions

We perform two-dimensional axisymmetric hydrodynamic simulations of matter mixing in aspherical core-collapse supernova explosions of a 16.3 $M_{\odot}$ star with a compact hydrogen envelope. Observations of SN 1987A have provided evidence that $^{56}$Ni synthesized by explosive nucleosynthesis is mixed into fast moving matter ($\gtrsim$ 3,500 km s$^{-1}$) in the exploding star. In order to clarify the key conditions for reproducing such high velocity of $^{56}$Ni, we revisit matter mixing in aspherical core-collapse supernova explosions. Explosions are initiated artificially by injecting thermal and kinetic energies around the interface between the iron core and the silicon-rich layer. Perturbations of 5% or 30% amplitude in the radial velocities are introduced at several points in time. We found that no high velocity $^{56}$Ni can be obtained if we consider bipolar explosions with perturbations (5% amplitude) of pre-supernova origins. If large perturbations (30% amplitude) are introduced or exist due to some unknown mechanism in a later phase just before the shock wave reaches the hydrogen envelope, $^{56}$Ni with a velocity of 3,000 km s$^{-1}$ can be obtained. Aspherical explosions that are asymmetric across the equatorial plane with clumpy structures in the initial shock waves are investigated. We found that the clump sizes affect the penetration of $^{56}$Ni. Finally, we report that an aspherical explosion model that is asymmetric across the equatorial plane with multiple perturbations of pre-supernova origins can cause the penetration of $^{56}$Ni clumps into fast moving matter of 3,000 km s$^{-1}$. We show that both aspherical explosion with clumpy structures and perturbations of pre-supernova origins may be necessary to reproduce the observed high velocity of $^{56}$Ni. To confirm this, more robust three-dimensional simulations are required.

preprint2013arXiv

Photospheric emission from stratified jets

We explore photospheric emissions from stratified two-component jets, wherein a highly relativistic spine outflow is surrounded by a wider and less relativistic sheath outflow. Thermal photons are injected in regions of high optical depth and propagated until they escape at the photosphere. Due to the presence of shear in velocity (Lorentz factor) at the boundary of the spine and sheath region, a fraction of the injected photons are accelerated via a Fermi-like acceleration mechanism such that a high energy power-law tail is formed in the resultant spectrum. We show, in particular, that if a velocity shear with a considerable variance in the bulk Lorentz factor is present, the high energy part of observed Gamma-ray Bursts (GRBs) photon spectrum can be explained by this photon acceleration mechanism. We also show that the accelerated photons may also account for the origin of the extra hard power-law component above the bump of the thermal-like peak seen in some peculiar bursts (e.g., GRB 090510, 090902B, 090926A). It is demonstrated that time-integrated spectra can also reproduce the low energy spectrum of GRBs consistently due to a multi-temperature effect when time evolution of the outflow is considered. Finally, we show that the empirical Ep-Lp relation can be explained by differences in the outflow properties of individual sources.

preprint2012arXiv

A Generalized Model of Nonlinear Diffusive Shock Acceleration Coupled to an Evolving Supernova Remnant

To better model the efficient production of cosmic rays (CRs) in supernova remnants (SNRs) with the associated coupling between CR production and SNR dynamics, we have generalized an existing cr-hydro-NEI code (i.e., Ellison et al. 2012) to include the following processes: (1) an explicit calculation of the upstream precursor structure including the position dependent flow speed, density, temperature, and magnetic field strength; (2) a momentum and space dependent CR diffusion coefficient; (3) an explicit calculation of magnetic field amplification (MFA); (4) calculation of the maximum CR momentum using the amplified magnetic field; (5) a finite Alfven speed for the particle scattering centers; and (6) the ability to accelerate a superthermal seed population of CRs as well as the ambient thermal plasma. While a great deal of work has been done modeling SNRs, most work has concentrated on either the continuum emission from relativistic electrons or ions, or the thermal emission from the shock heated plasma. Our generalized code combines these elements and describes the interplay between CR production and SNR evolution, including the nonlinear coupling of efficient diffusive shock acceleration (DSA), based mainly on the work of P. Blasi and co-workers, and a non-equilibrium ionization (NEI) calculation of thermal X-ray line emission. We believe our generalized model will provide a consistent modeling platform for SNRs, including those interacting with molecular clouds, and improve the interpretation of current and future observations, including the high-quality spectra expected from Astro-H. SNR RX J1713.7-3946 is modeled as an example.