Line Transfer Research Articles

Time-dependent modelling of the molecular line emission from shock waves in outflow sources

We have developed further the technique of time-dependent modelling of magnetohydrodynamic shock waves, with a view to interpreting the molecular line emission from outflow sources. The extensively observed source L1157 B1 was chosen as an exemplar of the application of this technique. The dynamical age of the shock wave model was varied in the range 500 ≤t≤ 5000 yr, with the best fit to the observed line intensities being obtained for t= 1000 yr; this is of the same order as the dynamical age derived by Gueth, Guilloteau & Bachiller from their observations of L1157 B1. The emission line spectra of H2, CO, SiO, ortho- and para-H2O, ortho- and para-NH3, and A- and E-type CH3OH were calculated in parallel with the dynamical and chemical parameters of the model, using the ‘large velocity gradient’ (LVG) approximation to the line transfer problem. We compared the predictions of the models with the observed intensities of emission lines of H2, CO, SiO, ortho-H2O, ortho-NH3 and CH3OH, which include recent Herschel satellite measurements. In the case of SiO, we show (in Appendix A) that extrapolations of the collisional rate coefficients beyond the range of kinetic temperature for which they were originally calculated lead to spurious rotational line intensities and profiles. The computed emission-line spectra of SiO, NH3 and CH3OH are shown to depend on the assumed initial composition of the grain mantles, from whence they are released, by sputtering in the shock wave, into the gas phase. The dependence of the model predictions on the adopted form of the grain-size distribution is investigated in Appendix B; the corresponding integral line intensities are given in tabular form, for a range of C-type shock speeds, in the online Supporting Information.

Rotating discs and non-kinematic double peaks

Double-peaked line profiles are commonly considered a hallmark of rotating disks, with the distance between the peaks a measure of the rotation velocity. However, double-peaks can arise also from radiative transfer effects in optically thick non-rotating sources. Utilizing exact solutions of the line transfer problem we present a detailed study of line emission from geometrically thin Keplerian disks. We derive the conditions for emergence of kinematic double peaks in optically thin and thick disks, and find that it is generally impossible to disentangle the effects of kinematics and line opacity in observed double-peaked profiles. Unless supplemented by additional information, a double-peaked profile alone is not a reliable indicator of a rotating disk. In certain circumstances, triple and quadruple profiles might be better indicators of rotation in optically thick disks.

CPHASH

CPHash is a concurrent hash table for multicore processors. CPHash partitions its table across the caches of cores and uses message passing to transfer lookups/inserts to a partition. CPHash's message passing avoids the need for locks, pipelines batches of asynchronous messages, and packs multiple messages into a single cache line transfer. Experiments on a 80-core machine with 2 hardware threads per core show that CPHash has ~1.6x higher throughput than a hash table implemented using fine-grained locks. An analysis shows that CPHash wins because it experiences fewer cache misses and its cache misses are less expensive, because of less contention for the on-chip interconnect and DRAM. CPServer, a key/value cache server using CPHash, achieves ~5% higher throughput than a key/value cache server that uses a hash table with fine-grained locks, but both achieve better throughput and scalability than memcached. The throughput of CPHash and CPServer also scale near-linearly with the number of cores.

POLARIZED LINE FORMATION IN MULTI-DIMENSIONAL MEDIA. V. EFFECTS OF ANGLE-DEPENDENT PARTIAL FREQUENCY REDISTRIBUTION

The solution of polarized radiative transfer equation with angle-dependent (AD) partial frequency redistribution (PRD) is a challenging problem. Modeling the observed, linearly polarized strong resonance lines in the solar spectrum often requires the solution of the AD line transfer problems in one-dimensional (1D) or multi-dimensional (multi-D) geometries. The purpose of this paper is to develop an understanding of the relative importance of the AD PRD effects and the multi-D transfer effects and particularly their combined influence on the line polarization. This would help in a quantitative analysis of the second solar spectrum (the linearly polarized spectrum of the Sun). We consider both non-magnetic and magnetic media. In this paper we reduce the Stokes vector transfer equation to a simpler form using a Fourier decomposition technique for multi-D media. A fast numerical method is also devised to solve the concerned multi-D transfer problem. The numerical results are presented for a two-dimensional medium with a moderate optical thickness (effectively thin), and are computed for a collisionless frequency redistribution. We show that the AD PRD effects are significant, and can not be ignored in a quantitative fine analysis of the line polarization. These effects are accentuated by the finite dimensionality of the medium (multi-D transfer). The presence of magnetic fields (Hanle effect) modifies the impact of these two effects to a considerable extent.

