Line formation in solar granulation: II. The photospheric Fe abundance

The solar photospheric Fe abundance has been determined using realistic ab initio 3D, time-dependent, hydrodynamical model atmospheres. The study is based on the excellent agreement between the predicted and observed line profiles directly rather than equivalent width, since the intrinsic Doppler broadening from the convective motions and oscillations provide the necessary non-thermal broadening. Thus, three of the four hotly debated parameters (equivalent widths, microturbulence and damping enhancement factors) in the center of the recent solar Fe abundance dispute regarding FeI lines no longer enter the analysis, leaving the transition probabilities as the main uncertainty. Both FeI (using the samples of lines of both the Oxford and Kiel studies) and FeII lines have been investigated, which give consistent results: log FeI = 7.44 +- 0.05 and log FeII = 7.45 +- 0.10. Also the wings of strong FeI lines return consistent abundances, log FeII = 7.42 +- 0.03, but due to the uncertainties inherent in analyses of strong lines we give this determination lower weight than the results from weak and intermediate strong lines. In view of the recent slight downward revision of the meteoritic Fe abundance log Fe = 7.46 +- 0.01, the agreement between the meteoritic and photospheric values is very good, thus appearingly settling the debate over the photospheric Fe abundance from FeI lines.Comment: Accepted for A&

Radon gas, useful for medical purposes, safely fixed in quartz

Radon gas is enclosed in quartz or glass ampules by subjecting the gas sealed at a low pressure in the ampules to an ionization process. This process is useful for preparing fixed radon sources for radiological treatment of malignancies, without the danger of releasing radioactive gases

An approximate buckling analysis for rectangular orthotropic plates with centrally located cutouts

An approximate analysis for predicting buckling of rectangular orthotropic composite plates with centrally located cutouts is presented. In this analysis, prebuckling and buckling problems are converted from a two-dimensional to a one-dimensional system of linear differential equations with variable coefficients. The conversion is accomplished by expressing the displacements as series with each element containing a trigonometric function of one coordinate and a coefficient that is an arbitrary function of the other coordinate. Ordinary differential equations are then obtained from a variational principle. Analytical results obtained from the approximate analysis are compared with finite element analyses for isotropic plates and for specially orthotropic plates with central circular cutouts of various sizes. Experimental results for the specially orthotropic plates are also presented. In nearly all cases, the approximate analysis predicts the buckling mode shapes correctly and predicts the buckling loads to within a few percent of the finite element and experimental results

Using rewards and penalties to obtain desired subject performance

Operant conditioning procedures, specifically the use of negative reinforcement, in achieving stable learning behavior is described. The critical tracking test (CTT) a method of detecting human operator impairment was tested. A pass level is set for each subject, based on that subject's asymptotic skill level while sober. It is critical that complete training take place before the individualized pass level is set in order that the impairment can be detected. The results provide a more general basis for the application of reward/penalty structures in manual control research

Line formation in solar granulation: I. Fe line shapes, shifts and asymmetries

Realistic ab-initio 3D, radiative-hydrodynamical convection simulations of the solar granulation have been applied to FeI and FeII line formation. In contrast to classical analyses based on hydrostatic 1D model atmospheres the procedure contains no adjustable free parameters but the treatment of the numerical viscosity in the construction of the 3D, time-dependent, inhomogeneous model atmosphere and the elemental abundance in the 3D spectral synthesis. However, the numerical viscosity is introduced purely for numerical stability purposes and is determined from standard hydrodynamical test cases with no adjustments allowed to improve the agreement with the observational constraints from the solar granulation. The non-thermal line broadening is mainly provided by the Doppler shifts arising from the convective flows in the solar photosphere and the solar oscillations. The almost perfect agreement between the predicted temporally and spatially averaged line profiles for weak Fe lines with the observed profiles and the absence of trends in derived abundances with line strengths, seem to imply that the micro- and macroturbulence concepts are obsolete in these 3D analyses. Furthermore, the theoretical line asymmetries and shifts show a very satisfactory agreement with observations with an accuracy of typically 50-100 m/s on an absolute velocity scale. The remaining minor discrepancies point to how the convection simulations can be refined further.Comment: Accepted for A&

