Open Access. Powered by Scholars. Published by Universities.®

Physical Sciences and Mathematics Commons

Open Access. Powered by Scholars. Published by Universities.®

Astrophysics and Astronomy

PDF

University of Kentucky

Physics and Astronomy Faculty Publications

Numerical methods

Articles 1 - 14 of 14

Full-Text Articles in Physical Sciences and Mathematics

Molpop-Cep: An Exact, Fast Code For Multi-Level Systems, Andrés Asensio Ramos, Moshe Elitzur Jul 2018

Molpop-Cep: An Exact, Fast Code For Multi-Level Systems, Andrés Asensio Ramos, Moshe Elitzur

Physics and Astronomy Faculty Publications

We present MOLPOP-CEP, a universal line transfer code that allows the exact calculation of multi-level line emission from a slab with variable physical conditions for any arbitrary atom or molecule for which atomic data exist. The code includes error control to achieve any desired level of accuracy, providing full confidence in its results. Publicly available, MOLPOP-CEP employs our recently developed coupled escape probability (CEP) technique, whose performance exceeds other exact methods by orders of magnitude. The program also offers the option of an approximate solution with different variants of the familiar escape probability method. As an illustration of the MOLPOP-CEP …


Direct Collapse To Supermassive Black Hole Seeds With Radiation Transfer: Cosmological Haloes, Kazem Ardaneh, Yang Luo, Isaac Shlosman, Kentaro Nagamine, John H. Wise, Michael C. Begelman Jun 2018

Direct Collapse To Supermassive Black Hole Seeds With Radiation Transfer: Cosmological Haloes, Kazem Ardaneh, Yang Luo, Isaac Shlosman, Kentaro Nagamine, John H. Wise, Michael C. Begelman

Physics and Astronomy Faculty Publications

We have modelled direct collapse of a primordial gas within dark matter haloes in the presence of radiative transfer, in high-resolution zoom-in simulations in a cosmological framework, down to the formation of the photosphere and the central object. Radiative transfer has been implemented in the flux-limited diffusion (FLD) approximation. Adiabatic models were run for comparison. We find that (a) the FLD flow forms an irregular central structure and does not exhibit fragmentation, contrary to adiabatic flow which forms a thick disc, driving a pair of spiral shocks, subject to Kelvin–Helmholtz shear instability forming fragments; (b) the …


Intermediate-Line Emission In Agns: The Effect Of Prescription Of The Gas Density, T. P. Adhikari, K. Hryniewicz, A. Różańska, B. Czerny, Gary J. Ferland Mar 2018

Intermediate-Line Emission In Agns: The Effect Of Prescription Of The Gas Density, T. P. Adhikari, K. Hryniewicz, A. Różańska, B. Czerny, Gary J. Ferland

Physics and Astronomy Faculty Publications

The requirement of an intermediate-line component in the recently observed spectra of several active galactic nuclei (AGNs) points to the possible existence of a physically separate region between the broad-line region (BLR) and narrow-line region (NLR). In this paper we explore the emission from the intermediate-line region (ILR) by using photoionization simulations of the gas clouds distributed radially from the center of the AGN. The gas clouds span distances typical for the BLR, ILR, and NLR, and the appearance of dust at the sublimation radius is fully taken into account in our model. The structure of a single cloud is …


Direct Collapse To Supermassive Black Hole Seeds With Radiative Transfer: Isolated Halos, Yang Luo, Kazem Ardaneh, Isaac Shlosman, Kentaro Nagamine, John H. Wise, Mitchell C. Begelman Feb 2018

Direct Collapse To Supermassive Black Hole Seeds With Radiative Transfer: Isolated Halos, Yang Luo, Kazem Ardaneh, Isaac Shlosman, Kentaro Nagamine, John H. Wise, Mitchell C. Begelman

Physics and Astronomy Faculty Publications

Direct collapse within dark matter haloes is a promising path to form supermassive black hole seeds at high redshifts. The outer part of this collapse remains optically thin. However, the innermost region of the collapse is expected to become optically thick and requires to follow the radiation field in order to understand its evolution. So far, the adiabatic approximation has been used exclusively for this purpose. We apply radiative transfer in the flux-limited diffusion (FLD) approximation to solve the evolution of coupled gas and radiation for isolated haloes. We find that (1) the photosphere forms at 10−6 pc and …


What Makes The Family Of Barred Disc Galaxies So Rich: Damping Stellar Bars In Spinning Haloes, Angela Collier, Isaac Shlosman, Clayton Heller Feb 2018

What Makes The Family Of Barred Disc Galaxies So Rich: Damping Stellar Bars In Spinning Haloes, Angela Collier, Isaac Shlosman, Clayton Heller

Physics and Astronomy Faculty Publications

