scholarly journals Resonant-line radiative transfer within power-law density profiles

2020 ◽  
Vol 497 (3) ◽  
pp. 3925-3942 ◽  
Author(s):  
Bing-Xin Lao ◽  
Aaron Smith

ABSTRACT Star-forming regions in galaxies are surrounded by vast reservoirs of gas capable of both emitting and absorbing Lyman α (Lyα) radiation. Observations of Lyα emitters and spatially extended Lyα haloes indeed provide insights into the formation and evolution of galaxies. However, due to the complexity of resonant scattering, only a few analytic solutions are known in the literature. We discuss several idealized but physically motivated scenarios to extend the existing formalism to new analytic solutions, enabling quantitative predictions about the transport and diffusion of Lyα photons. This includes a closed form solution for the radiation field and derived quantities including the emergent flux, peak locations, energy density, average internal spectrum, number of scatters, outward force multiplier, trapping time, and characteristic radius. To verify our predictions, we employ a robust gridless Monte Carlo radiative transfer (GMCRT) method, which is straightforward to incorporate into existing ray tracing codes but requires modifications to opacity-based calculations, including dynamical core-skipping acceleration schemes. We primarily focus on power-law density and emissivity profiles, however both the analytic and numerical methods can be generalized to other cases. Such studies provide additional intuition and understanding regarding the connection between the physical environments and observational signatures of galaxies throughout the Universe.

2017 ◽  
Vol 837 ◽  
pp. 210-229 ◽  
Author(s):  
E. V. Dontsov ◽  
O. Kresse

This study investigates the problem of a semi-infinite hydraulic fracture that propagates steadily in a permeable formation. The fracturing fluid rheology is assumed to follow a power-law behaviour, while the leak-off is modelled by Carter’s model. A non-singular formulation is employed to effectively analyse the problem and to construct a numerical solution. The problem under consideration features three limiting analytic solutions that are associated with dominance of either toughness, leak-off or viscosity. Transitions between all the limiting cases are analysed and the boundaries of applicability of all these limiting solutions are quantified. These bounds allow us to determine the regions in the parametric space, in which these limiting solutions can be used. The problem of a semi-infinite fracture, which is considered in this study, provides the solution for the tip region of a hydraulic fracture and can be used in hydraulic fracturing simulators to facilitate solving the moving fracture boundary problem. To cater for such applications, for which rapid evaluation of the solution is necessary, the last part of this paper constructs an approximate closed form solution for the problem and evaluates its accuracy against the numerical solution inside the parametric space.


2020 ◽  
Vol 494 (2) ◽  
pp. 2823-2838 ◽  
Author(s):  
Ana Trčka ◽  
Maarten Baes ◽  
Peter Camps ◽  
Sharon E Meidt ◽  
James Trayford ◽  
...  

ABSTRACT We compare the spectral energy distributions (SEDs) and inferred physical properties for simulated and observed galaxies at low redshift. We exploit UV-submillimetre mock fluxes of ∼7000 z = 0 galaxies from the EAGLE suite of cosmological simulations, derived using the radiative transfer code skirt. We compare these to ∼800 observed galaxies in the UV-submillimetre range, from the DustPedia sample of nearby galaxies. To derive global properties, we apply the SED fitting code cigale consistently to both data sets, using the same set of ∼80 million models. The results of this comparison reveal overall agreement between the simulations and observations, both in the SEDs and in the derived physical properties, with a number of discrepancies. The optical and far-infrared regimes, and the scaling relations based upon the global emission, diffuse dust, and stellar mass, show high levels of agreement. However, the mid-infrared fluxes of the EAGLE galaxies are overestimated while the far-UV domain is not attenuated enough, compared to the observations. We attribute these discrepancies to a combination of galaxy population differences between the samples and limitations in the subgrid treatment of star-forming regions in the EAGLE-skirt post-processing recipe. Our findings show the importance of detailed radiative transfer calculations and consistent comparison, and provide suggestions for improved numerical models.


2019 ◽  
Vol 629 ◽  
pp. A77
Author(s):  
A. I. Gómez-Ruiz ◽  
A. Gusdorf ◽  
S. Leurini ◽  
K. M. Menten ◽  
S. Takahashi ◽  
...  

Context. OMC-2/3 is one of the nearest embedded cluster-forming regions that includes intermediate-mass protostars at early stages of evolution. A previous CO (3–2) mapping survey towards this region revealed outflow activity related to sources at different evolutionary phases. Aims. The present work presents a study of the warm gas in the high-velocity emission from several outflows found in CO (3–2) emission by previous observations, determines their physical conditions, and makes a comparison with previous results in low-mass star-forming regions. Methods. We used the CHAMP+ heterodyne array on the APEX telescope to map the CO (6–5) and CO (7–6) emission in the OMC-2 FIR 6 and OMC-3 MMS 1-6 regions, and to observe 13CO (6–5) at selected positions. We analyzed these data together with previous CO (3–2) observations. In addition, we mapped the SiO (5–4) emission in OMC-2 FIR 6. Results. The CO (6–5) emission was detected in most of the outflow lobes in the mapped regions, while the CO (7–6) was found mostly in the OMC-3 outflows. In the OMC-3 MMS 5 outflow, a previously undetected extremely high-velocity gas was found in CO (6–5). This extremely high-velocity emission arises from the regions close to the central object MMS 5. Radiative transfer models revealed that the high-velocity gas from MMS 5 outflow consists of gas with nH2 = 104–105 cm−3 and T > 200 K, similar to what is observed in young Class 0 low-mass protostars. For the other outflows, values of nH2 > 104 cm−3 were found. Conclusions. The physical conditions and kinematic properties of the young intermediate-mass outflows presented here are similar to those found in outflows from Class 0 low-mass objects. Due to their excitation requirements, mid − J CO lines are good tracers of extremely high-velocity gas in young outflows likely related to jets.