Numerical constraints on the model of stochastic excitation of solar-type oscillations

Analyses of a 3D simulation of the upper layers of a solar convective envelope provide constraints on the physical quantities which enter the theoretical formulation of a stochastic excitation model of solar p modes, for instance the convective velocities and the turbulent kinetic energy spectrum. These constraints are then used to compute the acoustic excitation rate for solar p modes, P. The resulting values are found ~5 times larger than the values resulting from a computation in which convective velocities and entropy fluctuations are obtained with a 1D solar envelope model built with the time-dependent, nonlocal Gough (1977) extension of the mixing length formulation for convection (GMLT). This difference is mainly due to the assumed mean anisotropy properties of the velocity field in the excitation region. The 3D simulation suggests much larger horizontal velocities compared to vertical ones than in the 1D GMLT solar model. The values of P obtained with the 3D simulation constraints however are still too small compared with the values inferred from solar observations. Improvements in the description of the turbulent kinetic energy spectrum and its depth dependence yield further increased theoretical values of P which bring them closer to the observations. It is also found that the source of excitation arising from the advection of the turbulent fluctuations of entropy by the turbulent movements contributes ~ 65-75 % to the excitation and therefore remains dominant over the Reynolds stress contribution. The derived theoretical values of P obtained with the 3D simulation constraints remain smaller by a factor ~3 compared with the solar observations. This shows that the stochastic excitation model still needs to be improved.Comment: 11 pages, 9 figures, accepted for publication in A&

Realistic Magnetohydrodynamical Simulation of Solar Local Supergranulation

Three-dimensional numerical simulations of solar surface magnetoconvection using realistic model physics are conducted. The thermal structure of convective motions into the upper radiative layers of the photosphere, the main scales of convective cells and the penetration depths of convection are investigated. We take part of the solar photosphere with size of 60x60 Mm in horizontal direction and by depth 20 Mm from level of the visible solar surface. We use a realistic initial model of the Sun and apply equation of state and opacities of stellar matter. The equations of fully compressible radiation magnetohydrodynamics with dynamical viscosity and gravity are solved. We apply: 1) conservative TVD difference scheme for the magnetohydrodynamics, 2) the diffusion approximation for the radiative transfer, 3) dynamical viscosity from subgrid scale modeling. In simulation we take uniform two-dimesional grid in gorizontal plane and nonuniform grid in vertical direction with number of cells 600x600x204. We use 512 processors with distributed memory multiprocessors on supercomputer MVS-100k in the Joint Computational Centre of the Russian Academy of Sciences.Comment: 6 pages, 5 figures, submitted to the proceedings of the GONG 2008 / SOHO XXI conferenc

The effects of numerical resolution on hydrodynamical surface convection simulations and spectral line formation

The computationally demanding nature of radiative-hydrodynamical simulations of stellar surface convection warrants an investigation of the sensitivity of the convective structure and spectral synthesis to the numerical resolution and dimension of the simulations, which is presented here. With too coarse a resolution the predicted spectral lines tend to be too narrow, reflecting insufficient Doppler broadening from the convective motions, while at the currently highest affordable resolution the line shapes have converged essentially perfectly to the observed profiles. Similar conclusions are drawn from the line asymmetries and shifts. In terms of abundances, weak FeI and FeII lines show a very small dependence (~0.02 dex) while for intermediate strong lines with significant non-thermal broadening the sensitivity increases (~0.10 dex). Problems arise when using 2D convection simulations to describe an inherent 3D phenomenon, which translates to inaccurate atmospheric velocity fields and temperature and pressure structures. In 2D the theoretical line profiles tend to be too shallow and broad compared with the 3D calculations and observations, in particular for intermediate strong lines. In terms of abundances, the 2D results are systematically about 0.1 dex lower than for the 3D case for FeI lines. Furthermore, the predicted line asymmetries and shifts are much inferior in 2D. Given these shortcomings and computing time considerations it is better to use 3D simulations of even modest resolution than high-resolution 2D simulations.Comment: Accepted for A&