Source author record

T. W. Gorczyca

T. W. Gorczyca 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

11works
8topics
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

11 published item(s)

preprint2020arXiv

PSTGF : time-independent R-Matrix atomic electron-impact code

STGF is a community code employed for outer-region R-matrix calculations, describing electron-impact collisional processes. It is widely recognised that the original version of STGF was written by M. J. Seaton in 1983, but through constant refinement over the next decades by worldwide contributors has evolved into its current form that more reflects modern coding practice and current computer architectures. Despite its current wide acceptance, it was never formally published. Therefore, we present an updated high-performance parallel version of PSTGF, that balances the requirements of small university clusters, yet can exploit the computational power of cutting edge supercomputers. There are many improvements over the original STGF, but most noticeably, the full introduction of MQDT options that provide subsequent integration with ICFT (Intermediate Coupling Frame Transformation) codes, and for either Breit-Pauli / DARC (Dirac Atomic R-matrix Codes), better load balancing, high levels of vectorisation and simplified output. Semantically, the program is full fortran 90 in conjunction with MPI (Message Passing Interface) though has CUDA fortran options for the most numerically intensive code sections.

preprint2020arXiv

Silicon X-ray absorption in the Galactic ISM: the gaseous component

We present a detailed analysis of the gaseous component of the Si K edge using high-resolution Chandra spectra of low-mass X-ray binaries. We fit the spectra with amodified version of the ISMabs model, including new photoabsorption cross sectionscomputed for all Si ionic species. We estimate column densities for Si i, Si ii, Si iii, Si xii and Si xiii, which trace the warm, intermediate temperature and hot phases of the Galactic interstellar medium. We find that the ionic fractions of the first two phases are similar. This may be due to the physical state of the plasma determined by the temperature or to the presence of absorber material in the close vicinity of the sources. Our findings highlight the need for accurate modeling of the gaseous component before attempting to address the solid component.

preprint2015arXiv

Bootstrapping dielectronic recombination from second-row elements and the Orion Nebula

Dielectronic recombination (DR) is the dominant recombination process for most heavy elements in photoionized clouds. Accurate DR rates for a species can be predicted when the positions of autoionizing states are known. Unfortunately such data are not available for most third and higher-row elements. This introduces an uncertainty that is especially acute for photoionized clouds, where the low temperatures mean that DR occurs energetically through very low-lying autoionizing states. This paper discusses S$^{2+} \rightarrow$ S$^+$ DR, the process that is largely responsible for establishing the [S~III]/[S~II] ratio in nebulae. We derive an empirical rate coefficient using a novel method for second-row ions, which do have accurate data. Photoionization models are used to reproduce the [O~III] / [O~II] / [O~I] / [Ne~III] intensity ratios in central regions of the Orion Nebula. O and Ne have accurate atomic data and can be used to derive an empirical S$^{2+} \rightarrow$ S$^+$ DR rate coefficient at $\sim 10^{4}$~K. We present new calculations of the DR rate coefficient for S$^{2+} \rightarrow$ S$^+$ and quantify how uncertainties in the autoionizing level positions affect it. The empirical and theoretical results are combined and we derive a simple fit to the resulting rate coefficient at all temperatures for incorporation into spectral synthesis codes. This method can be used to derive empirical DR rates for other ions, provided that good observations of several stages of ionization of O and Ne are available.

preprint2014arXiv

ISMabs: a comprehensive X-ray absorption model for the interstellar medium

We present an X-ray absorption model for the interstellar medium, to be referred to as ISMabs, that takes into account both neutral and ionized species of cosmically abundant elements, and includes the most accurate atomic data available. Using high-resolution spectra from eight X-ray binaries obtained with the Chandra High Energy Transmission Grating Spectrometer, we proceed to benchmark the atomic data in the model particularly in the neon K-edge region. Compared with previous photoabsorption models, which solely rely on neutral species, the inclusion of ions leads to improved spectral fits. Fit parameters comprise the column densities of abundant contributors that allow direct estimates of the ionization states. ISMabs is provided in the appropriate format to be implemented in widely used X-ray spectral fitting packages such as XSPEC, ISIS and SHERPA.