We model and analyse the secular evolution of stellar bars in spinning dark matter (DM) haloes with the cosmological spin λ ∼ 0–0.09. Using high-resolution stellar and DM numerical simulations, we focus on angular momentum exchange between stellar discs and DM haloes of various axisymmetric shapes – spherical, oblate, and prolate. We find that stellar bars experience a diverse evolution that is guided by the ability of parent haloes to absorb angular momentum, J, lost by the disc through the action of gravitational torques, resonant and non-resonant. We confirm that dynamical bar instability is accelerated via resonant J-transfer …


The Abundance Discrepancy Factor And T2 In Nebulae: Are Non-Thermal Electrons The Culprits?, Gary J. Ferland, W. J. Henney, C. R. O'Dell, M. Peimbert Oct 2016

The Abundance Discrepancy Factor And T2 In Nebulae: Are Non-Thermal Electrons The Culprits?, Gary J. Ferland, W. J. Henney, C. R. O'Dell, M. Peimbert

Physics and Astronomy Faculty Publications

We discuss recent claims that the free electrons in ionized nebulae may not have a significantly Maxwellian velocity distribution. Supra-thermal electrons, electrons with much more energy than is encountered at electron temperatures found in nebulae, may solve the t2/ADF puzzle, the observations that abundances obtained from recombination and collisionally excited lines do not agree, and that different temperature indicators give different results. These non-Maxwellian electrons can be designated by the kappa formalism. We show that the distance over which heating rates change are much longer than the distance supra-thermal electrons can travel, and that the timescale to thermalize …


The Baryon Cycle At High Redshifts: Effects Of Galactic Winds On Galaxy Evolution In Overdense And Average Regions, Raphael Sadoun, Isaac Shlosman, Jun-Hwan Choi, Emilio Romano-Díaz Sep 2016

The Baryon Cycle At High Redshifts: Effects Of Galactic Winds On Galaxy Evolution In Overdense And Average Regions, Raphael Sadoun, Isaac Shlosman, Jun-Hwan Choi, Emilio Romano-Díaz

Physics and Astronomy Faculty Publications

We employ high-resolution cosmological zoom-in simulations focusing on a high-sigma peak and an average cosmological field at z ~ 6–12 in order to investigate the influence of environment and baryonic feedback on galaxy evolution in the reionization epoch. Strong feedback, e.g., galactic winds, caused by elevated star formation rates (SFRs) is expected to play an important role in this evolution. We compare different outflow prescriptions: (i) constant wind velocity (CW), (ii) variable wind scaling with galaxy properties (VW), and (iii) no outflows (NW). The overdensity leads to accelerated evolution of dark matter and baryonic structures, absent from the "normal" region, …


On The Observability Of Optically Thin Coronal Hyperfine Structure Lines, Marios Chatzikos, Gary J. Ferland, R. J. R. Williams, A. C. Fabian May 2014

On The Observability Of Optically Thin Coronal Hyperfine Structure Lines, Marios Chatzikos, Gary J. Ferland, R. J. R. Williams, A. C. Fabian

Physics and Astronomy Faculty Publications

We present CLOUDY calculations for the intensity of coronal hyperfine lines in various environments. We model indirect collisional and radiative transitions, and quantify the collisionally excited line emissivity in the density-temperature phase space. As an observational aid, we also express the emissivity in units of that in the 0.4-0.7 keV band. For most hyperfine lines, knowledge of the X-ray surface brightness and the plasma temperature is sufficient for rough estimates. We find that the radiation fields of both Perseus A and Virgo A can enhance the populations of highly ionized species within 1 kpc. They can also enhance line emissivity …


Molecular Hydrogen Regulated Star Formation In Cosmological Smoothed Particle Hydrodynamics Simulations, Robert Thompson, Kentaro Nagamine, Jason Jaacks, Jun-Hwan Choi Dec 2013

Molecular Hydrogen Regulated Star Formation In Cosmological Smoothed Particle Hydrodynamics Simulations, Robert Thompson, Kentaro Nagamine, Jason Jaacks, Jun-Hwan Choi

Physics and Astronomy Faculty Publications

Some observations have shown that star formation (SF) correlates tightly with the presence of molecular hydrogen (H2); therefore, it is important to investigate its implication on galaxy formation in a cosmological context. In the present work, we implement a sub-grid model (hereafter H2-SF model) that tracks the H2 mass fraction within our cosmological smoothed particle hydrodynamics code GADGET-3 by using an equilibrium analytic model of Krumholz et al. This model allows us to regulate the SF in our simulation by the local abundance of H2 rather than the total cold gas density, which naturally …


Effects Of External Radiation Fields On Line Emission—Application To Star-Forming Regions, Marios Chatzikos, Gary J. Ferland, R. J. R. Williams, Ryan Porter, P. A. M. Vanhoof Dec 2013

Effects Of External Radiation Fields On Line Emission—Application To Star-Forming Regions, Marios Chatzikos, Gary J. Ferland, R. J. R. Williams, Ryan Porter, P. A. M. Vanhoof

