A GPU-based finite-size pencil beam algorithm with 3D-density correction for radiotherapy dose calculation

Targeting at the development of an accurate and efficient dose calculation engine for online adaptive radiotherapy, we have implemented a finite size pencil beam (FSPB) algorithm with a 3D-density correction method on GPU. This new GPU-based dose engine is built on our previously published ultrafast FSPB computational framework [Gu et al. Phys. Med. Biol. 54 6287-97, 2009]. Dosimetric evaluations against Monte Carlo dose calculations are conducted on 10 IMRT treatment plans (5 head-and-neck cases and 5 lung cases). For all cases, there is improvement with the 3D-density correction over the conventional FSPB algorithm and for most cases the improvement is significant. Regarding the efficiency, because of the appropriate arrangement of memory access and the usage of GPU intrinsic functions, the dose calculation for an IMRT plan can be accomplished well within 1 second (except for one case) with this new GPU-based FSPB algorithm. Compared to the previous GPU-based FSPB algorithm without 3D-density correction, this new algorithm, though slightly sacrificing the computational efficiency (~5-15% lower), has significantly improved the dose calculation accuracy, making it more suitable for online IMRT replanning

The effects of leaflet material properties on the simulated function of regurgitant mitral valves

Advances in three-dimensional imaging provide the ability to construct and analyze finite element (FE) models to evaluate the biomechanical behavior and function of atrioventricular valves. However, while obtaining patient-specific valve geometry is now possible, non-invasive measurement of patient-specific leaflet material properties remains nearly impossible. Both valve geometry and tissue properties play a significant role in governing valve dynamics, leading to the central question of whether clinically relevant insights can be attained from FE analysis of atrioventricular valves without precise knowledge of tissue properties. As such we investigated 1) the influence of tissue extensibility and 2) the effects of constitutive model parameters and leaflet thickness on simulated valve function and mechanics. We compared metrics of valve function (e.g., leaflet coaptation and regurgitant orifice area) and mechanics (e.g., stress and strain) across one normal and three regurgitant mitral valve (MV) models with common mechanisms of regurgitation (annular dilation, leaflet prolapse, leaflet tethering) of both moderate and severe degree. We developed a novel fully-automated approach to accurately quantify regurgitant orifice areas of complex valve geometries. We found that the relative ordering of the mechanical and functional metrics was maintained across a group of valves using material properties up to 15% softer than the representative adult mitral constitutive model. Our findings suggest that FE simulations can be used to qualitatively compare how differences and alterations in valve structure affect relative atrioventricular valve function even in populations where material properties are not precisely known

GPU-based ultra-fast direct aperture optimization for online adaptive radiation therapy

Online adaptive radiation therapy (ART) has great promise to significantly reduce normal tissue toxicity and/or improve tumor control through real-time treatment adaptations based on the current patient anatomy. However, the major technical obstacle for clinical realization of online ART, namely the inability to achieve real-time efficiency in treatment re-planning, has yet to be solved. To overcome this challenge, this paper presents our work on the implementation of an intensity modulated radiation therapy (IMRT) direct aperture optimization (DAO) algorithm on graphics processing unit (GPU) based on our previous work on CPU. We formulate the DAO problem as a large-scale convex programming problem, and use an exact method called column generation approach to deal with its extremely large dimensionality on GPU. Five 9-field prostate and five 5-field head-and-neck IMRT clinical cases with 5\times5 mm2 beamlet size and 2.5\times2.5\times2.5 mm3 voxel size were used to evaluate our algorithm on GPU. It takes only 0.7~2.5 seconds for our implementation to generate optimal treatment plans using 50 MLC apertures on an NVIDIA Tesla C1060 GPU card. Our work has therefore solved a major problem in developing ultra-fast (re-)planning technologies for online ART

Monthly mean large-scale analyses of upper-tropospheric humidity and wind field divergence derived from three geostationary satellites

