Assumption Of Local Thermodynamic Equilibrium Research Articles

Titanium abundances in late-type stars. II. Grid of departure coefficients and application to a sample of 70 000 stars

Rapidly growing datasets from stellar spectroscopic surveys are providing unprecedented opportunities to analyse the chemical evolution history of our Galaxy. However, spectral analysis requires accurate modelling of synthetic stellar spectra for late-type stars, for which the assumption of local thermodynamic equilibrium (LTE) has been shown to be insufficient in many cases. Errors associated with LTE can be particularly large for Ti I which is susceptible to over-ionisation, particularly in metal-poor stars. The aims of this work are to study and quantify the 1D non-LTE effects on titanium abundances across the Hertzsprung-Russell diagram for a large sample of stars. A large grid of departure coefficients, $ were computed on standard MARCS model atmospheres. The grid extends from $3000\ K $ to $8000\ K $ in teff , $-0.5\ to $+5.5\ in logg , and $-5.0$ to $+1.0$ in feh , with non-LTE effects in this grid reaching up to 0.4 dex. This was used to compute abundance corrections that were subsequently applied to the LTE abundances of over 70$\,$000 stars selected from the GALAH survey in addition to a smaller sample of literature Keck data for metal-poor dwarfs. The non-LTE effects grow towards lower feh , lower logg , and higher teff , with a minimum and maximum $ A(Ti)_ Ti\,I $ of $0.02$ and $0.19$ in the GALAH sample. For metal-poor giants, the non-LTE modelling reduces the average ionisation imbalance ( from $-0.11\ to $-0.01\ at Fe/H = $-1.7$, and the enhancement in titanium abundances from Ti I lines results in a tife versus feh trend that more closely resembles the behaviour of Ti II at low metallicities. At higher metallicities, the results are limited by the precision of the GALAH DR3 LTE abundances and the effects are within the errors. For the most metal-poor dwarfs from the Keck sample, the average ionisation imbalance increases from $-0.1$ dex to $+0.2$ dex, a shortcoming that is consistent with previous 1D non-LTE studies and which we speculate could be related to 3D effects. Non-LTE effects on titanium abundances are significant. Neglecting them may alter our understanding of Galactic chemical evolution. We have made our grid of departure coefficients publicly available, with the caveat that the Ti abundances of metal-poor dwarfs need further study in 3D non-LTE.

Lithium Cepheid V708 Car with an unusual chemical composition

Aims. The purpose of this work is to spectroscopically analyse the classical Cepheid V708 Car. A preliminary check of the spectrum of V708 Car showed that this is a lithium-rich supergiant. We also found that V708 Car has an unusual chemical composition in that the abundances of various elements correlate with their condensation temperatures. We tried to find an explanation of this feature, which is unusual for classical Cepheids. Methods. For the spectroscopic analysis, we used methods based on the assumption of local and non-local thermodynamic equilibrium. Results. We determined the fundamental parameters of our program star V708 Car. This long-period Cepheid has a mass of about 12 M⊙. We derived the abundances of 27 chemical elements in this star. They are clearly correlated with their condensation temperature: the higher the condensation temperature, the lower the abundance (there are exceptions for sodium and barium, however). We explain this peculiar chemical composition of the V708 Car atmosphere by the gas–dust separation in the envelope of this star. A similar mechanism leads to the observed peculiarities of the chemical composition of λ Boo, W Vir, and asymptotic giant branch stars.

Thermoelastic effects in Bragg reflectors as a potential bottleneck for XFELs with megahertz repetition rate

Bragg reflectors are essential for beam transport in X-ray free-electron laser (XFEL) facilities. On interaction with Bragg reflectors, a part of the pulse energy will be absorbed, causing the propagation of displacement waves due to rapid thermal expansion. It is suspected that these waves may cause stability problems for XFELs operating with megahertz repetition rates. Here, we experimentally investigate the displacement of a diamond Bragg reflector induced by an optical ultra-violet laser pulse, simulating XFEL pulses with mJ pulse energy, both at room temperature and cryogenic temperatures. Our experiment shows negligible damping of the displacement waves on µs timescales, which could cause disruption for subsequent XFEL pulses. We compare our measurements to a simulation framework based on the assumptions of local thermodynamic equilibrium and classical mechanics, observing reasonable agreement. Our results show that thermoelastic effects are critical for a reliable stability assessment of Bragg reflectors, but are often overlooked.

Measurement of medium-voltage AC air arc temperature and particle number density based on dual-wavelength Moiré deflection technology

AC air arcs are generated in medium-voltage (MV) power systems under the effect of harsh weather conditions, equipment aging, and high penetration of distributed generation, threatening equipment and public safety. The arc current and temperature are low due to the wide application of arc suppression devices. In this scenario, the MV AC air arc does not satisfy the local thermodynamic equilibrium (LTE) condition. In addition, the repeated arcing and extinguishing processes further complicate the arc discharge mechanism, which bring challenges in the modeling and detection of MV AC air arcs. Experimental methods are a direct and efficient approach to determine the properties of arc plasmas. In this study, a dual-wavelength Moiré deflection diagnostic system was established to determine the time evolution of the particle density and radial distribution of the temperature in an MV AC air arc without relying on the LTE assumption. The electron number density and heavy particle number density change transiently during the arc discharge process and change gradient along the radial direction. The heavy particle temperature and electron temperature were then calculated based on the measured particle number density. During the arcing stage, the temperature of the electrons exceeded that of the heavy particles significantly, and the arc deviated from LTE. Finally, the limitations of the traditional single-wavelength Moiré deflection method are analyzed. The classic single-wavelength Moiré deflection method, while capable of estimating heavy particle temperature in plasma, exhibits a significant error in electron density estimation compared to the dual-wavelength Moiré deflection method.

Protonated acetylene in the z = 0.89 molecular absorber toward PKS 1830-211

We report the first interstellar identification of protonated acetylene, C2H3+, a fundamental hydrocarbon, in the z = 0.89 molecular absorber toward the gravitationally lensed quasar PKS 1830–211. The molecular species is identified from clear absorption features corresponding to the 212–101 (rest frequency 494.034 GHz) and 111–000 (431.316 GHz) ground-state transitions of ortho and para forms of C2H3+, respectively, in ALMA spectra toward the southwestern image of PKS 1830-211, where numerous molecules, including other hydrocarbons, have already been detected. From the simple assumption of local thermodynamic equilibrium (LTE) with cosmic microwave background photons and an ortho-to-para ratio of three, we estimate a total C2H3+ column density of 2 × 1012 cm−2 and an abundance of 10−10 compared to H2. However, formation pumping could affect the population of metastable states, yielding a C2H3+ column density higher than the LTE value by a factor of a few. We explore possible routes to the formation of C2H3+, mainly connected to acetylene and methane, and find that the methane route is more likely in PDR environment. As one of the initial hydrocarbon building blocks, C2H3+ is thought to play an important role in astrochemistry, in particular in the formation of more complex organic molecules.

Lithium-rich Cepheid V470 Cas

In this Letter, we report the discovery of a new lithium-rich yellow supergiant star -- the Cepheid V470 Cas -- that has a high lithium abundance, $ A(Li)$ = 3.29. This is highly unusual for supergiant stars. V470 Cas is joining a very select group of lithium-rich Cepheids, with only nine members known to date in our Galaxy. For the analysis of our high-resolution echelle spectrum obtained at the Canada-France-Hawaii Telescope, methods based on both local thermodynamic equilibrium (LTE) and non-local thermodynamic equilibrium (NLTE) assumptions were applied. In particular, the lithium abundance was derived by analysing the equivalent width of the 6707 line. Most lithium-rich Cepheids are located near the blue edge of the instability strip near the bottom of the Hertzsprung--Russell diagram. Their main sequence progenitors are distributed over a mass range of three to five solar masses. Of the nine lithium-rich Cepheids known, six are double-mode pulsators, including our programme Cepheid. It is very likely that the stars of this small group are entering the instability strip for the first time and have not yet passed the red giant phase. Therefore, they have not experienced the large-scale mixing event that could destroy lithium in their convective zones.

Irreversible Thermodynamics of Seawater Evaporation

Under typical marine conditions of about 80% relative humidity, evaporation of water from the ocean is an irreversible process accompanied by entropy production. In this article, equations are derived for the latent heat of irreversible evaporation and the related nonequilibrium entropy balance at the sea surface. To achieve this, linear irreversible thermodynamics is considered in a conceptual ocean evaporation model. The equilibrium thermodynamic standard TEOS-10, the International Thermodynamic Equation of Seawater—2010, is applied to irreversible evaporation under the assumption of local thermodynamic equilibrium. The relevance of local equilibrium conditions for irreversible thermodynamics is briefly explained. New equations are derived for the mass flux of evaporation and for the associated nonequilibrium enthalpies and entropies. The estimated entropy production rate of ocean evaporation amounts to 0.004 W m−2 K−1 as compared with the average terrestrial global entropy production of about 1 W m−2 K−1.

Shocking Sgr B2 (N1) with its own outflow

Aims. Because studies on complex organic molecules (COMs) in high-mass protostellar outflows are sparse, we want to investigate how a powerful outflow, such as that driven by the exciting source of the prominent hot core Sagittarius B2(N1), influences the gas molecular inventory of the surrounding medium with which it interacts. Identifying chemical differences to the hot core unaffected by the outflow and what causes them may help to better understand molecular segregation in other star-forming regions. Methods. We made use of the data taken as part of the 3 mm imaging spectral-line survey Re-exploring Molecular Complexity with ALMA (ReMoCA). We studied the morphology of the emission regions of simple and complex molecules in Sgr B2 (N1). For a selection of twelve COMs and four simpler species, spectra were modelled under the assumption of local thermodynamic equilibrium and population diagrams were derived at two positions, one in each lobe of the outflow. From this analysis, we obtained rotational temperatures and column densities. Abundances were subsequently compared to predictions of astrochemical models and to observations of L1157-B1, a position located in the well-studied outflow of the low-mass protostar L1157, and the source G+0.693-0.027 (G0.693), located in the Sgr B2 molecular cloud complex, which are other regions whose chemistry has been impacted by shocks. Results. Integrated intensity maps of SO and SiO emission reveal a bipolar structure with blue-shifted emission dominantly extending to the south-east from the centre of the hot core and red-shifted emission to the north-west. The morphology of both lobes is complex but can roughly be characterised by an emission component at a larger opening angle, containing most of the emission, and narrower features. The wider-angle component is also prominently observed in emission of S-bearing molecules and species that only contain N as a heavy element, including COMs, but also CH3OH, CH3CHO, HNCO, and NH2CHO. Rotational temperatures are found in the range of ~ 100–200 K. Abundances of N-bearing molecules with respect to CH3OH are enhanced in the outflow component compared to N1S, a position that is not impacted by the outflow. A comparison of molecular abundances with G+0.693–0.027 and L1157-B1 does not show any correlations, suggesting that a shock produced by the outflow impacts Sgr B2 (N1)’s material differently or that the initial conditions were different. Conclusions. The short distance of the analysed outflow positions to the centre of Sgr B2 (N1) lead us to propose a scenario in which a phase of hot-core chemistry (i.e. thermal desorption of ice species and high-temperature gas-phase chemistry) preceded a shock wave. The subsequent compression and further heating of the material resulted in the accelerated destruction of (mainly O-bearing) molecules. Gas-phase formation of cyanides seems to be able to compete with their destruction in the post-shock gas. The abundances of cyanopolyynes are enhanced in the outflow component pointing to (additional) gas-phase formation, possibly incorporating atomic N sourced from ammonia in the post-shock gas. To confirm such a scenario, chemical shock models need to be run that take into account the pre- and post-shock conditions of Sgr B2 (N1). In any case, the results provide new perspectives on shock chemistry and the importance of the environment in which it occurs.

Low-frequency Radio Recombination Lines Away from the Inner Galactic Plane

Diffuse radio recombination lines (RRLs) in the Galaxy are possible foregrounds for redshifted 21 cm experiments. We use EDGES drift scans centered at −26.°7 decl. to characterize diffuse RRLs across the southern sky. We find that RRLs averaged over the large antenna beam (72° × 110°) reach minimum amplitudes of R.A. = 2–6 hr. In this region, the Cα absorption amplitude is 33 ± 11 mK (1σ) averaged over 50–87 MHz (27 ≳ z ≳ 15 for the 21 cm line) and increases strongly as frequency decreases. Cβ and Hα lines are consistent with no detection with amplitudes of 13 ± 14 and 12 ± 10 mK (1σ), respectively. At 108–124.5 MHz (z ≈ 11) in the same region, we find no evidence for carbon or hydrogen lines at the noise level of 3.4 mK (1σ). Conservatively assuming that observed lines come broadly from the diffuse interstellar medium, as opposed to a few compact regions, these amplitudes provide upper limits on the intrinsic diffuse lines. The observations support expectations that Galactic RRLs can be neglected as significant foregrounds for a large region of sky until redshifted 21 cm experiments, particularly those targeting cosmic dawn, move beyond the detection phase. We fit models of the spectral dependence of the lines averaged over the large beam of EDGES, which may contain multiple line sources with possible line blending, and find that including degrees of freedom for expected smooth, frequency-dependent deviations from local thermodynamic equilibrium (LTE) is preferred over simple LTE assumptions for Cα and Hα lines. For Cα we estimate departure coefficients 0.79 < b n β n < 4.5 along the inner Galactic plane and 0 < b n β n < 2.3 away from the inner Galactic plane.

A Low-mass line-rich Core Found in Massive Star-forming Region IRAS 16351-4722

We present Atacama Large Millimeter/submillimeter Array subarcsecond-resolution observations of both continuum and molecular lines at 345 GHz toward the massive star-forming region IRAS 16351-4722 (hereafter I16351). A total of 12 dust cores were detected based on high-spatial-resolution observations of the continuum. Among them, a high-mass core (11.6 M ⊙) and a low-mass core (1.7 M ⊙) show abundant molecular line emissions. 164 molecular transitions from 29 species and 104 molecular transitions from 25 species are identified in the high-mass and low-mass cores, respectively. Complex organic molecules (COMs) such as CH3OH, CH3OCHO, CH3OCH3, C2H5OH, and C2H5CN are detected in the two cores. Under the assumption of local thermodynamic equilibrium, the rotational temperatures and column densities of the COMs are derived with the XCLASS software. The maximum rotation temperature values in the low-mass core and the high-mass core were found to be approximately 130 K and 198 K, respectively. Additionally, the line widths in the high-mass core are larger than those in the low-mass one. The abundant COM line transitions, high gas temperatures, and smaller line widths indicate the presence of a low-mass line-rich core in the massive star formation region for the first time, while the high-mass line-rich core shows hot core properties. When comparing the molecular abundances of CH3OH, CH3OCHO, CH3OCH3, and C2H5OH of the two cores with other hot cores and hot corinos reported in the literature, we further confirm that both a hot core and a low-mass line-rich core are simultaneously detected in I16351.

Non-LTE abundance analysis of A-B stars with low rotational velocities – II. Do A-B stars with normal abundances exist?

ABSTRACT We present chemical composition and fundamental parameters (the effective temperature, surface gravity, and radius) for four sharp-lined A-type stars γ Gem (HD 41705), o Peg (HD 214994), θ Vir (HD 114330), and ν Cap (HD 193432). Our analysis is based on a self-consistent model fitting of high-resolution spectra and spectrophotometric observations over a wide wavelength range. We refined the fundamental parameters of the stars with the sme (Spectroscopy Made Easy) package and verified their accuracy by comparing with the spectral energy distribution and hydrogen line profiles. We found Teff/log g = 9190 ± 130 K/3.56 ± 0.08, 9600 ± 50 K/3.81 ± 0.04, 9600 ± 140 K/3.61 ± 0.12, and 10200 ± 220 K/3.88 ± 0.08 for γ Gem, o Peg, θ Vir, and ν Cap, respectively. Our detailed abundance analysis employs a hybrid technique for spectrum synthesis based on classical model atmospheres calculated in local thermodynamic equilibrium (LTE) assumption together with the non-LTE (NLTE) line formation for 18 of 26 investigated species. Comparison of the abundance patterns observed in A stars of different types (normal A, Am, and Ap) with similar fundamental parameters reveals significant abundance diversity that cannot be explained by the current mechanisms of abundance peculiarity formation in stellar atmospheres. We found a rise of the heavy element (Zn, Sr, Y, Zr, and Ba) abundance excess up to +1 dex with Teff increasing from 7200 to 10 000 K, with a further decrease down to solar value at Teff = 13 000 K, indicating that stars with solar element abundances can be found among late B-type stars.

Non-LTE radio recombination line analysis towards H ii region G34.25+0.14

ABSTRACT The electron temperature $T^{*}_\mathrm{ e}$ of H ii regions in local thermodynamic equilibrium (LTE) is usually calculated based on the ratio of radio combination lines (RRLs) and continuum emission. In this work, the LTE $T^{*}_\mathrm{ e}$ in the H ii region G34.25+0.14, which is especially luminous and regarded as a prime target to investigate the physical properties of an H ii region, is obtained by analysing both RRLs and continuum in L band, and employing the correction factor $b_n^{*}$ to correct for non-LTE conditions. The LTE $T^{*}_\mathrm{ e}$ values in L band are calculated to be $2580\, {\rm K}$ for H i/OH/Recombination line survey (THOR) observations and $3690\, {\rm K}$ for Survey of Ionized Gas in our Galaxy, using the Arecibo telescope (SIGGMA), respectively. By employing a non-LTE correction, we derive Te of 8900 $\, {\rm K}$ and the corresponding $b^{*}_n$ values are calculated to be 1 and close to 1 in X band and C band, respectively. However, $b_n^{*}$ at L band (${\sim} 1.4\, {\rm GHz}$) is calculated as 4.2 and 2.4 for THOR and SIGGMA, suggesting that non-LTE corrections are necessary and significantly affect the results. This work demonstrates that the assumption of LTE to calculate $T^{*}_\mathrm{ e}$ is not valid for the observed transitions in L band, and hence, non-LTE corrections are needed. Furthermore, Ne for G34.25+0.14 is calculated from multiple bands independently, consistent with previous work. This study provides a method which could be used for other H ii regions to investigate their degree of departure from LTE.

Rotational spectrum and interstellar detection of the first torsionally excited state of methylamine

Context. Methylamine (CH3NH2) was first detected in the interstellar medium (ISM) toward Sgr B2 almost 50 years ago by observation of rotational transitions in its torsional ground state. Methylamine exhibits two large-amplitude motions (LAMs), the methyl torsion and amine wagging, which complicate the spectral analysis, especially in excited vibration states. The lack of an accurate model of the two coupled LAMs has also hampered the identification in the ISM of rotational transitions in excited vibrational states. Aims. The aim of this work is to study the terahertz and microwave rotational spectra of methylamine experimentally and theoretically in order to provide a reliable basis for the detection of its rotational transitions in the first torsionally excited state, υt = 1, in the ISM. Methods. The terahertz spectrum of methylamine was measured from 150 to 1520 GHz with the Lille fast scan spectrometer. Using a new “hybrid” Hamiltonian model, we were able to analyze the nuclear quadrupole hyperfine structure and to accurately fit the rotational spectrum of the υt = 1 state of methylamine. We used the imaging spectral line survey ReMoCA performed with the Atacama Large Millimeter/submillimeter Array (ALMA) to search for rotational transitions of methylamine in its first torsionally excited state toward the high-mass star forming region Sgr B2(N). The observed spectra are modeled under the assumption of local thermodynamic equilibrium (LTE). Results. Accurate spectral predictions were obtained for the ground and first excited states of CH3NH2. We report the first interstellar detection of methylamine in the υt = 1 state toward the offset position Sgr B2(N1S) in the hot molecular core Sgr B2(N1). The LTE parameters derived previously from the rotational emission of methylamine in its torsional ground state toward Sgr B2(N1S) yield synthetic spectra of methylamine in the υt = 1 state that are fully consistent with the ALMA spectra and allow us to identify five rotational lines of this state.

Molecular Abundance of the Circumnuclear Region Surrounding an Active Galactic Nucleus in NGC 1068 Based on an Imaging Line Survey in the 3 mm Band with ALMA

We present an imaging molecular line survey in the 3 mm band (85–114 GHz) focused on one of the nearest galaxies with an active galactic nucleus (AGN), NGC 1068, based on observations taken with the Atacama Large Millimeter/submillimeter Array. Distributions of 23 molecular transitions are obtained in the central ∼3 kpc region, including both the circumnuclear disk (CND) and starburst ring (SBR) with 60 and 350 pc resolution. The column densities and relative abundances of all the detected molecules are estimated under the assumption of local thermodynamic equilibrium in the CND and SBR. Then, we discuss the physical and chemical effects of the AGN on molecular abundance corresponding to the observation scale. We found that H13CN, SiO, HCN, and H13CO+ are abundant in the CND relative to the SBR. In contrast, 13CO is more abundant in the SBR. Based on the calculated column density ratios of N(HCN)/N(HCO+), N(HCN)/N(CN), and other molecular distributions, we conclude that the enhancement of HCN in the CND may be due to high-temperature environments resulting from strong shocks, which are traced by the SiO emission. Moreover, the abundance of CN in the CND is significantly lower than the expected value of the model calculations in the region affected by strong radiation. The expected strong X-ray irradiation from the AGN has a relatively lower impact on the molecular abundance in the CND than mechanical feedback.

The Molecular Clouds in a Section of the Third Galactic Quadrant: Observational Properties and Chemical Abundance Ratio between CO and its Isotopologues

We compare the observational properties between 12CO, 13CO, and C18O and summarize the observational parameters based on a 7069-cloud sample from the Milky Way Imaging Scroll Painting CO survey in a section of the third Galactic quadrant. We find that the 13CO angular area () generally increases with that of 12CO (), and the ratio of to is 0.38 by linear fitting. We find that the 12CO and 13CO flux are tightly correlated as with both fluxes calculated within the 13CO-bright region. This indicates that the abundance is a constant as for all samples under the assumption of local thermodynamic equilibrium (LTE). Additionally, we observed that the X-factor is approximately constant in a large sample of molecular clouds. Similarly, we find with both fluxes calculated within the C18O-bright region, which indicates that the abundance ratio stays at the same value of across the molecular clouds under the LTE assumption. The linear relationships of versus and versus hold not only for the 13CO-bright region and C18O-bright region, but also for the entire molecular cloud scale with a lower flux ratio. The abundance ratio inside the clouds shows a strong correlation with the column density and temperature. This indicates that the ratio is dominated by a combination of chemical fractionation, selective dissociation, and the self-shielding effect inside the clouds.

Fourth-order accurate numerical modeling of the multi-fluid plasma equations with adaptive mesh refinement

This work develops a multi-fluid plasma modeling capability in extension of a parallel, adaptive, fourth-order, finite-volume method, computational fluid dynamics algorithm. The governing equations of the multi-fluid plasma model are derived from moments of the Boltzmann equation, under the assumption of local thermodynamic equilibrium. Maxwell's equations describe the electromagnetic field evolution and are coupled to the fluid equations of charged species. The present work concentrates on two, charged, collisionless fluids as the first step, although the algorithm is designed to accommodate more fluids, both charged and neutral, and collisions between fluids; the extension to multiple fluids is straightforward. The numerical algorithm featuring parallel, high-order (≥4 for smooth flows), adaptive mesh refinement is expected to help cope with the inherent stiffness of the physics involved. The standard, fourth-order Runge-Kutta scheme is used to evolve the solution in time. The current study is primarily focused on verification and testing of the new, adaptive algorithm; we demonstrate 4th-order accuracy using problems with exact solutions and benchmark our solutions for common plasma simulations against references from literature. Using the 4th-order scheme, our multi-fluid model is shown to provide accurate solutions to plasma problems using half the cell count as traditional, 2nd-order methods. Furthermore, adaptive mesh refinement is successfully applied to three, plasma test problems, providing a significant reduction in mesh size. In subsequent work, an optimized version of our multi-fluid plasma algorithm will be applied to plasma test cases with more practical configurations, with a rigorous assessment of the improvements to computational efficiency afforded by our adaptive, higher-order scheme.

Local equilibrium in liquid phase shock waves.

We have assessed the assumption of local thermodynamic equilibrium in a shock wave by comparing local thermodynamic data generated with nonequilibrium molecular dynamics (NEMD) simulations with results from corresponding equilibrium simulations. The shock had a Mach number approximately equal to 2 in a Lennard-Jones spline liquid. We found that the local equilibrium assumption holds perfectly well behind the wave front, and is a very good approximation in the front itself. This was supported by calculations of the excess entropy production in the shock front with four different methods that use the local equilibrium assumption in different ways. Two of the methods assume local equilibrium between excess thermodynamic variables by treating the shock as an interface in Gibbs's sense. The other two methods are based on the local equilibrium assumption in a continuous description of the shock front. We show for the shock studied in this work that all four methods give excess entropy productions that are in excellent agreement, with an average variance of 3.5% for the nonequilibrium molecular dynamics (NEMD) simulations. In addition, we solved the Navier-Stokes (N-S) equationsnumerically for the same shock wave using an equilibrium equationof state (EoS) based on a recently developed perturbation theory. The results for the density, pressure, and temperature profiles agree well with the profiles from the NEMD simulations. For instance, the shock waves generated in the two simulations travel with almost the same speed; the average absolute Mach-number deviation of the N-S simulations relative to NEMD is 2.6% in the investigated time interval.

Droplet evaporation in inert gases

A general mixed kinetic-diffusion boundary condition is formulated to account for the out-of-equilibrium kinetics in the Knudsen layer. The mixed boundary condition is used to investigate the problem of quasi-steady evaporation of a droplet in an infinite domain containing inert gases. The widely adopted local thermodynamic equilibrium assumption is found to be the limiting case of infinitely large kinetic Péclet number ${{Pe}_k}$ , and it introduces significant error for ${{Pe}_k} \leqslant O(10)$ , which corresponds to a typical droplet radius $a$ of a few micrometres or smaller. When compared with experimental data, solutions based on the mixed boundary condition, which take into account the temperature jump across the Knudsen layer, better predict the time evolution of $a$ than the classical $D^2$ -law (i.e. $a^2 \propto t$ , where $t$ denotes time). In the slow evaporation limit, an analytical solution is obtained by linearising the full formulation about the equilibrium condition, which shows that the $D^2$ -law can be recovered only in the large ${{Pe}_k}$ limit. For small ${{Pe}_k}$ , where the process is dominated by kinetics, a linear relation, i.e. $a \propto t$ , emerges. When the gas phase density approaches the liquid density (e.g. at high-pressure or low-temperature conditions), the increase in the chemical potential of the liquid phase due to the presence of inert gases needs to be accounted for when formulating the mixed boundary condition, an effect largely ignored in the literature so far.

First detection of radio recombination lines of ions heavier than helium

We report the first detection of radio recombination lines (RRLs) of ions heavier than helium. In a highly sensitive multi-band (12–50 GHz) line survey toward Orion KL with the TianMa 65-m Radio Telescope (TMRT), we successfully detected more than fifteen unblendedαlines of RRLs of singly ionized species (XII) recombined from XIII. TheKa-band (26–35 GHz) spectrum also shows tentative signals ofβlines of ions. The detected lines can be successfully crossmatched with the rest frequencies of RRLs of C IIand/or O II. This finding greatly expands the scope of our understanding of ion RRLs, since prior to this work, only two blended lines (105αand 121α) of He IIhad been reported. Our detected lines can be fitted simultaneously under assumption of local thermodynamic equilibrium (LTE). An abundance of C IIIand O IIIof 8.8×10−4was obtained, avoiding the complexities of optical and infrared observations and the blending of RRLs of atoms. It is consistent with but approaches the upper bound of the value (10−4–10−3) estimated from optical and infrared observations. The effects of dielectronic recombination may contribute to enhancing the level populations even at largen. We expect future observations using radio interferometers could break the degeneracy between C and O, while also helping to reveal the ionization structure and dynamical evolution of various ionized regions.

Interpretation of temperature measurements by the Boltzmann plot method on spatially integrated plasma oxygen spectral lines

In thermal plasma spectroscopy, the Boltzmann plot method is particularly popular for measuring plasma temperature independently of pressure. In order to better understand the results of the Botzmann plot measurements performed with atomic oxygen lines on our thermal plasma, which has a non-negligible thickness, a Python code was developed based on the assumption of local thermodynamic equilibrium (LTE) and the calculated plasma composition and properties. This code allows us to simulate a measurement of the oxygen line intensity resulting from an integration over the plasma thickness in a chosen direction for a given temperature profile. From these a simulated Boltzmann plot gave us a simulated temperature. It resulted that this measurement is governed by the maximum of the temperature profile until the maximum temperature exceeds that of the maximum emissivity of the atomic oxygen lines. Above that temperature a limitation was observed, it is possible to measure higher temperatures but the interpretation of that measurement is difficult. Plasma pressure has a limited effect on this limitation. When no limitation is observed, the temperature measurement from the Boltzmann plot method that we simulated is always at least 90% of the maximum of the temperature profile.