Transfer Function Estimation of Telephone Lines from Input Impedance Measurements

The ability of a specific telephone line to support a certain digital subscriber line (DSL) service is determined by its downstream and upstream data rates, which are mainly dependent on the line's transfer function. In this way, methods for transfer function estimation play an important role on proper DSL deployment. Most of the existing methods derive the transfer function via line topology identification (LTI) processes. This paper proposes a method which directly estimates the transfer function of telephone lines without any previous LTI process. The results obtained from both simulations and experimental procedure using twisted-pair cables indicate that the proposed method achieves accurate estimations even for lines with bridged-taps.

Simulation of spectrum of doped Ar in indirectly driven implosion target

X-ray spectrum of tracer in ICF implosion target is usually used to infer electron temperature, density and mix of fuel. As the plasma in fuel is in non-local thermodynamic equilibrium (non-LTE), a line transfer code Alpha is developed. Taking the electron temperature and density provided by radiation hydrodynamic as input condition, atomic kinetics and radiation transfer equation are self-consistently solved with the detailed configuration atom (DCA) model. The opacity for specified frequency intervals is obtained, and X-ray spectrum in the image plane is also presented.As application of Alpha program, the spectrum of doped Ar in implosion target on SG Ⅱ laser facility is simulated. The effect of self-absorption of K α line is studied. And it is shown that self-absorption of K α line affects both the intensity and shape of the spectrum, and it should be considered in simulating X-ray spectrum of Ar. And as the spectrum of local thermodynamic equilibrium (LTE) simulation gives large intensity and different shapes compared with the non-LTE results, non-LTE simulation is necessary in such a simulation.

Modeling of Power Line Transfer of Data for Computer Simulation

,

IMPROVED log( gf ) VALUES OF SELECTED LINES IN Mn I AND Mn II FOR ABUNDANCE DETERMINATIONS IN FGK DWARFS AND GIANTS

The goal of the present work is to produce transition probabilities with very low uncertainties for a selected set of multiplets of Mn I and Mn II. Multiplets are chosen based upon their suitability for stellar abundance analysis. We report on new radiative lifetime measurements for 22 levels of Mn I from the e {sup 8} D, z {sup 6} P, z {sup 6} D, z {sup 4} F, e {sup 8} S, and e {sup 6} S terms and six levels of Mn II from the z {sup 5} P and z {sup 7} P terms using time-resolved laser-induced fluorescence on a slow atom/ion beam. New branching fractions for transitions from these levels, measured using a Fourier-transform spectrometer, are reported. When combined, these measurements yield transition probabilities for 47 transitions of Mn I and 15 transitions of Mn II. Comparisons are made to data from the literature and to Russell-Saunders (LS) theory. In keeping with the goal of producing a set of transition probabilities with the highest possible accuracy and precision, we recommend a weighted mean result incorporating our measurements on Mn I and II as well as independent measurements or calculations that we view as reliable and ofmore » a quality similar to ours. In a forthcoming paper, these Mn I/II transition probability data will be utilized to derive the Mn abundance in stars with spectra from both space-based and ground-based facilities over a 4000 A wavelength range. With the employment of a local thermodynamic equilibrium line transfer code, the Mn I/II ionization balance will be determined for stars of different evolutionary states.« less

Observations and radiative transfer modelling of a massive dense cold core in G333

Cold massive cores are one of the earliest manifestations of high mass star formation. Following the detection of SiO emission from G333.125-0.562, a cold massive core, further investigations of the physics, chemistry and dynamics of this object has been carried out. Mopra and NANTEN2 molecular line profile observations, Australia Telescope Compact Array (ATCA) line and continuum emission maps, and Spitzer 24 and 70 \mum images were obtained. These new data further constrain the properties of this prime example of the very early stages of high mass star formation. A model for the source was constructed and compared directly with the molecular line data using a 3D molecular line transfer code - MOLLIE. The ATCA data reveal that G333.125-0.562 is composed of two sources. One of the sources is responsible for the previously detected molecular outflow and is detected in the Spitzer 24 and 70 \mum band data. Turbulent velocity widths are lower than other more active regions of G333 which reflects the younger evolutionary stage and/or lower mass of this core. The molecular line modelling requires abundances of the CO isotopes that strongly imply heavy depletion due to freeze-out of this species onto dust grains. The principal cloud is cold, moderately turbulent and possesses an outflow which indicates the presence of a central driving source. The secondary source could be an even less evolved object as no apparent associations with continuum emissions at (far-)infrared wavelengths.