This paper describes the results from a collaborative study between the European Space Operations Center, the European Organization for the Exploitation of Meteorological Satellites, the National Oceanic and Atmospheric Administration, and the Cooperative Institute for Meteorological Satellite Studies investigating the relationship between satellite-derived monthly mean fields of wind and humidity in the upper troposphere for March 1994. Three geostationary meteorological satellites GOES-7, Meteosat-3, and Meteosat-5 are used to cover an area from roughly 160 deg W to 50 deg E. The wind fields are derived from tracking features in successive images of upper-tropospheric water vapor (WV) as depicted in the 6.5-micron absorption band. The upper-tropospheric relative humidity (UTH) is inferred from measured water vapor radiances with a physical retrieval scheme based on radiative forward calculations. Quantitative information on large-scale circulation patterns in the upper-troposphere is possible with the dense spatial coverage of the WV wind vectors. The monthly mean wind field is used to estimate the large-scale divergence; values range between about-5 x 10(exp -6) and 5 x 10(exp 6)/s when averaged over a scale length of about 1000-2000 km. The spatial patterns of the UTH field and the divergence of the wind field closely resemble one another, suggesting that UTH patterns are principally determined by the large-scale circulation. Since the upper-tropospheric humidity absorbs upwelling radiation from lower-tropospheric levels and therefore contributes significantly to the atmospheric greenhouse effect, this work implies that studies on the climate relevance of water vapor should include three-dimensional modeling of the atmospheric dynamics. The fields of UTH and WV winds are useful parameters for a climate-monitoring system based on satellite data. The results from this 1-month analysis suggest the desirability of further GOES and Meteosat studies to characterize the changes in the upper-tropospheric moisture sources and sinks over the past decade

The Three Dimensional Structure of EUV Accretion Regions in AM Herculis Stars: Modeling of EUV Photometric and Spectroscopic Observations

We have developed a model of the high-energy accretion region for magnetic cataclysmic variables and applied it to {\it Extreme Ultraviolet Explorer} observations of 10 AM Herculis type systems. The major features of the EUV light curves are well described by the model. The light curves exhibit a large variety of features such as eclipses of the accretion region by the secondary star and the accretion stream, and dips caused by material very close to the accretion region. While all the observed features of the light curves are highly dependent on viewing geometry, none of the light curves are consistent with a flat, circular accretion spot whose lightcurve would vary solely from projection effects. The accretion region immediately above the WD surface is a source of EUV radiation caused by either a vertical extent to the accretion spot, or Compton scattering off electrons in the accretion column, or, very likely, both. Our model yields spot sizes averaging 0.06 R$_{WD}$, or $f \sim 1 \times 10^{-3}$ the WD surface area, and average spot heights of 0.023 R$_{WD}$. Spectra extracted during broad dip phases are softer than spectra during the out-of-dip phases. This spectral ratio measurement leads to the conclusion that Compton scattering, some absorption by a warm absorber, geometric effects, an asymmetric temperature structure in the accretion region and an asymmetric density structure of the accretion columnare all important components needed to fully explain the data. Spectra extracted at phases where the accretion spot is hidden behind the limb of the WD, but with the accretion column immediately above the spot still visible, show no evidence of emission features characteristic of a hot plasma.Comment: 30 Pages, 11 Figure

A Survey of Air-to-Ground Propagation Channel Modeling for Unmanned Aerial Vehicles

In recent years, there has been a dramatic increase in the use of unmanned aerial vehicles (UAVs), particularly for small UAVs, due to their affordable prices, ease of availability, and ease of operability. Existing and future applications of UAVs include remote surveillance and monitoring, relief operations, package delivery, and communication backhaul infrastructure. Additionally, UAVs are envisioned as an important component of 5G wireless technology and beyond. The unique application scenarios for UAVs necessitate accurate air-to-ground (AG) propagation channel models for designing and evaluating UAV communication links for control/non-payload as well as payload data transmissions. These AG propagation models have not been investigated in detail when compared to terrestrial propagation models. In this paper, a comprehensive survey is provided on available AG channel measurement campaigns, large and small scale fading channel models, their limitations, and future research directions for UAV communication scenarios

Discovery and Rossiter-McLaughlin Effect of Exoplanet Kepler-8b

We report the discovery and the Rossiter-McLaughlin effect of Kepler-8b, a transiting planet identified by the NASA Kepler Mission. Kepler photometry and Keck-HIRES radial velocities yield the radius and mass of the planet around this F8IV subgiant host star. The planet has a radius RP = 1.419 RJ and a mass, MP = 0.60 MJ, yielding a density of 0.26 g cm^-3, among the lowest density planets known. The orbital period is P = 3.523 days and orbital semima jor axis is 0.0483+0.0006/-0.0012 AU. The star has a large rotational v sin i of 10.5 +/- 0.7 km s^-1 and is relatively faint (V = 13.89 mag), both properties deleterious to precise Doppler measurements. The velocities are indeed noisy, with scatter of 30 m s^-1, but exhibit a period and phase consistent with the planet implied by the photometry. We securely detect the Rossiter-McLaughlin effect, confirming the planet's existence and establishing its orbit as prograde. We measure an inclination between the projected planetary orbital axis and the projected stellar rotation axis of lambda = -26.9 +/- 4.6 deg, indicating a moderate inclination of the planetary orbit. Rossiter-McLaughlin measurements of a large sample of transiting planets from Kepler will provide a statistically robust measure of the true distribution of spin-orbit orientations for hot jupiters in general.Comment: 26 pages, 8 figures, 2 tables; In preparation for submission to the Astrophysical Journa