preprint2014arXiv

K-Shell Photoabsorption of Magnesium Ions

X-ray photoabsorption cross sections have been computed for all magnesium ions using the $R$-matrix method. A comparison with the other available data for Mg II-Mg X shows good qualitative agreement in the resultant resonance shapes. However, for the lower ionization stages, and for singly-ionized Mg II in particular, the previous $R$-matrix results (Witthoeft et al.2009; Witthoeft et al. 2011) overestimate the K-edge position due to the neglect of important orbital relaxation effects, and a global shift downward in photon energy of those cross sections is therefore warranted. We have found that the cross sections for Mg I and Mg II are further complicated by the M-shell ($n=3$) occupancy. As a result, the treatment of spectator Auger decay of $1s\rightarrow np$ resonances using a method based on multichannel quantum defect theory and an optical potential becomes problematic, making it necessary to implement an alternative, approximate treatment of Auger decay for neutral Mg. The new cross sections are used to fit the Mg K edge in XMM-Newton spectra of the low-mass X-ray binary GS 1826-238, where most of the interstellar Mg is found to be in ionized form.

preprint2014arXiv

Physical properties of the interstellar medium using high-resolution Chandra spectra: O K-edge absorption

Chandra high-resolution spectra toward eight low-mass Galactic binaries have been analyzed with a photoionization model that is capable of determining the physical state of the interstellar medium. Particular attention is given to the accuracy of the atomic data. Hydrogen column densities are derived with a broadband fit that takes into account pileup effects, and in general are in good agreement with previous results. The dominant features in the oxygen-edge region are O I and O II K$α$ absorption lines whose simultaneous fits lead to average values of the ionization parameter of $\logξ=-2.90$ and oxygen abundance of $A_{\rm O}=0.70$. The latter is relative to the standard by Grevesse & Sauval (1998), but a rescaling with the revision by Asplund et al. (2009) would lead to an average abundance value fairly close to solar. The low average oxygen column density ($N_{\rm O}=9.2 \times 10^{17}$ cm$^{-2}$) suggests a correlation with the low ionization parameters, the latter also being in evidence in the column density ratios OII/OI and OIII/OI that are estimated to be less than 0.1. We do not find conclusive evidence for absorption by any other compound but atomic oxygen.

preprint2013arXiv

A Comprehensive X-ray Absorption Model for Atomic Oxygen

An analytical formula is developed to represent accurately the photoabsorption cross section of O I for all energies of interest in X-ray spectral modeling. In the vicinity of the Kedge, a Rydberg series expression is used to fit R-matrix results, including important orbital relaxation effects, that accurately predict the absorption oscillator strengths below threshold and merge consistently and continuously to the above-threshold cross section. Further minor adjustments are made to the threshold energies in order to reliably align the atomic Rydberg resonances after consideration of both experimental and observed line positions. At energies far below or above the K-edge region, the formulation is based on both outer- and inner-shell direct photoionization, including significant shake-up and shake-off processes that result in photoionization-excitation and double photoionization contributions to the total cross section. The ultimate purpose for developing a definitive model for oxygen absorption is to resolve standing discrepancies between the astronomically observed and laboratory measured line positions, and between the inferred atomic and molecular oxygen abundances in the interstellar medium from XSTAR and SPEX spectral models.

preprint2013arXiv

Suppression of Dielectronic Recombination Due to Finite Density Effects