THE ABUNDANCES OF NEUTRON-CAPTURE SPECIES IN THE VERY METAL-POOR GLOBULAR CLUSTER M15: A UNIFORM ANALYSIS OF RED GIANT BRANCH AND RED HORIZONTAL BRANCH STARS

The globular cluster M15 is unique in its display of star-to-star variations in the neutron-capture elements. Comprehensive abundance surveys have been previously conducted for handfuls of M15 red giant branch (RGB) and red horizontal branch (RHB) stars. No attempt has been made to perform a single, self-consistent analysis of these stars, which exhibit a wide range in atmospheric parameters. In the current effort, a new comparative abundance derivation is presented for three RGB and six RHB members of the cluster. The analysis employs an updated version of the line transfer code MOOG, which now appropriately treats coherent, isotropic scattering. The apparent discrepancy in the previously reported values for the metallicity of M15 RGB and RHB stars is addressed and a resolute disparity of $\Delta(RHB-RGB) \approx 0.1$ dex in the iron abundance was found. The anti-correlative behavior of the light neutron capture elements (Sr, Y, Zr) is clearly demonstrated with both Ba and Eu, standard markers of the {\it s}- and {\it r}-process, respectively. No conclusive detection of Pb was made in the RGB targets. Consequently for the M15 cluster, this suggests that the main component of the {\it s}-process has made a negligible contribution to those elements normally dominated by this process in solar system material. Additionally for the M15 sample, a large Eu abundance spread is confirmed, which is comparable to that of the halo field at the same metallicity. These abundance results are considered in the discussion of the chemical inhomogeneity and nucleosynthetic history of M15.

Spectral line polarization with angle-dependent partial frequency redistribution

Context. The linear polarization of strong resonance lines observed in the solar spectrum is created by the scattering of the photospheric radiation field. This polarization is sensitive to the form of the partial frequency redistribution (PRD) function used in the line radiative transfer equation. Observations have been analyzed until now with angle-averaged PRD functions. With an increase in the polarimetric sensitivity and resolving power of the present-day telescopes, it will become possible to detect finer effects caused by the angle dependence of the PRD functions. Aims. We devise new efficient numerical methods to solve the polarized line transfer equation with angle-dependent PRD, in planeparallel cylindrically symmetrical media. We try to bring out the essential differences between the polarized spectra formed under angle-averaged and the more realistic case of angle-dependent PRD functions. Methods. We use a recently developed Stokes vector decomposition technique to formulate three different iterative methods tailored for angle-dependent PRD functions. Two of them are of the accelerated lambda iteration type, one is based on the core-wing approach, and the other one on the frequency by frequency approach suitably generalized to handle angle-dependent PRD. The third one is based on a series expansion in the mean number of scattering events (Neumann series expansion). Results. We show that all these methods work well on this difficult problem of polarized line formation with angle-dependent PRD. We present several benchmark solutions with isothermal atmospheres to show the performance of the three numerical methods and to analyze the role of the angle-dependent PRD effects. For weak lines, we find no significant effects when the angle-dependence of the PRD functions is taken into account. For strong lines, we find a significant decrease in the polarization, the largest effect occurring in the near wing maxima.

Acoustic waves in the solar atmosphere at high spatial resolution

Aims. We investigate the energy supply of the solar chromosphere by acoustic waves. Methods. A time sequence with high spatial and temporal resolution from the quiet Sun disc centre in Fe i 5434 A (Lande factor g = 0) is analysed. We used models from a numerical simulation of granular convection and apply NLTE spectral line transfer to determine the height of formation. For estimates of acoustic energy flux, we adopted wave propagation with inclinations of the wave vector with respect to the vertical of 0 ◦ ,3 0 ◦ , and 45 ◦ . For a granular and an intergranular model, the transmissions of the atmosphere to highfrequency waves were determined for the three inclination angles. Wavelet and Fourier analyses were performed and the resulting power spectra were corrected for atmospheric transmission. Results. We find waves with periods down to ∼40 s. They occur intermittently in space and time. The velocity signal is formed at a height of 500 km in the granular model and at 620 km in the intergranule. At periods shorter than the acoustic cutoff (∼190 s), ∼40% of the waves occur above granules and ∼60% above intergranules. By adopting vertical propagation, we estimate total fluxes above granules of 2750–3360 W m −2 , and of 910–1 000 W m −2 above intergranules. The weighted average is 1730–2 060 W m −2 .T he estimates of the total fluxes increase by 15% when inclined wave propagation of 45 ◦ is assumed.