Physics and Astronomy Faculty Publications

A variety of astronomical environments contain clouds irradiated by a combination of isotropic and beamed radiation fields. For example, molecular clouds may be irradiated by the isotropic cosmic microwave background, as well as by a nearby active galactic nucleus. These radiation fields excite atoms and molecules and produce emission in different ways. We revisit the escape probability theorem and derive a novel expression that accounts for the presence of external radiation fields. We show that when the field is isotropic the escape probability is reduced relative to that in the absence of external radiation. This is in agreement with previous …


Detecting The Rapidly Expanding Outer Shell Of The Crab Nebula: Where To Look, Xiang Wang, Gary J. Ferland, J. A. Baldwin, E. D. Loh, C. T. Richardson Sep 2013

Detecting The Rapidly Expanding Outer Shell Of The Crab Nebula: Where To Look, Xiang Wang, Gary J. Ferland, J. A. Baldwin, E. D. Loh, C. T. Richardson

Physics and Astronomy Faculty Publications

We present a range of steady-state photoionization simulations, corresponding to different assumed shell geometries and compositions, of the unseen postulated rapidly expanding outer shell to the Crab Nebula. The properties of the shell are constrained by the mass that must lie within it, and by limits to the intensities of hydrogen recombination lines. In all cases the photoionization models predict very strong emissions from high ionization lines that will not be emitted by the Crab's filaments, alleviating problems with detecting these lines in the presence of light scattered from brighter parts of the Crab. The near-NIR [Ne VI] λ7.652 μ …


Supermassive Black Hole Formation At High Redshifts Via Direct Collapse: Physical Processes In The Early Stage, Jun-Hwan Choi, Isaac Shlosman, Mitchell C. Begelman Aug 2013

Supermassive Black Hole Formation At High Redshifts Via Direct Collapse: Physical Processes In The Early Stage, Jun-Hwan Choi, Isaac Shlosman, Mitchell C. Begelman

Physics and Astronomy Faculty Publications

We use numerical simulations to explore whether direct collapse can lead to the formation of supermassive black hole (SMBH) seeds at high redshifts. Using the adaptive mesh refinement code ENZO, we follow the evolution of gas within slowly tumbling dark matter (DM) halos of Mvir ~ 2 × 108 M and Rvir ~ 1 kpc. For our idealized simulations, we adopt cosmologically motivated DM and baryon density profiles and angular momentum distributions. Our principal goal is to understand how the collapsing flow overcomes the centrifugal barrier and whether it is subject to fragmentation which can potentially …


Expanded Iron Uta Spectra—Probing The Thermal Stability Limits In Agn Clouds, Gary J. Ferland, R. Kisielius, F. P. Keenan, P. A. M. Van Hoof, V. Jonauskas, Matt L. Lykins, R. L. Porter, R. J. R. Williams Apr 2013

Expanded Iron Uta Spectra—Probing The Thermal Stability Limits In Agn Clouds, Gary J. Ferland, R. Kisielius, F. P. Keenan, P. A. M. Van Hoof, V. Jonauskas, Matt L. Lykins, R. L. Porter, R. J. R. Williams

Physics and Astronomy Faculty Publications

The Fe unresolved transition arrays (UTAs) produce prominent features in the ~15-17 Å wavelength range in the spectra of active galactic nuclei (AGNs). Here, we present new calculations of the energies and oscillator strengths of inner-shell lines from Fe XIV, Fe XV, and Fe XVI. These are crucial ions since they are dominant at inflection points in the gas thermal stability curve, and UTA excitation followed by autoionization is an important ionization mechanism for these species. We incorporate these, and data reported in previous papers, into the plasma simulation code Cloudy. This updated physics is subsequently employed to reconsider the …


Two-Component Structure Of The HΒ Broad-Line Region In Quasars. I. Evidence From Spectral Principal Component Analysis, Chen Hu, Jian-Min Wang, Luis C. Ho, Gary J. Ferland, Jack A. Baldwin, Ye Wang Nov 2012

Two-Component Structure Of The HΒ Broad-Line Region In Quasars. I. Evidence From Spectral Principal Component Analysis, Chen Hu, Jian-Min Wang, Luis C. Ho, Gary J. Ferland, Jack A. Baldwin, Ye Wang

Physics and Astronomy Faculty Publications

We report on a spectral principal component analysis (SPCA) of a sample of 816 quasars, selected to have small Fe II velocity shifts with spectral coverage in the rest wavelength range 3500-5500 Å. The sample is explicitly designed to mitigate spurious effects on SPCA induced by Fe II velocity shifts. We improve the algorithm of SPCA in the literature and introduce a new quantity, the fractional-contribution spectrum, that effectively identifies the emission features encoded in each eigenspectrum. The first eigenspectrum clearly records the power-law continuum and very broad Balmer emission lines. Narrow emission lines dominate the second eigenspectrum. The …