We have developed a general model for determining density-dependent effective dielectronic recombination (DR) rate coefficients in order to explore finite-density effects on the ionization balance of plasmas. Our model consists of multiplying by a suppression factor those highly-accurate total zero-density DR rate coefficients which have been produced from state-of-the-art theoretical calculations and which have been benchmarked by experiment. The suppression factor is based-upon earlier detailed collision-radiative calculations which were made for a wide range of ions at various densities and temperatures, but used a simplified treatment of DR. A general suppression formula is then developed as a function of isoelectronic sequence, charge, density, and temperature. These density-dependent effective DR rate coefficients are then used in the plasma simulation code Cloudy to compute ionization balance curves for both collisionally ionized and photoionized plasmas at very low (ne = 1 cm^-3) and finite (ne=10^10 cm^-3) densities. We find that the denser case is significantly more ionized due to suppression of DR, warranting further studies of density effects on DR by detailed collisional-radiative calculations which utilize state-of-the-art partial DR rate coefficients. This is expected to impact the predictions of the ionization balance in denser cosmic gases such as those found in nova and supernova shells, accretion disks, and the broad emission line regions in active galactic nuclei.

preprint2012arXiv

Radiation Damping in the Photoionization of Fe^{14+}

A theoretical investigation of photoabsorption and photoionization of Fe^{14+} extending beyond an earlier frame transformation R-matrix implementation is performed using a fully-correlated, Breit-Pauli R-matrix formulation including both fine-structure splitting of strongly-bound resonances and radiation damping. The radiation damping of $2p\rightarrow nd$ resonances gives rise to a resonant photoionization cross section that is significantly lower than the total photoabsorption cross section. Furthermore, the radiation-damped photoionization cross section is found to be in good agreement with recent experimental results once a global shift in energy of $\approx -3.5$ eV is applied. These findings have important implications. Firstly, the presently available synchrotron experimental data are applicable only to photoionization processes and not to photoabsorption; the latter is required in opacity calculations. Secondly, our computed cross section, for which the L-shell ionization threshold is aligned with the NIST value, shows a series of $2p \rightarrow nd$ Rydberg resonances that are uniformly 3-4 eV higher in energy than the corresponding experimental profiles, indicating that the L-shell threshold energy values currently recommended by NIST are likely in error.

preprint2010arXiv

Dielectronic Recombination of Argon-Like Ions

We present a theoretical investigation of dielectronic recombination (DR) of Ar-like ions that sheds new light on the behavior of the rate coefficient at low-temperatures where these ions form in photoionized plasmas. We provide results for the total and partial Maxwellian-averaged DR rate coefficients from the initial ground level of K II -- Zn XIII ions. It is expected that these new results will advance the accuracy of the ionization balance for Ar-like M-shell ions and pave the way towards a detailed modeling of astrophysically relevant X-ray absorption features. We utilize the AUTOSTRUCTURE computer code to obtain the accurate core-excitation thresholds in target ions and carry out multiconfiguration Breit-Pauli (MCBP) calculations of the DR cross section in the independent-processes, isolated-resonance, distorted-wave (IPIRDW) approximation. Our results mediate the complete absence of direct DR calculations for certain Ar-like ions and question the reliability of the existing empirical rate formulas, often inferred from renormalized data within this isoelectronic sequence.

preprint2010arXiv

K-Shell Photoabsorption Studies of the Carbon Isonuclear Sequence

K-shell photoabsorption cross sections for the isonuclear C I - C IV ions have been computed using the R-matrix method. Above the K-shell threshold, the present results are in good agreement with the independent-particle results of Reilman & Manson (1979). Below threshold, we also compute the strong 1s -> np absorption resonances with the inclusion of important spectator Auger broadening effects. For the lowest 1s -> 2p, 3p resonances, comparisons to available C II, C III, and C IV experimental results show good agreement in general for the resonance strengths and positions, but unexplained discrepancies exist. Our results also provide detailed information on the C I K-shell photoabsorption cross section including the strong resonance features, since very limited laboratory experimental data exist. The resultant R-matrix cross sections are then used to model the Chandra X-ray absorption spectrum of the blazar Mkn 421.