Three-dimensional molecular line transfer: a simulated star-forming region

We present the first non-LTE, co-moving frame molecular line calculations of a star-forming cluster simulated using smoothed particle hydrodynamics (SPH), from which we derive high-resolution synthetic observations. We have resampled a particle representation onto an adaptive mesh and self-consistently solved the equations of statistical equilibrium in the co-moving frame, using TORUS, a three-dimensional adaptive mesh refined (AMR) radiative transfer (RT) code. We verified the applicability of the code to the conditions of the SPH simulation by testing its output against other codes. We find that the level populations obtained for optically thick and thin scenarios closely match the ensemble average of the other codes. We have used the code to obtain non-LTE level populations of multiple molecular species throughout the cluster and have created three-dimensional velocity-resolved spatial maps of the emergent intensity. Line profiles of cores traced by N2H+ (1-0) are compared to probes of low density gas, 13CO (1-0) and C18O (1-0), surrounding the cores along the line of sight. The relative differences of the line-centre velocities are shown to be small compared to the velocity dispersion, matching recent observations. We conclude that one cannot reject competitive accretion as a viable theory of star formation based on observed velocity profiles.

FREQUENCY PROPERTIES OF THE NON-HOMOGENEOUS PERIODICAL LINE

Frequency properties of the non-homogeneous periodical line are analyzed. The model of the signal path containing the non-homogeneous periodical line is presented. The programs for calculation of line input impedance and transfer coefficient were developed. Calculation results are presented and analyzed. They show that the stop–band appears due to periodical inhomogeneities. Its central frequency depends on the length of the period of the line. The width of the stop–band and the attenuation in it depend on intensities of periodical inhomogeneities.

Iteration factors in the line formation problem with frequency dependent source function

In this paper some iteration factors families introduced previously to solve the pure line transfer problem are generalized to the case when the back- ground continuum is taken into account. The convergence properties of these factors are discussed when they are applied to the solution of the two-level atom line transfer problem in a constant and variable property media.

The use of iteration factors in the solution of the NLTE line transfer problem—II. Multilevel atom

The iteration factors method (IFM) developed in Paper I [1] (Atanacković-Vukmanović and Simonneau, 1994) to solve the NLTE line transfer problem for a two-level atom model, is extended here to deal with a multilevel atom case. At the beginning of each iteration step, for each line transition, angle and frequency averaged depth-dependent iteration factors are computed from the formal solution of radiative transfer (RT) equation and used to close the system of the RT equation moments, non-linearly coupled with the statistical equilibrium (SE) equations. Non-linear coupling of the atomic level populations and the corresponding line radiation field intensities is tackled in two ways. One is based on the linearization of the equations with respect to the relevant variables, and the other on the use of the old (known from the previous iteration) level populations in the line-opacity-like terms of the SE equations. In both cases the use of quasi-invariant iteration factors provided very fast and accurate solution. The properties of the proposed procedures are investigated in detail by applying them to the solution of the prototype multilevel RT problem of Avrett and Loeser [2], and compared with the properties of some other methods.

Critical line-set configurations to epipolar geometry determination and application to image line transfer

We address the issue of what linear line structures in space constitute critical configuration to three generally positioned cameras. By critical configuration we mean a set of lines in space whose image projections to the cameras do not allow unique determination of the cameras’ relative positions. The work is focused on the domain of a particular but common line structure group – the linear line structures – which include linear ruled surface, linear line congruence, and linear line complex, and more specifically line pencil, point star, and ruled plane. We tackle the issue by looking into the rank of a matrix that is related to the estimation of the trifocal tensor – a quantity that captures the relative geometry of the cameras. Our result is a summary of the following: which of the linear line structure families are critical structures, how critical each of them is, and how many lines need be present minimally in each of the families for the full information of the particular family to be revealed in their image projections. Real image results are presented to illustrate the findings. The findings are important to the validity and stability of algorithms related to structure from motion and projective reconstruction using line correspondences. This article is an extended version of the report ( Zhao and Chung, 2008) that received IAPR’s Piero Zamperoni Best Paper Award in 2008.

