-
A Photoionization model for the Infrared Coronal Line Emission in the Classical Nova V1716 Scorpii
Authors:
C. E. Woodward,
G. Shaw,
S. Starrfield,
A. Evans,
K. L. Page
Abstract:
A near-infrared spectrum of nova V1716 Scorpii (PNV J17224490-4137160), a recent bright (V_max = 7.3 mag), Fermi-LAT detected gamma-ray source, was modeled using the photoionization code CLOUDY. Abundances were estimated for He, C, N, O, Si, Al, Mg, Fe, Ne, S, Ca, and P. Notably, P (a factor of 120) and N (a factor of 248) are highly overabundant. It was necessary to assume the ejecta consist of t…
▽ More
A near-infrared spectrum of nova V1716 Scorpii (PNV J17224490-4137160), a recent bright (V_max = 7.3 mag), Fermi-LAT detected gamma-ray source, was modeled using the photoionization code CLOUDY. Abundances were estimated for He, C, N, O, Si, Al, Mg, Fe, Ne, S, Ca, and P. Notably, P (a factor of 120) and N (a factor of 248) are highly overabundant. It was necessary to assume the ejecta consist of two components (with a cylindrical geometry): a dense component from which the bulk of the H, He, and neutral O~I and N emission arises and a more diffuse component from which most of the coronal lines arise. Some of the coronal lines are found to originate from both the dense and diffuse components. The mass of the ejecta, including neutral and ionized gas, is ~ 4.19e-4 solar masses. Our analysis indicates that in the case of V1716 Sco (which has a carbon-oxygen white dwarf), a fraction of 25% white dwarf matter rather than 50% is favored for the mixing between white dwarf and the accreted envelope before the outburst. This mixing ratio is like that found for Oxygen-Neon novae where a 25% mixing fraction is also indicated. Helium hydride -- the first molecule to form after the Big Bang -- may have formed in the ejecta of V1716 Sco based on photoionization modeling. This prediction suggests that novae may be potential formation sites of this important molecular ion.
△ Less
Submitted 18 April, 2024;
originally announced April 2024.
-
Revealing Physical Properties of a Tidal Disruption Event: iPTF16fnl
Authors:
T. Mageshwaran,
Gargi Shaw,
Sudip Bhattacharyya,
Kimitake Hayasaki
Abstract:
Tidal disruption event (TDE) iPTF16fnl shows a relatively low optical flare with observationally very weak X-ray emission and the spectroscopic property that the helium emission line from the source dominates over the hydrogen emission line at early times. We explore these observed signatures by calculating spectral emission lines with the publicly available code, CLOUDY. We estimate five physical…
▽ More
Tidal disruption event (TDE) iPTF16fnl shows a relatively low optical flare with observationally very weak X-ray emission and the spectroscopic property that the helium emission line from the source dominates over the hydrogen emission line at early times. We explore these observed signatures by calculating spectral emission lines with the publicly available code, CLOUDY. We estimate five physical parameters by fitting the observed optical UV spectra on multiple days to a theoretical model of a steady-state, slim disk with a spherical outflow. The resultant key parameters among them are black hole mass $M_{\bullet} = (6.73 \pm 0.44) \times 10^5 M_{\odot}$, stellar mass $M_{\star} = (2.59 \pm 0.17) M_{\odot}$, and wind velocity $v_{\rm w} = 7447.43 \pm 183.9~{\rm km~s^{-1}}$. The disk-wind model also estimates the radiative efficiency to be $0.01\lesssimη\lesssim0.02$ over the observational time, resulting in the disk being radiatively inefficient, and the disk X-ray luminosity is consistent with the observed low luminosity. In our CLOUDY model, the filling factor of the wind is also estimated to be 0.8, suggesting that the wind is moderately clumpy. We reveal that the helium-to-hydrogen number density ratio of the wind lies between 0.1 and 0.15, which is nearly the same as the solar case, suggesting the tidally disrupted star is originally a main sequence star. Because the optical depth of the helium line is lower than the hydrogen line by two orders of magnitude, the helium line is significantly optically thinner than the hydrogen line. Consequently, our results indicate that the helium line luminosity dominates the hydrogen line luminosity due to the optical depth effect despite a small helium-to-hydrogen number density ratio value.
△ Less
Submitted 28 December, 2023;
originally announced December 2023.
-
The 2023 release of Cloudy
Authors:
Marios Chatzikos,
Stefano Bianchi,
Francesco Camilloni,
Priyanka Chakraborty,
Chamani M. Gunasekera,
Francisco Guzmán,
Jonathan S. Milby,
Arnab Sarkar,
Gargi Shaw,
Peter A. M. van Hoof,
Gary J. Ferland
Abstract:
We describe the 2023 release of the spectral synthesis code Cloudy. Since the previous major release, migrations of our online services motivated us to adopt git as our version control system. This change alone led us to adopt an annual release scheme, accompanied by a short release paper, the present being the inaugural. Significant changes to our atomic and molecular data have improved the accur…
▽ More
We describe the 2023 release of the spectral synthesis code Cloudy. Since the previous major release, migrations of our online services motivated us to adopt git as our version control system. This change alone led us to adopt an annual release scheme, accompanied by a short release paper, the present being the inaugural. Significant changes to our atomic and molecular data have improved the accuracy of Cloudy predictions: we have upgraded our instance of the Chianti database from version 7 to 10; our H- and He-like collisional rates to improved theoretical values; our molecular data to the most recent LAMDA database, and several chemical reaction rates to their most recent UDfA and KiDA values. Finally, we describe our progress on upgrading Cloudy's capabilities to meet the requirements of the X-ray microcalorimeters aboard the upcoming XRISM and Athena missions, and outline future development that will make Cloudy of use to the X-ray community.
△ Less
Submitted 11 August, 2023;
originally announced August 2023.
-
Revisiting the gas-phase chemical rate coefficients at high temperatures in CLOUDY
Authors:
Gargi Shaw,
Gary Ferland,
M. Chatzikos
Abstract:
A two-body gas-phase reaction rate coefficient can be given by the usual Arrhenius-type formula which depends on temperature. The UMIST Database for Astrochemistry is a widely used database for reaction rate coefficients. They provide fittings for coefficients valid over a particular range of temperatures. The permissible upper-temperature limits vary over a wide range: from 100 K to 41000K. A wid…
▽ More
A two-body gas-phase reaction rate coefficient can be given by the usual Arrhenius-type formula which depends on temperature. The UMIST Database for Astrochemistry is a widely used database for reaction rate coefficients. They provide fittings for coefficients valid over a particular range of temperatures. The permissible upper-temperature limits vary over a wide range: from 100 K to 41000K. A wide range of temperatures occurs in nature; thus, it requires evaluating the rate coefficients at temperatures outside the range of validity. As a result, a simple extrapolation of the rate coefficients can lead to unphysically large values at high temperatures. These result in unrealistic predictions. Here we present a solution to prevent the gas-phase reaction coefficients from diverging at a very high temperature. We implement this into the spectral synthesis code CLOUDY which operates over a wide range of temperatures from CMB to 10$^{10}$ K subject to different astrophysical environments.
△ Less
Submitted 20 July, 2023;
originally announced August 2023.
-
Recent updates of gas-phase chemical reactions and molecular lines of SiS in CLOUDY
Authors:
Gargi Shaw,
Gary Ferland,
M. Chatzikos
Abstract:
Here we present our current update of CLOUDY on gas-phase chemical reactions for the formation and destruction of the SiS molecule, its energy levels, and collisional rate coefficients with H$_2$, H, and He over a wide range of temperatures. As a result, henceforth the spectral synthesis code CLOUDY predicts SiS line intensities and column densities for various astrophysical environments.
Here we present our current update of CLOUDY on gas-phase chemical reactions for the formation and destruction of the SiS molecule, its energy levels, and collisional rate coefficients with H$_2$, H, and He over a wide range of temperatures. As a result, henceforth the spectral synthesis code CLOUDY predicts SiS line intensities and column densities for various astrophysical environments.
△ Less
Submitted 8 March, 2023;
originally announced March 2023.
-
Probing the tidal disruption event iPTF16axa with CLOUDY and disc-wind models
Authors:
T. Mageshwaran,
Gargi Shaw,
Sudip Bhattacharyya
Abstract:
We present both a disc-wind model on the optical/UV emission continuum and CLOUDY modelling on the spectral lines of the tidal disruption event (TDE) iPTF16axa to understand the disc-wind emission and the properties of the atmosphere that impacts the line luminosity of the TDE. Assuming the optical/UV emission from the wind due to the disc super-Eddington phase, we use the steady structured disc-w…
▽ More
We present both a disc-wind model on the optical/UV emission continuum and CLOUDY modelling on the spectral lines of the tidal disruption event (TDE) iPTF16axa to understand the disc-wind emission and the properties of the atmosphere that impacts the line luminosity of the TDE. Assuming the optical/UV emission from the wind due to the disc super-Eddington phase, we use the steady structured disc-wind model with a spherical wind with constant velocity to fit the observations on multiple days. The extracted parameters are stellar-mass $M_{\star} = 6.20 \pm 1.19 M_{\odot}$, disc radiative efficiency $\log_{\rm 10}(η) = -1.22 \pm 1.327$, wind inner radius $r_l = (2.013 \pm 0.551) \times 10^{14}~{\rm cm}$ and velocity $v_w = 18999.4 \pm 1785.1 ~{\rm km~s^{-1}}$. The photosphere temperature for wind emission is $ \sim 2 \times 10^4~{\rm K}$ and the disc single blackbody temperature is $\sim 0.995 \times 10^5~{\rm K}$. We also perform CLOUDY modelling to explain the observed He and H line luminosities that estimate a wind inner radius $r_l = 7.07 \times 10^{14}~{\rm cm}$ and velocity $v_w = 1.3 \times 10^4~{\rm km~s^{-1}}$. The independent analyses of iPTF16axa using CLOUDY and disc-wind models show comparable results that agree with observations. The CLOUDY modelling finds that both the super solar abundance of He and a smaller He II line optical depth is responsible for the enhancement of He II line luminosity over the H$α$ line luminosity. The super-solar abundance of He II agrees with a relatively large stellar mass and suggests that the disrupted star might have been a red giant.
△ Less
Submitted 25 November, 2022;
originally announced November 2022.
-
A recent update of gas-phase chemical reactions and molecular lines in CLOUDY: its effects on millimeter and sub-millimeter molecular line predictions
Authors:
Gargi Shaw,
G. J. Ferland,
M. Chatzikos
Abstract:
Here we present our current updates of the gas-phase chemical reaction rates and molecular lines in the spectral synthesis code CLOUDY, and its implications in spectroscopic modelling of various astrophysical environments. We include energy levels, radiative and collisional rates for HF, CF$^+$, HC$_3$N, ArH$^+$, HCl, HCN, CN, CH, and CH$_2$. Simultaneously, we expand our molecular network involvi…
▽ More
Here we present our current updates of the gas-phase chemical reaction rates and molecular lines in the spectral synthesis code CLOUDY, and its implications in spectroscopic modelling of various astrophysical environments. We include energy levels, radiative and collisional rates for HF, CF$^+$, HC$_3$N, ArH$^+$, HCl, HCN, CN, CH, and CH$_2$. Simultaneously, we expand our molecular network involving these molecules. For this purpose, we have added 561 new reactions and have updated the existing 165 molecular reaction rates involving these molecules. As a result, CLOUDY now predicts all the lines arising from these nine molecules. In addition, we also update H$_2$--H$_2$ collisional data up to rotational levels $J$=31 for $v$=0. We demonstrate spectroscopic simulations of these molecules for a few astrophysical environments. Our existing model for globules in the Crab nebula successfully predicts the observed column density of ArH$^+$. Our model predicts a detectable amount of HeH$^+$, OH$^+$, and CH$^+$ for the Crab nebula. We also model the ISM towards HD185418, W31C, NGC 253, and our predictions match with most of the observed column densities within the observed error bars. Very often molecular lines trace various physical conditions. Hence, this update will be very supportive for spectroscopic modelling of various astrophysical environments, particularly involving sub-millimeter and mid-infrared observations using ALMA and JWST, respectively.
△ Less
Submitted 9 June, 2022;
originally announced June 2022.
-
Detecting PAHs in high-z galaxies in proxy: Modelling physical conditions in an extremely strong damped Lyman-alpha absorber towards QSO SDSS J1143+1420 at z=2.323
Authors:
Gargi Shaw,
A. Ranjan
Abstract:
We explore indirect methods to detect Polycyclic Aromatic Hydrocarbons (PAHs) in gas-rich, absorption-selected galaxies at high redshift. We look at the optical VLT/X-shooter observations of an intervening, extremely strong damped Lyman-alpha absorber (or ESDLA, with log(N(HI))>~21.7)) towards QSO SDSS J1143+1420 at redshift, z(ESDLA)=2.323. Literature studies have shown that this ESDLA contains s…
▽ More
We explore indirect methods to detect Polycyclic Aromatic Hydrocarbons (PAHs) in gas-rich, absorption-selected galaxies at high redshift. We look at the optical VLT/X-shooter observations of an intervening, extremely strong damped Lyman-alpha absorber (or ESDLA, with log(N(HI))>~21.7)) towards QSO SDSS J1143+1420 at redshift, z(ESDLA)=2.323. Literature studies have shown that this ESDLA contains signatures of dust and diffuse molecular hydrogen and it was specifically chosen for our study due to its close spatial proximity (impact parameter, rho=0.6+/-0.3 kpc) with its associated galaxy. There is no direct detection of PAHs emission in the limited observations of infrared(IR)-spectra along this sight-line. Hence, we use CLOUDY numerical simulation modelling to indirectly probe the presence of PAH in the ESDLA. We note that PAHs need to be included in the models to reproduce the observed column densities of warm H2 and CI. Thus, we infer the presence of PAHs indirectly in our ESDLA, with an abundance of PAH/H = 10^(-7.046). We also measure a low 2175 A bump strength (E(bump)~0.03-0.19 mag) relative to star-forming galaxies by modelling extinction of QSO spectra by dust at the absorber rest-frame. This is consistent with the low PAH abundance obtained indirectly using CLOUDY modelling. Our study highlights the usage of CLOUDY modelling to indirectly detect PAH in high-redshift gas-rich absorption-selected galaxies.
△ Less
Submitted 6 May, 2022;
originally announced May 2022.
-
Multi-phase gas properties of extremely strong intervening DLAs towards quasars
Authors:
A. Ranjan,
R. Srianand,
P. Petitjean,
G. Shaw,
Y. -K. Sheen,
S. A. Balashev,
N. Gupta,
C. Ledoux,
K. N. Telikova
Abstract:
We present the results of a spectroscopic analysis of extremely strong damped Lyman-α absorbers (ESDLAs, log N(Hi)>=21.7) observed with VLT-XShooter. ESDLAs probe gas from within the star-forming disk of the associated galaxies and thus ESDLAs provide a unique opportunity to study the interstellar medium of galaxies at high-redshift. We report column densities (N), equivalent widths (w), and the k…
▽ More
We present the results of a spectroscopic analysis of extremely strong damped Lyman-α absorbers (ESDLAs, log N(Hi)>=21.7) observed with VLT-XShooter. ESDLAs probe gas from within the star-forming disk of the associated galaxies and thus ESDLAs provide a unique opportunity to study the interstellar medium of galaxies at high-redshift. We report column densities (N), equivalent widths (w), and the kinematic spread (Δ v90) of species from neutral, singly ionised, and higher ionisation species. We find that, using the dust correction prescription, the measured metallicities are consistent for singly ionised gas species such as Pii, S ii, Si ii, Mnii and Crii, and Znii in all ESDLAs within 3-sigma uncertainty. We find that the distributions of N(Ari)/N(Hi) ratio in DLAs and ESDLAs are similar. We further report that ESDLAs do not show a strong deficiency of Ari relative to other α-capture elements as is seen in DLAs. This supports the idea that the mentioned under-abundance of Ari in DLAs is possibly caused by the presence of background UV photons that penetrate the low N(Hi) clouds to ionise Ari, but they cannot penetrate deep enough in the high N(Hi) ESDLA environment. The w(Mgii lambda2796) distribution in ESDLAs is found to be similar to that of metal-rich Ci-selected absorbers, but the velocity spread of their Mgii profile is different. For higher ionisation species (such as C iv and Si iv), Δ v90 is similar in the two populations, while the Δ v90 of singly ionised species is smaller for ESDLAs. This suggests that the ESDLAs sample a different Hi region of their associated galaxy compared to the general DLA population. We further study the N(Cl i) distribution in high-redshift DLA and ESDLA sightlines, as Cl i is a good tracer of H2 gas. The N(Cl i)-N(H2) correlation is followed by all the clouds (ESDLAs and otherwise) having log N(H2)<22.
△ Less
Submitted 6 May, 2022; v1 submitted 17 January, 2022;
originally announced January 2022.
-
Photoionization modeling of the dusty nova V1280 Scorpii
Authors:
Ruchi Pandey,
Ramkrishna Das,
Gargi Shaw,
Soumen Mondal
Abstract:
We perform photoionization modeling of the dusty nova V1280 Scorpii (V1280 Sco) with an aim to study the changes in the physical and chemical parameters. We model pre and post dust phase, optical and near-Infrared (NIR), spectra using the photoionization code \textsc{cloudy}, v.17.02, considering a two-component (low density and high density region) model. From the best-fit model, we find that the…
▽ More
We perform photoionization modeling of the dusty nova V1280 Scorpii (V1280 Sco) with an aim to study the changes in the physical and chemical parameters. We model pre and post dust phase, optical and near-Infrared (NIR), spectra using the photoionization code \textsc{cloudy}, v.17.02, considering a two-component (low density and high density region) model. From the best-fit model, we find that the temperature and luminosity of the central ionizing source in the pre-dust phase are in the range 1.32 - 1.50 $\times 10^4$ K and 2.95 - 3.16 $\times 10^{36}$ ergs$^{-1}$, respectively, which increase to 1.58 - 1.62 $\times 10^4$ K and 3.23 - 3.31 $\times 10^{36}$ ergs$^{-1}$, respectively, in the post-dust phase. It is found that a very high hydrogen density ($\sim 10^{13} - 10^{14}$ cm$^{-3}$) is required for the generation of spectra properly. Dust condensation conditions are achieved at high ejecta density ($\sim 3.16 \times 10^{8}$cm$^{-3}$) and low temperature ($\sim$2000 K) in the outer region of the ejecta. It is found that a mixture of small (0.005 - 0.25$μ$m) amorphous carbon dust grains and large (0.03 - 3.0$μ$m) astrophysical silicate dust grains iis present n the ejecta in the post-dust phase. Our model yields very high elemental abundance values as C/H = 13.5 - 20, N/H = 250, O/H = 27 - 35, by number, relative to solar in the ejecta, during the pre-dust phase, which decrease in the post-dust phase.
△ Less
Submitted 3 November, 2021;
originally announced November 2021.
-
Spectroscopic modelling of two high-mass X-ray binaries, Cyg X--3 and 4U 1538--522
Authors:
Gargi Shaw,
Sudip Bhattacharyya
Abstract:
We report a detailed modelling of soft X-ray emission lines from two stellar-wind fed Galactic high mass X-ray binary (HMXB) systems, Cyg X-3 and 4U 1538-522, and estimate physical parameters, e.g., hydrogen density, radiation field, chemical abundances, wind velocity, etc. The spectral synthesis code CLOUDY is utilized for this modelling. We model highly ionised X-ray spectral lines such as Fe XX…
▽ More
We report a detailed modelling of soft X-ray emission lines from two stellar-wind fed Galactic high mass X-ray binary (HMXB) systems, Cyg X-3 and 4U 1538-522, and estimate physical parameters, e.g., hydrogen density, radiation field, chemical abundances, wind velocity, etc. The spectral synthesis code CLOUDY is utilized for this modelling. We model highly ionised X-ray spectral lines such as Fe XXV (6.700 keV), Fe XXVI (6.966 keV), and reproduce the observed line flux values. We find that for Cyg X--3 and 4U 1538-522, the inner radius of the ionised gas is at a distance of 10$^{12.25}$ cm and 10$^{10.43}$ cm respectively from the primary star, which is the main source of ionisation. The densities of the ionised gas for Cyg X--3 and 4U 1538--522 are found to be $\sim$ 10$^{11.35}$ cm$^{-3}$ and 10$^{11.99}$ cm$^{-3}$, respectively. The corresponding wind velocities are 2000 km s$^{-1}$ and 1500 km s$^{-1}$. The respective predicted hydrogen column densities for Cyg X--3 and 4U 1538--522 are $10^{23.2}$ cm$^{-2}$ and 10$^{22.25}$ cm$^{-2}$. In addition, we find that a magnetic field affects the strength of the spectral lines through cyclotron cooling. Hence, we perform separate model comparisons including a magnetic field for both the sources. Most of the parameters, except the hydrogen column density, have similar values with and without a magnetic field. We estimate that the most probable strength of the magnetic field for Cyg X--3 and 4U 1538--522, where the Fe XXV and Fe XXVI lines originate, is $\sim$ 10$^{2.5}$G.
△ Less
Submitted 30 July, 2021;
originally announced July 2021.
-
Role of Polycyclic Aromatic Hydrocarbons on the Cosmic-Ray ionization rate in the Galaxy
Authors:
Gargi Shaw,
G. J. Ferland
Abstract:
The cosmic-ray ionization rate ($ζ$, s$^{-1}$) plays an important role in the interstellar medium. It controls ion-molecular chemistry and provides a source of heating. Here we perform a grid of calculations using the spectral synthesis code CLOUDY along nine sightlines towards, HD 169454, HD 110432, HD 204827, $λ$ Cep, X Per, HD 73882, HD 154368, Cyg OB2 5, Cyg OB2 12. The value of $ζ$ is determi…
▽ More
The cosmic-ray ionization rate ($ζ$, s$^{-1}$) plays an important role in the interstellar medium. It controls ion-molecular chemistry and provides a source of heating. Here we perform a grid of calculations using the spectral synthesis code CLOUDY along nine sightlines towards, HD 169454, HD 110432, HD 204827, $λ$ Cep, X Per, HD 73882, HD 154368, Cyg OB2 5, Cyg OB2 12. The value of $ζ$ is determined by matching the observed column densities of H$_3^+$ and H$_2$. The presence of polycyclic aromatic hydrocarbons (PAHs) affects the free electron density, which changes the H$_3^+$ density and the derived ionization rate. PAHs are ubiquitous in the Galaxy, but there are also regions where PAHs do not exist. Hence, we consider clouds with a range of PAH abundances and show their effects on the H$_3^+$ abundance. We predict an average cosmic-ray ionization rate for H$_2$ ($ζ$(H$_2$))= (7.88 $\pm$ 2.89) $\times$ 10$^{-16}$ s$^{-1}$ for models with average Galactic PAHs abundances, (PAH/H =10$^{-6.52}$), except Cyg OB2 5 and Cyg OB2 12. The value of $ζ$ is nearly 1 dex smaller for sightlines toward Cyg OB2 12. We estimate the average value of $ζ$(H$_2$)= (95.69 $\pm$ 46.56) $\times$ 10$^{-16}$ s$^{-1}$ for models without PAHs.
△ Less
Submitted 12 January, 2021; v1 submitted 11 January, 2021;
originally announced January 2021.
-
Cosmic ray dissociation of molecular hydrogen and dense cloud chemistry
Authors:
Gargi Shaw,
G. J. Ferland,
S. Ploeckinger
Abstract:
Dissociation of molecular hydrogen by secondary electrons produced by cosmic ray or X-ray ionization plays a crucial role in the chemistry of the densest part of molecular clouds. Here we study the effect of the mean kinetic energy of secondary electrons on this process. We compare predictions using a range of secondary electron energies and predictions of the cross-sections with the values in the…
▽ More
Dissociation of molecular hydrogen by secondary electrons produced by cosmic ray or X-ray ionization plays a crucial role in the chemistry of the densest part of molecular clouds. Here we study the effect of the mean kinetic energy of secondary electrons on this process. We compare predictions using a range of secondary electron energies and predictions of the cross-sections with the values in the UMIST database. We find that the predicted column densities change by nearly one dex.
△ Less
Submitted 30 October, 2020; v1 submitted 19 October, 2020;
originally announced October 2020.
-
Physical conditions in two high-redshift H$_2$-bearing GRB-DLAs, 120815A and 121024A
Authors:
Gargi Shaw,
G. J. Ferland
Abstract:
The gamma-ray burst (GRB) afterglows provide an unique opportunity to study the interstellar medium (ISM) of star-forming galaxies at high-redshift. The GRB-DLAs (damped Lyman-$α$ absorbers) contain a large neutral hydrogen column density, N(H I), and are observed against the GRB afterglow. A large fraction of GRB-DLAs show presence of molecular hydrogen (H$_2$) which is an indicator of star-forma…
▽ More
The gamma-ray burst (GRB) afterglows provide an unique opportunity to study the interstellar medium (ISM) of star-forming galaxies at high-redshift. The GRB-DLAs (damped Lyman-$α$ absorbers) contain a large neutral hydrogen column density, N(H I), and are observed against the GRB afterglow. A large fraction of GRB-DLAs show presence of molecular hydrogen (H$_2$) which is an indicator of star-formation. Hence it is important to study those GRB-DLAs which have H$_2$ lines to decipher and understand their physical conditions. The GRB-DLAs 121024A and 120815A, situated at redshift 2.30 and 2.36 respectively, are two such important H$_2$-bearing GRB-DLAs. Besides H$_2$, these two GRB-DLAs also show many metal lines. In this work we have carried out a detail numerical study on the H$_2$ lines, as well as on those metal lines, in GRB-DLAs 121024A and 120815A self-consistently. We use the spectral synthesis code CLOUDY for this study. This modeling helps us to determine the underlying physical conditions which give rise such lines and hence to understand these two GRB-DLAs in much more detail than any other previous investigation. We find that the hydrogen densities for these two H$_2$-bearing DLAs are $\geq$ 60 cm$^{-3}$. Moreover our study infers that the linear sizes are $\leq 17.7$ pc for these two GRB-DLAs, and the mean gas temperatures averaged over the cloud thickness, are $\leq$ 140 K.
△ Less
Submitted 1 March, 2020;
originally announced March 2020.
-
Spectroscopic modelling of four neutron star low-mass X-ray binaries using CLOUDY
Authors:
Gargi Shaw,
Sudip Bhattacharyya
Abstract:
Low-mass X-ray binaries (LMXBs) have a wide range of X-ray properties which can be utilised to reveal many physical conditions of the associated accretion discs. We use the spectral synthesis code CLOUDY to perform a detailed modelling of neutron star LMXBs GX 13+1, MXB 1659--298, 4U 1323--62 and XB 1916--053; and characterise the underlying physical conditions, such as density, radiation field, m…
▽ More
Low-mass X-ray binaries (LMXBs) have a wide range of X-ray properties which can be utilised to reveal many physical conditions of the associated accretion discs. We use the spectral synthesis code CLOUDY to perform a detailed modelling of neutron star LMXBs GX 13+1, MXB 1659--298, 4U 1323--62 and XB 1916--053; and characterise the underlying physical conditions, such as density, radiation field, metallicity, wind velocity, etc. For this purpose we model highly ionised spectra of Fe, Ca, S, Si, Mg, Al in the soft X-ray band, and compare the predicted line flux ratios with the observed values. We also find that the strength and profile of these spectral lines get modified in the presence of magnetic field in the accretion disc. Using this, we estimate an upper limit of the existing magnetic field to be about a few hundred to a few thousand G in the accretion discs of these four LMXBs.
△ Less
Submitted 13 March, 2019;
originally announced March 2019.
-
A photoionization model grid for novae: estimation of physical parameters
Authors:
Anindita Mondal,
Ramkrishna Das,
Gargi Shaw,
Soumen Mondal
Abstract:
We present here a method to estimate physical parameters of novae systems using an extensive grid of photoionization models for novae. We use the photoionization code CLOUDY to construct grid of models covering a wide range of different parameters, e.g. total hydrogen density ($n_H$), source temperature ($T_{BB}$) and luminosity ($L$), inner radius ($R_{in}$) and thickness of ejecta ($ΔR$), keepin…
▽ More
We present here a method to estimate physical parameters of novae systems using an extensive grid of photoionization models for novae. We use the photoionization code CLOUDY to construct grid of models covering a wide range of different parameters, e.g. total hydrogen density ($n_H$), source temperature ($T_{BB}$) and luminosity ($L$), inner radius ($R_{in}$) and thickness of ejecta ($ΔR$), keeping other elements at solar metallicity. In this way, a total of 1792 models have been generated. From the model generated spectra which cover a wide wavelength region from ultra-violet to infrared, we calculate ratios of hydrogen and helium emission lines fluxes which are generally strong in novae spectra. We show that physical parameters associated with novae system could be estimated by comparing these line ratios with those obtained from observed spectra. We elaborate the idea with examples and estimate the parameter values in case of few other novae. The results of the grid model are available online.
△ Less
Submitted 16 November, 2018;
originally announced November 2018.
-
Multicomponent H2 in DLA at zabs = 2.05: physical conditions through observations and numerical models
Authors:
Katherine Rawlins,
Raghunathan Srianand,
Gargi Shaw,
Hadi Rahmani,
Rajeshwari Dutta,
Sajeev Chacko
Abstract:
We perform detailed spectroscopic analysis and numerical modelling of an H2-bearing damped Lyman-alpha absorber (DLA) at zabs = 2.05 towards the quasar FBQS J2340-0053. Metal absorption features arise from fourteen components spread over $Δv_{90}$ = 114 km s$^{-1}$, seven of which harbour H2. Column densities of atomic and molecular species are derived through Voigt profile analysis of their absor…
▽ More
We perform detailed spectroscopic analysis and numerical modelling of an H2-bearing damped Lyman-alpha absorber (DLA) at zabs = 2.05 towards the quasar FBQS J2340-0053. Metal absorption features arise from fourteen components spread over $Δv_{90}$ = 114 km s$^{-1}$, seven of which harbour H2. Column densities of atomic and molecular species are derived through Voigt profile analysis of their absorption lines. We measure total N(H I), N(H2) and N(HD) to be 20.35+/-0.05, 17.99+/-0.05 and 14.28+/-0.08 (log cm$^{-2}$) respectively. H2 is detected in the lowest six rotational levels of the ground vibrational state. The DLA has metallicity, Z = 0.3 Z$_\sun$ ([S/H] = -0.52+/-0.06) and dust-to-gas ratio, $κ$ = 0.34+/-0.07. Numerical models of the H2 components are constrained individually to understand the physical structure of the DLA. We conclude that the DLA is subjected to the metagalactic background radiation and cosmic ray ionization rate of $\sim$ 10$^{-15.37}$ s$^{-1}$. Dust grains in this DLA are smaller than grains in the Galactic interstellar medium. The inner molecular regions of the H2 components have density, temperature and gas pressure in the range 30-120 cm$^{-3}$, 140-360 K and 7,000-23,000 cm$^{-3}$ K respectively. Micro-turbulent pressure is a significant constituent of the total pressure, and can play an important role in these innermost regions. Our H2 component models enable us to constrain component-wise N(H I), and elemental abundances of sulphur, silicon, iron and carbon. We deduce the line-of-sight thickness of the H2-bearing parts of the DLA to be 7.2 pc.
△ Less
Submitted 23 August, 2018;
originally announced August 2018.
-
The validity of 21 cm spin temperature as a kinetic temperature indicator in atomic and molecular gas
Authors:
Gargi Shaw,
G. J. Ferland,
I. Hubeny
Abstract:
The gas kinetic temperature ($T_K$) of various interstellar environments is often inferred from observations that can deduce level populations of atoms, ions, or molecules using spectral line observations; H I 21 cm is perhaps the most widely used with a long history. Usually the H I 21 cm line is assumed to be in thermal equilibrium and the populations are given by the Boltzmann distribution. A v…
▽ More
The gas kinetic temperature ($T_K$) of various interstellar environments is often inferred from observations that can deduce level populations of atoms, ions, or molecules using spectral line observations; H I 21 cm is perhaps the most widely used with a long history. Usually the H I 21 cm line is assumed to be in thermal equilibrium and the populations are given by the Boltzmann distribution. A variety of processes, many involving Lyman alpha ($Lyα$), can affect the 21 cm line. Here we show how this is treated in the spectral simulation code Cloudy, and present numerical simulations of environments where this temperature indicator is used, with a detailed treatment of the physical processes that determine level populations within $H^0$. We discuss situations where this temperature indicator traces $T_K$, cases where they fail, as well as the effects of $Lyα$ pumping on the 21 cm spin temperature. We also show that the $Lyα$ excitation temperature rarely traces the gas kinetic temperature.
△ Less
Submitted 5 June, 2017;
originally announced June 2017.
-
Physical conditions in three high-z H2-bearing DLAs: implications for grain size
Authors:
Katherine Rawlins,
Gargi Shaw,
Raghunathan Srianand
Abstract:
We present results of our numerical simulation of three H2-bearing damped Lyman alpha absorbers (DLAs) in the redshift interval ~ 2-3. The systems we have modelled are the DLAs at zabs = 2.3377 towards the quasar LBQS 1232+0815, at zabs = 2.41837 towards SDSS J143912.04+111740.5 and at zabs = 2.6265 towards FBQS J081240.6+320808. We have used the spectral synthesis code CLOUDY to simulate the phys…
▽ More
We present results of our numerical simulation of three H2-bearing damped Lyman alpha absorbers (DLAs) in the redshift interval ~ 2-3. The systems we have modelled are the DLAs at zabs = 2.3377 towards the quasar LBQS 1232+0815, at zabs = 2.41837 towards SDSS J143912.04+111740.5 and at zabs = 2.6265 towards FBQS J081240.6+320808. We have used the spectral synthesis code CLOUDY to simulate the physical environment of these DLAs, and constrain the density, radiation field, geometry and dust-grain properties of the DLAs self-consistently based on the observed column densities of various atomic and molecular species such as H I, fine structure lines of C I and the rotational level population of H2. In our models, we explore the effect of grain size distribution on the predicted column densities of different species. Within the allowed uncertainties in the inferred dust-to-gas ratio, both models with standard ISM grains and smaller-sized grains reproduce the observations equally well. Improved constraints on dust-to-gas ratio and line-of-sight extinction are important for probing the grain size distribution in high-z DLAs. We find the H2-bearing clouds to have line-of-sight thickness in the range 1-6 pc, consistent with what has been found using partial coverage or 21-cm observations in some high-z DLAs.
△ Less
Submitted 6 April, 2016;
originally announced April 2016.
-
Crustal Failure on Icy Moons from a Strong Tidal Encounter
Authors:
Alice C. Quillen,
David Giannella,
John G. Shaw,
Cynthia Ebinger
Abstract:
Close tidal encounters among large planetesimals and moons should have been more common than grazing or normal impacts. Using a mass spring model within an N-body simulation, we simulate the deformation of the surface of an elastic spherical body caused by a close parabolic tidal encounter with a body that has similar mass as that of the primary body. Such an encounter can induce sufficient stress…
▽ More
Close tidal encounters among large planetesimals and moons should have been more common than grazing or normal impacts. Using a mass spring model within an N-body simulation, we simulate the deformation of the surface of an elastic spherical body caused by a close parabolic tidal encounter with a body that has similar mass as that of the primary body. Such an encounter can induce sufficient stress on the surface to cause brittle failure of an icy crust and simulated fractures can extend a large fraction of the radius of body. Strong tidal encounters may be responsible for the formation of long graben complexes and chasmata in ancient terrain of icy moons such as Dione, Tethys, Ariel and Charon.
△ Less
Submitted 5 April, 2016; v1 submitted 7 December, 2015;
originally announced December 2015.
-
The 2013 Release of Cloudy
Authors:
G. J. Ferland,
R. L. Porter,
P. A. M. van Hoof,
R. J. R. Williams,
N. P. Abel,
M. L. Lykins,
Gargi Shaw,
W. J. Henney,
P. C. Stancil
Abstract:
This is a summary of the 2013 release of the plasma simulation code Cloudy. Cloudy models the ionization, chemical, and thermal state of material that may be exposed to an external radiation field or other source of heating, and predicts observables such as emission and absorption spectra. It works in terms of elementary processes, so is not limited to any particular temperature or density regime.…
▽ More
This is a summary of the 2013 release of the plasma simulation code Cloudy. Cloudy models the ionization, chemical, and thermal state of material that may be exposed to an external radiation field or other source of heating, and predicts observables such as emission and absorption spectra. It works in terms of elementary processes, so is not limited to any particular temperature or density regime. This paper summarizes advances made since the last major review in 1998. Much of the recent development has emphasized dusty molecular environments, improvements to the ionization / chemistry solvers, and how atomic and molecular data are used. We present two types of simulations to demonstrate the capability of the code. We consider a molecular cloud irradiated by an X-ray source such as an Active Nucleus and show how treating EUV recombination lines and the full SED affects the observed spectrum. A second example illustrates the very wide range of particle and radiation density that can be considered.
△ Less
Submitted 18 February, 2013;
originally announced February 2013.
-
Rotationally Warm Molecular Hydrogen in the Orion Bar
Authors:
Gargi Shaw,
G. J. Ferland,
W. J. Henney,
P. C. Stancil,
N. P. Abel,
E. W. Pellegrini,
J. A. Baldwin,
P. A. M. van Hoof
Abstract:
The Orion Bar is one of the nearest and best-studied photodissociation or photon-dominated regions (PDRs). Observations reveal the presence of H2 lines from vibrationally or rotationally excited upper levels that suggest warm gas temperatures (400 to 700 K). However, standard models of PDRs are unable to reproduce such warm rotational temperatures. In this paper we attempt to explain these obser…
▽ More
The Orion Bar is one of the nearest and best-studied photodissociation or photon-dominated regions (PDRs). Observations reveal the presence of H2 lines from vibrationally or rotationally excited upper levels that suggest warm gas temperatures (400 to 700 K). However, standard models of PDRs are unable to reproduce such warm rotational temperatures. In this paper we attempt to explain these observations with new comprehensive models which extend from the H+ region through the Bar and include the magnetic field in the equation of state. We adopt the model parameters from our previous paper which successfully reproduced a wide variety of spectral observations across the Bar. In this model the local cosmic-ray density is enhanced above the galactic background, as is the magnetic field, and which increases the cosmic-ray heating elevating the temperature in the molecular region. The pressure is further enhanced above the gas pressure in the H+ region by the momentum transferred from the absorbed starlight. Here we investigate whether the observed H2 lines can be reproduced with standard assumptions concerning the grain photoelectric emission. We also explore the effects due to the inclusion of recently computed H2 + H2, H2 + H and H2 + He collisional rate coefficients.
△ Less
Submitted 12 June, 2009;
originally announced June 2009.
-
Orion's Bar: Physical Conditions across the Definitive H+ / H0 / H2 Interface
Authors:
E. W. Pellegrini,
J. A. Baldwin,
G. J. Ferland,
Gargi Shaw,
S. Heathcote
Abstract:
Previous work has shown the Orion Bar to be an interface between ionized and molecular gas, viewed roughly edge on, which is excited by the light from the Trapezium cluster. Much of the emission from any star-forming region will originate from such interfaces, so the Bar serves as a foundation test of any emission model. Here we combine X-ray, optical, IR and radio data sets to derive emission s…
▽ More
Previous work has shown the Orion Bar to be an interface between ionized and molecular gas, viewed roughly edge on, which is excited by the light from the Trapezium cluster. Much of the emission from any star-forming region will originate from such interfaces, so the Bar serves as a foundation test of any emission model. Here we combine X-ray, optical, IR and radio data sets to derive emission spectra along the transition from H+ to H0 to H2 regions. We then reproduce the spectra of these layers with a simulation that simultaneously accounts for the detailed microphysics of the gas, the grains, and molecules, especially H2 and CO. The magnetic field, observed to be the dominant pressure in another region of the Orion Nebula, is treated as a free parameter, along with the density of cosmic rays. Our model successfully accounts for the optical, IR and radio observations across the Bar by including a significant magnetic pressure and also heating by an excess density of cosmic rays, which we suggest is due to cosmic rays being trapped in the compressed magnetic field. In the Orion Bar, as we had previously found in M17, momentum carried by radiation and winds from the newly formed stars pushes back and compresses the surrounding gas. There is a rough balance between outward momentum in starlight and the total pressure in atomic and molecular gas surrounding the H+ region. If the gas starts out with a weak magnetic field, the starlight from a newly formed cluster will push back the gas and compress the gas, magnetic field, and cosmic rays until magnetic pressure becomes an important factor.
△ Less
Submitted 14 November, 2008; v1 submitted 7 November, 2008;
originally announced November 2008.
-
Sensitivity of PDR Calculations to Microphysical Details
Authors:
N. P. Abel,
P. A. M. van Hoof,
G. Shaw,
G. J. Ferland,
T. Elwert
Abstract:
Our understanding of physical processes in Photodissociation regions or Photon Dominated Regions (PDRs) largely depends on the ability of spectral synthesis codes to reproduce the observed infrared emission-line spectrum. In this paper, we explore the sensitivity of a single PDR model to microphysical details. Our calculations use the Cloudy spectral synthesis code, recently modified to include…
▽ More
Our understanding of physical processes in Photodissociation regions or Photon Dominated Regions (PDRs) largely depends on the ability of spectral synthesis codes to reproduce the observed infrared emission-line spectrum. In this paper, we explore the sensitivity of a single PDR model to microphysical details. Our calculations use the Cloudy spectral synthesis code, recently modified to include a wealth of PDR physical processes. We show how the chemical/thermal structure of a PDR, along with the calculated spectrum, changes when the treatment of physical processes such as grain physics and atomic/molecular rates are varied. We find a significant variation in the intensities of PDR emission lines, depending on different treatments of the grain physics. We also show how different combinations of the cosmic-ray ionization rate, inclusion of grain-atom/ion charge transfer, and the grain size distribution can lead to very similar results for the chemical structure. Additionally, our results show the utility of Cloudy for the spectral modeling of molecular environments.
△ Less
Submitted 19 August, 2008;
originally announced August 2008.
-
Ultraviolet Survey of CO and H_2 in Diffuse Molecular Clouds: The Reflection of Two Photochemistry Regimes in Abundance Relationships
Authors:
Y. Sheffer,
M. Rogers,
S. R. Federman,
N. P. Abel,
R. Gredel,
D. L. Lambert,
G. Shaw
Abstract:
(Abridged) We carried out a comprehensive far-ultraviolet (UV) survey of ^12CO and H_2 column densities along diffuse molecular Galactic sight lines in order to explore in detail the relationship between CO and H_2. We measured new CO abundances from HST spectra, new H_2 abundances from FUSE data, and new CH, CH^+, and CN abundances from the McDonald and European Southern Observatories. A plot o…
▽ More
(Abridged) We carried out a comprehensive far-ultraviolet (UV) survey of ^12CO and H_2 column densities along diffuse molecular Galactic sight lines in order to explore in detail the relationship between CO and H_2. We measured new CO abundances from HST spectra, new H_2 abundances from FUSE data, and new CH, CH^+, and CN abundances from the McDonald and European Southern Observatories. A plot of log N(CO) versus log N(H_2) shows that two power-law relationships are needed for a good fit of the entire sample, with a break located at log N(CO, cm^-2) = 14.1 and log N(H_2) = 20.4, corresponding to a change in production route for CO in higher-density gas. Similar logarithmic plots among all five diatomic molecules allow us to probe their relationships, revealing additional examples of dual slopes in the cases of CO versus CH (break at log N = 14.1, 13.0), CH^+ versus H_2 (13.1, 20.3), and CH^+ versus CO (13.2, 14.1). These breaks are all in excellent agreement with each other, confirming the break in the CO versus H_2 relationship, as well as the one-to-one correspondence between CH and H_2 abundances. Our new sight lines were selected according to detectable amounts of CO in their spectra and they provide information on both lower-density (< 100 cm^-3) and higher-density diffuse clouds. The CO versus H_2 correlation and its intrinsic width are shown to be empirically related to the changing total gas density among the sight lines of the sample. We employ both analytical and numerical chemical schemes in order to derive details of the molecular environments. In the low-density gas, where equilibrium-chemistry studies have failed to reproduce the abundance of CH^+, our numerical analysis shows that nonequilibrium chemistry must be employed for correctly predicting the abundances of both CH^+ and CO.
△ Less
Submitted 7 July, 2008;
originally announced July 2008.
-
Rotational quenching rate coefficients for H_2 in collisions with H_2 from 2 to 10,000 K
Authors:
T. -G. Lee,
N. Balakrishnan,
R. C. Forrey,
P. C. Stancil,
G. Shaw,
D. R. Schultz,
G. J. Ferland
Abstract:
Rate coefficients for rotational transitions in H_2 induced by H_2 impact are presented. Extensive quantum mechanical coupled-channel calculations based on a recently published (H_2)_2 potential energy surface were performed. The potential energy surface used here is presumed to be more reliable than surfaces used in previous work. Rotational transition cross sections with initial levels J <= 8…
▽ More
Rate coefficients for rotational transitions in H_2 induced by H_2 impact are presented. Extensive quantum mechanical coupled-channel calculations based on a recently published (H_2)_2 potential energy surface were performed. The potential energy surface used here is presumed to be more reliable than surfaces used in previous work. Rotational transition cross sections with initial levels J <= 8 were computed for collision energies ranging between 0.0001 and 2.5 eV, and the corresponding rate coefficients were calculated for the temperature range 2 < T <10,000 K. In general, agreement with earlier calculations, which were limited to 100-6000 K, is good though discrepancies are found at the lowest and highest temperatures. Low-density-limit cooling functions due to para- and ortho-H_2 collisions are obtained from the collisional rate coefficients. Implications of the new results for non-thermal H_2 rotational distributions in molecular regions are also investigated.
△ Less
Submitted 12 May, 2008;
originally announced May 2008.
-
Merged ionization/dissociation fronts in planetary nebulae
Authors:
William J. Henney,
R. J. R. Williams,
Gary J. Ferland,
Gargi Shaw,
C. R. O'Dell
Abstract:
The hydrogen ionization and dissociation front around an ultraviolet radiation source should merge when the ratio of ionizing photon flux to gas density is sufficiently low and the spectrum is sufficiently hard. This regime is particularly relevant to the molecular knots that are commonly found in evolved planetary nebulae, such as the Helix Nebula, where traditional models of photodissociation…
▽ More
The hydrogen ionization and dissociation front around an ultraviolet radiation source should merge when the ratio of ionizing photon flux to gas density is sufficiently low and the spectrum is sufficiently hard. This regime is particularly relevant to the molecular knots that are commonly found in evolved planetary nebulae, such as the Helix Nebula, where traditional models of photodissociation regions have proved unable to explain the high observed luminosity in H_2 lines. In this paper we present results for the structure and steady-state dynamics of such advection-dominated merged fronts, calculated using the Cloudy plasma/molecular physics code. We find that the principal destruction processes for H_2 are photoionization by extreme ultraviolet radiation and charge exchange reactions with protons, both of which form H_2^+, which rapidly combines with free electrons to undergo dissociative recombination. Advection moves the dissociation front to lower column densities than in the static case, which vastly increases the heating in the partially molecular gas due to photoionization of He^0, H_2, and H^0. This causes a significant fraction of the incident bolometric flux to be re-radiated as thermally excited infrared H_2 lines, with the lower excitation pure rotational lines arising in 1000 K gas and higher excitation H_2 lines arising in 2000 K gas, as is required to explain the H_2 spectrum of the Helix cometary knots.
△ Less
Submitted 27 November, 2007; v1 submitted 27 November, 2007;
originally announced November 2007.
-
On the enhanced cosmic-ray ionization rate in the diffuse cloud towards Zeta Persei
Authors:
Gargi Shaw,
G. J. Ferland,
R. Srianand,
N. P. Abel,
P. A. M. van Hoof,
P. C. Stancil
Abstract:
The spatial distribution of the cosmic-ray flux is important in understanding the Interstellar Medium (ISM) of the Galaxy. This distribution can be analyzed by studying different molecular species along different sight lines whose abundances are sensitive to the cosmic-ray ionization rate. Recently several groups have reported an enhanced cosmic-ray ionization rate in diffuse clouds compared to…
▽ More
The spatial distribution of the cosmic-ray flux is important in understanding the Interstellar Medium (ISM) of the Galaxy. This distribution can be analyzed by studying different molecular species along different sight lines whose abundances are sensitive to the cosmic-ray ionization rate. Recently several groups have reported an enhanced cosmic-ray ionization rate in diffuse clouds compared to the standard value, zeta(average)=2.5e-17 s^-1, measured toward dense molecular clouds. In an earlier work we reported an enhancement in the cosmic ray rate of 20 towards HD185418. McCall et al. have reported enhancements of 48 towards zeta Persei based on the observed abundance of H$_{3}^+ while Le Petit et al. found a cosmic ray enhancement of ~10 to be consistent with their models for this same sight line. Here we revisit zeta Persei and perform a detailed calculation using a self-consistent treatment of the hydrogen chemistry, grain physics, energy and ionization balance, and excitation physics. We show that the value of zeta deduced from the H$_{3}^+$ column density in the diffuse region of the sightline depends strongly on the properties of the grains because they remove free electrons and change the hydrogen chemistry. The observations are largely consistent with a cosmic ray enhancement of 40, with several diagnostics indicating higher values. This underscores the importance of a full treatment of grain physics in studies of interstellar chemistry.
△ Less
Submitted 2 November, 2007;
originally announced November 2007.
-
A PDR-Code Comparison Study
Authors:
M. Roellig,
N. P. Abel,
T. Bell,
F. Bensch,
J. Black,
G. J. Ferland,
B. Jonkheid,
I. Kamp,
M. J. Kaufman,
J. Le Bourlot,
F. Le Petit,
R. Meijerink,
O. Morata,
V. Ossenkopf,
E. Roueff,
G. Shaw,
M. Spaans,
A. Sternberg,
J. Stutzki,
W. -F. Thi,
E. F. van Dishoeck,
P. A. M. van Hoof,
S. Viti,
M. G. Wolfire
Abstract:
We present a comparison between independent computer codes, modeling the physics and chemistry of interstellar photon dominated regions (PDRs). Our goal was to understand the mutual differences in the PDR codes and their effects on the physical and chemical structure of the model clouds, and to converge the output of different codes to a common solution. A number of benchmark models have been cr…
▽ More
We present a comparison between independent computer codes, modeling the physics and chemistry of interstellar photon dominated regions (PDRs). Our goal was to understand the mutual differences in the PDR codes and their effects on the physical and chemical structure of the model clouds, and to converge the output of different codes to a common solution. A number of benchmark models have been created, covering low and high gas densities and far ultraviolet intensities. The benchmark models were computed in two ways: one set assuming constant temperatures, thus testing the consistency of the chemical network and photo-processes, and a second set determining the temperature selfconsistently. We investigated the impact of PDR geometry and agreed on the comparison of results from spherical and plane-parallel PDR models. We identified a number of key processes governing the chemical network which have been treated differently in the various codes, and defined a proper common treatment. We established a comprehensive set of reference models for ongoing and future PDR model bench-marking and were able to increase the agreement in model predictions for all benchmark models significantly.
△ Less
Submitted 8 February, 2007;
originally announced February 2007.
-
A Magnetically-Supported Photodissociation Region in M17
Authors:
E. W. Pellegrini,
J. A. Baldwin,
C. L. Brogan,
M. M. Hanson,
N. P. Abel,
G. J. Ferland,
H. B. Nemala,
G. Shaw,
T. H. Troland
Abstract:
The southwestern (SW) part of the Galactic H II region M17 contains an obscured ionization front that is most easily seen at infrared and radio wavelengths. It is nearly edge-on, thus offering an excellent opportunity to study the way in which the gas changes from fully ionized to molecular as radiation from the ionizing stars penetrates into the gas. M17 is also one of the very few H II regions f…
▽ More
The southwestern (SW) part of the Galactic H II region M17 contains an obscured ionization front that is most easily seen at infrared and radio wavelengths. It is nearly edge-on, thus offering an excellent opportunity to study the way in which the gas changes from fully ionized to molecular as radiation from the ionizing stars penetrates into the gas. M17 is also one of the very few H II regions for which the magnetic field strength can be measured in the photodissociation region ( PDR) that forms the interface between the ionized and molecular gas. Here we model an observed line of sight through the gas cloud, including the H+, H0 (PDR), and molecular layers, in a fully self-consistent single calculation. An interesting aspect of the M17 SW bar is that the PDR is very extended. We show that the strong magnetic field that is observed to be present inevitably leads to a very deep PDR, because the structure of the neutral and molecular gas is dominated by magnetic pressure, rather than by gas pressure, as previously had been supposed.We also show that a wide variety of observed facts can be explained if a hydrostatic geometry prevails, in which the gas pressure from an inner X-ray hot bubble and the outward momentum of the stellar radiation field compress the gas and its associated magnetic field in the PDR, as has already been shown to occur in the Orion Nebula. The magnetic field compression may also amplify the local cosmic-ray density. The pressure in the observed magnetic field balances the outward forces, suggesting that the observed geometry is a natural consequence of the formation of a star cluster within a molecular cloud.
△ Less
Submitted 1 August, 2017; v1 submitted 27 November, 2006;
originally announced November 2006.
-
Physical Conditoins in Orion's Veil II: A Multi-Component Study of the Line of Sight Toward the Trapezium
Authors:
N. P. Abel,
G. J. Ferland,
C. R. O'Dell,
G. Shaw,
T. H. Troland
Abstract:
Orion's Veil is an absorbing screen that lies along the line of sight to the Orion H II region. It consists of two or more layers of gas that must lie within a few parsecs of the Trapezium cluster. Our previous work considered the Veil as a whole and found that the magnetic field dominates the energetics of the gas in at least one component. Here we use high-resolution STIS UV spectra that resol…
▽ More
Orion's Veil is an absorbing screen that lies along the line of sight to the Orion H II region. It consists of two or more layers of gas that must lie within a few parsecs of the Trapezium cluster. Our previous work considered the Veil as a whole and found that the magnetic field dominates the energetics of the gas in at least one component. Here we use high-resolution STIS UV spectra that resolve the two velocity components in absorption and determine the conditions in each. We derive a volume hydrogen density, 21 cm spin temperature, turbulent velocity, and kinetic temperature, for each. We combine these estimates with magnetic field measurements to find that magnetic energy significantly dominates turbulent and thermal energies in one component, while the other component is close to equipartition between turbulent and magnetic energies. We observe molecular hydrogen absorption for highly excited v, J levels that are photoexcited by the stellar continuum, and detect blueshifted S III and P III. These ions must arise from ionized gas between the mostly neutral portions of the Veil and the Trapezium and shields the Veil from ionizing radiation. We find that this layer of ionized gas is also responsible for He I absorption in the Veil, which resolves a 40-year-old debate on the origin of He I absorption towards the Trapezium. Finally, we determine that the ionized and mostly atomic layers of the Veil will collide in less than 85,000 years.
△ Less
Submitted 21 December, 2005;
originally announced December 2005.
-
Physical conditions in the ISM towards HD185418
Authors:
Gargi Shaw,
G. J. Ferland,
R. Srianand,
N. P. Abel
Abstract:
We have developed a complete model of the hydrogen molecule as part of the spectral simulation code Cloudy. Our goal is to apply this to spectra of high-redshift star-forming regions where H2 absorption is seen, but where few other details are known, to understand its implication for star formation. The microphysics of H2 is intricate, and it is important to validate these numerical simulations…
▽ More
We have developed a complete model of the hydrogen molecule as part of the spectral simulation code Cloudy. Our goal is to apply this to spectra of high-redshift star-forming regions where H2 absorption is seen, but where few other details are known, to understand its implication for star formation. The microphysics of H2 is intricate, and it is important to validate these numerical simulations in better-understood environments. This paper studies a well-defined line-of-sight through the Galactic interstellar medium (ISM) as a test of the microphysics and methods we use. We present a self-consistent calculation of the observed absorption-line spectrum to derive the physical conditions in the ISM towards HD185418, a line-of-sight with many observables. We deduce density, temperature, local radiation field, cosmic ray ionization rate, chemical composition and compare these conclusions with conditions deduced from analytical calculations. We find a higher density, similar abundances, and require a cosmic ray flux enhanced over the Galactic background value, consistent with enhancements predicted by MHD simulations.
△ Less
Submitted 23 October, 2005;
originally announced October 2005.
-
Molecular hydrogen in the diffuse interstellar medium at high redshift
Authors:
Raghunathan Srianand,
Gargi Shaw,
Gary Ferland,
Patrick Petitjean,
Cedric Ledoux
Abstract:
Physical conditions within DLAs can reveal the star formation history, determine the chemical composition of the associated ISM, and hence document the first steps in the formation of present day galaxies. Here we present calculations that self-consistently determine the gas ionization, level populations (atomic fine-structure levels and rotational levels of H_2), grain physics, and chemistry. W…
▽ More
Physical conditions within DLAs can reveal the star formation history, determine the chemical composition of the associated ISM, and hence document the first steps in the formation of present day galaxies. Here we present calculations that self-consistently determine the gas ionization, level populations (atomic fine-structure levels and rotational levels of H_2), grain physics, and chemistry. We show that for a low-density gas (n< 0.1 cm^-3) the meta-galactic UV background due to quasars is sufficient to maintain H_2 column densities below the detection limit (i.e N(H_2)< 10^14 cm^-2) irrespective of the metallicity and dust content in the gas. Such a gas will have a 21 cm spin temperature in excess of 7000 K and very low C I and C II* column densities for H I column densities typically observed in DLAs. We show that the observed properties of the ~15% per cent of the DLAs that do show detectable H_2 absorption cannot be reproduced with only the quasar dominated meta-galactic UV radiation field. Gas with higher densities (n>10 cm^-3), a moderate radiation field (flux density one to ten times that of the background radiation of the Galactic ISM), the observed range of metallicity and dust-to-gas ratio reproduce all the observed properties of the DLAs that show H_2 absorption lines. This favors the presence of ongoing star formation in DLAs with H_2. The absence of detectable H_2 and C I absorption in a large fraction of DLAs can be explained if they originate either in a low-density gas or in a high-density gas with a large ambient radiation field. The absence of 21 cm absorption and C II* absorption will be consistent with the first possibility. The presence of 21 cm absorption and strong C II* without H_2 and C I absorption will suggest the second alternative. (Abridged)
△ Less
Submitted 23 June, 2005;
originally announced June 2005.
-
The VLT-UVES survey for molecular hydrogen in high-redshift damped Lyman-alpha systems: Physical conditions in the neutral gas
Authors:
Raghunathan Srianand,
Patrick Petitjean,
Cedric Ledoux,
Gary Ferland,
Gargi Shaw
Abstract:
We study the physical conditions in DLAs, using a sample of 33 systems toward 26 QSOs acquired for a recently completed survey of H_2 by Ledoux et al. (2003). We use the column densities of H_2 in different rotational levels, together with those of C I, C I*, C I**, C II* and singly ionized atomic species to discuss the kinetic temperature, the density of hydrogen and the electronic density in t…
▽ More
We study the physical conditions in DLAs, using a sample of 33 systems toward 26 QSOs acquired for a recently completed survey of H_2 by Ledoux et al. (2003). We use the column densities of H_2 in different rotational levels, together with those of C I, C I*, C I**, C II* and singly ionized atomic species to discuss the kinetic temperature, the density of hydrogen and the electronic density in the gas together with the ambient UV radiation field. Detailed comparisons are made between the observed properties in DLAs, the interstellar medium (ISM) of the Galaxy, LMC and SMC. The mean kinetic temperature of the gas corresponding to DLA subcomponents in which H_2 absorption line is detected, derived from the ortho-to-para ratio (153+/-78 K), is higher than that measured in the ISM (77+/-17 K) and the Magellanic clouds (82+/-21 K). Typical pressure in these components (corresponding to T = 100-300 K and n = 10-200 cm^-3), measured using C I fine-structure excitation, are higher than what is measured along ISM sightlines. This is consistent with the corresponding higher values for N(H_2,J=2)/N(H_2,J=0) seen in DLAs. From the column densities of the high-J rotational levels, we derive that the typical radiation field in the H_2 bearing components is of the order of or slightly higher than the mean UV field in the Galactic ISM. Determination of electron density in the gas with H_2 and C I show the ionization rate is similar to that of a cold neutral medium (CNM) in a moderate radiation field. This, together with the fact that we see H_2 in 13-20% of the DLAs, can be used to conclude that DLAs at z>1.9 could contribute as much as 50% star formation rate density seen in Lyman break galaxies (LBGs). (Abridged)
△ Less
Submitted 23 June, 2005;
originally announced June 2005.
-
The H II Region/PDR Connection: Self-Consistent Calculations of Physical Conditions in Star-Forming Regions
Authors:
N. P. Abel,
G. J. Ferland,
G. Shaw,
P. A. M. van Hoof
Abstract:
We have performed a series of calculations designed to reproduce infrared diagnostics used to determine physical conditions in star forming regions. We self-consistently calculate the thermal and chemical structure of an H II region and photodissociation region (PDR) that are in pressure equilibrium. This differs from previous work, which used separate calculations for each gas phase. Our calcul…
▽ More
We have performed a series of calculations designed to reproduce infrared diagnostics used to determine physical conditions in star forming regions. We self-consistently calculate the thermal and chemical structure of an H II region and photodissociation region (PDR) that are in pressure equilibrium. This differs from previous work, which used separate calculations for each gas phase. Our calculations span a wide range of stellar temperatures, gas densities, and ionization parameters. We describe improvements made to the spectral synthesis code Cloudy that made these calculations possible. These include the addition of a molecular network with ~1000 reactions involving 68 molecular species and improved treatment of the grain physics. Data from the Spitzer First Look Survey, along with other archives, are used to derive important physical characteristics of the H II region and PDR. These include stellar temperatures, electron densities, ionization parameters, UV radiation flux, and PDR density. Finally, we calculate the contribution of the H II region to PDR emission line diagnostics, which allows for a more accurate determination of physical conditions in the PDR.
△ Less
Submitted 21 June, 2005;
originally announced June 2005.
-
Molecular Hydrogen in Star-forming regions: implementation of its micro-physics in Cloudy
Authors:
G. Shaw,
G. J. Ferland,
N. P. Abel,
P. C. Stancil,
P. A. M. van Hoof
Abstract:
Much of the baryonic matter in the Universe is in the form of H2 which includes most of the gas in galactic and extragalactic interstellar clouds. Molecular hydrogen plays a significant role in establishing the thermal balance in many astrophysical environments and can be important as a spectral diagnostic of the gas. Modeling and interpretation of observations of such environments requires a qu…
▽ More
Much of the baryonic matter in the Universe is in the form of H2 which includes most of the gas in galactic and extragalactic interstellar clouds. Molecular hydrogen plays a significant role in establishing the thermal balance in many astrophysical environments and can be important as a spectral diagnostic of the gas. Modeling and interpretation of observations of such environments requires a quantitatively complete and accurate treatment of H2. Using this micro-physical model of H2, illustrative calculations of prototypical astrophysical environments are presented. This work forms the foundation for future investigations of these and other environments where H2 is an important constituent.
△ Less
Submitted 22 January, 2005;
originally announced January 2005.
-
Physical Conditions in Orion's Veil
Authors:
N. P. Abel,
C. L. Brogan,
G. J. Ferland,
C. R. O'Dell,
G. Shaw,
T. H. Troland
Abstract:
Orion's veil consists of several layers of largely neutral gas lying between us and the main ionizing stars of the Orion nebula. It is visible in 21cm H I absorption and in optical and UV absorption lines of H I and other species. Toward the Trapezium, the veil has two remarkable properties, high magnetic field (~100 microGauss) and a surprising lack of molecular hydrogen given its total hydroge…
▽ More
Orion's veil consists of several layers of largely neutral gas lying between us and the main ionizing stars of the Orion nebula. It is visible in 21cm H I absorption and in optical and UV absorption lines of H I and other species. Toward the Trapezium, the veil has two remarkable properties, high magnetic field (~100 microGauss) and a surprising lack of molecular hydrogen given its total hydrogen column density. Here we compute photoionization models of the veil to establish its gas density and its distance from the Trapezium. We use a greatly improved model of the hydrogen molecule that determines level populations in ~1e5 rotational/vibrational levels and provides improved estimates of molecular hydrogen destruction via the Lyman-Werner bands. Our best fit photoionization models place the veil 1-3 pc in front of the star at a density of 1e3-1e4 cubic centimeters. Magnetic energy dominates the energy of non-thermal motions in at least one of the 21cm H I velocity components. Therefore, the veil is the first interstellar environment where magnetic dominance appears to exist. We find that the low ratio of molecular to atomic hydrogen (< 1e-4) is a consequence of high UV flux incident upon the veil due to its proximity to the Trapezium stars and the absence of small grains in the region.
△ Less
Submitted 20 April, 2004;
originally announced April 2004.
-
Observational consequences of fine structure line optical depths on infrared spectral diagnostics
Authors:
Nicholas Abel,
Adam Bryant,
Prabodh Dhakal,
Ashley Gale,
Alva Gibson,
William Goddard,
Chad Howard,
Ameya Kolarkar,
Pey Lian Lim,
Gargi Shaw,
Gary Ferland
Abstract:
It has long been known that infrared fine structure lines of abundant ions, like the [O III] 88 micron line, can become optically thick in H II regions under certain high luminosity conditions. This could mitigate their potential as diagnostic tools, especially if the source is too dusty for optical spectroscopy to otherwise determine the system's parameters. We examined a series of photoionizat…
▽ More
It has long been known that infrared fine structure lines of abundant ions, like the [O III] 88 micron line, can become optically thick in H II regions under certain high luminosity conditions. This could mitigate their potential as diagnostic tools, especially if the source is too dusty for optical spectroscopy to otherwise determine the system's parameters. We examined a series of photoionization calculations which were designed to push the nebulae into the limit where many IR lines should be quite optically thick. We find that radiative transfer effects do not significantly change the observed emission line spectrum. This is due to a combination of grain absorption of the hydrogen ionizing continuum and the fact that the correction for stimulated emission in these lines is large. Given these results, and the likelihood that real objects have non-thermal line broadening, it seems unlikely that line optical depth presents a problem in using these lines as diagnostics of the physical conditions or chemical composition.
△ Less
Submitted 20 November, 2002;
originally announced November 2002.