Further development of discrete computational techniques for calculation of restricted diffusion propagators in porous media

Magnetic resonance is a well-established tool for structural characterisation of porous media. Features of pore-space morphology can be inferred from NMR diffusion-diffraction plots or the time-dependence of the apparent diffusion coefficient. Diffusion NMR signal attenuation can be computed from the restricted diffusion propagator, which describes the distribution of diffusing particles for a given starting position and diffusion time. We present two techniques for efficient evaluation of restricted diffusion propagators for use in NMR porous-media characterisation. The first is the Lattice Path Count (LPC). Its physical essence is that the restricted diffusion propagator connecting points A and B in time t is proportional to the number of distinct length-t paths from A to B. By using a discrete lattice, the number of such paths can be counted exactly. The second technique is the Markov transition matrix (MTM). The matrix represents the probabilities of jumps between every pair of lattice nodes within a single timestep. The propagator for an arbitrary diffusion time can be calculated as the appropriate matrix power. For periodic geometries, the transition matrix needs to be defined only for a single unit cell. This makes MTM ideally suited for periodic systems. Both LPC and MTM are closely related to existing computational techniques: LPC, to combinatorial techniques; and MTM, to the Fokker-Planck master equation. The relationship between LPC, MTM and other computational techniques is briefly discussed in the paper. Both LPC and MTM perform favourably compared to Monte Carlo sampling, yielding highly accurate and almost noiseless restricted diffusion propagators. Initial tests indicate that their computational performance is comparable to that of finite element methods. Both LPC and MTM can be applied to complicated pore-space geometries with no analytic solution. We discuss the new methods in the context of diffusion propagator calculation in porous materials and model biological tissues.

Estimation of inhaled ultrafine particle surface area dose for urban environments

There is considerable scientific interest in personal exposure to ultrafine particles. Owing to their small size, these particles are able to penetrate deep into the lungs, where they may cause adverse respiratory, pulmonary and cardiovascular health effects. This article presents Bayesian hierarchical models for estimating and comparing inhaled particle surface area in the lung.

Calculation of the speed of walking of individuals with transfemoral amputation fitted with bone-anchored prosthesis

Typically, the walking ability of individuals with a transfemoral amputation (TFA) can be represented by the speed of walking (SofW) obtained in experimental settings. Recent developments in portable kinetic systems allow assessing the level of activity of TFA during actual daily living outside the confined space of a gait lab. Unfortunately, only minimal spatio-temporal characteristics could be extracted from the kinetic data including the cadence and the duration on gait cycles. Therefore, there is a need for a way to use some of these characteristics to assess the instantaneous speed of walking during daily living. The purpose of the study was to compare several methods to determine SofW using minimal spatial gait characteristics.

Enhanced cortisol suppression following administration of low-dose dexamethasone in first-episode psychosis patients

OBJECTIVE Impaired regulation of the hypothalamic-pituitary-adrenal (HPA) axis and hyper-activity of this system have been described in patients with psychosis. Conversely, some psychiatric disorders such as post-traumatic stress disorder (PTSD) are characterised by HPA hypo-activity, which could be related to prior exposure to trauma. This study examined the cortisol response to the administration of low-dose dexamethasone in first-episode psychosis (FEP) patients and its relationship to childhood trauma. METHOD The low-dose (0.25 mg) Dexamethasone Suppression Test (DST) was performed in 21 neuroleptic-naive or minimally treated FEP patients and 20 healthy control participants. Childhood traumatic events were assessed in all participants using the Childhood Trauma Questionnaire (CTQ) and psychiatric symptoms were assessed in patients using standard rating scales. RESULTS FEP patients reported significantly higher rates of childhood trauma compared to controls (p = 0.001) and exhibited lower basal (a.m.) cortisol (p = 0.04) and an increased rate of cortisol hyper-suppression following dexamethasone administration compared to controls (33% (7/21) vs 5% (1/20), respectively; p = 0.04). There were no significant group differences in mean cortisol decline or percent cortisol suppression following the 0.25 mg DST. This study shows for the first time that a subset of patients experiencing their first episode of psychosis display enhanced cortisol suppression. CONCLUSIONS These findings suggest there may be distinct profiles of HPA axis dysfunction in psychosis which should be further explored.

Inhaled particle surface area dose: Stationary vs personal monitoring

The aim of this paper is to determine the suitability of solely stationary measurements for exposure assessment and management applications. For this purpose, quantified inhaled particle surface area (IPSA) doses using both stationary and personal particle exposure monitors were evaluated and compared.

Anatomical modelling of the pregnant radiotherapy patient

