An extension of implicit Monte Carlo diffusion: Multigroup and the difference formulation

Implicit Monte Carlo (IMC) and Implicit Monte Carlo Diffusion (IMD) are approaches to the numerical solution of the equations of radiative transfer. IMD was previously derived and numerically tested on grey, or frequency-integrated problems [1]. In this research, we extend Implicit Monte Carlo Diffusion (IMD) to account for frequency dependence, and we implement the difference formulation [2] as a source manipulation variance reduction technique. We derive the relevant probability distributions and present the frequency dependent IMD algorithm, with and without the difference formulation. The IMD code with and without the difference formulation was tested using both grey and frequency dependent benchmark problems. The Su and Olson semi-analytic Marshak wave benchmark was used to demonstrate the validity of the code for grey problems [3]. The Su and Olson semi-analytic picket fence benchmark was used for the frequency dependent problems [4]. The frequency dependent IMD algorithm reproduces the results of both Su and Olson benchmark problems. Frequency group refinement studies indicate that the computational cost of refining the group structure is likely less than that of group refinement in deterministic solutions of the radiation diffusion methods. Our results show that applying the difference formulation to the IMD algorithm can result in an overall increase in the figure of merit for frequency dependent problems. However, the creation of negatively weighted particles from the difference formulation can cause significant numerical instabilities in regions of the problem with sharp spatial gradients in the solution. An adaptive implementation of the difference formulation may be necessary to focus its use in regions that are at or near thermal equilibrium.     

Voice source characteristics in six premier country singers

Voice source characteristics as derived from inverse filtering were analyzed in 6 country singers' speech and singing. Results showed that the closed quotient varied systematically with vocal loudness, and that glottal compliance (the ratio between transglottal AC volume displacement and subglottal pressure) decreased with increases in fundamental frequency but remained unaffected by vocal loudness. No striking differences were found in source characteristics between speech and singing within subjects. The degree of phonatory press, as judged by a panel of 19 expert listeners, appeared related to the range in which the singer was singing and to the sound pressure level gain from a doubling of subglottal pressure.     

Detection of tobacco smoke deposition by hyperpolarized krypton-83 MRI

Despite the importance of the tobacco smoke particulate matter in the lungs to the etiology of pulmonary disease in cigarette smokers, little is currently known about the spatial distribution of particle deposition or the persistence of the resulting deposits in humans, and no satisfactory technique currently exists to directly observe tobacco smoke condensate in airways. In this proof-of-principle work, hyperpolarized (hp) 83Kr MRI and NMR spectroscopy are introduced as probes for tobacco smoke deposition in porous media. A reduction in the hp-83Kr longitudinal (T1) relaxation of up to 95% under near-ambient humidity, pressure and temperature conditions was observed when the krypton gas was brought into contact with surfaces that had been exposed to cigarette smoke. This smoke-induced acceleration of the 83Kr self-relaxation was observed for model glass surfaces that, in some experiments, were coated with bovine lung surfactant extract. However, a similar effect was not observed with hp-(129)Xe indicating that the 83Kr sensitivity to smoke deposition was not caused by paramagnetic species but rather by quadrupolar relaxation due to high adsorption affinity for the smoke deposits. The 83Kr T1 differences between smoke-treated and untreated surfaces were sufficient to produce a strong contrast in variable flip angle FLASH hp-83Kr MRI, suggesting that hp-83Kr may be a promising contrast agent for in vivo pulmonary MRI.     

Independent measurement of the total active 8B solar neutrino flux using an array of 3He proportional counters at the Sudbury Neutrino Observatory

The Sudbury Neutrino Observatory (SNO) used an array of 3He proportional counters to measure the rate of neutral-current interactions in heavy water and precisely determined the total active (nu_x) 8B solar neutrino flux. This technique is independent of previous methods employed by SNO. The total flux is found to be 5.54_-0.31;+0.33(stat)-0.34+0.36(syst)x10(6) cm(-2) s(-1), in agreement with previous measurements and standard solar models. A global analysis of solar and reactor neutrino results yields Deltam2=7.59_-0.21;+0.19x10(-5) eV2 and theta=34.4_-1.2;+1.3 degrees. The uncertainty on the mixing angle has been reduced from SNO's previous results.     