2020 ◽  
Vol 635 ◽  
pp. A154 ◽  
Author(s):  
L. Michel-Dansac ◽  
J. Blaizot ◽  
T. Garel ◽  
A. Verhamme ◽  
T. Kimm ◽  
...  

Context. Resonant lines are powerful probes of the interstellar and circumgalactic medium of galaxies. Their transfer in gas being a complex process, the interpretation of their observational signatures, either in absorption or in emission, is often not straightforward. Numerical radiative transfer simulations are needed to accurately describe the travel of resonant line photons in real and in frequency space, and to produce realistic mock observations. Aims. This paper introduces RASCAS, a new public 3D radiative transfer code developed to perform the propagation of any resonant line in numerical simulations of astrophysical objects. RASCAS was designed to be easily customisable and to process simulations of arbitrarily large sizes on large supercomputers. Methods. RASCAS performs radiative transfer on an adaptive mesh with an octree structure using the Monte Carlo technique. RASCAS features full MPI parallelisation, domain decomposition, adaptive load-balancing, and a standard peeling algorithm to construct mock observations. The radiative transport of resonant line photons through different mixes of species (e.g. H I, Si II, Mg II, Fe II), including their interaction with dust, is implemented in a modular fashion to allow new transitions to be easily added to the code. Results. RASCAS is very accurate and efficient. It shows perfect scaling up to a minimum of a thousand cores. It has been fully tested against radiative transfer problems with analytic solutions and against various test cases proposed in the literature. Although it was designed to describe accurately the many scatterings of line photons, RASCAS may also be used to propagate photons at any wavelength (e.g. stellar continuum or fluorescent lines), or to cast millions of rays to integrate the optical depths of ionising photons, making it highly versatile.


1984 ◽  
Vol 51 (2) ◽  
pp. 289-293 ◽  
Author(s):  
J. R. Walton

In a previous paper, the title problem was solved for a homogeneous power-law linearly viscoelastic half-plane. Such material has a constant Poisson’s ratio and a shear modulus with a power-law dependence on time. In this paper, the shear modulus is assumed also to have a power-law dependence on depth from the half-plane boundary. As in the earlier paper, only a quasi-static analysis is presented, that is, the enertial terms in the equations of motion are not retained and the indentor is assumed to slide with constant speed. The resulting boundary value problem is reduced to a generalized Abel integral equation. A simple closed-form solution is obtained from which all relevant physical parameters are easily computed.


2014 ◽  
Vol 10 (S309) ◽  
pp. 310-310
Author(s):  
Ilse De Looze ◽  
Jacopo Fritz ◽  
Maarten Baes ◽  

AbstractWe present a new technique developed to model the radiative transfer (RT) effects in nearby face-on galaxies. The face-on perspective provides insight into the star-forming regions and clumpy structure, imposing the need for high-resolution 3D models to recover the asymmetric stellar and dust geometries observed in galaxies. RT modeling of the continuum emission of stars and its interaction with the embedding dust in a galaxy's interstellar medium enables a self-consistent study of the main dust heating mechanisms in galaxies. The main advantage of RT calculations is the non-local character of dust heating that can be addressed by tracing the propagation of stellar radiation through the dusty galaxy medium.


2020 ◽  
Vol 644 ◽  
pp. A122
Author(s):  
D. Dall’Olio ◽  
W. H. T. Vlemmings ◽  
B. Lankhaar ◽  
G. Surcis

Context. Astronomical masers have been effective tools in the study of magnetic fields for years. Observations of the linear and circular polarisation of different maser species allow for the determination of magnetic field properties, such as morphology and strength. In particular, methanol can be used to probe different parts of protostars, such as accretion discs and outflows, since it produces one of the strongest and the most commonly observed masers in massive star-forming regions. Aims. We investigate the polarisation properties of selected methanol maser transitions in light of newly calculated methanol Landé g-factors and in consideration of hyperfine components. We compare our results with previous observations and evaluate the effect of preferred hyperfine pumping and non-Zeeman effects. Methods. We ran simulations using the radiative transfer code, CHAMP, for different magnetic field values, hyperfine components, and pumping efficiencies. Results. We find a dependence between the linear polarisation fraction and the magnetic field strength as well as the hyperfine transitions. The circular polarisation fraction also shows a dependence on the hyperfine transitions. Preferred hyperfine pumping can explain some high levels of linear and circular polarisation and some of the peculiar features seen in the S-shape of observed V-profiles. By comparing a number of methanol maser observations taken from the literature with our simulations, we find that the observed methanol masers are not significantly affected by non-Zeeman effects related to the competition between stimulated emission rates and Zeeman rates, such as the rotation of the symmetry axis. We also consider the relevance of other non-Zeeman effects that are likely to be at work for modest saturation levels, such as the effect of magnetic field changes along the maser path and anisotropic resonant scattering. Conclusions. Our models show that for methanol maser emission, both the linear and circular polarisation percentages depend on which hyperfine transition is masing and the degree to which it is being pumped. Since non-Zeeman effects become more relevant at high values of brightness temperatures, it is important to obtain good estimates of these quantities and the maser beaming angles. Better constraints on the brightness temperature will help improve our understanding of the extent to which non-Zeeman effects contribute to the observed polarisation percentages. In order to detect separate hyperfine components, an intrinsic thermal line width that is significantly smaller than the hyperfine separation is required.


Sign in / Sign up

Export Citation Format

Share Document