This study aimed to take existing anatomical models of pregnant women, currently used for radiation pro-tection and nuclear medicine dose calculations, and adapt them for use in the calculation of fetal dose from external beam radiotherapy (EBRT). The models investigated were ‘KATJA’, which was provided as an MCNPX geometry file, and ‘RPI-P6’, which was provided in a simple, voxelized bina-ry format. In-house code was developed, to convert both mod-els into an `egsphant’ format, suitable for use with DOSXYZnrc. The geometries and densities of the resulting phantoms were evaluated and found to accurately represent the source data. As an example of the use of the phantoms, the delivery of a cranial EBRT treatment was simulated using the BEAMnrc and DOSXYZnrc Monte Carlo codes and the likely out-of-field doses to the fetus in each model was calculated. The results of these calculations showed good agreement (with-in one standard deviation) between the doses calculated in KATJA and PRI-P6, despite substantial anatomical differ-ences between the two models. For a 36 Gy prescription dose to a 233.2 cm3 target in the right brain, the mean doses calcu-lated in a region of interest covering the entire uterus were 1.0 +/- 0.6 mSv for KATJA and 1.3 +/- 0.9 mSv for RPI-P6. This work is expected to lead to more comprehensive studies of EBRT treatment plan design and its effects on fetal dose in the future.

Suitability of diodes for point dose measurements in IMRT/VMAT beams

This study investigated a potential source of inaccuracy for diode measurements in modulated beams; the effect of diode housing asymmetry on measurement results. The possible effects of diode housing asymmetry on the measurement of steep dose gradients were evaluated by measuring 5x5 cm2 beam profiles, with three cylindrical diodes and two commonly used ionization chambers, with each dosimeter positioned in a 3D scanning water tank with its stem perpendicular to the beam axis (horizontal) and parallel to the direction of scanning. The resulting profiles were used to compare the penumbrae measured with the diode stem pointing into (equivalent to a “stem-first” setup) and out of the field (equivalent to a “stem-last” setup) in order to evaluate the effects of dosimeter alignment and thereby identify the effects of dosimeter asymmetry. The stem-first and stem-last orientations resulted in differences of up to 0.2 mm in the measured 20-80% penumbra widths and differences of up to 0.4 mm in the off axis position of the 90% isodose. These differences, which are smaller than previously reported for older model dosimeters, were apparent in the profile results for both diodes and small volume ionization chambers. As an extension to this study, the practical use of all five dosimeters was exemplified by measuring point doses in IMRT test beams. These measurements showed good agreement (within 2%) between the diodes and the small volume ionization chamber, with all of these dosimeters being able to identify a region 3% under-dosage which was not identified by a larger volume (6 mm diameter) ionization chamber. The results of this work should help to remove some of the barriers to the use of diodes for modulated radiotherapy dosimetry in the future.

Morphometric analysis of the thoracic intervertebral foramen osseous anatomy in adolescent idiopathic scoliosis using low dose computerized tomography

Introduction. The dimensions of the thoracic intervertebral foramen in adolescent idiopathic scoliosis (AIS) have not previously been quantified. During posterior approach scoliosis correction surgery pedicle screws may occasionally breach into the foramen. Better understanding of the dimensions of the foramen may be useful in surgical planning. This study describes a reproducible method for measurement of the thoracic foramen in AIS using computerized tomography (CT). Methods. In 23 pre-operative female patients with Lenke 1 type AIS with right side convexity major curves confined to the thoracic spine the foraminal height (FH), foraminal width (FW), pedicle to superior articular process distance (P-SAP) and cross sectional foraminal area (FA) were measured using multiplanar reconstructed CT. Measurements were made at entrance, midpoint and exit of the thoracic foramina from T1/T2 to T11/T12. Results were correlated with potential dependent variables of major curve Cobb Angle measured on X-ray and CT, Age, Weight, Lenke classification subtype, Risser Grade and number of spinal levels in the major curve. Results. The FH, FW, P-SAP and FA dimensions and ratios are all significantly larger on the convexity of the major curve and maximal at or close to the apex. Mean thoracic foraminal dimensions change in a predictable manner relative to position on the major thoracic curve. There was no significant correlation with the measured foraminal dimensions or ratios and the potential dependent variables. The average ratio of convexity to concavity dimensions at the apex foramina for entrance, midpoint and exit respectively are FH (1.50, 1.38, 1.25), FW (1.28, 1.30, 0.98), FA (2.06, 1.84, 1.32), P-SAP (1.61, 1.47, 1.30). Conclusion. Foraminal dimensions of the thoracic spine are significantly affected by AIS. Foraminal dimensions have a predictable convexity to concavity ratio relative to the proximity to the major curve apex. Surgeons should be aware of these anatomical differences during scoliosis correction surgery.

Using low dose computerized tomography to morphometrically analyse the thoracic intervertebral foramen osseous anatomy in adolescent idiopathic scoliosis

INTRODUCTION The dimensions of the thoracic intervertebral foramen in adolescent idiopathic scoliosis (AIS) have not previously been quantified. During posterior approach scoliosis correction surgery pedicle screws may occasionally breach into the foramen. Better understanding of the dimensions of the foramen may be useful in surgical planning. This study describes a reproducible method for measurement of the thoracic foramen in AIS using computerized tomography (CT). METHODS In 23 pre-operative female patients with Lenke 1 type AIS with right side convexity major curves confined to the thoracic spine the foraminal height (FH), foraminal width (FW), pedicle to superior articular process distance (P-SAP) and cross sectional foraminal area (FA) were measured using multiplanar reconstructed CT. Measurements were made at entrance, midpoint and exit of the thoracic foramina from T1/T2 to T11/T12. Results were correlated with potential dependent variables of major curve Cobb Angle measured on X-ray and CT, Age, Weight, Lenke classification subtype, Risser Grade and number of spinal levels in the major curve. RESULTS The FH, FW, P-SAP and FA dimensions and ratios are all significantly larger on the convexity of the major curve and maximal at or close to the apex. Mean thoracic foraminal dimensions change in a predictable manner relative to position on the major thoracic curve. There was no significant correlation with the measured foraminal dimensions or ratios and the potential dependent variables. The average ratio of convexity to concavity dimensions at the apex foramina for entrance, midpoint and exit respectively are FH (1.50, 1.38, 1.25), FW (1.28, 1.30, 0.98), FA (2.06, 1.84, 1.32), P-SAP (1.61, 1.47, 1.30). CONCLUSION Foraminal dimensions of the thoracic spine are significantly affected by AIS. Foraminal dimensions have a predictable convexity to concavity ratio relative to the proximity to the major curve apex. Surgeons should be aware of these anatomical differences during scoliosis correction surgery.

Ultrasound attenuation computed tomography assessment of PAGAT gel dose

Ultrasound has been previously investigated as an alternative readout method for irradiated polymer gel dosimeters, with authors reporting varying dose responses. We extend previous work utilizing a new computed tomography ultrasound scanner comprising of two identical 5 MHz, 128-element linear-array ultrasound transducers, co-axially aligned and submerged in water as a coupling agent, with rotational of the gel dosimeter between the transducers facilitated by a robotic arm. We have investigated the dose-dependence of both ultrasound bulk attenuation and broadband ultrasound attenuation (BUA) for the PAGAT gel dosimeter. The ultrasound bulk attenuation dose sensitivity was found to be 1.46  ±  0.04 dB m −1 Gy −1, being in agreement with previously published results for PAG and MAGIC gels. BUA was also found to be dose dependent and was measured to be 0.024  ±  0.003 dB MHz −1 Gy −1; the advantage of BUA being its insensitivity to frequency-independent attenuation mechanisms including reflection and refraction, thereby minimizing image reconstruction artefacts.

Low-dose curcumin stimulates proliferation, migration and phagocytic activity of olfactory ensheathing cells

One of the promising strategies for neural repair therapies is the transplantation of olfactory ensheathing cells (OECs) which are the glial cells of the olfactory system. We evaluated the effects of curcumin on the behaviour of mouse OECs to determine if it could be of use to further enhance the therapeutic potential of OECs. Curcumin, a natural polyphenol compound found in the spice turmeric, is known for its anti-cancer properties at doses over 10 µM, and often at 50 µM, and it exerts its effects on cancer cells in part by activation of MAP kinases. In contrast, we found that low-dose curcumin (0.5 µM) applied to OECs strikingly modulated the dynamic morphology, increased the rate of migration by up to 4-fold, and promoted significant proliferation of the OECs. Most dramatically, low-dose curcumin stimulated a 10-fold increase in the phagocytic activity of OECs. All of these potently stimulated behavioural characteristics of OECs are favourable for neural repair therapies. Importantly, low-dose curcumin gave a transient activation of p38 kinases, which is in contrast to the high dose curcumin effects on cancer cells in which these MAP kinases tend to undergo prolonged activation. Low-dose curcumin mediated effects on OECs demonstrate cell-type specific stimulation of p38 and ERK kinases. These results constitute the first evidence that low-dose curcumin can modulate the behaviour of olfactory glia into a phenotype potentially more favourable for neural repair and thereby improve the therapeutic use of OECs for neural repair therapies

trans-Dichloro-bis-(arylazoimidazole)palladium(II): Synthesis, Structure, Photoisomerization, and DFT Calculation

Reaction between PdCl2 and 1-alkyl-2-(arylazo)imidazole (RaaiR') or 1-alkyl-2-(naphthyl-alpha/beta-azo)imidazole (alpha/beta-NaiR') under reflux in ethanol has isolated complexes of compositions Pd(RaaiR')(2)Cl-2 (5, 6) and Pd(alpha/beta-NaiR')(2)Cl-2 (7, 8). The X-ray structure determination of one of the molecules, Pd(alpha-NaiBz)(2)Cl-2 (7c), has reported a trans-PdCl2 configuration, and alpha-NaiBz acts as monodentate N(imidazole) donor ligand. The spectral (IR, UV-vis, H-1 NMR) data support the structure. UV light irradiation (light source: Perkin-Elmer LS 55 spectrofluorimeter, Xenon discharge lamp, lambda = 360-396 nm) in a MeCN solution of the complexes shows E-to-Z isomerization of the coordinated azoimidazole unit. The reverse transformation, Z-to-E, is very slow with visible light irradiation. Quantum yields (phi(E-Z)) of E-to-Z isomerization are calculated, and phi is lower than that of the free ligand but comparable with those of Cd(II) and Hg(II) complexes of the same ligand. The Z-to-E isomerization is a thermally induced process. The activation energy (E-a) of Z-to-E isomerization is calculated by controlled-temperature experimentation. cis-Pd(azoimidazole)Cl-2 complexes (azomidazole acts as N(imidazole) and N(azo) Chelating ligand) do not respond upon light irradiation, which supports the idea that the presence of noncoordinated azo-N to make free azo (-N=N-) function is important to reveal photochromic activity. DFT calculation of Pd(alpha-NaiBz)(2)Cl-2 (7c) has suggested that the HOMO of the molecule is constituted of Pd (32%) and Cl (66%), and hence photo excitation may use the energy of Pd and Cl instead of that of the photofunctional -N=N-Ar motif; thus, the rate of photoisomerization and quantum yield decrease versus the free ligand values.

A simple Bayesian decision-theoretic design for dose-finding trials

A flexible and simple Bayesian decision-theoretic design for dose-finding trials is proposed in this paper. In order to reduce the computational burden, we adopt a working model with conjugate priors, which is flexible to fit all monotonic dose-toxicity curves and produces analytic posterior distributions. We also discuss how to use a proper utility function to reflect the interest of the trial. Patients are allocated based on not only the utility function but also the chosen dose selection rule. The most popular dose selection rule is the one-step-look-ahead (OSLA), which selects the best-so-far dose. A more complicated rule, such as the two-step-look-ahead, is theoretically more efficient than the OSLA only when the required distributional assumptions are met, which is, however, often not the case in practice. We carried out extensive simulation studies to evaluate these two dose selection rules and found that OSLA was often more efficient than two-step-look-ahead under the proposed Bayesian structure. Moreover, our simulation results show that the proposed Bayesian method's performance is superior to several popular Bayesian methods and that the negative impact of prior misspecification can be managed in the design stage.

Statistical power calculation and sample size determination for environmental studies with data below detection limits

Power calculation and sample size determination are critical in designing environmental monitoring programs. The traditional approach based on comparing the mean values may become statistically inappropriate and even invalid when substantial proportions of the response values are below the detection limits or censored because strong distributional assumptions have to be made on the censored observations when implementing the traditional procedures. In this paper, we propose a quantile methodology that is robust to outliers and can also handle data with a substantial proportion of below-detection-limit observations without the need of imputing the censored values. As a demonstration, we applied the methods to a nutrient monitoring project, which is a part of the Perth Long-Term Ocean Outlet Monitoring Program. In this example, the sample size required by our quantile methodology is, in fact, smaller than that by the traditional t-test, illustrating the merit of our method.

Decision-theoretic designs for dose-finding clinical trials with multiple outcomes

A decision-theoretic framework is proposed for designing sequential dose-finding trials with multiple outcomes. The optimal strategy is solvable theoretically via backward induction. However, for dose-finding studies involving k doses, the computational complexity is the same as the bandit problem with k-dependent arms, which is computationally prohibitive. We therefore provide two computationally compromised strategies, which is of practical interest as the computational complexity is greatly reduced: one is closely related to the continual reassessment method (CRM), and the other improves CRM and approximates to the optimal strategy better. In particular, we present the framework for phase I/II trials with multiple outcomes. Applications to a pediatric HIV trial and a cancer chemotherapy trial are given to illustrate the proposed approach. Simulation results for the two trials show that the computationally compromised strategy can perform well and appear to be ethical for allocating patients. The proposed framework can provide better approximation to the optimal strategy if more extensive computing is available.