Intercalation of trioxatriangulenium ion in DNA: binding,electron transfer,x-ray crystallography,and electronic structure

Trioxatriangulenium ion (TOTA(+)) is a flat, somewhat hydrophobic compound that has a low-energy unoccupied molecular orbital. It binds to duplex DNA by intercalation with a preference for G-C base pairs. Irradiation of intercalated TOTA(+) causes charge (radical cation) injection that results in strand cleavage (after piperidine treatment) primarily at GG steps. The X-ray crystal structure of TOTA(+) intercalated in the hexameric duplex d[CGATCG](2) described here reveals that intercalation of TOTA(+) results in an unusually large extension of the helical rise of the DNA and that the orientation of TOTA(+) is sensitive to hydrogen-bonding interactions with backbone atoms of the DNA. Electronic structure calculations reveal no meaningful charge transfer from DNA to TOTA(+) because the lowest unoccupied molecular orbital of TOTA(+), (LUMO)(T), falls in the gap between the highest occupied molecular orbital, (HOMO)(D), and the (LUMO)(D) of the DNA bases. These calculations reveal the importance of backbone, water, and counterion interactions, which shift the energy levels of the bases and the intercalated TOTA(+) orbitals significantly. The calculations also show that the inserted TOTA(+) strongly polarizes the intercalation cavity where a sheet of excess electron density surrounds the TOTA(+).     

0.6~T人体核磁共振成象及其图象的彩色化处理

在0.6~τ超导磁体人体核磁共振成象设备上获得了横断、矢状、冠状各种断层的成象。采用了直接两维傅里叶变换成象方法,应用了自旋回波脉冲序列。由于人体心脏跳动和呼吸运动会造成图象模糊而使用心电延迟触发及呼吸门控方法。作者用一种相位调制光密度分割的方法将原始黑白灰度表示的核磁共振成象进行彩色化处理。提高了原图象低密度梯度区的清晰度,为进行医学图象诊断及判读解译提供了方便度。山了光学系统图,解释了光学信息处理方法。由所处理的人体核磁共振成象图片可以看出,在低密度梯给区更易辨认。参考文献10种。     

How to measure energy dissipation in dynamic mode atomic force microscopy

When studying a mechanical system like an atomic force microscope (AFM) in dynamic mode it is intuitive and instructive to analyse the forces involved in tip–sample interaction. A different but complementary approach is based on analysing the energy that is dissipated when the tip periodically interacts with the sample surface. This method does not require solving the differential equation of motion for the oscillating cantilever, but is based entirely on the analysis of the energy flow in and out of the dynamic system. Therefore the problem of finding a realistic model to describe the tip–sample interaction in terms of non-linear force–distance dependencies and damping effects is omitted. Instead, it is possible to determine the energy dissipated by the tip–sample interaction directly by measuring such quantities as oscillation amplitude, frequency, phase shift and drive amplitude. The method proved to be important when interpreting phase data obtained in tapping mode, but is also applicable to a variety of scanning probe microscopes operating in different dynamic modes. Additional electronics were designed to allow a direct mapping of local energy dissipation while scanning a sample surface. By applying this technique to the cross-section of a polymer blend a material specific contrast was observed.     

A weight-function approach for determining watershed initial conditions for combustion waves

In many generic combustion models, one finds that a combustionwave will develop with a specific wave speed. However, thereare possible initial temperature profiles which do not evolveinto such waves, but rather die out to the ambient temperature.There can exist, in some models, a clear distinction betweenthose initial conditions that do evolve into combustion wavesand those that do not; this is sometimes referred to as thewatershed initial condition. When fuel consumption is consideredto be negligible, analytical methods can be used to obtain theexact watershed. In this paper, we consider the problem of determiningpseudo-watersheds and ascertaining the relationship betweenthese pseudo-watersheds and the exact watersheds. In the processa novel weight-function approach for infinite spatial domainsis developed.