Estimating survival in patients with gastrointestinal cancers and brain metastases: An update of the graded prognostic assessment for gastrointestinal cancers (GI-GPA).

BackgroundPatients with gastrointestinal cancers and brain metastases (BM) represent a unique and heterogeneous population. Our group previously published the Diagnosis-Specific Graded Prognostic Assessment (DS-GPA) for patients with GI cancers (GI-GPA) (1985-2007, n = 209). The purpose of this study is to update the GI-GPA based on a larger contemporary database.MethodsAn IRB-approved consortium database analysis was performed using a multi-institutional (18), multi-national (3) cohort of 792 patients with gastrointestinal (GI) cancers, with newly-diagnosed BM diagnosed between 1/1/2006 and 12/31/2017. Survival was measured from date of first treatment for BM. Multiple Cox regression was used to select and weight prognostic factors in proportion to their hazard ratios. These factors were incorporated into the updated GI-GPA.ResultsMedian survival (MS) varied widely by primary site and other prognostic factors. Four significant factors (KPS, age, extracranial metastases and number of BM) were used to formulate the updated GI-GPA. Overall MS for this cohort remains poor; 8 months. MS by GPA was 3, 7, 11 and 17 months for GPA 0-1, 1.5-2, 2.5-3.0 and 3.5-4.0, respectively. >30% present in the worst prognostic group (GI-GPA of ≤1.0).ConclusionsBrain metastases are not uncommon in GI cancer patients and MS varies widely among them. This updated GI-GPA index improves our ability to estimate survival for these patients and will be useful for therapy selection, end-of-life decision-making and stratification for future clinical trials. A user-friendly, free, on-line app to calculate the GPA score and estimate survival for an individual patient is available at brainmetgpa.com

On the finite-time splash and splat singularities for the 3-D free-surface Euler equations

We prove that the 3-D free-surface incompressible Euler equations with regular initial geometries and velocity fields have solutions which can form a finite-time "splash" (or "splat") singularity first introduced in [9], wherein the evolving 2-D hypersurface, the moving boundary of the fluid domain, self-intersects at a point (or on surface). Such singularities can occur when the crest of a breaking wave falls unto its trough, or in the study of drop impact upon liquid surfaces. Our approach is founded upon the Lagrangian description of the free-boundary problem, combined with a novel approximation scheme of a finite collection of local coordinate charts; as such we are able to analyze a rather general set of geometries for the evolving 2-D free-surface of the fluid. We do not assume the fluid is irrotational, and as such, our method can be used for a number of other fluid interface problems, including compressible flows, plasmas, as well as the inclusion of surface tension effects.Comment: 40 pages, 5 figures, to appear in Comm. Math. Phys, abstract added for UK RE

TRPV4 receptor as a functional sensory molecule in bladder urothelium: Stretch‐independent, tissue‐specific actions and pathological implications

The newly recognized sensory role of bladder urothelium has generated intense interest in identifying its novel sensory molecules. Sensory receptor TRPV4 may serve such function. However, specific and physiologically relevant tissue actions of TRPV4, stretch‐independent responses, and underlying mechanisms are unknown and its role in human conditions has not been examined. Here we showed TRPV4 expression in guinea‐pig urothelium, suburothelium, and bladder smooth muscle, with urothelial predominance. Selective TRPV4 activation without stretch evoked significant ATP release—key urothelial sensory process, from live mucosa tissue, full‐thickness bladder but not smooth muscle, and sustained muscle contractions. ATP release was mediated by Ca2+‐dependent, pannexin/connexin‐conductive pathway involving protein tyrosine kinase, but independent from vesicular transport and chloride channels. TRPV4 activation generated greater Ca2+ rise than purinergic activation in urothelial cells. There was intrinsic TRPV4 activity without exogeneous stimulus, causing ATP release. TRPV4 contributed to 50% stretch‐induced ATP release. TRPV4 activation also triggered superoxide release. TRPV4 expression was increased with aging. Human bladder mucosa presented similarities to guinea pigs. Overactive bladders exhibited greater TRPV4‐induced ATP release with age dependence. These data provide the first evidence in humans for the key functional role of TRPV4 in urothelium with specific mechanisms and identify TRPV4 up‐regulation in aging and overactive bladders