EMISSION FROM A YOUNG PROTOSTELLAR OBJECT. I. SIGNATURES OF YOUNG EMBEDDED OUTFLOWS

We examine emission from a young protostellar object (YPO) with three-dimensional ideal MHD simulations and three-dimensional non-local thermodynamic equilibrium (non-LTE) line transfer calculations, and show the first results. To calculate the emission field, we employed a snapshot result of an MHD simulation having young bipolar outflows and a dense protostellar disk (a young circumstellar disk) embedded in an infalling envelope. Synthesized line emission of two molecular species (CO and SiO) show that subthermally excited SiO lines as a high density tracer can provide a better probe of the complex velocity field of a YPO, compared to fully thermalized CO lines. In a YPO at the earliest stage when the outflows are still embedded in the collapsing envelope, infall, rotation and outflow motions have similar speeds. We find that the combined velocity field of these components introduces a great complexity in the line emissions through varying optical thickness and emissivity, such as asymmetric double-horn profiles. We show that the rotation of the outflows, one of the features that characterizes an outflow driven by magneto-centrifugal forces, appears clearly in velocity channel maps and intensity-weighted mean velocity (first moment of velocity) maps. The somewhat irregular morphology of the line emission at this youngest stage is dissimilar to a more evolved object such as young Class 0. High angular resolution observation by e.g., the Atacama Large Millimeter/submillimeter Array (ALMA) telescope can reveal these features.

Hot and cool water in Herbig Ae protoplanetary disks

The spatial origin and detectability of rotational H2O emission lines from Herbig Ae type protoplanetary disks beyond 70 micron is discussed. We use the recently developed disk code ProDiMo to calculate the thermo-chemical structure of a Herbig Ae type disk and apply the non-LTE line radiative transfer code Ratran to predict water line profiles and intensity maps. The model shows three spatially distinct regions in the disk where water concentrations are high, related to different chemical pathways to form the water: (1) a big water reservoir in the deep midplane behind the inner rim, (2) a belt of cold water around the distant icy midplane beyond the snow-line r>20AU, and (3) a layer of irradiated hot water at high altitudes z/r=0.1...0.3, extending from about 1AU to 30AU, where the kinetic gas temperature ranges from 200K to 1500K. Although region 3 contains only little amounts of water vapour (~3x10^-5 M_Earth), we find this warm layer to be almost entirely responsible for the rotational water emission lines, execpt for the 3 lowest excitation lines. Thus, Herschel will probe first and foremost the conditions and radial extension of region 3, where water is predominantly formed via neutral-neutral reactions and the gas is thermally decoupled from the dust T_gas>T_dust. The observations do not allow for a determination of the snow-line, because the snow-line truncates the radial extension of region 1, whereas the lines originate from region 3. Different line transfer approximations (LTE, escape probability, Monte Carlo) are discussed. A non-LTE treatment is required in most cases, and the results obtained with the escape probability method are found to underestimate the Monte Carlo results by 2%...45%.

NLTE models of line-driven stellar winds - III. Influence of X-ray radiation on wind structure of O stars

We study the influence of X-rays on the wind structure of selected O stars. For this purpose we use our non-local thermodynamic equilibrium (NLTE) wind code with inclusion of additional artificial source of X-rays, assumed to originate in the wind shocks. We show that the influence of shock X-ray emission on wind mass-loss rate is relatively small. Wind terminal velocity may be slightly influenced by the presence of strong X-ray sources, especially for stars cooler than Teff < 35 000 K. We discuss the origin of the Lx/L \sim 10^-7 relation. For stars with thick wind this relation can be explained assuming that the cooling time depends on wind density. Stars with optically thin winds exhibiting the weak wind problem display enhanced X-ray emission which may be connected with large shock cooling length. We propose that this effect can explain the weak wind Inclusion of X-rays leads to a better agreement of the model ionization structure with observations. However, we do not found any significant influence of X-rays on Pv ionization fraction implying that the presence of X-rays cannot explain the Pv problem. We study the implications of modified ionization equilibrium due to shock emission on the line transfer in the X-ray region. We conclude that the X-ray line profiles of helium-like ions may be affected by the line absorption within the cool wind.