Bailey, Daniel W; Wang, Iris Z; Lakeman, Tara; Hales, Lee D; Singh, Anurag K; Podgorsak, Matthew B
2015-01-01
This study compares lung dose distributions for two common techniques of total body photon irradiation (TBI) at extended source-to-surface distance calculated with, and without, tissue density correction (TDC). Lung dose correction factors as a function of lateral thorax separation are approximated for bilateral opposed TBI (supine), similar to those published for anteroposterior-posteroanterior (AP-PA) techniques in AAPM Report 17 (i.e., Task Group 29). 3D treatment plans were created retrospectively for 24 patients treated with bilateral TBI, and for whom CT data had been acquired from the head to the lower leg. These plans included bilateral opposed and AP-PA techniques- each with and without - TDC, using source-to-axis distance of 377 cm and largest possible field size. On average, bilateral TBI requires 40% more monitor units than AP-PA TBI due to increased separation (26% more for 23 MV). Calculation of midline thorax dose without TDC leads to dose underestimation of 17% on average (standard deviation, 4%) for bilateral 6 MV TBI, and 11% on average (standard deviation, 3%) for 23 MV. Lung dose correction factors (CF) are calculated as the ratio of midlung dose (with TDC) to midline thorax dose (without TDC). Bilateral CF generally increases with patient separation, though with high variability due to individual uniqueness of anatomy. Bilateral CF are 5% (standard deviation, 4%) higher than the same corrections calculated for AP-PA TBI in the 6 MV case, and 4% higher (standard deviation, 2%) for 23 MV. The maximum lung dose is much higher with bilateral TBI (up to 40% higher than prescribed, depending on patient anatomy) due to the absence of arm tissue blocking the anterior chest. Dose calculations for bilateral TBI without TDC are incorrect by up to 24% in the thorax for 6 MV and up to 16% for 23 MV. Bilateral lung CF may be calculated as 1.05 times the values published in Table 6 of AAPM Report 17, though a larger patient pool is necessary to better
NASA Astrophysics Data System (ADS)
Disher, Brandon; Hajdok, George; Wang, An; Craig, Jeff; Gaede, Stewart; Battista, Jerry J.
2013-06-01
Cone-beam computed tomography (CBCT) has rapidly become a clinically useful imaging modality for image-guided radiation therapy. Unfortunately, CBCT images of the thorax are susceptible to artefacts due to scattered photons, beam hardening, lag in data acquisition, and respiratory motion during a slow scan. These limitations cause dose errors when CBCT image data are used directly in dose computations for on-line, dose adaptive radiation therapy (DART). The purpose of this work is to assess the magnitude of errors in CBCT numbers (HU), and determine the resultant effects on derived tissue density and computed dose accuracy for stereotactic body radiation therapy (SBRT) of lung cancer. Planning CT (PCT) images of three lung patients were acquired using a Philips multi-slice helical CT simulator, while CBCT images were obtained with a Varian On-Board Imaging system. To account for erroneous CBCT data, three practical correction techniques were tested: (1) conversion of CBCT numbers to electron density using phantoms, (2) replacement of individual CBCT pixel values with bulk CT numbers, averaged from PCT images for tissue regions, and (3) limited replacement of CBCT lung pixels values (LCT) likely to produce artificial lateral electron disequilibrium. For each corrected CBCT data set, lung SBRT dose distributions were computed for a 6 MV volume modulated arc therapy (VMAT) technique within the Philips Pinnacle treatment planning system. The reference prescription dose was set such that 95% of the planning target volume (PTV) received at least 54 Gy (i.e. D95). Further, we used the relative depth dose factor as an a priori index to predict the effects of incorrect low tissue density on computed lung dose in regions of severe electron disequilibrium. CT number profiles from co-registered CBCT and PCT patient lung images revealed many reduced lung pixel values in CBCT data, with some pixels corresponding to vacuum (-1000 HU). Similarly, CBCT data in a plastic lung
Disher, Brandon; Hajdok, George; Wang, An; Craig, Jeff; Gaede, Stewart; Battista, Jerry J
2013-06-21
Cone-beam computed tomography (CBCT) has rapidly become a clinically useful imaging modality for image-guided radiation therapy. Unfortunately, CBCT images of the thorax are susceptible to artefacts due to scattered photons, beam hardening, lag in data acquisition, and respiratory motion during a slow scan. These limitations cause dose errors when CBCT image data are used directly in dose computations for on-line, dose adaptive radiation therapy (DART). The purpose of this work is to assess the magnitude of errors in CBCT numbers (HU), and determine the resultant effects on derived tissue density and computed dose accuracy for stereotactic body radiation therapy (SBRT) of lung cancer. Planning CT (PCT) images of three lung patients were acquired using a Philips multi-slice helical CT simulator, while CBCT images were obtained with a Varian On-Board Imaging system. To account for erroneous CBCT data, three practical correction techniques were tested: (1) conversion of CBCT numbers to electron density using phantoms, (2) replacement of individual CBCT pixel values with bulk CT numbers, averaged from PCT images for tissue regions, and (3) limited replacement of CBCT lung pixels values (LCT) likely to produce artificial lateral electron disequilibrium. For each corrected CBCT data set, lung SBRT dose distributions were computed for a 6 MV volume modulated arc therapy (VMAT) technique within the Philips Pinnacle treatment planning system. The reference prescription dose was set such that 95% of the planning target volume (PTV) received at least 54 Gy (i.e. D95). Further, we used the relative depth dose factor as an a priori index to predict the effects of incorrect low tissue density on computed lung dose in regions of severe electron disequilibrium. CT number profiles from co-registered CBCT and PCT patient lung images revealed many reduced lung pixel values in CBCT data, with some pixels corresponding to vacuum (-1000 HU). Similarly, CBCT data in a plastic lung
Vertical gradients of lung density in healthy supine men.
Millar, A B; Denison, D M
1989-01-01
Computed tomography was used to determine the vertical gradient of physical density in peripheral lung tissue of 12 healthy supine subjects, at total lung capacity and residual volume. At total lung capacity the mean (SD) density of peripheral lung tissue at the level of the mid right atrium was 0.0715 (0.017) g/cm3 and the vertical gradient of density was slight. At residual volume the density of peripheral tissue at the same level was 0.272 (0.067) g/cm3 and the vertical density gradient was curvilinear and more pronounced. Predictions of the gradient at residual volume were made on the basis of the known compliance of the lung and measured effects were attributed to the action of gravity on blood vessel distensibility at total lung capacity. These predictions agreed closely with the actual density gradient measured at residual volume and provide a basis for forecasting the vertical density gradient that would exist in healthy lungs at any degree of inflation. Departure from these gradients would imply local abnormalities of lung compliance, distribution of mechanical stress, or distensibility of vessels. Images PMID:2763259
Improved correction for the tissue fraction effect in lung PET/CT imaging
NASA Astrophysics Data System (ADS)
Holman, Beverley F.; Cuplov, Vesna; Millner, Lynn; Hutton, Brian F.; Maher, Toby M.; Groves, Ashley M.; Thielemans, Kris
2015-09-01
Recently, there has been an increased interest in imaging different pulmonary disorders using PET techniques. Previous work has shown, for static PET/CT, that air content in the lung influences reconstructed image values and that it is vital to correct for this ‘tissue fraction effect’ (TFE). In this paper, we extend this work to include the blood component and also investigate the TFE in dynamic imaging. CT imaging and PET kinetic modelling are used to determine fractional air and blood voxel volumes in six patients with idiopathic pulmonary fibrosis. These values are used to illustrate best and worst case scenarios when interpreting images without correcting for the TFE. In addition, the fractional volumes were used to determine correction factors for the SUV and the kinetic parameters. These were then applied to the patient images. The kinetic parameters K1 and Ki along with the static parameter SUV were all found to be affected by the TFE with both air and blood providing a significant contribution to the errors. Without corrections, errors range from 34-80% in the best case and 29-96% in the worst case. In the patient data, without correcting for the TFE, regions of high density (fibrosis) appeared to have a higher uptake than lower density (normal appearing tissue), however this was reversed after air and blood correction. The proposed correction methods are vital for quantitative and relative accuracy. Without these corrections, images may be misinterpreted.
Method and apparatus for measuring lung density by Compton backscattering
Loo, Billy W.; Goulding, Frederick S.
1991-01-01
The density of the lung of a patient suffering from pulmonary edema is monitored by irradiating the lung by a single collimated beam of monochromatic photons and measuring the energies of photons Compton backscattered from the lung by a single high-resolution, high-purity germanium detector. A compact system geometry and a unique data extraction scheme are utilized to monimize systematic errors due to the presence of the chestwall and multiple scattering.
Method and apparatus for measuring lung density by Compton backscattering
Loo, B.W.; Goulding, F.S.
1988-03-11
The density of the lung of a patient suffering from pulmonary edema is monitored by irradiating the lung by a single collimated beam of monochromatic photons and measuring the energies of photons compton back-scattered from the lung by a single high-resolution, high-purity germanium detector. A compact system geometry and a unique data extraction scheme are utilized to minimize systematic errors due to the presence of the chestwall and multiple scattering. 11 figs., 1 tab.
Regional differences in alveolar density in the human lung are related to lung height.
McDonough, John E; Knudsen, Lars; Wright, Alexander C; Elliott, W Mark; Ochs, Matthias; Hogg, James C
2015-06-01
The gravity-dependent pleural pressure gradient within the thorax produces regional differences in lung inflation that have a profound effect on the distribution of ventilation within the lung. This study examines the hypothesis that gravitationally induced differences in stress within the thorax also influence alveolar density in terms of the number of alveoli contained per unit volume of lung. To test this hypothesis, we measured the number of alveoli within known volumes of lung located at regular intervals between the apex and base of four normal adult human lungs that were rapidly frozen at a constant transpulmonary pressure, and used microcomputed tomographic imaging to measure alveolar density (number alveoli/mm3) at regular intervals between the lung apex and base. These results show that at total lung capacity, alveolar density in the lung apex is 31.6 ± 3.4 alveoli/mm3, with 15 ± 6% of parenchymal tissue consisting of alveolar duct. The base of the lung had an alveolar density of 21.2 ± 1.6 alveoli/mm3 and alveolar duct volume fraction of 29 ± 6%. The difference in alveolar density can be negated by factoring in the effects of alveolar compression due to the pleural pressure gradient at the base of the lung in vivo and at functional residual capacity. PMID:25882386
Association Between RT-Induced Changes in Lung Tissue Density and Global Lung Function
Ma Jinli; Zhang Junan; Zhou Sumin; Hubbs, Jessica L.; Foltz, Rodney J.; Hollis, Donna R.; Light, Kim L.; Wong, Terence Z.; Kelsey, Christopher R.; Marks, Lawrence B.
2009-07-01
Purpose: To assess the association between radiotherapy (RT)-induced changes in computed tomography (CT)-defined lung tissue density and pulmonary function tests (PFTs). Methods and Materials: Patients undergoing incidental partial lung RT were prospectively assessed for global (PFTs) and regional (CT and single photon emission CT [SPECT]) lung function before and, serially, after RT. The percent reductions in the PFT and the average changes in lung density were compared (Pearson correlations) in the overall group and subgroups stratified according to various clinical factors. Comparisons were also made between the CT- and SPECT-based computations using the Mann-Whitney U test. Results: Between 1991 and 2004, 343 patients were enrolled in this study. Of these, 111 patients had a total of 203 concurrent post-RT evaluations of changes in lung density and PFTs available for the analyses, and 81 patients had a total of 141 concurrent post-RT SPECT images. The average increases in lung density were related to the percent reductions in the PFTs, albeit with modest correlation coefficients (range, 0.20-0.43). The analyses also indicated that the association between lung density and PFT changes is essentially equivalent to the corresponding association with SPECT-defined lung perfusion. Conclusion: We found a weak quantitative association between the degree of increase in lung density as defined by CT and the percent reduction in the PFTs.
Adiabatic corrections to density functional theory energies and wave functions.
Mohallem, José R; Coura, Thiago de O; Diniz, Leonardo G; de Castro, Gustavo; Assafrão, Denise; Heine, Thomas
2008-09-25
The adiabatic finite-nuclear-mass-correction (FNMC) to the electronic energies and wave functions of atoms and molecules is formulated for density-functional theory and implemented in the deMon code. The approach is tested for a series of local and gradient corrected density functionals, using MP2 results and diagonal-Born-Oppenheimer corrections from the literature for comparison. In the evaluation of absolute energy corrections of nonorganic molecules the LDA PZ81 functional works surprisingly better than the others. For organic molecules the GGA BLYP functional has the best performance. FNMC with GGA functionals, mainly BLYP, show a good performance in the evaluation of relative corrections, except for nonorganic molecules containing H atoms. The PW86 functional stands out with the best evaluation of the barrier of linearity of H2O and the isotopic dipole moment of HDO. In general, DFT functionals display an accuracy superior than the common belief and because the corrections are based on a change of the electronic kinetic energy they are here ranked in a new appropriate way. The approach is applied to obtain the adiabatic correction for full atomization of alcanes C(n)H(2n+2), n = 4-10. The barrier of 1 mHartree is approached for adiabatic corrections, justifying its insertion into DFT. PMID:18537228
Ions in solution: Density corrected density functional theory (DC-DFT)
Kim, Min-Cheol; Sim, Eunji; Burke, Kieron
2014-05-14
Standard density functional approximations often give questionable results for odd-electron radical complexes, with the error typically attributed to self-interaction. In density corrected density functional theory (DC-DFT), certain classes of density functional theory calculations are significantly improved by using densities more accurate than the self-consistent densities. We discuss how to identify such cases, and how DC-DFT applies more generally. To illustrate, we calculate potential energy surfaces of HO·Cl{sup −} and HO·H{sub 2}O complexes using various common approximate functionals, with and without this density correction. Commonly used approximations yield wrongly shaped surfaces and/or incorrect minima when calculated self consistently, while yielding almost identical shapes and minima when density corrected. This improvement is retained even in the presence of implicit solvent.
Self-interaction corrections in density functional theory
Tsuneda, Takao; Hirao, Kimihiko
2014-05-14
Self-interaction corrections for Kohn-Sham density functional theory are reviewed for their physical meanings, formulations, and applications. The self-interaction corrections get rid of the self-interaction error, which is the sum of the Coulomb and exchange self-interactions that remains because of the use of an approximate exchange functional. The most frequently used self-interaction correction is the Perdew-Zunger correction. However, this correction leads to instabilities in the electronic state calculations of molecules. To avoid these instabilities, several self-interaction corrections have been developed on the basis of the characteristic behaviors of self-interacting electrons, which have no two-electron interactions. These include the von Weizsäcker kinetic energy and long-range (far-from-nucleus) asymptotic correction. Applications of self-interaction corrections have shown that the self-interaction error has a serious effect on the states of core electrons, but it has a smaller than expected effect on valence electrons. This finding is supported by the fact that the distribution of self-interacting electrons indicates that they are near atomic nuclei rather than in chemical bonds.
Empirical corrections for atmospheric neutral density derived from thermospheric models
NASA Astrophysics Data System (ADS)
Forootan, Ehsan; Kusche, Jürgen; Börger, Klaus; Henze, Christina; Löcher, Anno; Eickmans, Marius; Agena, Jens
2016-04-01
Accurately predicting satellite positions is a prerequisite for various applications from space situational awareness to precise orbit determination (POD). Given the fact that atmospheric drag represents a dominant influence on the position of low-Earth orbit objects, an accurate evaluation of thermospheric mass density is of great importance to low Earth orbital prediction. Over decades, various empirical atmospheric models have been developed to support computation of density changes within the atmosphere. The quality of these models is, however, restricted mainly due to the complexity of atmospheric density changes and the limited resolution of indices used to account for atmospheric temperature and neutral density changes caused by solar and geomagnetic activity. Satellite missions, such as Challenging Mini-Satellite Payload (CHAMP) and Gravity Recovery and Climate Experiment (GRACE), provide a direct measurement of non-conservative accelerations, acting on the surface of satellites. These measurements provide valuable data for improving our knowledge of thermosphere density and winds. In this paper we present two empirical frameworks to correct model-derived neutral density simulations by the along-track thermospheric density measurements of CHAMP and GRACE. First, empirical scale factors are estimated by analyzing daily CHAMP and GRACE acceleration measurements and are used to correct the density simulation of Jacchia and MSIS (Mass-Spectrometer-Incoherent-Scatter) thermospheric models. The evolution of daily scale factors is then related to solar and magnetic activity enabling their prediction in time. In the second approach, principal component analysis (PCA) is applied to extract the dominant modes of differences between CHAMP/GRACE observations and thermospheric model simulations. Afterwards an adaptive correction procedure is used to account for long-term and high-frequency differences. We conclude the study by providing recommendations on possible
Level densities and shell corrections of superheavy nuclei
NASA Astrophysics Data System (ADS)
Bezbakh, A. N.; Shneidman, T. M.; Adamian, G. G.; Antonenko, N. V.
2015-02-01
The intrinsic level densities of superheavy nuclei in the α-decay chains of 296;298;300120 nuclei are calculated using the single-particle spectra obtained with the modifed two-center shell model. The level density parameters are extracted and compared with their phenomenological values used in the calculations of the survival of excited heavy nuclei. The dependences of the level density parameters on the mass and charge numbers as well as on the ground-state shell corrections are studied.
Atmospheric Density Corrections Estimated from Fitted Drag Coefficients
NASA Astrophysics Data System (ADS)
McLaughlin, C. A.; Lechtenberg, T. F.; Mance, S. R.; Mehta, P.
2010-12-01
Fitted drag coefficients estimated using GEODYN, the NASA Goddard Space Flight Center Precision Orbit Determination and Geodetic Parameter Estimation Program, are used to create density corrections. The drag coefficients were estimated for Stella, Starlette and GFZ using satellite laser ranging (SLR) measurements; and for GEOSAT Follow-On (GFO) using SLR, Doppler, and altimeter crossover measurements. The data analyzed covers years ranging from 2000 to 2004 for Stella and Starlette, 2000 to 2002 and 2005 for GFO, and 1995 to 1997 for GFZ. The drag coefficient was estimated every eight hours. The drag coefficients over the course of a year show a consistent variation about the theoretical and yearly average values that primarily represents a semi-annual/seasonal error in the atmospheric density models used. The atmospheric density models examined were NRLMSISE-00 and MSIS-86. The annual structure of the major variations was consistent among all the satellites for a given year and consistent among all the years examined. The fitted drag coefficients can be converted into density corrections every eight hours along the orbit of the satellites. In addition, drag coefficients estimated more frequently can provide a higher frequency of density correction.
Long-Range Corrected Hybrid Density Functionals with Damped Atom-Atom Dispersion Corrections
Chai, Jeng-Da; Head-Gordon, Martin
2008-06-14
We report re-optimization of a recently proposed long-range corrected (LC) hybrid density functionals [J.-D. Chai and M. Head-Gordon, J. Chem. Phys. 128, 084106 (2008)] to include empirical atom-atom dispersion corrections. The resulting functional, {omega}B97X-D yields satisfactory accuracy for thermochemistry, kinetics, and non-covalent interactions. Tests show that for non-covalent systems, {omega}B97X-D shows slight improvement over other empirical dispersion-corrected density functionals, while for covalent systems and kinetics, it performs noticeably better. Relative to our previous functionals, such as {omega}B97X, the new functional is significantly superior for non-bonded interactions, and very similar in performance for bonded interactions.
Association between lung function and airway wall density
NASA Astrophysics Data System (ADS)
Leader, J. Ken; Zheng, Bin; Fuhrman, Carl R.; Tedrow, John; Park, Sang C.; Tan, Jun; Pu, Jiantao; Drescher, John M.; Gur, David; Sciurba, Frank C.
2009-02-01
Computed tomography (CT) examination is often used to quantify the relation between lung function and airway remodeling in chronic obstructive pulmonary disease (COPD). In this preliminary study, we examined the association between lung function and airway wall computed attenuation ("density") in 200 COPD screening subjects. Percent predicted FVC (FVC%), percent predicted FEV1 (FEV1%), and the ratio of FEV1 to FVC as a percentage (FEV1/FVC%) were measured post-bronchodilator. The apical bronchus of the right upper lobe was manually selected from CT examinations for evaluation. Total airway area, lumen area, wall area, lumen perimeter and wall area as fraction of the total airway area were computed. Mean HU (meanHU) and maximum HU (maxHU) values were computed across pixels assigned membership in the wall and with a HU value greater than -550. The Pearson correlation coefficients (PCC) between FVC%, FEV1%, and FEV1/FVC% and meanHU were -0.221 (p = 0.002), -0.175 (p = 0.014), and -0.110 (p = 0.123), respectively. The PCCs for maxHU were only significant for FVC%. The correlations between lung function and the airway morphometry parameters were slightly stronger compared to airway wall density. MeanHU was significantly correlated with wall area (PCC = 0.720), airway area (0.498) and wall area percent (0.611). This preliminary work demonstrates that airway wall density is associated with lung function. Although the correlations in our study were weaker than a recent study, airway wall density initially appears to be an important parameter in quantitative CT analysis of COPD.
Mehranian, Abolfazl; Zaidi, Habib
2015-06-21
In standard segmentation-based MRI-guided attenuation correction (MRAC) of PET data on hybrid PET/MRI systems, the inter/intra-patient variability of linear attenuation coefficients (LACs) is ignored owing to the assignment of a constant LAC to each tissue class. This can lead to PET quantification errors, especially in the lung regions. In this work, we aim to derive continuous and patient-specific lung LACs from time-of-flight (TOF) PET emission data using the maximum likelihood reconstruction of activity and attenuation (MLAA) algorithm. The MLAA algorithm was constrained for estimation of lung LACs only in the standard 4-class MR attenuation map using Gaussian lung tissue preference and Markov random field smoothness priors. MRAC maps were derived from segmentation of CT images of 19 TOF-PET/CT clinical studies into background air, lung, soft tissue and fat tissue classes, followed by assignment of predefined LACs of 0, 0.0224, 0.0864 and 0.0975 cm(-1), respectively. The lung LACs of the resulting attenuation maps were then estimated from emission data using the proposed MLAA algorithm. PET quantification accuracy of MRAC and MLAA methods was evaluated against the reference CT-based AC method in the lungs, lesions located in/near the lungs and neighbouring tissues. The results show that the proposed MLAA algorithm is capable of retrieving lung density gradients and compensate fairly for respiratory-phase mismatch between PET and corresponding attenuation maps. It was found that the mean of the estimated lung LACs generally follow the trend of the reference CT-based attenuation correction (CTAC) method. Quantitative analysis revealed that the MRAC method resulted in average relative errors of -5.2 ± 7.1% and -6.1 ± 6.7% in the lungs and lesions, respectively. These were reduced by the MLAA algorithm to -0.8 ± 6.3% and -3.3 ± 4.7%, respectively. In conclusion, we demonstrated the potential and capability of emission-based methods in deriving patient
Semilocal density functional theory with correct surface asymptotics
NASA Astrophysics Data System (ADS)
Constantin, Lucian A.; Fabiano, Eduardo; Pitarke, J. M.; Della Sala, Fabio
2016-03-01
Semilocal density functional theory is the most used computational method for electronic structure calculations in theoretical solid-state physics and quantum chemistry of large systems, providing good accuracy with a very attractive computational cost. Nevertheless, because of the nonlocality of the exchange-correlation hole outside a metal surface, it was always considered inappropriate to describe the correct surface asymptotics. Here, we derive, within the semilocal density functional theory formalism, an exact condition for the imagelike surface asymptotics of both the exchange-correlation energy per particle and potential. We show that this condition can be easily incorporated into a practical computational tool, at the simple meta-generalized-gradient approximation level of theory. Using this tool, we also show that the Airy-gas model exhibits asymptotic properties that are closely related to those at metal surfaces. This result highlights the relevance of the linear effective potential model to the metal surface asymptotics.
NASA Astrophysics Data System (ADS)
Maintas, Dimitris; Houzard, Claire; Ksyar, Rachid; Mognetti, Thomas; Maintas, Catherine; Scheiber, Christian; Itti, Roland
2006-12-01
It is considered that one of the great strengths of PET imaging is the ability to correct for body attenuation. This enables better lesion uptake quantification and quality of PET images. The aim of this work is to compare the sensitivity of non-attenuation-corrected (NAC) PET images, the gamma photons (GPAC) and CT attenuation-corrected (CTAC) images in detecting and staging of lung cancer. We have studied 66 patients undergoing PET/CT examinations for detecting and staging NSC lung cancer. The patients were injected with 18-FDG; 5 MBq/kg under fasting conditions and examination was started 60 min later. Transmission data were acquired by a spiral CT X-ray tube and by gamma photons emitting Cs-137l source and were used for the patient body attenuation correction without correction for respiratory motion. In 55 of 66 patients we performed both attenuation correction procedures and in 11 patients only CT attenuation correction. In seven patients with solitary nodules PET was negative and in 59 patients with lung cancer PET/CT was positive for pulmonary or other localization. In the group of 55 patients we found 165 areas of focal increased 18-FDG uptake in NAC, 165 in CTAC and 164 in GPAC PET images.In the patients with only CTAC we found 58 areas of increased 18-FDG uptake on NAC and 58 areas lesions on CTAC. In the patients with positive PET we found 223 areas of focal increased uptake in NAC and 223 areas in CTAC images. The sensitivity of NAC was equal to the sensitivity of CTAC and GPAC images. The visualization of peripheral lesions was better in NAC images and the lesions were better localized in attenuation-corrected images. In three lesions of the thorax the localization was better in GPAC and fused images than in CTAC images.
NASA Astrophysics Data System (ADS)
Lederman, Dror; Leader, Joseph K.; Zheng, Bin; Sciurba, Frank C.; Tan, Jun; Gur, David
2011-03-01
Quantitative computed tomography (CT) has been widely used to detect and evaluate the presence (or absence) of emphysema applying the density masks at specific thresholds, e.g., -910 or -950 Hounsfield Unit (HU). However, it has also been observed that subjects with similar density-mask based emphysema scores could have varying lung function, possibly indicating differences of disease severity. To assess this possible discrepancy, we investigated whether density distribution of "viable" lung parenchyma regions with pixel values > -910 HU correlates with lung function. A dataset of 38 subjects, who underwent both pulmonary function testing and CT examinations in a COPD SCCOR study, was assembled. After the lung regions depicted on CT images were automatically segmented by a computerized scheme, we systematically divided the lung parenchyma into different density groups (bins) and computed a number of statistical features (i.e., mean, standard deviation (STD), skewness of the pixel value distributions) in these density bins. We then analyzed the correlations between each feature and lung function. The correlation between diffusion lung capacity (DLCO) and STD of pixel values in the bin of -910HU <= PV < -750HU was -0.43, as compared with a correlation of -0.49 obtained between the post-bronchodilator ratio (FEV1/FVC) measured by the forced expiratory volume in 1 second (FEV1) dividing the forced vital capacity (FVC) and the STD of pixel values in the bin of -1024HU <= PV < -910HU. The results showed an association between the distribution of pixel values in "viable" lung parenchyma and lung function, which indicates that similar to the conventional density mask method, the pixel value distribution features in "viable" lung parenchyma areas may also provide clinically useful information to improve assessments of lung disease severity as measured by lung functional tests.
Metallophilic interactions from dispersion-corrected density-functional theory
Otero-de-la-Roza, Alberto Mallory, Joel D.; Johnson, Erin R.
2014-05-14
In this article, we present the first comprehensive study of metallophilic (aurophilic) interactions using dispersion-corrected density-functional theory. Dispersion interactions (an essential component of metallophilicity) are treated using the exchange-hole dipole moment (XDM) model. By comparing against coupled-cluster benchmark calculations on simple dimers, we show that LC-ωPBE-XDM is a viable functional to study interactions between closed-shell transition metals and that it performs uniformly better than second-order Møller-Plesset theory, the basic computational technique used in previous works. We apply LC-ωPBE-XDM to address several open questions regarding metallophilicity, such as the interplay between dispersion and relativistic effects, the interaction strength along group 11, the additivity of homo- and hetero-metallophilic effects, the stability of [E(AuPH{sub 3}){sub 4}]{sup +} cations (E = N, P, As, Sb), and the role of metallophilic effects in crystal packing. We find that relativistic effects explain the prevalence of aurophilicity not by stabilizing metal-metal contacts, but by preventing gold from forming ionic structures involving bridge anions (which are otherwise common for Ag and Cu) as a result of the increased electron affinity of the metal. Dispersion effects are less important than previously assumed and their stabilization contribution is relatively independent of the metal.
Metallophilic interactions from dispersion-corrected density-functional theory.
Otero-de-la-Roza, Alberto; Mallory, Joel D; Johnson, Erin R
2014-05-14
In this article, we present the first comprehensive study of metallophilic (aurophilic) interactions using dispersion-corrected density-functional theory. Dispersion interactions (an essential component of metallophilicity) are treated using the exchange-hole dipole moment (XDM) model. By comparing against coupled-cluster benchmark calculations on simple dimers, we show that LC-ωPBE-XDM is a viable functional to study interactions between closed-shell transition metals and that it performs uniformly better than second-order Møller-Plesset theory, the basic computational technique used in previous works. We apply LC-ωPBE-XDM to address several open questions regarding metallophilicity, such as the interplay between dispersion and relativistic effects, the interaction strength along group 11, the additivity of homo- and hetero-metallophilic effects, the stability of [E(AuPH3)4](+) cations (E = N, P, As, Sb), and the role of metallophilic effects in crystal packing. We find that relativistic effects explain the prevalence of aurophilicity not by stabilizing metal-metal contacts, but by preventing gold from forming ionic structures involving bridge anions (which are otherwise common for Ag and Cu) as a result of the increased electron affinity of the metal. Dispersion effects are less important than previously assumed and their stabilization contribution is relatively independent of the metal. PMID:24832312
Relative density measurements in a simple lung phantom by Compton backscatter.
Wolf, E A; Munro, T R
1985-02-01
Compton backscatter of 60 keV gamma radiation from a simple lung phantom has been used to measure changes in "lung" density. It was shown how introduction of a small volume of air can increase as well as decrease the count. Radiation scattered from the "chest wall" was prevented from entering the detector by careful choice of geometry. The remaining count increased linearly with "lung" density. The relative increase of count rate with density was entirely independent of "chest wall" thickness. With our apparatus a change of 0.01 kg/L in "lung" density produced a change in count rate of 2.2%. PMID:3980122
Firth, Amy L; Menon, Tushar; Parker, Gregory S; Qualls, Susan J; Lewis, Benjamin M; Ke, Eugene; Dargitz, Carl T; Wright, Rebecca; Khanna, Ajai; Gage, Fred H; Verma, Inder M
2015-09-01
Lung disease is a major cause of death in the United States, with current therapeutic approaches serving only to manage symptoms. The most common chronic and life-threatening genetic disease of the lung is cystic fibrosis (CF) caused by mutations in the cystic fibrosis transmembrane regulator (CFTR). We have generated induced pluripotent stem cells (iPSCs) from CF patients carrying a homozygous deletion of F508 in the CFTR gene, which results in defective processing of CFTR to the cell membrane. This mutation was precisely corrected using CRISPR to target corrective sequences to the endogenous CFTR genomic locus, in combination with a completely excisable selection system, which significantly improved the efficiency of this correction. The corrected iPSCs were subsequently differentiated to mature airway epithelial cells where recovery of normal CFTR expression and function was demonstrated. This isogenic iPSC-based model system for CF could be adapted for the development of new therapeutic approaches. PMID:26299960
Corrections and improvements of lung imaging under Optical Coherence Tomography (OCT)
NASA Astrophysics Data System (ADS)
Golabchi, Ali
Visualization and correct assessment of alveolar volume via intact lung imaging is important to study and assess respiratory mechanics. Optical Coherence Tomography (OCT), a real time imaging technique based on near-infrared interferometry, can image several layers of distal alveoli in intact, ex-vivo lung tissue. However optical effects associated with heterogeneity of lung tissue, including the refraction caused by air-tissue interfaces along alveoli and duct walls, and changes in speed of light as it travels through the tissue, result in inaccurate measurement of alveolar volume. Experimentally such errors have been difficult to analyze because of lack of ''ground truth,'' as the lung has a unique microstructure of liquid-coated thin walls surrounding relatively large airspaces, which is difficult to model with synthetic foams. In addition, both lung and foams contain airspaces of highly irregular shape, further complicating quantitative measurement of optical artifacts and correction. To address this we have adapted the Bragg-Nye bubble raft, a crystalline two-dimensional arrangement of elements similar in geometry to alveoli (up to several hundred um in diameter with thin walls) as an inflated lung phantom in order to understand, analyze and correct these errors. By applying exact optical ray tracing on OCT images of the bubble raft, the errors are predicted and corrected. The results are validated by imaging the bubble raft with OCT from one edge and with a charged coupled device (CCD) camera in transillumination from top, providing ground truth for the OCT. We also developed a tomographic technique based on incoherent summation of multiple angle-diverse images by utilizing image registration to increase our depth of imaging and our results were validated by utilizing the inflated lung phantom. In this thesis also, an experimental apparatus for macro-scale mechanical probing of lung with in-situ micro-scale imaging of alveolar deformation was analyzed
A Method for Lung Boundary Correction Using Split Bregman Method and Geometric Active Contour Model.
Feng, Changli; Zhang, Jianxun; Liang, Rui
2015-01-01
In order to get the extracted lung region from CT images more accurately, a model that contains lung region extraction and edge boundary correction is proposed. Firstly, a new edge detection function is presented with the help of the classic structure tensor theory. Secondly, the initial lung mask is automatically extracted by an improved active contour model which combines the global intensity information, local intensity information, the new edge information, and an adaptive weight. It is worth noting that the objective function of the improved model is converted to a convex model, which makes the proposed model get the global minimum. Then, the central airway was excluded according to the spatial context messages and the position relationship between every segmented region and the rib. Thirdly, a mesh and the fractal theory are used to detect the boundary that surrounds the juxtapleural nodule. Finally, the geometric active contour model is employed to correct the detected boundary and reinclude juxtapleural nodules. We also evaluated the performance of the proposed segmentation and correction model by comparing with their popular counterparts. Efficient computing capability and robustness property prove that our model can correct the lung boundary reliably and reproducibly. PMID:26089976
Refractive errors and corrections for OCT images in an inflated lung phantom
Golabchi, Ali; Faust, J.; Golabchi, F. N.; Brooks, D. H.; Gouldstone, A.; DiMarzio, C. A.
2012-01-01
Visualization and correct assessment of alveolar volume via intact lung imaging is important to study and assess respiratory mechanics. Optical Coherence Tomography (OCT), a real-time imaging technique based on near-infrared interferometry, can image several layers of distal alveoli in intact, ex vivo lung tissue. However optical effects associated with heterogeneity of lung tissue, including the refraction caused by air-tissue interfaces along alveoli and duct walls, and changes in speed of light as it travels through the tissue, result in inaccurate measurement of alveolar volume. Experimentally such errors have been difficult to analyze because of lack of ’ground truth,’ as the lung has a unique microstructure of liquid-coated thin walls surrounding relatively large airspaces, which is difficult to model with cellular foams. In addition, both lung and foams contain airspaces of highly irregular shape, further complicating quantitative measurement of optical artifacts and correction. To address this we have adapted the Bragg-Nye bubble raft, a crystalline two-dimensional arrangement of elements similar in geometry to alveoli (up to several hundred μm in diameter with thin walls) as an inflated lung phantom in order to understand, analyze and correct these errors. By applying exact optical ray tracing on OCT images of the bubble raft, the errors are predicted and corrected. The results are validated by imaging the bubble raft with OCT from one edge and with a charged coupled device (CCD) camera in transillumination from top, providing ground truth for the OCT. PMID:22567599
Trends in corrected lung cancer mortality rates in Brazil and regions
Malta, Deborah Carvalho; de Abreu, Daisy Maria Xavier; de Moura, Lenildo; Lana, Gustavo C; Azevedo, Gulnar; França, Elisabeth
2016-01-01
ABSTRACT OBJECTIVE To describe the trend in cancer mortality rates in Brazil and regions before and after correction for underreporting of deaths and redistribution of ill-defined and nonspecific causes. METHODS The study used data of deaths from lung cancer among the population aged from 30 to 69 years, notified to the Mortality Information System between 1996 and 2011, corrected for underreporting of deaths, non-registered sex and age , and causes with ill-defined or garbage codes according to sex, age, and region. Standardized rates were calculated by age for raw and corrected data. An analysis of time trend in lung cancer mortality was carried out using the regression model with autoregressive errors. RESULTS Lung cancer in Brazil presented higher rates among men compared to women, and the South region showed the highest death risk in 1996 and 2011. Mortality showed a trend of reduction for males and increase for women. CONCLUSIONS Lung cancer in Brazil presented different distribution patterns according to sex, with higher rates among men and a reduction in the mortality trend for men and increase for women. PMID:27355467
Motion correction for synthesis and analysis of respiratory-gated lung SPECT image
NASA Astrophysics Data System (ADS)
Ue, Hidenori; Haneishi, Hideaki; Iwanaga, Hideyuki; Suga, Kazuyoshi
2005-04-01
A conventional SPECT image of lung is obtained by accumulating the detected count of gamma rays over long acquisition time that contains many respiratory cycles. The lung motion due to respiration during the acquisition makes reconstructed image blurred and may lead to a misdiagnosis. If a respiratory-gated SPECT is used, reconstructed images at various phase of respiration are obtained and the blur in a image can be avoided. However, the respiratory-gated SPECT requires long time to accumulate sufficient number of counts at each phase. If the acquisition time is not long enough, the detected count becomes inadequately small and hence the reconstructed image becomes noisy. We propose a method for correcting the motion between different phase images obtained with the respiratory-gated SPECT. In this method, an objective function consisting of both the degree of similarity between a reference and a deformed image and the smoothness of deformation is defined and optimized. The expansion ratio defined as a ratio of the change of the local volume due to the deformation is introduced to preserve the total activity during the motion correction process. By summing each phase images corrected by this method, a less noisy and less blurred SPECT image can be obtained. Furthermore, this method allows us to analyze the local movement of lung. This method was applied to the computer phantom, the real phantom and some clinical data and the motion correction and visualization of local movements between inspiration and expiration phase images were successfully achieved.
High density, optically corrected, micro-channel cooled, v-groove monolithic laser diode array
Freitas, Barry L.
1998-01-01
An optically corrected, micro-channel cooled, high density laser diode array achieves stacking pitches to 33 bars/cm by mounting laser diodes into V-shaped grooves. This design will deliver>4kW/cm2 of directional pulsed laser power. This optically corrected, micro-channel cooled, high density laser is usable in all solid state laser systems which require efficient, directional, narrow bandwidth, high optical power density pump sources.
High density, optically corrected, micro-channel cooled, v-groove monolithic laser diode array
Freitas, B.L.
1998-10-27
An optically corrected, micro-channel cooled, high density laser diode array achieves stacking pitches to 33 bars/cm by mounting laser diodes into V-shaped grooves. This design will deliver > 4kW/cm{sup 2} of directional pulsed laser power. This optically corrected, micro-channel cooled, high density laser is usable in all solid state laser systems which require efficient, directional, narrow bandwidth, high optical power density pump sources. 13 figs.
Gu, Xuejun; Jelen, Urszula; Li, Jinsheng; Jia, Xun; Jiang, Steve B
2011-06-01
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 graphics processing unit (GPU). This new GPU-based dose engine is built on our previously published ultrafast FSPB computational framework (Gu et al 2009 Phys. Med. Biol. 54 6287-97). Dosimetric evaluations against Monte Carlo dose calculations are conducted on ten IMRT treatment plans (five head-and-neck cases and five 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 s (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. PMID:21558589
NASA Astrophysics Data System (ADS)
Gu, Xuejun; Jelen, Urszula; Li, Jinsheng; Jia, Xun; Jiang, Steve B.
2011-06-01
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 graphics processing unit (GPU). This new GPU-based dose engine is built on our previously published ultrafast FSPB computational framework (Gu et al 2009 Phys. Med. Biol. 54 6287-97). Dosimetric evaluations against Monte Carlo dose calculations are conducted on ten IMRT treatment plans (five head-and-neck cases and five 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 s (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.
NASA Astrophysics Data System (ADS)
Disher, Brandon; Hajdok, George; Gaede, Stewart; Battista, Jerry J.
2012-03-01
Modern radiation therapy techniques such as intensity-modulated radiation therapy (IMRT) and stereotactic body radiation therapy (SBRT) use tightly conformed megavoltage x-ray fields to irradiate a tumour within lung tissue. For these conditions, lateral electron disequilibrium (LED) may occur, which systematically perturbs the dose distribution within tumour and nearby lung tissues. The goal of this work is to determine the combination of beam and lung density parameters that cause significant LED within and near the tumour. The Monte Carlo code DOSXYZnrc (National Research Council of Canada, Ottawa, ON) was used to simulate four 20 × 20 × 25 cm3 water-lung-water slab phantoms, which contained lung tissue only, or one of three different centrally located small tumours (sizes: 1 × 1 × 1, 3 × 3 × 3, 5 × 5 × 5 cm3). Dose calculations were performed using combinations of six beam energies (Co-60 up to 18 MV), five field sizes (1 × 1 cm2 up to 15 × 15 cm2), and 12 lung densities (0.001 g cm-3 up to 1 g cm-3) for a total of 1440 simulations. We developed the relative depth-dose factor (RDDF), which can be used to characterize the extent of LED (RDDF <1.0). For RDDF <0.7 severe LED occurred, and both lung and tumour dose were drastically reduced. For example, a 6 MV (3 × 3 cm2) field was used to irradiate a 1 cm3 tumour embedded in lung with ultra-low density of 0.001 g cm-3 (RDDF = 0.2). Dose in up-stream lung and tumour centre were reduced by as much as 80% with respect to the water density calculation. These reductions were worse for smaller tumours irradiated with high energy beams, small field sizes, and low lung density. In conclusion, SBRT trials based on dose calculations in homogeneous tissue are misleading as they do not reflect the actual dosimetric effects due to LED. Future clinical trials should only use dose calculation engines that can account for electron scatter, with special attention given to patients with low lung density (i.e. emphysema
Soh, R; Lee, J; Harianto, F
2014-06-01
Purpose: To determine and compare the correction factors obtained for TLDs in 2 × 2cm{sup 2} small field in lung heterogenous phantom using Acuros XB (AXB) and EGSnrc. Methods: This study will simulate the correction factors due to the perturbation of TLD-100 chips (Harshaw/Thermoscientific, 3 × 3 × 0.9mm{sup 3}, 2.64g/cm{sup 3}) in small field lung medium for Stereotactic Body Radiation Therapy (SBRT). A physical lung phantom was simulated by a 14cm thick composite cork phantom (0.27g/cm{sup 3}, HU:-743 ± 11) sandwiched between 4cm thick Plastic Water (CIRS,Norfolk). Composite cork has been shown to be a good lung substitute material for dosimetric studies. 6MV photon beam from Varian Clinac iX (Varian Medical Systems, Palo Alto, CA) with field size 2 × 2cm{sup 2} was simulated. Depth dose profiles were obtained from the Eclipse treatment planning system Acuros XB (AXB) and independently from DOSxyznrc, EGSnrc. Correction factors was calculated by the ratio of unperturbed to perturbed dose. Since AXB has limitations in simulating actual material compositions, EGSnrc will also simulate the AXB-based material composition for comparison to the actual lung phantom. Results: TLD-100, with its finite size and relatively high density, causes significant perturbation in 2 × 2cm{sup 2} small field in a low lung density phantom. Correction factors calculated by both EGSnrc and AXB was found to be as low as 0.9. It is expected that the correction factor obtained by EGSnrc wlll be more accurate as it is able to simulate the actual phantom material compositions. AXB have a limited material library, therefore it only approximates the composition of TLD, Composite cork and Plastic water, contributing to uncertainties in TLD correction factors. Conclusion: It is expected that the correction factors obtained by EGSnrc will be more accurate. Studies will be done to investigate the correction factors for higher energies where perturbation may be more pronounced.
Thermal Corrections to Density Functional Simulations of Warm Dense Matter
NASA Astrophysics Data System (ADS)
Smith, Justin; Pribram-Jones, Aurora; Burke, Kieron
Present density functional calculations of warm dense matter often use the Mermin-Kohn-Sham (MKS) scheme at finite temperature, but employ ground-state approximations to the exchange-correlation (XC) free energy. In the simplest solvable non-trivial model, an asymmetric Hubbard dimer, we calculate the exact many-body energies, the exact Mermin-Kohn-Sham functionals for this system, and extract the exact XC free energy. For moderate temperatures and weak correlation, we show this approximation is excellent, but fails for stronger correlations. Additionally, we use this system to test various conditions that must be satisfied.
Effect of point sampling density in quantifying mouse lung emphysema
Limjunyawong, Nathachit; Kearson, Alexandra; Das, Sandhya; Mitzner, Wayne
2016-01-01
In the official joint policy document of the American Thoracic Society and European Respiratory Society (Hsia et al., 2010), the need for proper stereologic assessment of lungs was emphasized. In this document it was emphasized that for the quantitative analysis of lung histologic sections, one of the most robust and reliable methods is point and intercept counting (Knudsen et al., 2010). One of the practical aspects of this method is how many points or intercepts are needed. The answer to this question has been considered from a theoretical perspective, and it depends on the relative magnitudes of the methodological and biologic variabilities. Although it is generally accepted that in a normal lung, one needs only 100–200 points to sufficiently lower the methodological variability, given the increased variability often seen in experimental emphysematous lung injury, the requisite number of points of intercepts has not been evaluated. In this study, we examined this question by focusing on some of the relevant sampling levels in mice with extensive elastase-induced emphysema. Using fixed samples of tissue blocks, we varied the number of sampling points or intercepts from about 25 to 1000 in control and emphysematous lungs. Our results show that, at the sampling levels investigated, even with the increased heterogeneity in the lung tissue damage caused by elastase, the number of sampling points needed to detect changes is similar to what is needed for control mice. PMID:25371008
Altunbas, Cem Kavanagh, Brian; Dzingle, Wayne; Stuhr, Kelly; Gaspar, Laurie; Miften, Moyed
2013-01-01
Early experience with stereotactic body radiation therapy (SBRT) of centrally located lung tumors indicated increased rate of high-grade toxicity in the lungs. These clinical results were based on treatment plans that were computed using pencil beam–like algorithms and without tissue inhomogeneity corrections. In this study, we evaluated the dosimetric errors in plans with and without inhomogeneity corrections and with planning target volumes (PTVs) that were within the zone of the proximal bronchial tree (BT). For 10 patients, the PTV, lungs, and sections of the BT either inside or within 2 cm of the PTV were delineated. Two treatment plans were generated for each patient using the following dose-calculation methods: (1) pencil beam (PB) algorithm without inhomogeneity correction (IC) (PB − IC) and (2) PB with inhomogeneity correction (PB + IC). Both plans had identical beam geometry but different beam segment shapes and monitor units (MU) to achieve similar conformal dose coverage of PTV. To obtain the baseline dose distributions, each plan was recalculated using a Monte Carlo (MC) algorithm by keeping MUs the same in the respective plans. The median maximum dose to the proximal BT and PTV dose coverage in the PB + IC plans were overestimated by 8% and 11%, respectively. However, the median maximum dose to the proximal BT and PTV dose coverage in PB − IC plans were underestimated by 15% and 9%. Similar trends were observed in low-dose regions of the lung within the irradiated volume. Our study indicates that dosimetric bias introduced by unit tissue density plans cannot be characterized as underestimation or overestimation of dose without taking the tumor location into account. This issue should be considered when analyzing clinical toxicity data from early lung SBRT trials that utilized unit tissue density for dose calculations.
Communication: Self-interaction correction with unitary invariance in density functional theory
Pederson, Mark R.; Ruzsinszky, Adrienn; Perdew, John P.; Department of Chemistry, Temple University, Philadelphia, Pennsylvania 19122
2014-03-28
Standard spin-density functionals for the exchange-correlation energy of a many-electron ground state make serious self-interaction errors which can be corrected by the Perdew-Zunger self-interaction correction (SIC). We propose a size-extensive construction of SIC orbitals which, unlike earlier constructions, makes SIC computationally efficient, and a true spin-density functional. The SIC orbitals are constructed from a unitary transformation that is explicitly dependent on the non-interacting one-particle density matrix. When this SIC is applied to the local spin-density approximation, improvements are found for the atomization energies of molecules.
Why Density-Gradient Corrections Improve Atomization Energies and Barrier Heights
NASA Astrophysics Data System (ADS)
Perdew, John P.; Ernzerhof, Matthias; Zupan, Ales; Burke, Kieron
While the Hartree-Fock (HF) approximation typically underestimates the strength of the chemical bond, the local spin density (LSD) approximation overestimates it. Thus LSD overbinds atoms in molecules, and underestimates the heights of energy barriers when the transition state is more highly bonded than the initial state. Generalized gradient approximations (GGA's), which incorporate density-gradient corrections to LSD, improve the agreement between calculated and measured energetics. This has been previously understood as a consequence of the fact that gradient corrections favor density inhomogeneity, which increases when a bond is stretched or broken. We show that gradient corrections also favor high density, which increases when a bond is compressed or formed, but that the inhomogeneity effect usually prevails. To quantify the discussion, we present a thermodynamic-like inequality which is satisfied when gradient corrections favor a process.
NASA Astrophysics Data System (ADS)
Zheng, Bin; Leader, J. K.; Coxson, Harvey O.; Scuirba, Frank C.; Fuhrman, Carl R.; Balkan, Arzu; Weissfeld, Joel L.; Maitz, Glenn S.; Gur, David
2006-03-01
The fraction of lung voxels below a pixel value "cut-off" has been correlated with pathologic estimates of emphysema. We performed a "standard" quantitative CT (QCT) lung analysis using a -950 HU cut-off to determine the volume fraction of emphysema (below the cut-off) and a "corrected" QCT analysis after removing small group (5 and 10 pixels) of connected pixels ("blobs") below the cut-off. CT examinations two dataset of 15 subjects each with a range of visible emphysema and pulmonary obstruction were acquired at "low-dose and conventional dose reconstructed using a high-spatial frequency kernel at 2.5 mm section thickness for the same subject. The "blob" size (i.e., connected-pixels) removed was inversely related to the computed fraction of emphysema. The slopes of emphysema fraction versus blob size were 0.013, 0.009, and 0.005 for subjects with both no emphysema and no pulmonary obstruction, moderate emphysema and pulmonary obstruction, and severe emphysema and severe pulmonary obstruction, respectively. The slopes of emphysema fraction versus blob size were 0.008 and 0.006 for low-dose and conventional CT examinations, respectively. The small blobs of pixels removed are most likely CT image artifacts and do not represent actual emphysema. The magnitude of the blob correction was appropriately associated with COPD severity. The blob correction appears to be applicable to QCT analysis in low-dose and conventional CT exams.
Diot, Quentin; Kavanagh, Brian; Schefter, Tracey; Gaspar, Laurie; Stuhr, Kelly; Miften, Moyed
2012-11-15
Purpose: To describe regional lung tissue density changes in normal lung tissue of patients with primary and metastatic lung tumors who received stereotactic body radiation therapy (SBRT). Methods and Materials: A total of 179 post-SBRT follow-up computed tomography (CT) scans of 62 patients who received SBRT between 2003 and 2009 were studied. Median prescription dose was 54 Gy (range, 30-60 Gy) in 3 to 5 fractions. SBRT-induced lung density changes on post-SBRT follow-up CT were evaluated at approximately 3, 6, 12, 18, 24, and 30 months after treatment. Dose-response curves (DRC) were generated for SBRT-induced lung damage by averaging CT number (HU) changes for regions of the lungs receiving the same dose at 5-Gy intervals. Results: For all follow-up interval periods, CT numbers linearly increased with dose until 35 Gy and were constant thereafter. For 3, 18, 24, and 30 months, the rate of relative electron density increase with dose was approximately 0.24% per Gy. At 6 months, the rate was also similar below 20 Gy but then rose to 0.6% per Gy above this threshold. After 6 months, DRCs were mostly time-independent. When split between patients treated with 3 fractions of 12 to 20 Gy (median, 20 Gy; average tumor volume, 12 {+-} 16 cm{sup 3}) and with >3 fractions of 6 to 12.5 Gy (median, 9 Gy; average tumor volume, 30 {+-} 40 cm{sup 3}), DRCs differed significantly. In both cases, CT changes at 3, 18, 24, and 30 months were identical to those of the population DRC; however, patients who received >3 fractions showed 6-month CT changes that were more than twice those for the group that received 3 fractions. Conclusions: This analysis of SBRT-induced normal lung density changes indicates that lung normal tissue has more pronounced self-limited acute effects than late effects. Differences in acute CT changes following treatments in 3 fractions were considerably less than for treatments in >3 fractions.
Self-Energy Correction to Momentum-Density Distribution of Positron-Electron Pairs
NASA Astrophysics Data System (ADS)
Tang, Z.; Nagai, Y.; Inoue, K.; Toyama, T.; Chiba, T.; Saito, M.; Hasegawa, M.
2005-03-01
Positron two-dimensional angular correlation of annihilation radiation (2D ACAR), i.e., the 2D projection of the electron momentum densities sampled by positron, in Si is employed to verify the prediction of the density functional theory within the local-density approximation (LDA). Carefully conducted test shows that the LDA introduces small but definite discrepancies to the 2D-ACAR anisotropies. Self-energy calculation using the GW method indicates that density-fluctuation contributes anisotropic momentum-density correction and thus improves the agreement between theory and experiment. These results provide valuable annotations to the arguments concerning the accuracy and validity of the LDA and GW schemes.
Du, Kaifang; Reinhardt, Joseph M.; Christensen, Gary E.; Ding, Kai; Bayouth, John E.
2013-12-15
Purpose: Four-dimensional computed tomography (4DCT) can be used to make measurements of pulmonary function longitudinally. The sensitivity of such measurements to identify change depends on measurement uncertainty. Previously, intrasubject reproducibility of Jacobian-based measures of lung tissue expansion was studied in two repeat prior-RT 4DCT human acquisitions. Difference in respiratory effort such as breathing amplitude and frequency may affect longitudinal function assessment. In this study, the authors present normalization schemes that correct ventilation images for variations in respiratory effort and assess the reproducibility improvement after effort correction.Methods: Repeat 4DCT image data acquired within a short time interval from 24 patients prior to radiation therapy (RT) were used for this analysis. Using a tissue volume preserving deformable image registration algorithm, Jacobian ventilation maps in two scanning sessions were computed and compared on the same coordinate for reproducibility analysis. In addition to computing the ventilation maps from end expiration to end inspiration, the authors investigated the effort normalization strategies using other intermediated inspiration phases upon the principles of equivalent tidal volume (ETV) and equivalent lung volume (ELV). Scatter plots and mean square error of the repeat ventilation maps and the Jacobian ratio map were generated for four conditions: no effort correction, global normalization, ETV, and ELV. In addition, gamma pass rate was calculated from a modified gamma index evaluation between two ventilation maps, using acceptance criterions of 2 mm distance-to-agreement and 5% ventilation difference.Results: The pattern of regional pulmonary ventilation changes as lung volume changes. All effort correction strategies improved reproducibility when changes in respiratory effort were greater than 150 cc (p < 0.005 with regard to the gamma pass rate). Improvement of reproducibility was
Study of density and stability of a lung-equivalent gel
NASA Astrophysics Data System (ADS)
Claeys, Carolien; de Deene, Yves; Truyens, Bart; de Wagter, Carlos
2006-12-01
Gel dosimetry is a useful tool for the verification of radiation treatments in water-equivalent tissues. In order to extend the application of gel dosimetry to the lung, the density of the dosimeter should be reduced. Some methods have been proposed for the fabrication of low-density gels. Major challenges in the fabrication of these gel dosimeters are to achieve a density that equals the electron-density of lung tissue and to obtain an acceptable homogeneity. Both polymer and Fricke gel formulations have been used as basic chemical compositions for low-density gel dosimeters. To reduce the density, two approaches have been suggested: (1) Styrofoam beads can be added to the gel or (2) the gel can be beaten until a foam is obtained. In this study we followed the latter method and added sodium-dodecyl-sulphate (SDS) as a surfactant to increase the surface tension of the gel.
Li, J; Yu, Y
2014-06-01
Purpose: RTOG 1005 requires density correction in the dose calculation of breast cancer radiation treatment. The aim of the study was to evaluate the impact of density correction on the dose calculation. Methods: Eight cases were studied, which were planned on an XiO treatment planning system with pixel-by-pixel density correction using a superposition algorithm, following RTOG 1005 protocol requirements. Four were protocol Arm 1 (standard whole breast irradiation with sequential boost) cases and four were Arm 2 (hypofractionated whole breast irradiation with concurrent boost) cases. The plans were recalculated with the same monitor units without density correction. Dose calculations with and without density correction were compared. Results: Results of Arm 1 and Arm 2 cases showed similar trends in the comparison. The average differences between the calculations with and without density correction (difference = Without - With) among all the cases were: -0.82 Gy (range: -2.65∼−0.18 Gy) in breast PTV Eval D95, −0.75 Gy (range: −1.23∼0.26 Gy) in breast PTV Eval D90, −1.00 Gy (range: −2.46∼−0.29 Gy) in lumpectomy PTV Eval D95, −0.78 Gy (range: −1.30∼0.11 Gy) in lumpectomy PTV Eval D90, −0.43% (range: −0.95∼−0.14%) in ipsilateral lung V20, −0.81% (range: −1.62∼−0.26%) in V16, −1.95% (range: −4.13∼−0.84%) in V10, −2.64% (−5.55∼−1.04%) in V8, −4.19% (range: −6.92∼−1.81%) in V5, and −4.95% (range: −7.49∼−2.01%) in V4, respectively. The differences in other normal tissues were minimal. Conclusion: The effect of density correction was observed in breast target doses (an average increase of ∼1 Gy in D95 and D90, compared to the calculation without density correction) and exposed ipsilateral lung volumes in low dose region (average increases of ∼4% and ∼5% in V5 and V4, respectively)
Lung Density Changes After Stereotactic Radiotherapy: A Quantitative Analysis in 50 Patients
Palma, David A.; Soernsen de Koste, John van; Verbakel, Wilko F.A.R.; Vincent, Andrew; Senan, Suresh
2011-11-15
Purpose: Radiologic lung density changes are observed in more than 50% of patients after stereotactic body radiotherapy (SBRT) for lung cancer. We studied the relationship between SBRT dose and posttreatment computed tomography (CT) density changes, a surrogate for lung injury. Methods and Materials: The SBRT fractionation schemes used to treat Stage I lung cancer with RapidArc were three fractions of 18 Gy, five fractions of 11 Gy, or eight fractions of 7.5 Gy, prescribed at the 80% isodose. Follow-up CT scans performed at less than 6 months (n = 50) and between 6 and 9 months (n = 30) after SBRT were reviewed. Posttreatment scans were coregistered with baseline scans using a B-spline deformable registration algorithm. Voxel-Hounsfield unit histograms were created for doses between 0.5 and 50 Gy. Linear mixed effects models were used to assess the effects of SBRT dose on CT density, and the influence of possible confounders was tested. Results: Increased CT density was associated with higher dose, increasing planning target volume size, and increasing time after SBRT (all p < 0.0001). Density increases were apparent in areas receiving >6 Gy, were most prominent in areas receiving >20 Gy, and seemed to plateau above 40 Gy. In regions receiving >36 Gy, the reduction in air-filled fraction of lung after treatment was up to 18%. No increase in CT density was observed in the contralateral lung receiving {>=}3 Gy. Conclusions: A dose-response relationship exists for quantitative CT density changes after SBRT. A threshold of effect is seen at low doses, and a plateau at highest doses.
Ruggiero, Michael T; Gooch, Jonathan; Zubieta, Jon; Korter, Timothy M
2016-02-18
The problem of nonlocal interactions in density functional theory calculations has in part been mitigated by the introduction of range-corrected functional methods. While promising solutions, the continued evaluation of range corrections in the structural simulations of complex molecular crystals is required to judge their efficacy in challenging chemical environments. Here, three pyridinium-based crystals, exhibiting a wide range of intramolecular and intermolecular interactions, are used as benchmark systems for gauging the accuracy of several range-corrected density functional techniques. The computational results are compared to low-temperature experimental single-crystal X-ray diffraction and terahertz spectroscopic measurements, enabling the direct assessment of range correction in the accurate simulation of the potential energy surface minima and curvatures. Ultimately, the simultaneous treatment of both short- and long-range effects by the ωB97-X functional was found to be central to its rank as the top performer in reproducing the complex array of forces that occur in the studied pyridinium solids. These results demonstrate that while long-range corrections are the most commonly implemented range-dependent improvements to density functionals, short-range corrections are vital for the accurate reproduction of forces that rapidly diminish with distance, such as quadrupole-quadrupole interactions. PMID:26814572
Guidez, Emilie B; Gordon, Mark S
2015-03-12
The modeling of dispersion interactions in density functional theory (DFT) is commonly performed using an energy correction that involves empirically fitted parameters for all atom pairs of the system investigated. In this study, the first-principles-derived dispersion energy from the effective fragment potential (EFP) method is implemented for the density functional theory (DFT-D(EFP)) and Hartree-Fock (HF-D(EFP)) energies. Overall, DFT-D(EFP) performs similarly to the semiempirical DFT-D corrections for the test cases investigated in this work. HF-D(EFP) tends to underestimate binding energies and overestimate intermolecular equilibrium distances, relative to coupled cluster theory, most likely due to incomplete accounting for electron correlation. Overall, this first-principles dispersion correction yields results that are in good agreement with coupled-cluster calculations at a low computational cost. PMID:25651435
Effect of lung and target density on small-field dose coverage and PTV definition
Higgins, Patrick D. Ehler, Eric D.; Cho, Lawrence C.; Dusenbery, Kathryn E.
2015-04-01
We have studied the effect of target and lung density on block margin for small stereotactic body radiotherapy (SBRT) targets. A phantom (50 × 50 × 50 cm{sup 3}) was created in the Pinnacle (V9.2) planning system with a 23-cm diameter lung region of interest insert. Diameter targets of 1.6, 2.0, 3.0, and 4.0 cm were placed in the lung region of interest and centered at a physical depth of 15 cm. Target densities evaluated were 0.1 to 1.0 g/cm{sup 3}, whereas the surrounding lung density was varied between 0.05 and 0.6 g/cm{sup 3}. A dose of 100 cGy was delivered to the isocenter via a single 6-MV field, and the ratio of the average dose to points defining the lateral edges of the target to the isocenter dose was recorded for each combination. Field margins were varied from none to 1.5 cm in 0.25-cm steps. Data obtained in the phantom study were used to predict planning treatment volume (PTV) margins that would match the clinical PTV and isodose prescription for a clinical set of 39 SBRT cases. The average internal target volume (ITV) density was 0.73 ± 0.17, average local lung density was 0.33 ± 0.16, and average ITV diameter was 2.16 ± 0.8 cm. The phantom results initially underpredicted PTV margins by 0.35 cm. With this offset included in the model, the ratio of predicted-to-clinical PTVs was 1.05 ± 0.32. For a given target and lung density, it was found that treatment margin was insensitive to target diameter, except for the smallest (1.6-cm diameter) target, for which the treatment margin was more sensitive to density changes than the larger targets. We have developed a graphical relationship for block margin as a function of target and lung density, which should save time in the planning phase by shortening the design of PTV margins that can satisfy Radiation Therapy Oncology Group mandated treatment volume ratios.
Goel, Himanshu; Butler, Charles L; Windom, Zachary W; Rai, Neeraj
2016-07-12
Recent developments in dispersion corrected and nonlocal density functionals are aimed at accurately capturing dispersion interactions, a key shortcoming of local and semilocal approximations of density functional theory. These functionals have shown significant promise for dimers and small clusters of molecules as well as crystalline materials. However, their efficacy for predicting vapor liquid equilibria is largely unexplored. In this work, we examine the accuracy of dispersion-corrected and nonlocal van der Waals functionals by computing the vapor liquid coexistence curves (VLCCs) of hydrofluoromethanes. Our results indicate that the PBE-D3 functional performs significantly better in predicting saturated liquid densities than the rVV10 functional. With the PBE-D3 functional, we also find that as the number of fluorine atoms increase in the molecule, the accuracy of saturated liquid density prediction improves as well. All the functionals significantly underpredict the saturated vapor densities, which also result in an underprediction of saturated vapor pressure of all compounds. Despite the differences in the bulk liquid densities, the local microstructures of the liquid CFH3 and CF2H2 are relatively insensitive to the density functional employed. For CF3H, however, rVV10 predicts slightly more structured liquid than the PBE-D3 functional. PMID:27295451
Hoffman, Eric A; Jiang, Rui; Baumhauer, Heather; Brooks, Michael A; Carr, J Jeffrey; Detrano, Robert; Reinhardt, Joseph; Rodriguez, Josanna; Stukovsky, Karen; Wong, Nathan; Barr, R Graham
2009-01-01
Purpose Cardiac CT scans for the assessment of coronary calcium scores include approximately 70% of the lung volume and may be useful for the quantitative assessment of emphysema. The reproducibility of lung density measures from cardiac CTs and their validity compared to lung density measures from full-lung scans is unknown. Methods and Methods The Multi-Ethnic Study of Atherosclerosis (MESA) performed paired cardiac CT scans for 6,814 participants at baseline and at follow-up. The MESA-Lung Study assessed lung density measures in the lung fields of these cardiac scans, counting voxels below -910 HU as moderate-to-severe emphysema-like lung regions. We evaluated: 1) the reproducibility of lung density measures among 120 randomly selected participants, 2) the comparability of measures acquired on electron-beam CT (EBT) and multidetector CT (MDCT) scanners among 10 participants; and 3) the validity of these measures compared to full-lung scans among 42 participants. Limits of agreement were determined using Bland-Altman approaches. Results Percent emphysema measures from paired cardiac scans were highly correlated (r=0.92-0.95) with mean difference of -0.05% (95% limits of agreement: -8.3, 8.4%). Measures from EBT and MDCT scanners were comparable (mean difference -0.9%; 95% limits of agreement: -5.1, 3.3%). Percent emphysema measures from MDCT cardiac and MDCT full-lung scans were highly correlated (r=0.93) and demonstrated reasonable agreement (mean difference 2.2%; 95% limits of agreement: -9.2, 13.8%). Conclusion While full-lung imaging is preferred for the quantification of emphysema, the lung imaging from paired cardiac CTs provided a reproducible and valid quantitative assessment of emphysema in a population-based sample. PMID:19427979
Orthometric corrections from leveling, gravity, density and elevation data: a case study in Taiwan
NASA Astrophysics Data System (ADS)
Hwang, C.; Hsiao, Y.-S.
2003-08-01
A new orthometric correction (OC) formula is presented and tested with various mean gravity reduction methods using leveling, gravity, elevation, and density data. For mean gravity computations, the Helmert method, a modified Helmert method with variable density and gravity anomaly gradient, and a modified Mader method were used. An improved method of terrain correction computation based on Gaussian quadrature is used in the modified Mader method. These methods produce different results and yield OCs that are greater than 10 cm between adjacent benchmarks (separated by sim2 km) at elevations over 3000 m. Applying OC reduces misclosures at closed leveling circuits and improves the results of leveling network adjustments. Variable density yields variation of OC at millimeter level everywhere, while gravity anomaly gradient introduces variation of OC of greater than 10 cm at higher elevations, suggesting that these quantities must be considered in OC. The modified Mader method is recommended for computing OC.
Lee, J.Y.; Shank, B.; Bonfiglio, P.; Reid, A.
1984-10-01
Sequential changes in lung density measured by CT are potentially sensitive and convenient monitors of lung abnormalities following total body irradiation (TBI). Methods have been developed to compare pre- and post-TBI CT of lung. The average local features of a cross-sectional lung slice are extracted from three peripheral regions of interest in the anterior, posterior, and lateral portions of the CT image. Also, density profiles across a specific region may be obtained. These may be compared first for verification of patient position and breathing status and then for changes between pre- and post-TBI. These may also be compared with radiation dose profiles through the lung. A preliminary study on 21 leukemia patients undergoing total body irradiation indicates the following: (a) Density gradients of patients' lungs in the antero-posterior direction show a marked heterogeneity before and after transplantation compared with normal lungs. The patients with departures from normal density gradients pre-TBI correlate with later pulmonary complications. (b) Measurements of average peripheral lung densities have demonstrated that the average lung density in the younger age group is substantially higher: pre-TBI, the average CT number (1,000 scale) is -638 +/- 39 Hounsfield unit (HU) for 0-10 years old and -739 +/- 53 HU for 21-40 years old. (c) Density profiles showed no post-TBI regional changes in lung density corresponding to the dose profile across the lung, so no differentiation of a radiation-specific effect has yet been possible. Computed tomographic density profiles in the antero-posterior direction are successfully used to verify positioning of the CT slice and the breathing level of the lung.
Ding, Huanjun; Johnson, Travis; Lin, Muqing; Le, Huy Q.; Ducote, Justin L.; Su, Min-Ying; Molloi, Sabee
2013-12-15
Purpose: Quantification of breast density based on three-dimensional breast MRI may provide useful information for the early detection of breast cancer. However, the field inhomogeneity can severely challenge the computerized image segmentation process. In this work, the effect of the bias field in breast density quantification has been investigated with a postmortem study. Methods: T1-weighted images of 20 pairs of postmortem breasts were acquired on a 1.5 T breast MRI scanner. Two computer-assisted algorithms were used to quantify the volumetric breast density. First, standard fuzzy c-means (FCM) clustering was used on raw images with the bias field present. Then, the coherent local intensity clustering (CLIC) method estimated and corrected the bias field during the iterative tissue segmentation process. Finally, FCM clustering was performed on the bias-field-corrected images produced by CLIC method. The left–right correlation for breasts in the same pair was studied for both segmentation algorithms to evaluate the precision of the tissue classification. Finally, the breast densities measured with the three methods were compared to the gold standard tissue compositions obtained from chemical analysis. The linear correlation coefficient, Pearson'sr, was used to evaluate the two image segmentation algorithms and the effect of bias field. Results: The CLIC method successfully corrected the intensity inhomogeneity induced by the bias field. In left–right comparisons, the CLIC method significantly improved the slope and the correlation coefficient of the linear fitting for the glandular volume estimation. The left–right breast density correlation was also increased from 0.93 to 0.98. When compared with the percent fibroglandular volume (%FGV) from chemical analysis, results after bias field correction from both the CLIC the FCM algorithms showed improved linear correlation. As a result, the Pearson'sr increased from 0.86 to 0.92 with the bias field correction
Dispersion- and Exchange-Corrected Density Functional Theory for Sodium Ion Hydration.
Soniat, Marielle; Rogers, David M; Rempe, Susan B
2015-07-14
A challenge in density functional theory is developing functionals that simultaneously describe intermolecular electron correlation and electron delocalization. Recent exchange-correlation functionals address those two issues by adding corrections important at long ranges: an atom-centered pairwise dispersion term to account for correlation and a modified long-range component of the electron exchange term to correct for delocalization. Here we investigate how those corrections influence the accuracy of binding free energy predictions for sodium-water clusters. We find that the dual-corrected ωB97X-D functional gives cluster binding energies closest to high-level ab initio methods (CCSD(T)). Binding energy decomposition shows that the ωB97X-D functional predicts the smallest ion-water (pairwise) interaction energy and larger multibody contributions for a four-water cluster than most other functionals - a trend consistent with CCSD(T) results. Also, ωB97X-D produces the smallest amounts of charge transfer and the least polarizable waters of the density functionals studied, which mimics the lower polarizability of CCSD. When compared with experimental binding free energies, however, the exchange-corrected CAM-B3LYP functional performs best (error <1 kcal/mol), possibly because of its parametrization to experimental formation enthalpies. For clusters containing more than four waters, "split-shell" coordination must be considered to obtain accurate free energies in comparison with experiment. PMID:26575733
Local respiratory motion correction for PET/CT imaging: Application to lung cancer
Lamare, F. Fernandez, P.; Fayad, H.; Visvikis, D.
2015-10-15
Purpose: Despite multiple methodologies already proposed to correct respiratory motion in the whole PET imaging field of view (FOV), such approaches have not found wide acceptance in clinical routine. An alternative can be the local respiratory motion correction (LRMC) of data corresponding to a given volume of interest (VOI: organ or tumor). Advantages of LRMC include the use of a simple motion model, faster execution times, and organ specific motion correction. The purpose of this study was to evaluate the performance of LMRC using various motion models for oncology (lung lesion) applications. Methods: Both simulated (NURBS based 4D cardiac-torso phantom) and clinical studies (six patients) were used in the evaluation of the proposed LRMC approach. PET data were acquired in list-mode and synchronized with respiration. The implemented approach consists first in defining a VOI on the reconstructed motion average image. Gated PET images of the VOI are subsequently reconstructed using only lines of response passing through the selected VOI and are used in combination with a center of gravity or an affine/elastic registration algorithm to derive the transformation maps corresponding to the respiration effects. Those are finally integrated in the reconstruction process to produce a motion free image over the lesion regions. Results: Although the center of gravity or affine algorithm achieved similar performance for individual lesion motion correction, the elastic model, applied either locally or to the whole FOV, led to an overall superior performance. The spatial tumor location was altered by 89% and 81% for the elastic model applied locally or to the whole FOV, respectively (compared to 44% and 39% for the center of gravity and affine models, respectively). This resulted in similar associated overall tumor volume changes of 84% and 80%, respectively (compared to 75% and 71% for the center of gravity and affine models, respectively). The application of the nonrigid
NASA Astrophysics Data System (ADS)
Zupan, Ales; Burke, Kieron; Ernzerhof, Matthias; Perdew, John P.
1997-06-01
We analyze the electron densities n(r) of atoms, molecules, solids, and surfaces. The distributions of values of the Seitz radius rs=(3/4πn)1/3 and the reduced density gradient s=|∇n|/(2(3π2)1/3n4/3) in an electron density indicate which ranges of these variables are significant for physical processes. We also define energy-weighted averages of these variables, , from which local spin density (LSD) and generalized gradient approximation (GGA) exchange-correlation energies may be estimated. The changes in these averages upon rearrangement of the nuclei (atomization of molecules or solids, stretching of bond lengths or lattice parameters, change of crystal structure, etc.) are used to explain why GGA corrects LSD in the way it does. A thermodynamic-like inequality (essentially d/>d
Generalized-exchange local-spin-density-functional theory: Self-interaction correction
NASA Astrophysics Data System (ADS)
Manoli, S.; Whitehead, M. A.
1988-07-01
The local-spin-density (LSD) generalized-exchange (GX) theory is corrected for self-interaction by splitting the single-particle Fermi hole into pure-exchange and self-interaction holes. An analysis of these components shows that the non-self-interaction-corrected GX scheme overestimates the pure exchange while underestimating the self-interaction. This self-interaction-corrected scheme is called the GX-SI scheme. Using this method of correcting for self-interaction, two other approximate self-interaction-corrected (SIC) GX schemes can be derived in which (1) the GX-LSD-SI total exchange does not include the nonlocal, self-interaction potential and (2) the GX-SIX exchange is very similar to the exchange derived by Gopinathan [Phys. Rev. A 15, 2135 (1977)]. Neither of these exchanges obeys the sum rule. The GX-SI scheme contains correction terms to the LSD GX which are smaller than the corresponding ones derived in the SIC of Perdew and Zunger [Phys. Rev. B 23, 5048 (1981)]. This shows that the LSD-GX exchange is a better approximation to the true exchange of an inhomogeneous electron gas around an atom than the LSD free-electron gas exchange.
Ensemble density variational methods with self- and ghost-interaction-corrected functionals
Pastorczak, Ewa; Pernal, Katarzyna
2014-05-14
Ensemble density functional theory (DFT) offers a way of predicting excited-states energies of atomic and molecular systems without referring to a density response function. Despite a significant theoretical work, practical applications of the proposed approximations have been scarce and they do not allow for a fair judgement of the potential usefulness of ensemble DFT with available functionals. In the paper, we investigate two forms of ensemble density functionals formulated within ensemble DFT framework: the Gross, Oliveira, and Kohn (GOK) functional proposed by Gross et al. [Phys. Rev. A 37, 2809 (1988)] alongside the orbital-dependent eDFT form of the functional introduced by Nagy [J. Phys. B 34, 2363 (2001)] (the acronym eDFT proposed in analogy to eHF – ensemble Hartree-Fock method). Local and semi-local ground-state density functionals are employed in both approaches. Approximate ensemble density functionals contain not only spurious self-interaction but also the so-called ghost-interaction which has no counterpart in the ground-state DFT. We propose how to correct the GOK functional for both kinds of interactions in approximations that go beyond the exact-exchange functional. Numerical applications lead to a conclusion that functionals free of the ghost-interaction by construction, i.e., eDFT, yield much more reliable results than approximate self- and ghost-interaction-corrected GOK functional. Additionally, local density functional corrected for self-interaction employed in the eDFT framework yields excitations energies of the accuracy comparable to that of the uncorrected semi-local eDFT functional.
Ensemble density variational methods with self- and ghost-interaction-corrected functionals.
Pastorczak, Ewa; Pernal, Katarzyna
2014-05-14
Ensemble density functional theory (DFT) offers a way of predicting excited-states energies of atomic and molecular systems without referring to a density response function. Despite a significant theoretical work, practical applications of the proposed approximations have been scarce and they do not allow for a fair judgement of the potential usefulness of ensemble DFT with available functionals. In the paper, we investigate two forms of ensemble density functionals formulated within ensemble DFT framework: the Gross, Oliveira, and Kohn (GOK) functional proposed by Gross et al. [Phys. Rev. A 37, 2809 (1988)] alongside the orbital-dependent eDFT form of the functional introduced by Nagy [J. Phys. B 34, 2363 (2001)] (the acronym eDFT proposed in analogy to eHF--ensemble Hartree-Fock method). Local and semi-local ground-state density functionals are employed in both approaches. Approximate ensemble density functionals contain not only spurious self-interaction but also the so-called ghost-interaction which has no counterpart in the ground-state DFT. We propose how to correct the GOK functional for both kinds of interactions in approximations that go beyond the exact-exchange functional. Numerical applications lead to a conclusion that functionals free of the ghost-interaction by construction, i.e., eDFT, yield much more reliable results than approximate self- and ghost-interaction-corrected GOK functional. Additionally, local density functional corrected for self-interaction employed in the eDFT framework yields excitations energies of the accuracy comparable to that of the uncorrected semi-local eDFT functional. PMID:24832322
NASA Astrophysics Data System (ADS)
de La Fuente Herman, Tania
Early stage lung cancer is found with increasing frequency by screening high risk patients. Recently, the use of Stereotactic Body Radiation Therapy (SBRT) has been found to be highly successful. The hypothesis being tested here is that the use of respiratory gating and tissue heterogeneity corrections are necessary to optimize tumor and normal tissue dose distributions for SBRT.
The effect of respiratory induced density variations on non-TOF PET quantitation in the lung
NASA Astrophysics Data System (ADS)
Holman, Beverley F.; Cuplov, Vesna; Hutton, Brian F.; Groves, Ashley M.; Thielemans, Kris
2016-04-01
Accurate PET quantitation requires a matched attenuation map. Obtaining matched CT attenuation maps in the thorax is difficult due to the respiratory cycle which causes both motion and density changes. Unlike with motion, little attention has been given to the effects of density changes in the lung on PET quantitation. This work aims to explore the extent of the errors caused by pulmonary density attenuation map mismatch on dynamic and static parameter estimates. Dynamic XCAT phantoms were utilised using clinically relevant 18F-FDG and 18F-FMISO time activity curves for all organs within the thorax to estimate the expected parameter errors. The simulations were then validated with PET data from 5 patients suffering from idiopathic pulmonary fibrosis who underwent PET/Cine-CT. The PET data were reconstructed with three gates obtained from the Cine-CT and the average Cine-CT. The lung TACs clearly displayed differences between true and measured curves with error depending on global activity distribution at the time of measurement. The density errors from using a mismatched attenuation map were found to have a considerable impact on PET quantitative accuracy. Maximum errors due to density mismatch were found to be as high as 25% in the XCAT simulation. Differences in patient derived kinetic parameter estimates and static concentration between the extreme gates were found to be as high as 31% and 14%, respectively. Overall our results show that respiratory associated density errors in the attenuation map affect quantitation throughout the lung, not just regions near boundaries. The extent of this error is dependent on the activity distribution in the thorax and hence on the tracer and time of acquisition. Consequently there may be a significant impact on estimated kinetic parameters throughout the lung.
Pederson, Mark R
2015-02-14
A recent modification of the Perdew-Zunger self-interaction-correction to the density-functional formalism has provided a framework for explicitly restoring unitary invariance to the expression for the total energy. The formalism depends upon construction of Löwdin orthonormalized Fermi-orbitals which parametrically depend on variational quasi-classical electronic positions. Derivatives of these quasi-classical electronic positions, required for efficient minimization of the self-interaction corrected energy, are derived and tested, here, on atoms. Total energies and ionization energies in closed-shell singlet atoms, where correlation is less important, using the Perdew-Wang 1992 Local Density Approximation (PW92) functional, are in good agreement with experiment and non-relativistic quantum-Monte-Carlo results albeit slightly too low. PMID:25681892
Galerani, Ana Paula; Grills, Inga; Hugo, Geoffrey; Kestin, Larry; Mohammed, Nasiruddin; Chao, K. Kenneth; Suen, Andrew; Martinez, Alvaro; Yan, Di
2010-12-01
Purpose: To evaluate the dosimetric impact of online cone-beam computed tomography (CBCT) guided correction in lung stereotactic body radiation therapy (SBRT). Methods and Materials: Twenty planning and 162 CBCT images from 20 patients undergoing lung SBRT were analyzed. The precorrection CBCT (CBCT after patient setup, no couch correction) was registered to planning CT using soft tissue; couch shift was applied, with a second CBCT for verification (postcorrection CBCT). Targets and normal structures were delineated on CBCTs: gross tumor volume (GTV), clinical target volume (CTV), cord, esophagus, lung, proximal bronchial tree, and aorta. Dose distributions on all organs manifested on each CBCT were compared with those planned on the CT. Results: Without CBCT guided target position correction, target dose reduced with respect to treatment plan. Mean and standard deviation of treatment dose discrepancy from the plan were -3.2% (4.9%), -2.1% (4.4%), -6.1% (10.7%), and -3.5% (7%) for GTV D{sub 99%}, GTV D{sub 95%}, CTV D{sub 99%}, and CTV D{sub 95%}, respectively. With CBCT correction, the results were -0.4% (2.6%), 0.1% (1.7%), -0.3% (4.2%), and 0.5% (3%). Mean and standard deviation of the difference in normal organ maximum dose were 2.2% (6.5%) before correction and 2.4% (5.9%) after correction for esophagus; 6.1% (14.1%) and 3.8% (8.1%) for cord; 3.1% (17.5%) and 6.2% (9.8%) for proximal bronchial tree; and 17.7% (19.5%) and 14.1% (17%) for aorta. Conclusion: Online CBCT guidance improves the accuracy of target dose delivery for lung SBRT. However, treatment dose to normal tissue can vary regardless of the correction. Normal tissues should be considered during target registration, according to target proximity.
A neural network based error correction method for radio occultation electron density retrieval
NASA Astrophysics Data System (ADS)
Pham, Viet-Cuong; Juang, Jyh-Ching
2015-12-01
Abel inversion techniques have been widely employed to retrieve electron density profiles (EDPs) from radio occultation (RO) measurements, which are available by observing Global Navigation Satellite System (GNSS) satellites from low-earth-orbit (LEO) satellites. It is well known that the ordinary Abel inversion might introduce errors in the retrieval of EDPs when the spherical symmetry assumption is violated. The error, however, is case-dependent; therefore it is desirable to associate an error index or correction coefficient with respect to each retrieved EDP. Several error indices have been proposed but they only deal with electron density at the F2 peak and suffer from some drawbacks. In this paper we propose an artificial neural network (ANN) based error correction method for EDPs obtained by the ordinary Abel inversion. The ANN is first trained to learn the relationship between vertical total electron content (TEC) measurements and retrieval errors at the F2 peak, 220 km and 110 km altitudes; correction coefficients are then estimated to correct the retrieved EDPs at these three altitudes. Experiments using the NeQuick2 model and real FORMOSAT-3/COSMIC RO geometry show that the proposed method outperforms existing ones. Real incoherent scatter radar (ISR) measurements at the Jicamarca Radio Observatory and the global TEC map provided by the International GNSS Service (IGS) are also used to valid the proposed method.
Assessment of dispersion-corrected density functional approaches for extended systems
NASA Astrophysics Data System (ADS)
Al-Saidi, Wissam; Voora, Vamsee; Jordan, Ken
2011-03-01
Standard density functional (DFT) methods are known to fail in describing the long range van der Waals interactions, and currently, there is a great interest in incorporating dispersion corrections in density functionals. Recently, Tkatchenko and Scheffler introduced a new scheme where dispersion corrections are included by a summation of damped interatomic C6 / R6 terms. However, contrary to the DFT-D2 approach of Grimme, the C6 coefficients depend on the electron density through a Hirshfeld atom-in-a-molecule decomposition scheme. We have implemented the vdW-TS approach in VASP and applied it to the study of a series of prototype dispersion-dominated systems including layered materials, noble-gas solids and molecular crystals. Full optimization of all degrees of freedom is possible in our implementation because dispersion corrections are computed for the forces acting on the atoms, and also the stresses on the unitcell. Our results show that the vdW-TS method yield good structure, bulk moduli, and cohesive energies of weakly bonded systems in much better agreement with experiment than those obtained with standard DFT approaches.
How important is self-consistency for the dDsC density dependent dispersion correction?
Brémond, Éric; Corminboeuf, Clémence; Golubev, Nikolay; Steinmann, Stephan N.
2014-05-14
The treatment of dispersion interactions is ubiquitous but computationally demanding for seamless ab initio approaches. A highly popular and simple remedy consists in correcting for the missing interactions a posteriori by adding an attractive energy term summed over all atom pairs to standard density functional approximations. These corrections were originally based on atom pairwise parameters and, hence, had a strong touch of empiricism. To overcome such limitations, we recently proposed a robust system-dependent dispersion correction, dDsC, that is computed from the electron density and that provides a balanced description of both weak inter- and intramolecular interactions. From the theoretical point of view and for the sake of increasing reliability, we here verify if the self-consistent implementation of dDsC impacts ground-state properties such as interaction energies, electron density, dipole moments, geometries, and harmonic frequencies. In addition, we investigate the suitability of the a posteriori scheme for molecular dynamics simulations, for which the analysis of the energy conservation constitutes a challenging tests. Our study demonstrates that the post-SCF approach in an excellent approximation.
How important is self-consistency for the dDsC density dependent dispersion correction?
Brémond, Éric; Golubev, Nikolay; Steinmann, Stephan N; Corminboeuf, Clémence
2014-05-14
The treatment of dispersion interactions is ubiquitous but computationally demanding for seamless ab initio approaches. A highly popular and simple remedy consists in correcting for the missing interactions a posteriori by adding an attractive energy term summed over all atom pairs to standard density functional approximations. These corrections were originally based on atom pairwise parameters and, hence, had a strong touch of empiricism. To overcome such limitations, we recently proposed a robust system-dependent dispersion correction, dDsC, that is computed from the electron density and that provides a balanced description of both weak inter- and intramolecular interactions. From the theoretical point of view and for the sake of increasing reliability, we here verify if the self-consistent implementation of dDsC impacts ground-state properties such as interaction energies, electron density, dipole moments, geometries, and harmonic frequencies. In addition, we investigate the suitability of the a posteriori scheme for molecular dynamics simulations, for which the analysis of the energy conservation constitutes a challenging tests. Our study demonstrates that the post-SCF approach in an excellent approximation. PMID:24832324
Correcting for dispersion interaction and beyond in density functional theory through force matching
NASA Astrophysics Data System (ADS)
Song, Yang; Akin-Ojo, Omololu; Wang, Feng
2010-11-01
The force matching method is used to improve density functional theory (DFT) by designing a supplemental potential to capture the difference in atomic forces between a DFT functional and a high-quality post Hartree-Fock method. The supplemental potential has two-body terms designed to correct for dispersion and hydrogen bond interactions. The potential also has one-body terms to improve the description of the intramolecular potential energy surface. Our procedure is tested by providing corrections to the Becke-Lee-Yang-Parr exchange-correlation functional for water and is found to perform significantly better than the standard DFT-D approach, giving QCISD quality predictions for relative cluster energies, atomic forces, and molecular structures. It is found that a simple Lennard-Jones term does a good job at correcting for van der Waals interactions and possibly also providing corrections to exchange repulsion. The one-body corrections, while contributing only slightly to improving relative cluster energies, significantly reduce the errors in binding energies and atomic forces for the systems studied.
2010-01-01
Background Calibrating mammograms to produce a standardized breast density measurement for breast cancer risk analysis requires an accurate spatial measure of the compressed breast thickness. Thickness inaccuracies due to the nominal system readout value and compression paddle orientation induce unacceptable errors in the calibration. Method A thickness correction was developed and evaluated using a fully specified two-component surrogate breast model. A previously developed calibration approach based on effective radiation attenuation coefficient measurements was used in the analysis. Water and oil were used to construct phantoms to replicate the deformable properties of the breast. Phantoms consisting of measured proportions of water and oil were used to estimate calibration errors without correction, evaluate the thickness correction, and investigate the reproducibility of the various calibration representations under compression thickness variations. Results The average thickness uncertainty due to compression paddle warp was characterized to within 0.5 mm. The relative calibration error was reduced to 7% from 48-68% with the correction. The normalized effective radiation attenuation coefficient (planar) representation was reproducible under intra-sample compression thickness variations compared with calibrated volume measures. Conclusion Incorporating this thickness correction into the rigid breast tissue equivalent calibration method should improve the calibration accuracy of mammograms for risk assessments using the reproducible planar calibration measure. PMID:21080916
Goldbogen, J A; Calambokidis, J; Oleson, E; Potvin, J; Pyenson, N D; Schorr, G; Shadwick, R E
2011-01-01
Lunge feeding by rorqual whales (Balaenopteridae) is associated with a high energetic cost that decreases diving capacity, thereby limiting access to dense prey patches at depth. Despite this cost, rorquals exhibit high rates of lipid deposition and extremely large maximum body size. To address this paradox, we integrated kinematic data from digital tags with unsteady hydrodynamic models to estimate the energy budget for lunges and foraging dives of blue whales (Balaenoptera musculus), the largest rorqual and living mammal. Our analysis suggests that, despite the large amount of mechanical work required to lunge feed, a large amount of prey and, therefore, energy is obtained during engulfment. Furthermore, we suggest that foraging efficiency for blue whales is significantly higher than for other marine mammals by nearly an order of magnitude, but only if lunges target extremely high densities of krill. The high predicted efficiency is attributed to the enhanced engulfment capacity, rapid filter rate and low mass-specific metabolic rate associated with large body size in blue whales. These results highlight the importance of high prey density, regardless of prey patch depth, for efficient bulk filter feeding in baleen whales and may explain some diel changes in foraging behavior in rorqual whales. PMID:21147977
On the validity of density overrides for VMAT lung SBRT planning
Wiant, David Vanderstraeten, Caroline; Maurer, Jacqueline; Pursley, Jan; Terrell, Jonathon; Sintay, Benjamin J.
2014-08-15
Purpose: Modeling dose to a moving target in lung is a very difficult task. Current approaches to planning lung stereotactic body radiotherapy (SBRT) generally calculate dose on either free breathing or average computed tomography (CT) scans, which do not always accurately predict dose to parts of the target volume not occupied by tumor on the planning scan. In this work, the authors look at using density overrides of the target volumes to more accurately predict dose for lung SBRT using the analytic anisotropic algorithm (AAA). Methods: Volumetric modulated arc therapy plans were created on free breathing scans (FBP), time average scans (AVGP), free breathing scans with the internal target volume overridden to tumor density (ITVP), free breathing scans with the planning target volume overridden to tumor density (PTVP), and free breathing scan using a hybrid scheme with the internal target volume set to tumor density and the planning target volume minus the internal target volume set to a density intermediate between lung and tumor (HP) for the case of a 4D motion phantom and five patient cases. Radiochromic film measurements were made for the phantom plans, with gamma analysis used to compare the planned to delivered dose. The patient plans were recalculated on each of the phases of a 4DCT to evaluate tumor coverage and conformity index (CI). A modified modulation complexity score (MCSv) and average open area per control point (AA) metrics were used to evaluate multileaf collimator (MLC) modulation for each of the plans. Results: The HP plans showed significantly higher gamma passing rates (p < 0.05) than the FBP, AVGP, and ITVP for criteria of 2 mm/2% and 1 mm/1%. No significant correlation was observed between gamma values and AA or MCSv. The tumor volume was covered by the prescription dose on all phases of the 4DCT for all patient plans. The PTVP and HP yielded lower mean CI than the other plans for all five patients, with three of the cases showing
Molecular density functional theory for water with liquid-gas coexistence and correct pressure
Jeanmairet, Guillaume Levesque, Maximilien; Sergiievskyi, Volodymyr; Borgis, Daniel
2015-04-21
The solvation of hydrophobic solutes in water is special because liquid and gas are almost at coexistence. In the common hypernetted chain approximation to integral equations, or equivalently in the homogenous reference fluid of molecular density functional theory, coexistence is not taken into account. Hydration structures and energies of nanometer-scale hydrophobic solutes are thus incorrect. In this article, we propose a bridge functional that corrects this thermodynamic inconsistency by introducing a metastable gas phase for the homogeneous solvent. We show how this can be done by a third order expansion of the functional around the bulk liquid density that imposes the right pressure and the correct second order derivatives. Although this theory is not limited to water, we apply it to study hydrophobic solvation in water at room temperature and pressure and compare the results to all-atom simulations. The solvation free energy of small molecular solutes like n-alkanes and hard sphere solutes whose radii range from angstroms to nanometers is now in quantitative agreement with reference all atom simulations. The macroscopic liquid-gas surface tension predicted by the theory is comparable to experiments. This theory gives an alternative to the empirical hard sphere bridge correction used so far by several authors.
2015-11-01
In the article by Heuslein et al, which published online ahead of print on September 3, 2015 (DOI: 10.1161/ATVBAHA.115.305775), a correction was needed. Brett R. Blackman was added as the penultimate author of the article. The article has been corrected for publication in the November 2015 issue. PMID:26490278
NASA Astrophysics Data System (ADS)
Hahn, Torsten; Liebing, Simon; Kortus, Jens; Pederson, Mark
The correction of the self-interaction error that is inherent to all standard density functional theory (DFT) calculations is an object of increasing interest. We present our results on the application of the recently developed Fermi-orbital based approach for the self-interaction correction (FO-SIC) to a set of different molecular systems. Our study covers systems ranging from simple diatomic to large organic molecules. Our focus lies on the direct estimation of the ionization potential from orbital eigenvalues and on the ordering of electronic levels in metal-organic molecules. Further, we show that the Fermi orbital positions in structurally similar molecules appear to be transferable. Support by DFG FOR1154 is greatly acknowledged.
NASA Astrophysics Data System (ADS)
Hahn, T.; Liebing, S.; Kortus, J.; Pederson, Mark R.
2015-12-01
The correction of the self-interaction error that is inherent to all standard density functional theory calculations is an object of increasing interest. In this article, we apply the very recently developed Fermi-orbital based approach for the self-interaction correction [M. R. Pederson et al., J. Chem. Phys. 140, 121103 (2014) and M. R. Pederson, J. Chem. Phys. 142, 064112 (2015)] to a set of different molecular systems. Our study covers systems ranging from simple diatomic to large organic molecules. We focus our analysis on the direct estimation of the ionization potential from orbital eigenvalues. Further, we show that the Fermi orbital positions in structurally similar molecules appear to be transferable.
Andersson, M P
2016-07-28
We have performed density functional theory calculations using our modified DFT-D2 dispersion correction for metals to investigate adsorption of a range of molecules on Pt(111). The agreement between our calculations and experimental adsorption energies ranging from 0 to 3 eV was excellent with a mean absolute deviation of 0.19 eV and a maximum deviation of 0.37 eV. Our results show that the DFT-D2 semiempirical dispersion correction can provide accurate results also for describing adsorption on metals, provided that relevant physical properties of the system are taken into account, such as shorter ranged dispersion because of screening by the conducting electrons and a lower polarizability of the core electrons in metals compared to isolated atoms. PMID:27357643
Hahn, T. Liebing, S.; Kortus, J.; Pederson, Mark R.
2015-12-14
The correction of the self-interaction error that is inherent to all standard density functional theory calculations is an object of increasing interest. In this article, we apply the very recently developed Fermi-orbital based approach for the self-interaction correction [M. R. Pederson et al., J. Chem. Phys. 140, 121103 (2014) and M. R. Pederson, J. Chem. Phys. 142, 064112 (2015)] to a set of different molecular systems. Our study covers systems ranging from simple diatomic to large organic molecules. We focus our analysis on the direct estimation of the ionization potential from orbital eigenvalues. Further, we show that the Fermi orbital positions in structurally similar molecules appear to be transferable.
Long-range corrected density functional theory with linearly-scaled HF exchange
NASA Astrophysics Data System (ADS)
Song, Jong-Won; Hirao, Kimihiko
2015-12-01
Long-range corrected density functional theory (LC-DFT) attracts many chemists' attentions as a quantum chemical method to be applied to large molecular system and its property calculations. However, the expensive time cost to evaluate the long-range HF exchange is a big obstacle to be overcome to be applied to the large molecular systems and the solid state materials. Upon this problem, we propose a linear-scaling method of the HF exchange integration, in particular, for the LC-DFT hybrid functional.
Long-range corrected density functional theory with linearly-scaled HF exchange
Song, Jong-Won; Hirao, Kimihiko
2015-12-31
Long-range corrected density functional theory (LC-DFT) attracts many chemists’ attentions as a quantum chemical method to be applied to large molecular system and its property calculations. However, the expensive time cost to evaluate the long-range HF exchange is a big obstacle to be overcome to be applied to the large molecular systems and the solid state materials. Upon this problem, we propose a linear-scaling method of the HF exchange integration, in particular, for the LC-DFT hybrid functional.
Deconvolution-based correction of alkali beam emission spectroscopy density profile measurements
Pusztai, I.; Pokol, G.; Refy, D.; Por, G.; Dunai, D.; Anda, G.; Zoletnik, S.; Schweinzer, J.
2009-08-15
A deconvolution-based correction method of the beam emission spectroscopy (BES) density profile measurement is demonstrated by its application to simulated measurements of the COMPASS and TEXTOR tokamaks. If the line of sight is far from tangential to the flux surfaces, and the beam width is comparable to the scale length on which the light profile varies, the observation may cause an undesired smoothing of the light profile, resulting in a non-negligible underestimation of the calculated density profile. This effect can be reduced significantly by the emission reconstruction method, which gives an estimate of the emissivity along the beam axis from the measured light profile, taking the finite beam width and the properties of the measurement into account in terms of the transfer function of the observation. Characteristics and magnitude of the mentioned systematic error and its reduction by the introduced method are studied by means of the comprehensive alkali BES simulation code RENATE.
NASA Astrophysics Data System (ADS)
Landerville, Aaron; Oleynik, Ivan
2015-06-01
Dispersion Corrected Density Functional Theory (DFT+vdW) calculations are performed to predict vibrational and thermal properties of the bulk energetic materials (EMs) β-octahydrocyclotetramethylene-tetranitramine (β-HMX) and triaminotrinitrobenzene (TATB). DFT+vdW calculations of optimized unit cells along the hydrostatic equation of state are followed by frozen-phonon calculations of their respective vibration spectra. These are then used under the quasi-harmonic approximation to obtain zero-point and thermal free energy contributions to the pressure, resulting in PVT equations of state for each material that is in excellent agreement with experiment. Further, heat capacities, thermal expansion coefficients, and Gruneissen parameters as functions of temperature are calculated and compared with experiment. The vibrational properties, including phonon densities of states and pressure dependencies of individual modes, are also analyzed and compared with experiment.
Masoomi, Michael A; McLean, Anne H; Bouchareb, Yassine; Ryder, Will; Robinson, Andy
2013-01-01
Objective(s): To investigate the impact of respiratory motion on localization, and quantification of lung lesions for the Gross Tumor Volume utilizing a fully automated Auto3Dreg program and dynamic NURBS-based cardiac-torso digitized phantom (NCAT). Methods: Respiratory motion may result in more than 30% underestimation of the SUV values of lung, liver and kidney tumor lesions. The motion correction technique adopted in this study was an image-based motion correction approach using, a voxel-intensity-based and a multi-resolution multi-optimization (MRMO) algorithm. The NCAT phantom was used to generate CT attenuation maps and activity distribution volumes for the lung regions. All the generated frames were co-registered to a reference frame using a time efficient scheme. Quantitative assessment including Region of Interest (ROI), image fidelity and image correlation techniques, as well as semi-quantitative line profile analysis and qualitatively overlaying non-motion and motion corrected image frames were performed. Results: The largest motion was observed in the Z-direction. The greatest translation was for the frame 3, end inspiration, and the smallest for the frame 5 which was closet frame to the reference frame at 67% expiration. Visual assessment of the lesion sizes, 20-60mm at 3 different locations, apex, mid and base of lung showed noticeable improvement for all the foci and their locations. The maximum improvements for the image fidelity were from 0.395 to 0.930 within the lesion volume of interest. The greatest improvement in activity concentration underestimation was 7.7% below the true activity for the 20 mm lesion in comparison to 34.4% below, prior to correction. The discrepancies in activity underestimation were reduced with increasing the lesion sizes. Overlaying activity distribution on the attenuation map showed improved localization of the PET metabolic information to the anatomical CT images. Conclusion: The respiratory motion correction for the
Dzyubak, Oleksandr; Kincaid, Russell; Hertanto, Agung; Hu, Yu-Chi; Pham, Hai; Yorke, Ellen; Zhang, Qinghui; Mageras, Gig S.; Rimner, Andreas
2014-10-15
Purpose: Target localization accuracy of cone-beam CT (CBCT) images used in radiation treatment of respiratory disease sites is affected by motion artifacts (blurring and streaking). The authors have previously reported on a method of respiratory motion correction in thoracic CBCT at end expiration (EE). The previous retrospective study was limited to examination of reducing motion artifacts in a small number of patient cases. They report here on a prospective study in a larger group of lung cancer patients to evaluate respiratory motion-corrected (RMC)-CBCT ability to improve lung tumor localization accuracy and reduce motion artifacts in Linac-mounted CBCT images. A second study goal examines whether the motion correction derived from a respiration-correlated CT (RCCT) at simulation yields similar tumor localization accuracy at treatment. Methods: In an IRB-approved study, 19 lung cancer patients (22 tumors) received a RCCT at simulation, and on one treatment day received a RCCT, a respiratory-gated CBCT at end expiration, and a 1-min CBCT. A respiration monitor of abdominal displacement was used during all scans. In addition to a CBCT reconstruction without motion correction, the motion correction method was applied to the same 1-min scan. Projection images were sorted into ten bins based on abdominal displacement, and each bin was reconstructed to produce ten intermediate CBCT images. Each intermediate CBCT was deformed to the end expiration state using a motion model derived from RCCT. The deformed intermediate CBCT images were then added to produce a final RMC-CBCT. In order to evaluate the second study goal, the CBCT was corrected in two ways, one using a model derived from the RCCT at simulation [RMC-CBCT(sim)], the other from the RCCT at treatment [RMC-CBCT(tx)]. Image evaluation compared uncorrected CBCT, RMC-CBCT(sim), and RMC-CBCT(tx). The gated CBCT at end expiration served as the criterion standard for comparison. Using automatic rigid image
2015-12-01
In the article by Narayan et al (Narayan O, Davies JE, Hughes AD, Dart AM, Parker KH, Reid C, Cameron JD. Central aortic reservoir-wave analysis improves prediction of cardiovascular events in elderly hypertensives. Hypertension. 2015;65:629–635. doi: 10.1161/HYPERTENSIONAHA.114.04824), which published online ahead of print December 22, 2014, and appeared in the March 2015 issue of the journal, some corrections were needed.On page 632, Figure, panel A, the label PRI has been corrected to read RPI. In panel B, the text by the upward arrow, "10% increase in kd,” has been corrected to read, "10% decrease in kd." The corrected figure is shown below.The authors apologize for these errors. PMID:26558821
Error Correction using Quantum Quasi-Cyclic Low-Density Parity-Check(LDPC) Codes
NASA Astrophysics Data System (ADS)
Jing, Lin; Brun, Todd; Quantum Research Team
Quasi-cyclic LDPC codes can approach the Shannon capacity and have efficient decoders. Manabu Hagiwara et al., 2007 presented a method to calculate parity check matrices with high girth. Two distinct, orthogonal matrices Hc and Hd are used. Using submatrices obtained from Hc and Hd by deleting rows, we can alter the code rate. The submatrix of Hc is used to correct Pauli X errors, and the submatrix of Hd to correct Pauli Z errors. We simulated this system for depolarizing noise on USC's High Performance Computing Cluster, and obtained the block error rate (BER) as a function of the error weight and code rate. From the rates of uncorrectable errors under different error weights we can extrapolate the BER to any small error probability. Our results show that this code family can perform reasonably well even at high code rates, thus considerably reducing the overhead compared to concatenated and surface codes. This makes these codes promising as storage blocks in fault-tolerant quantum computation. Error Correction using Quantum Quasi-Cyclic Low-Density Parity-Check(LDPC) Codes.
Dipietro, Laura; Poizner, Howard; Krebs, Hermano I.
2015-01-01
The ability to control online motor corrections is key to dealing with unexpected changes arising in the environment with which we interact. How the CNS controls online motor corrections is poorly understood, but evidence has accumulated in favor of a submovement-based model in which apparently continuous movement is segmented into distinct submovements. Although most studies have focused on submovements’ kinematic features, direct links with the underlying neural dynamics have not been extensively explored. This study sought to identify an electroencephalographic signature of submovements. We elicited kinematic submovements using a double-step displacement paradigm. Participants moved their wrist toward a target whose direction could shift mid-movement with a 50% probability. Movement kinematics and cortical activity were concurrently recorded with a low-friction robotic device and high-density electroencephalography. Analysis of spatiotemporal dynamics of brain activation and its correlation with movement kinematics showed that the production of each kinematic submovement was accompanied by (1) stereotyped topographic scalp maps and (2) frontoparietal ERPs time-locked to submovements. Positive ERP peaks from frontocentral areas contralateral to the moving wrist preceded kinematic submovement peaks by 220–250 msec and were followed by positive ERP peaks from contralateral parietal areas (140–250 msec latency, 0–80 msec before submovement peaks). Moreover, individual subject variability in the latency of frontoparietal ERP components following the target shift significantly predicted variability in the latency of the corrective submovement. Our results are in concordance with evidence for the intermittent nature of continuous movement and elucidate the timing and role of frontoparietal activations in the generation and control of corrective submovements. PMID:24564462
NASA Astrophysics Data System (ADS)
1995-04-01
Seismic images of the Brooks Range, Arctic Alaska, reveal crustal-scale duplexing: Correction Geology, v. 23, p. 65 68 (January 1995) The correct Figure 4A, for the loose insert, is given here. See Figure 4A below. Corrected inserts will be available to those requesting copies of the article from the senior author, Gary S. Fuis, U.S. Geological Survey, 345 Middlefield Road, Menlo Park, CA 94025. Figure 4A. P-wave velocity model of Brooks Range region (thin gray contours) with migrated wide-angle reflections (heavy red lines) and migreated vertical-incidence reflections (short black lines) superimposed. Velocity contour interval is 0.25 km/s; 4,5, and 6 km/s contours are labeled. Estimated error in velocities is one contour interval. Symbols on faults shown at top are as in Figure 2 caption.
Otero-de-la-Roza, A; Johnson, Erin R; DiLabio, Gino A
2014-12-01
Halogen bonds are formed when a Lewis base interacts with a halogen atom in a different molecule, which acts as an electron acceptor. Due to its charge transfer component, halogen bonding is difficult to model using many common density-functional approximations because they spuriously overstabilize halogen-bonded dimers. It has been suggested that dispersion-corrected density functionals are inadequate to describe halogen bonding. In this work, we show that the exchange-hole dipole moment (XDM) dispersion correction coupled with functionals that minimize delocalization error (for instance, BH&HLYP, but also other half-and-half functionals) accurately model halogen-bonded interactions, with average errors similar to other noncovalent dimers with less charge-transfer effects. The performance of XDM is evaluated for three previously proposed benchmarks (XB18 and XB51 by Kozuch and Martin, and the set proposed by Bauzá et al.) spanning a range of binding energies up to ∼50 kcal/mol. The good performance of BH&HLYP-XDM is comparable to M06-2X, and extends to the "extreme" cases in the Bauzá set. This set contains anionic electron donors where charge transfer occurs even at infinite separation, as well as other charge transfer dimers belonging to the pnictogen and chalcogen bonding classes. We also show that functional delocalization error results in an overly delocalized electron density and exact-exchange hole. We propose intermolecular Bader delocalization indices as an indicator of both the donor-acceptor character of an intermolecular interaction and the delocalization error coming from the underlying functional. PMID:26583227
Kim, Won June; Kim, Minho; Lee, Eok Kyun; Lebègue, Sébastien; Kim, Hyungjun
2016-08-18
Previous density functional dispersion corrections to density functional theory lead to an unphysical description of metallic systems, as exemplified by alkali and alkaline earth compounds. We demonstrate that it is possible to remedy this limitation by including screening effects into the form of interacting smeared-out dipoles in the many-body expansion of the interaction. Our new approach, called the coupled fluctuating smeared dipole model, describes equally well noncovalent systems, such as molecular pairs and crystals, and metallic systems. PMID:27487413
2016-02-01
Neogi T, Jansen TLTA, Dalbeth N, et al. 2015 Gout classification criteria: an American College of Rheumatology/European League Against Rheumatism collaborative initiative. Ann Rheum Dis 2015;74:1789–98. The name of the 20th author was misspelled. The correct spelling is Janitzia Vazquez-Mellado. We regret the error. PMID:26881284
NASA Technical Reports Server (NTRS)
Chadwick, C.
1984-01-01
This paper describes the development and use of an algorithm to compute approximate statistics of the magnitude of a single random trajectory correction maneuver (TCM) Delta v vector. The TCM Delta v vector is modeled as a three component Cartesian vector each of whose components is a random variable having a normal (Gaussian) distribution with zero mean and possibly unequal standard deviations. The algorithm uses these standard deviations as input to produce approximations to (1) the mean and standard deviation of the magnitude of Delta v, (2) points of the probability density function of the magnitude of Delta v, and (3) points of the cumulative and inverse cumulative distribution functions of Delta v. The approximates are based on Monte Carlo techniques developed in a previous paper by the author and extended here. The algorithm described is expected to be useful in both pre-flight planning and in-flight analysis of maneuver propellant requirements for space missions.
Self-interaction-corrected local-spin-density calculations for rare earth materials
Svane, A.; Temmerman, W.M.; Szotek, Z.; Laegsgaard, J.; Winter, H.
2000-04-20
The ab initio self-interaction-corrected (SIC) local-spin-density (LSD) approximation is discussed with emphasis on the ability to describe localized f-electron states in rare earth solids. Two methods for minimizing the SIC-LSD total energy functional are discussed, one using a unified Hamiltonian for all electron states, thus having the advantages of Bloch's theorem, the other one employing an iterative scheme in real space. Results for cerium and cerium compounds as well as other rare earths are presented. For the cerium compounds the onset of f-electron delocalization can be accurately described, including the intricate isostructural phase transitions in elemental cerium and CeP. In Pr and Sm the equilibrium lattice constant and zero temperature equation of state is greatly improved in comparison with the LSD results.
NASA Astrophysics Data System (ADS)
Sloan, J.; Bichoutskaia, E.; Liu, Z.; Kuganathan, N.; Faulques, E.; Suenaga, K.; Shannon, I. J.
2012-07-01
80 kV aberration-corrected transmission electron microscopy (AC-TEM) of discrete [W6O19]2- polyoxometalate ions mounted within double walled carbon nanotubes (DWNTs) allow high precision structural studies to be performed. W atom column separations within the octahedral W6 tungsten template can be visualized with sufficient clarity that correlation with full-scale density functional theory (DFT) can be achieved. Calculations performed on the gas phase and DWNT-mounted [W6O19]2- anions show good agreement, in the latter case, with measured separations between pairs of W2 atom columns imaged within equatorial WO6 octahedra and single W atoms within axial WO6 octahedra. Structural data from the tilted chiral encapsulating DWNT was also determined simultaneously with the anion structural measurements, allowing the nanotube conformation to be addressed in the DFT calculations.
Statistical mechanics of low-density parity check error-correcting codes over Galois fields
NASA Astrophysics Data System (ADS)
Nakamura, K.; Kabashima, Y.; Saad, D.
2001-11-01
A variation of low-density parity check (LDPC) error-correcting codes defined over Galois fields (GF(q)) is investigated using statistical physics. A code of this type is characterised by a sparse random parity check matrix composed of C non-zero elements per column. We examine the dependence of the code performance on the value of q, for finite and infinite C values, both in terms of the thermodynamical transition point and the practical decoding phase characterised by the existence of a unique (ferromagnetic) solution. We find different q-dependence in the cases of C = 2 and C >= 3; the analytical solutions are in agreement with simulation results, providing a quantitative measure to the improvement in performance obtained using non-binary alphabets.
2016-02-01
In the article by Guessous et al (Guessous I, Pruijm M, Ponte B, Ackermann D, Ehret G, Ansermot N, Vuistiner P, Staessen J, Gu Y, Paccaud F, Mohaupt M, Vogt B, Pechère-Bertschi A, Martin PY, Burnier M, Eap CB, Bochud M. Associations of ambulatory blood pressure with urinary caffeine and caffeine metabolite excretions. Hypertension. 2015;65:691–696. doi: 10.1161/HYPERTENSIONAHA.114.04512), which published online ahead of print December 8, 2014, and appeared in the March 2015 issue of the journal, a correction was needed.One of the author surnames was misspelled. Antoinette Pechère-Berstchi has been corrected to read Antoinette Pechère-Bertschi.The authors apologize for this error. PMID:26763012
Filatov, Michael; Huix-Rotllant, Miquel; Burghardt, Irene
2015-05-14
State-averaged (SA) variants of the spin-restricted ensemble-referenced Kohn-Sham (REKS) method, SA-REKS and state-interaction (SI)-SA-REKS, implement ensemble density functional theory for variationally obtaining excitation energies of molecular systems. In this work, the currently existing version of the SA-REKS method, which included only one excited state into the ensemble averaging, is extended by adding more excited states to the averaged energy functional. A general strategy for extension of the REKS-type methods to larger ensembles of ground and excited states is outlined and implemented in extended versions of the SA-REKS and SI-SA-REKS methods. The newly developed methods are tested in the calculation of several excited states of ground-state multi-reference systems, such as dissociating hydrogen molecule, and excited states of donor–acceptor molecular systems. For hydrogen molecule, the new method correctly reproduces the distance dependence of the lowest excited state energies and describes an avoided crossing between the doubly excited and singly excited states. For bithiophene–perylenediimide stacked complex, the SI-SA-REKS method correctly describes crossing between the locally excited state and the charge transfer excited state and yields vertical excitation energies in good agreement with the ab initio wavefunction methods.
Calibrated and completeness-corrected optical stellar density maps of the northern Galactic plane
NASA Astrophysics Data System (ADS)
Farnhill, H. J.; Drew, J. E.; Barentsen, G.; González-Solares, E. A.
2016-03-01
Following on from the second release of calibrated photometry from IPHAS, the INT/WFC Photometric Hα Survey of the Northern Galactic Plane, we present incompleteness-corrected stellar density maps in the r and i photometric bands. These have been computed to a range of limiting magnitudes reaching to 20th magnitude in r and 19th in i (Vega system), and with different angular resolutions - the highest resolution available being 1 arcmin2. The maps obtained cover 94 per cent of the 1800 square degree IPHAS footprint, spanning the Galactic latitude range, -5° < b < +5°, north of the celestial equator. The corrections for incompleteness, due to confusion and sensitivity loss at the faint limit, have been deduced by the method of artificial source injection. The presentation of this method is preceded by a discussion of other more approximate methods of determining completeness. Our method takes full account of position-dependent seeing and source ellipticity in the survey data base. The application of the star counts to testing reddened Galactic disc models is previewed by a comparison with predicted counts along three constant-longitude cuts at ℓ ≃ 30°, 90° and 175°: some overprediction of the most heavily reddened ℓ ≃ 30° counts is found, alongside good agreement at ℓ ≃ 90° and 175°.
Taylor, John S.; Folta, James A.; Montcalm, Claude
2005-01-18
Figure errors are corrected on optical or other precision surfaces by changing the local density of material in a zone at or near the surface. Optical surface height is correlated with the localized density of the material within the same region. A change in the height of the optical surface can then be caused by a change in the localized density of the material at or near the surface.
Self-interaction corrected density functional calculations of molecular Rydberg states
Gudmundsdóttir, Hildur; Zhang, Yao; Weber, Peter M.; Jónsson, Hannes
2013-11-21
A method is presented for calculating the wave function and energy of Rydberg excited states of molecules. A good estimate of the Rydberg state orbital is obtained using ground state density functional theory including Perdew-Zunger self-interaction correction and an optimized effective potential. The total energy of the excited molecule is obtained using the Delta Self-Consistent Field method where an electron is removed from the highest occupied orbital and placed in the Rydberg orbital. Results are presented for the first few Rydberg states of NH{sub 3}, H{sub 2}O, H{sub 2}CO, C{sub 2}H{sub 4}, and N(CH{sub 3}){sub 3}. The mean absolute error in the energy of the 33 molecular Rydberg states presented here is 0.18 eV. The orbitals are represented on a real space grid, avoiding the dependence on diffuse atomic basis sets. As in standard density functional theory calculations, the computational effort scales as NM{sup 2} where N is the number of orbitals and M is the number of grid points included in the calculation. Due to the slow scaling of the computational effort with system size and the high level of parallelism in the real space grid approach, the method presented here makes it possible to estimate Rydberg electron binding energy in large molecules.
Deng, Wang; Deng, Yue; Deng, Jia; Wang, Dao-Xin; Zhang, Ting
2015-01-01
Introduction: Recent study has shown that renin-angiotensin system plays an important role in the development of acute lung injury (ALI) with high level of angiotensin II (AngII) generated form AngI catalyzed by angiotensin-converting enzyme. AngII plays a major effect mainly through AT1 receptor. Therefore, we speculate inhibition of AT1 receptor may possibly attenuate the lung injury. Losartan, an antagonist of AT1 receptor for angiotensin II, attenuated lung injury by alleviation of the inflammation response in ALI, but the mechanism of losartan in ALI still remains unclear. Methods: Thirty male Sprague-Dawley rats were randomly divided into Control group, ALI group (LPS), and Losartan group (LPS + Losartan). Bronchoalveolar lavage fluid (BALF) and lung tissue were obtained for analysis. The expressions of lectin-like oxidized low-density lipoprotein receptor-1 (LOX-1), intercellular adhesion molecule-1 (ICAM-1) and caspase-3 were detected by reverse transcriptase polymerase chain reaction (RT-PCR) and western blotting. Results: In ALI group, TNF-α and protein level in BALF, MPO activity in lung tissue, pulmonary edema and lung injury were significantly increased. Losartan significantly reduced LPS-induced increase in TNF-α and protein level in BALF, MPO activity, pulmonary edema and lung injury in LPS-induced lung injury. The mRNA and protein expression levels of LOX-1 were significantly decreased with the administration of losartan in LPS-induced lung injury. Also, losartan blocked the protein levels of caspase-3 and ICAM-1 mediated by LOX-1 in LPS-induced lung injury. Conclusions: Losartan attenuated lung injury by alleviation of the inflammation and cell apoptosis by inhibition of LOX-1 in LPS-induced lung injury. PMID:26884836
2015-05-22
The Circulation Research article by Keith and Bolli (“String Theory” of c-kitpos Cardiac Cells: A New Paradigm Regarding the Nature of These Cells That May Reconcile Apparently Discrepant Results. Circ Res. 2015:116:1216-1230. doi: 10.1161/CIRCRESAHA.116.305557) states that van Berlo et al (2014) observed that large numbers of fibroblasts and adventitial cells, some smooth muscle and endothelial cells, and rare cardiomyocytes originated from c-kit positive progenitors. However, van Berlo et al reported that only occasional fibroblasts and adventitial cells derived from c-kit positive progenitors in their studies. Accordingly, the review has been corrected to indicate that van Berlo et al (2014) observed that large numbers of endothelial cells, with some smooth muscle cells and fibroblasts, and more rarely cardiomyocytes, originated from c-kit positive progenitors in their murine model. The authors apologize for this error, and the error has been noted and corrected in the online version of the article, which is available at http://circres.ahajournals.org/content/116/7/1216.full ( PMID:25999426
NASA Astrophysics Data System (ADS)
1998-12-01
Alleged mosasaur bite marks on Late Cretaceous ammonites are limpet (patellogastropod) home scars Geology, v. 26, p. 947 950 (October 1998) This article had the following printing errors: p. 947, Abstract, line 11, “sepia” should be “septa” p. 947, 1st paragraph under Introduction, line 2, “creep” should be “deep” p. 948, column 1, 2nd paragraph, line 7, “creep” should be “deep” p. 949, column 1, 1st paragraph, line 1, “creep” should be “deep” p. 949, column 1, 1st paragraph, line 5, “19774” should be “1977)” p. 949, column 1, 4th paragraph, line 7, “in particular” should be “In particular” CORRECTION Mammalian community response to the latest Paleocene thermal maximum: An isotaphonomic study in the northern Bighorn Basin, Wyoming Geology, v. 26, p. 1011 1014 (November 1998) An error appeared in the References Cited. The correct reference appears below: Fricke, H. C., Clyde, W. C., O'Neil, J. R., and Gingerich, P. D., 1998, Evidence for rapid climate change in North America during the latest Paleocene thermal maximum: Oxygen isotope compositions of biogenic phosphate from the Bighorn Basin (Wyoming): Earth and Planetary Science Letters, v. 160, p. 193 208.
Stram, D O; Langholz, B; Huberman, M; Thomas, D C
1999-09-01
The exposure estimates used to date for the analysis of lung cancer mortality in the Colorado Plateau Uranium Miners cohort were developed from radon progeny measurements taken in mines beginning in 1951. Since uranium miners were often exposed over long periods of time and since mines were not continuously monitored, much extrapolation and/or interpolation of measured dose-rates was needed in order to develop estimates of exposure for each of the miners in the cohort. We have recently re-examined the interpolation scheme used to create the histories in the light of the fit of a statistical model for the radon progeny measurements taken in mines within the Plateau, and we have computed revised exposure estimates for the large majority of miners in the cohort. This report describes the use of these new model-based revised exposure estimates in the analysis of lung cancer mortality, using follow-up data current through 1990. Specific issues addressed here are (1) the strength of the association between exposure and risk of lung cancer mortality; (2) effects of attained age and time since exposure upon risk of lung cancer mortality; and (3) exposure-rate effects upon risk. Results using the revised exposure estimates are compared to those obtained fitting the same models using the original Public Health Service (PHS) exposure estimates. We found evidence that the new exposure histories provide a better fit to the lung cancer mortality data than do the histories based upon the original PHS dose-rate estimates. In general, the new results show a stronger overall relationship (larger slope estimate) between lung cancer mortality and exposure per unit exposure compared to those obtained with the original estimates, while displaying similar age at exposure and time since exposure effects. In the reanalysis the impact of low dose-rate exposure is found to be relatively unchanged before and after exposure error correction, while the estimate of the effect of high dose
Kronik, Leeor; Tkatchenko, Alexandre
2014-11-18
CONSPECTUS: Molecular crystals are ubiquitous in many areas of science and engineering, including biology and medicine. Until recently, our ability to understand and predict their structure and properties using density functional theory was severely limited by the lack of approximate exchange-correlation functionals able to achieve sufficient accuracy. Here we show that there are many cases where the simple, minimally empirical pairwise correction scheme of Tkatchenko and Scheffler provides a useful prediction of the structure and properties of molecular crystals. After a brief introduction of the approach, we demonstrate its strength through some examples taken from our recent work. First, we show the accuracy of the approach using benchmark data sets of molecular complexes. Then we show its efficacy for structural determination using the hemozoin crystal, a challenging system possessing a wide range of strong and weak binding scenarios. Next, we show that it is equally useful for response properties by considering the elastic constants exhibited by the supramolecular diphenylalanine peptide solid and the infrared signature of water libration movements in brushite. Throughout, we emphasize lessons learned not only for the methodology but also for the chemistry and physics of the crystals in question. We further show that in many other scenarios where the simple pairwise correction scheme is not sufficiently accurate, one can go beyond it by employing a computationally inexpensive many-body dispersive approach that results in useful, quantitative accuracy, even in the presence of significant screening and/or multibody contributions to the dispersive energy. We explain the principles of the many-body approach and demonstrate its accuracy for benchmark data sets of small and large molecular complexes and molecular solids. PMID:24901508
The low density lipoprotein receptor-related protein (LRP) 1 and its function in lung diseases.
Wujak, L; Markart, P; Wygrecka, M
2016-07-01
The low density lipoprotein receptor-related protein (LRP) 1 is a ubiquitously expressed, versatile cell surface transmembrane receptor involved in embryonic development and adult tissue homeostasis. LRP1 binds and endocytoses a broad spectrum of over 40 ligands identified thus far, including lipoproteins, extracellular matrix proteins, proteases and protease/inhibitor complexes and growth factors. Interactions with other membrane receptors and intracellular adaptors/scaffolding proteins allow LRP1 to modulate cell migration, survival, proliferation and (trans) differentiation. Because LRP1 displays a wide-range of interactions and activities, its expression and function is temporally and spatially tightly controlled. It is not, therefore, surprising that deregulation of LRP1 production and/or activity is observed in several diseases. In this review, we will systematically examine the evidence for the role of LRP1 in human pathologies placing special emphasis on LRP1-mediated pathogenesis of the lung. PMID:26926950
Chan, Poh Yin; Tong, Chi Ming; Durrant, Marcus C
2011-09-01
An empirical method for estimation of the boiling points of organic molecules based on density functional theory (DFT) calculations with polarized continuum model (PCM) solvent corrections has been developed. The boiling points are calculated as the sum of three contributions. The first term is calculated directly from the structural formula of the molecule, and is related to its effective surface area. The second is a measure of the electronic interactions between molecules, based on the DFT-PCM solvation energy, and the third is employed only for planar aromatic molecules. The method is applicable to a very diverse range of organic molecules, with normal boiling points in the range of -50 to 500 °C, and includes ten different elements (C, H, Br, Cl, F, N, O, P, S and Si). Plots of observed versus calculated boiling points gave R²=0.980 for a training set of 317 molecules, and R²=0.979 for a test set of 74 molecules. The role of intramolecular hydrogen bonding in lowering the boiling points of certain molecules is quantitatively discussed. PMID:21798775
Prediction of d^0 magnetism in self-interaction corrected density functional theory
NASA Astrophysics Data System (ADS)
Das Pemmaraju, Chaitanya
2010-03-01
Over the past couple of years, the phenomenon of ``d^0 magnetism'' has greatly intrigued the magnetism community [1]. Unlike conventional magnetic materials, ``d^0 magnets'' lack any magnetic ions with open d or f shells but surprisingly, exhibit signatures of ferromagnetism often with a Curie temperature exceeding 300 K. Current research in the field is geared towards trying to understand the mechanism underlying this observed ferromagnetism which is difficult to explain within the conventional m-J paradigm [1]. The most widely studied class of d^0 materials are un-doped and light element doped wide gap Oxides such as HfO2, MgO, ZnO, TiO2 all of which have been put forward as possible d0 ferromagnets. General experimental trends suggest that the magnetism is a feature of highly defective samples leading to the expectation that the phenomenon must be defect related. In particular, based on density functional theory (DFT) calculations acceptor defects formed from the O-2p states in these Oxides have been proposed as being responsible for the ferromagnetism [2,3]. However. predicting magnetism originating from 2p orbitals is a delicate problem, which depends on the subtle interplay between covalency and Hund's coupling. DFT calculations based on semi-local functionals such as the local spin-density approximation (LSDA) can lead to qualitative failures on several fronts. On one hand the excessive delocalization of spin-polarized holes leads to half-metallic ground states and the expectation of room-temperature ferromagnetism. On the other hand, in some cases a magnetic ground state may not be predicted at all as the Hund's coupling might be under estimated. Furthermore, polaronic distortions which are often a feature of acceptor defects in Oxides are not predicted [4,5]. In this presentation, we argue that the self interaction error (SIE) inherent to semi-local functionals is responsible for the failures of LSDA and demonstrate through various examples that beyond
NASA Astrophysics Data System (ADS)
Szewczyk, J.; Gientka, D.
2003-04-01
Paleoclimatic ground surface temperature (GST) changes in last 100 ka years are a major factor causing vertical variation of terrestrial heat flow density (HFD). The value of this parameter important for thermal and rheological modelling may be considerably influenced by paleoclimatic factor and should be corrected for this reason. Very important criteria for studying paleoclimatic events on boreholes is the knowledge of depth distribution of thermal conductivity. However, core samples from majority of deep boreholes are hardly available and laboratory measurements of thermal conductivity are very scarce and sometimes not confident. We used a method of estimating the thermal conductivity from well logging data interpretation. The thermal conductivity was calculated using volumetric model of rock with mean geometric formula. The synthetic temperature logs (T_s) based on this data are an "active" method of investigation of vertical variation of HFD and GST determination. For a majority of deep boreholes in Polish Lowlands in uppermost part (<2000 m) we have observed dramatic disagreement between measured temperature (T) and synthetic data (T_s). We consider that the observed vertical variations of HFD in shallow part of profiles are mainly due to Holocene warming. The lower parts of profiles are still in thermal regime of the Weichselian glaciation. Presented results of GST in the Late Pleistocene for the representative data for 59 deep boreholes for the N of Poland. The GST was -5.17 +/- 5.45^oC. The observed big scatter of presented results seems to be consequence of unstable thermal conditions and bad calibration of old temperature logs. The amplitude of post glacial warming (ΔGST) is not less then +13.1^oC. The history of climate for the last 500 ka years shows that this time was spent mainly in ice age and this is "normal" state of HFD. The presented method of investigations seems to be very effective for determination of HFD for this condition.
Ong, A.; Berengut, J. C.; Flambaum, V. V.
2010-07-15
In this paper we consider the contribution of the anomalous magnetic moments of protons and neutrons to the nuclear charge density. We show that the spin-orbit contribution to the mean-square charge radius, which has been neglected in recent nuclear calculations, can be important in light halonuclei. We estimate the size of the effect in helium, lithium, and beryllium nuclei. It is found that the spin-orbit contribution represents a approx2% correction to the charge density at the center of the {sup 7}Be nucleus. We derive a simple expression for the correction to the mean-square charge radius due to the spin-orbit term and find that in light halonuclei it may be larger than the Darwin-Foldy term and comparable to finite size corrections. A comparison of experimental and theoretical mean-square radii including the spin-orbit contribution is presented.
Nie, K; Yue, N; Chen, T; Millevoi, R; Qin, S; Guo, J
2014-06-15
Purpose: In lung radiation treatment, PTV is formed with a margin around GTV (or CTV/ITV). Although GTV is most likely of water equivalent density, the PTV margin may be formed with the surrounding low-density tissues, which may lead to unreal dosimetric plan. This study is to evaluate whether the concern of dose calculation inside the PTV with only low density margin could be justified in lung treatment. Methods: Three SBRT cases were analyzed. The PTV from the original plan (Plan-O) was created with a 5–10 mm margin outside the ITV to incorporate setup errors and all mobility from 10 respiratory phases. Test plans were generated with the GTV shifted to the PTV edge to simulate the extreme situations with maximum setup uncertainties. Two representative positions as the very posterior-superior (Plan-PS) and anterior-inferior (Plan-AI) edge were considered. The virtual GTV was assigned a density of 1.0 g.cm−3 and surrounding lung, including the PTV margin, was defined as 0.25 g.cm−3. Also, additional plan with a 1mm tissue-margin instead of full lung-margin was created to evaluate whether a composite-margin (Plan-Comp) has a better approximation for dose calculation. All plans were generated on the average CT using Analytical Anisotropic Algorithm with heterogeneity correction on and all planning parameters/monitor unites remained unchanged. DVH analyses were performed for comparisons. Results: Despite the non-static dose distribution, the high-dose region synchronized with tumor positions. This might due to scatter conditions as greater doses were absorbed in the solid-tumor than in the surrounding low-density lungtissue. However, it still showed missing target coverage in general. Certain level of composite-margin might give better approximation for the dosecalculation. Conclusion: Our exploratory results suggest that with the lungmargin only, the planning dose of PTV might overestimate the coverage of the target during treatment. The significance of this
NASA Astrophysics Data System (ADS)
Rak, Michal Bartosz; Wozniak, Adam; Mayer, J. R. R.
2016-06-01
Coordinate measuring techniques rely on computer processing of coordinate values of points gathered from physical surfaces using contact or non-contact methods. Contact measurements are characterized by low density and high accuracy. On the other hand optical methods gather high density data of the whole object in a short time but with accuracy at least one order of magnitude lower than for contact measurements. Thus the drawback of contact methods is low density of data, while for non-contact methods it is low accuracy. In this paper a method for fusion of data from two measurements of fundamentally different nature: high density low accuracy (HDLA) and low density high accuracy (LDHA) is presented to overcome the limitations of both measuring methods. In the proposed method the concept of virtual markers is used to find a representation of pairs of corresponding characteristic points in both sets of data. In each pair the coordinates of the point from contact measurements is treated as a reference for the corresponding point from non-contact measurement. Transformation enabling displacement of characteristic points from optical measurement to their match from contact measurements is determined and applied to the whole point cloud. The efficiency of the proposed algorithm was evaluated by comparison with data from a coordinate measuring machine (CMM). Three surfaces were used for this evaluation: plane, turbine blade and engine cover. For the planar surface the achieved improvement was of around 200 μm. Similar results were obtained for the turbine blade but for the engine cover the improvement was smaller. For both freeform surfaces the improvement was higher for raw data than for data after creation of mesh of triangles.
NASA Astrophysics Data System (ADS)
Garza, Jorge; Nichols, Jeffrey A.; Dixon, David A.
2000-05-01
The Krieger, Li, and Iafrate approximation to the optimized effective potential including the self-interaction correction for density functional theory has been implemented in a molecular code, NWChem, that uses Gaussian functions to represent the Kohn and Sham spin-orbitals. The differences between the implementation of the self-interaction correction in codes where planewaves are used with an optimized effective potential are discussed. The importance of the localization of the spin-orbitals to maximize the exchange-correlation of the self-interaction correction is discussed. We carried out exchange-only calculations to compare the results obtained with these approximations, and those obtained with the local spin density approximation, the generalized gradient approximation and Hartree-Fock theory. Interesting results for the energy difference (GAP) between the highest occupied molecular orbital, HOMO, and the lowest unoccupied molecular orbital, LUMO, (spin-orbital energies of closed shell atoms and molecules) using the optimized effective potential and the self-interaction correction have been obtained. The effect of the diffuse character of the basis set on the HOMO and LUMO eigenvalues at the various levels is discussed. Total energies obtained with the optimized effective potential and the self-interaction correction show that the exchange energy with these approximations is overestimated and this will be an important topic for future work.
Li, Shaohong L; Marenich, Aleksandr V; Xu, Xuefei; Truhlar, Donald G
2014-01-16
Linear response (LR) Kohn-Sham (KS) time-dependent density functional theory (TDDFT), or KS-LR, has been widely used to study electronically excited states of molecules and is the method of choice for large and complex systems. The Tamm-Dancoff approximation to TDDFT (TDDFT-TDA or KS-TDA) gives results similar to KS-LR and alleviates the instability problem of TDDFT near state intersections. However, KS-LR and KS-TDA share a debilitating feature; conical intersections of the reference state and a response state occur in F - 1 instead of the correct F - 2 dimensions, where F is the number of internal degrees of freedom. Here, we propose a new method, named the configuration interaction-corrected Tamm-Dancoff approximation (CIC-TDA), that eliminates this problem. It calculates the coupling between the reference state and an intersecting response state by interpreting the KS reference-state Slater determinant and linear response as if they were wave functions. Both formal analysis and test results show that CIC-TDA gives similar results to KS-TDA far from a conical intersection, but the intersection occurs with the correct dimensionality. We anticipate that this will allow more realistic application of TDDFT to photochemistry. PMID:26270707
Lehtola, Susi; Head-Gordon, Martin; Jónsson, Hannes
2016-07-12
Implentation of seminumerical stability analysis for calculations using the Perdew-Zunger self-interaction correction is described. It is shown that real-valued solutions of the Perdew-Zunger equations for gas phase atoms are unstable with respect to imaginary orbital rotations, confirming that a proper implementation of the correction requires complex-valued orbitals. The orbital density dependence of the self-interaction corrected functional is found to lead to multiple local minima in the case of the acrylic acid, H6, and benzene molecules. In the case of benzene, symmetry breaking that results in incorrect ground state geometry is found to occur, erroneously leading to alternating bond lengths in the molecule. PMID:27232582
Sergiievskyi, Volodymyr P; Jeanmairet, Guillaume; Levesque, Maximilien; Borgis, Daniel
2014-06-01
Molecular density functional theory (MDFT) offers an efficient implicit-solvent method to estimate molecule solvation free-energies, whereas conserving a fully molecular representation of the solvent. Even within a second-order approximation for the free-energy functional, the so-called homogeneous reference fluid approximation, we show that the hydration free-energies computed for a data set of 500 organic compounds are of similar quality as those obtained from molecular dynamics free-energy perturbation simulations, with a computer cost reduced by 2-3 orders of magnitude. This requires to introduce the proper partial volume correction to transform the results from the grand canonical to the isobaric-isotherm ensemble that is pertinent to experiments. We show that this correction can be extended to 3D-RISM calculations, giving a sound theoretical justification to empirical partial molar volume corrections that have been proposed recently. PMID:26273876
Level density of a Fermi gas and integer partitions: A Gumbel-like finite-size correction
Roccia, Jerome; Leboeuf, Patricio
2010-04-15
We investigate the many-body level density of a gas of noninteracting fermions. We determine its behavior as a function of the temperature and the number of particles. As the temperature increases, and beyond the usual Sommerfeld expansion that describes the degenerate gas behavior, corrections due to a finite number of particles lead to Gumbel-like contributions. We discuss connections with the partition problem in number theory, extreme value statistics, and differences with respect to the Bose gas.
VandeVondele, Joost; Sprik, Michiel
2005-04-01
We have performed density functional theory based molecular dynamics (MD) simulations of the *OH radical in solution using self-interaction corrected (SIC) methods. We use a scheme recently proposed by M. d'Avezac, M. Calandra and F. Mauri [arXiv:cond-mat/0407750] in which a correction is only applied to the spin density within a restricted open shell formulation. In addition to two correction formulas employed within this scheme by M. d'Avezac, M. Calandra and F. Mauri, we propose and test an new empirical form which only introduces a scaled Coulomb term. This new functional leads to good agreement with reference calculations on radical cation dimers and on the hydroxyl water dimer in the gas phase. Applied in ab initio MD simulations, these three SIC methods provide a picture of the *OH solvation that differs qualitatively from the one obtained using the standard generalised gradient approximation (GGA). Hemibonded water, observed in GGA simulations and believed to be an artefact due to self-interaction error, is not present. We find that the *OH acts as a good hydrogen bond donor, but accepts less than two hydrogen bonds on average. These hydrogen bonds are part of a mobile, otherwise quasi-hydrophobic solvation cage. Our results show the potential of this computationally expedient scheme, which might extend the range of problems that can be modelled adequately with density functional theory. PMID:19787955
Calbo, Joaquín; Ortí, Enrique; Sancho-García, Juan C; Aragó, Juan
2015-03-10
In this work, we present a thorough assessment of the performance of some representative double-hybrid density functionals (revPBE0-DH-NL and B2PLYP-NL) as well as their parent hybrid and GGA counterparts, in combination with the most modern version of the nonlocal (NL) van der Waals correction to describe very large weakly interacting molecular systems dominated by noncovalent interactions. Prior to the assessment, an accurate and homogeneous set of reference interaction energies was computed for the supramolecular complexes constituting the L7 and S12L data sets by using the novel, precise, and efficient DLPNO-CCSD(T) method at the complete basis set limit (CBS). The correction of the basis set superposition error and the inclusion of the deformation energies (for the S12L set) have been crucial for obtaining precise DLPNO-CCSD(T)/CBS interaction energies. Among the density functionals evaluated, the double-hybrid revPBE0-DH-NL and B2PLYP-NL with the three-body dispersion correction provide remarkably accurate association energies very close to the chemical accuracy. Overall, the NL van der Waals approach combined with proper density functionals can be seen as an accurate and affordable computational tool for the modeling of large weakly bonded supramolecular systems. PMID:26579747
NASA Astrophysics Data System (ADS)
Woo, K. M.; Yu, S. S.; Barnard, J. J.
2013-06-01
It is well known that the imperfection of pulse power sources that drive the linear induction accelerators can lead to time-varying fluctuation in the accelerating voltages, which in turn leads to longitudinal emittance growth. We show that this source of emittance growth is correctable, even in space-charge dominated beams with significant transients induced by space-charge waves. Two correction methods are proposed, and their efficacy in reducing longitudinal emittance is demonstrated with three-dimensional particle-in-cell simulations.
Kim, Leonard H.; Zhang Miao; Howell, Roger W.; Yue, Ning J.; Khan, Atif J.
2013-01-15
Purpose: Recent recommendations by the American Association of Physicists in Medicine Task Group 186 emphasize the importance of understanding material properties and their effect on inhomogeneity-corrected dose calculation for brachytherapy. Radiographic contrast is normally injected into breast brachytherapy balloons. In this study, the authors independently estimate properties of contrast solution that were expected to be incorrectly specified in a commercial brachytherapy dose calculation algorithm. Methods: The mass density and atomic weight fractions of a clinical formulation of radiographic contrast solution were determined using manufacturers' data. The mass density was verified through measurement and compared with the density obtained by the treatment planning system's CT calibration. The atomic weight fractions were used to determine the photon interaction cross section of the contrast solution for a commercial high-dose-rate (HDR) brachytherapy source and compared with that of muscle. Results: The density of contrast solution was 10% less than that obtained from the CT calibration. The cross section of the contrast solution for the HDR source was 1.2% greater than that of muscle. Both errors could be addressed by overriding the density of the contrast solution in the treatment planning system. Conclusions: The authors estimate the error in mass density and cross section parameters used by a commercial brachytherapy dose calculation algorithm for radiographic contrast used in a clinical breast brachytherapy practice. This approach is adaptable to other clinics seeking to evaluate dose calculation errors and determine appropriate density override values if desired.
NASA Astrophysics Data System (ADS)
Stisen, Simon; Lajer Højberg, Anker; Refsgaard, Jens Christian; Troldborg, Lars
2015-04-01
Precipitation data of good quality and with sufficient temporal and spatial resolution is paramount for integrated hydrological modelling. In Denmark precipitation has traditionally been collected in a network consisting of automated rain gauge stations supplemented by a large number of manual stations providing daily accumulated measures. Prior to its use for hydrological modelling, precipitation data have been corrected for under catch using historic mean monthly correction factors that were uniform for the entire country. Problems on closing the water balance in hydrological modelling have questioned this correction approach, leading to a detailed national water balance study. The backbone of the analysis was the national water resources model (DK-model), which is a physically based, coupled and fully distributed model for the entire Denmark, constructed using the MIKE SHE/MIKE 11 code. The results suggested that a time-space variable approach for rain gauge catch correction based on gridded daily wind speed and temperature data is superior to the correction approach historically used. The new correction approach enabled a far better model performance on simulated discharge throughout the country and is now used in all hydrological modelling in Denmark. The results illustrates the importance of choosing an appropriate rain gauge catch correction method, especially in mid-high latitudes where solid precipitation is common. The study was carried out utilizing data from the period 1990 - 2003. Following this period the network of automated rain gauge stations have been expanded, but all manual stations have been shut down, resulting in a significant reduction in the total number of stations, with only around one fifth remaining in 2010. In 2014 the national water resources model was updated, which included a new model calibration. The model was setup for 1990 - 2010 and due to the transfer from manual to automated rain gauge station, the number of stations varies
Carrier Density Profiling of Ultra-Shallow Junction Layers Through Corrected C-V Plotting
Chen, James; Dimitrov, Dimitar; Dimitrova, Tatiana; Timans, Paul; Gelpey, Jeff; McCoy, Steve; Lerch, Wilfried; Paul, Silke; Bolze, Detlef
2008-11-03
The aim of this report is to present and justify a new approach for carrier density profiling in ultra-shallow junction (USJ) layer. This new approach is based on a capacitance measurement model, which takes series impedance, shunt resistance and the presence of a boron skin on the USJ layer into account. It allows us to extract the depletion layer capacitances in the USJ layer from C-V plotting more accurately and hence to obtain better carrier density profiles. Based on this new approach the carrier density profiles of different USJ layers with and without halo-style implants are obtained and discussed.
Burkart, Kristin M; Manichaikul, Ani; Wilk, Jemma B; Ahmed, Firas S; Burke, Gregory L; Enright, Paul; Hansel, Nadia N; Haynes, Demondes; Heckbert, Susan R; Hoffman, Eric A; Kaufman, Joel D; Kurai, Jun; Loehr, Laura; London, Stephanie J; Meng, Yang; O'Connor, George T; Oelsner, Elizabeth; Petrini, Marcy; Pottinger, Tess D; Powell, Charles A; Redline, Susan; Rotter, Jerome I; Smith, Lewis J; Soler Artigas, María; Tobin, Martin D; Tsai, Michael Y; Watson, Karol; White, Wendy; Young, Taylor R; Rich, Stephen S; Barr, R Graham
2014-04-01
Chronic obstructive pulmonary disease (COPD) is linked to cardiovascular disease; however, there are few studies on the associations of cardiovascular genes with COPD. We assessed the association of lung function with 2100 genes selected for cardiovascular diseases among 20 077 European-Americans and 6900 African-Americans. We performed replication of significant loci in the other racial group and an independent consortium of Europeans, tested the associations of significant loci with per cent emphysema and examined gene expression in an independent sample. We then tested the association of a related lipid biomarker with forced expiratory volume in 1 s (FEV1)/forced vital capacity (FVC) ratio and per cent emphysema. We identified one new polymorphism for FEV1/FVC (rs805301) in European-Americans (p=1.3×10(-6)) and a second (rs707974) in the combined European-American and African-American analysis (p=1.38×10(-7)). Both single-nucleotide polymorphisms (SNPs) flank the gene for apolipoprotein M (APOM), a component of high-density lipoprotein (HDL) cholesterol. Both were replicated in an independent cohort. SNPs in a second gene related to apolipoprotein M and HDL, PCSK9, were associated with FEV1/FVC ratio among African-Americans. rs707974 was associated with per cent emphysema among European-Americans and African-Americans and APOM expression was related to FEV1/FVC ratio and per cent emphysema. Higher HDL levels were associated with lower FEV1/FVC ratio and greater per cent emphysema. These findings suggest a novel role for the apolipoprotein M/HDL pathway in the pathogenesis of COPD and emphysema. PMID:23900982
Hwang, Seho; Shin, Jehyun; Kim, Jongman; Won, Byeongho
2015-03-10
Density log is widely applied for a variety of fields such as the petroleum exploration, mineral exploration, and geotechnical survey. The logging condition of density log is normally open holes but there are frequently cased boreholes. The primary calibration curve by slim hole logging manufacturer is normally the calibration curves for the variation of borehole diameter. In this study, we have performed the correction of steel casing effects using numerical and experimental methods. We have performed numerical modeling using the Monte Carlo N-Particle (MCNP) code based on Monte Carlo method, and field experimental method from open and cased hole log. In this study, we used the FDGS (Formation Density Gamma Sonde) for slim borehole with a 100 mCi 137Cs source, three inch borehole and steel casing. The casing effect between numerical and experimental method is well matched.
Larrue, A; Kaster, F; Kadir, T; Gooding, M; Elmpt, W van
2014-06-01
Purpose: Deformable Image Registration (DIR) is gaining wider clinical acceptance in radiation oncology. The aim of this work is to characterise a DIR algorithm on publically available 4DCT lung images, such that comparison can be performed against other algorithms. We propose an evaluation method of registration accuracy that takes into account the initial misregistration of the datasets. Methods: The “DIR Validation dataset” ( http://www.creatis.insa-lyon.fr/rio/dir{sub v}alidation{sub d}ata ) provides benchmark data for evaluating 3D CT registration algorithms. It consists of six 4DCT lung datasets (1x1x2mm resolution) with 100 landmarks identified on the end-exhalation and end-inhalation phases. Images were registered to end-inhalation using proprietary form of optical flow in commercial software (Mirada RTx, Mirada Medical, UK). Target registration error was measured before and after DIR, referred to as Initial Registration Error (IRE) and Final Registration Error (FRE). Results: The mean FRE over all landmarks was 1.37±1.81mm. FRE increased with IRE. Mean FRE of 0.86, 0.86, 1.53, 3.38, 4.45, 7.58mm was observed for IRE in the ranges 0–5, 5–10, 10–15, 15–20, 20–25, >25 mm. Higher FRE was observed at the inferior lung, where IRE was greater. Out-of-plane motion contributed more to IRE, and therefore to FRE. Maximum FRE of 20.6mm was observed for IRE of 32.1mm, located at the posterior of the middle lobe for dataset 2. Sub-voxel registration accuracy was achieved for up to 10mm IRE, and increased linearly at 0.3mm FRE/mm IRE thereafter. Conclusion: Publicly available clinical datasets enable algorithms to be compared objectively between publications. However, only reporting average TRE after registration can be misleading as the ability of an algorithm to correct for displacements varies with the IRE or position within the patient. Consequently, algorithms should be characterized using the entire range of initial displacements. For the algorithm
Lutsker, V; Aradi, B; Niehaus, T A
2015-11-14
Bridging the gap between first principles methods and empirical schemes, the density functional based tight-binding method (DFTB) has become a versatile tool in predictive atomistic simulations over the past years. One of the major restrictions of this method is the limitation to local or gradient corrected exchange-correlation functionals. This excludes the important class of hybrid or long-range corrected functionals, which are advantageous in thermochemistry, as well as in the computation of vibrational, photoelectron, and optical spectra. The present work provides a detailed account of the implementation of DFTB for a long-range corrected functional in generalized Kohn-Sham theory. We apply the method to a set of organic molecules and compare ionization potentials and electron affinities with the original DFTB method and higher level theory. The new scheme cures the significant overpolarization in electric fields found for local DFTB, which parallels the functional dependence in first principles density functional theory (DFT). At the same time, the computational savings with respect to full DFT calculations are not compromised as evidenced by numerical benchmark data. PMID:26567646
Lutsker, V.; Niehaus, T. A.; Aradi, B.
2015-11-14
Bridging the gap between first principles methods and empirical schemes, the density functional based tight-binding method (DFTB) has become a versatile tool in predictive atomistic simulations over the past years. One of the major restrictions of this method is the limitation to local or gradient corrected exchange-correlation functionals. This excludes the important class of hybrid or long-range corrected functionals, which are advantageous in thermochemistry, as well as in the computation of vibrational, photoelectron, and optical spectra. The present work provides a detailed account of the implementation of DFTB for a long-range corrected functional in generalized Kohn-Sham theory. We apply the method to a set of organic molecules and compare ionization potentials and electron affinities with the original DFTB method and higher level theory. The new scheme cures the significant overpolarization in electric fields found for local DFTB, which parallels the functional dependence in first principles density functional theory (DFT). At the same time, the computational savings with respect to full DFT calculations are not compromised as evidenced by numerical benchmark data.
NASA Astrophysics Data System (ADS)
Lutsker, V.; Aradi, B.; Niehaus, T. A.
2015-11-01
Bridging the gap between first principles methods and empirical schemes, the density functional based tight-binding method (DFTB) has become a versatile tool in predictive atomistic simulations over the past years. One of the major restrictions of this method is the limitation to local or gradient corrected exchange-correlation functionals. This excludes the important class of hybrid or long-range corrected functionals, which are advantageous in thermochemistry, as well as in the computation of vibrational, photoelectron, and optical spectra. The present work provides a detailed account of the implementation of DFTB for a long-range corrected functional in generalized Kohn-Sham theory. We apply the method to a set of organic molecules and compare ionization potentials and electron affinities with the original DFTB method and higher level theory. The new scheme cures the significant overpolarization in electric fields found for local DFTB, which parallels the functional dependence in first principles density functional theory (DFT). At the same time, the computational savings with respect to full DFT calculations are not compromised as evidenced by numerical benchmark data.
NASA Technical Reports Server (NTRS)
Ingels, F. M.; Schoggen, W. O.
1982-01-01
The design to achieve the required bit transition density for the Space Shuttle high rate multiplexes (HRM) data stream of the Space Laboratory Vehicle is reviewed. It contained a recommended circuit approach, specified the pseudo random (PN) sequence to be used and detailed the properties of the sequence. Calculations showing the probability of failing to meet the required transition density were included. A computer simulation of the data stream and PN cover sequence was provided. All worst case situations were simulated and the bit transition density exceeded that required. The Preliminary Design Review and the critical Design Review are documented. The Cover Sequence Generator (CSG) Encoder/Decoder design was constructed and demonstrated. The demonstrations were successful. All HRM and HRDM units incorporate the CSG encoder or CSG decoder as appropriate.
NASA Technical Reports Server (NTRS)
Ingels, F.; Schoggen, W. O.
1981-01-01
The various methods of high bit transition density encoding are presented, their relative performance is compared in so far as error propagation characteristics, transition properties and system constraints are concerned. A computer simulation of the system using the specific PN code recommended, is included.
Statistical physics of regular low-density parity-check error-correcting codes
NASA Astrophysics Data System (ADS)
Murayama, Tatsuto; Kabashima, Yoshiyuki; Saad, David; Vicente, Renato
2000-08-01
A variation of Gallager error-correcting codes is investigated using statistical mechanics. In codes of this type, a given message is encoded into a codeword that comprises Boolean sums of message bits selected by two randomly constructed sparse matrices. The similarity of these codes to Ising spin systems with random interaction makes it possible to assess their typical performance by analytical methods developed in the study of disordered systems. The typical case solutions obtained via the replica method are consistent with those obtained in simulations using belief propagation decoding. We discuss the practical implications of the results obtained and suggest a computationally efficient construction for one of the more practical configurations.
NASA Astrophysics Data System (ADS)
Ferullo, Ricardo M.; Domancich, Nicolás F.; Castellani, Norberto J.
2010-11-01
The atomic hydrogen physisorption on graphite was studied using the hydrogen-coronene model system and the van der Waals corrected-density functional theory (DFT + vdW). The results show that H preferentially occupies the hollow site. The adsorption energy at this site is calculated as 38.1 meV, in very good agreement with the available experimental measurements on a single graphite layer (39.2 ± 0.5 meV) and with reported MP2/aug-cc-pVDZ calculations (39.7 meV). The results suggest that, in DFT simulations, dispersion corrections should be considered in order to obtain accurate distances, adsorption energies and diffusion barriers in physisorption processes such as those occurring in the cold interstellar medium.
Importance of far-field Topographic and Isostatic corrections for regional density modeling
NASA Astrophysics Data System (ADS)
Szwillus, Ebbing, Holzrichter
2016-07-01
The long-wavelength gravity field contains information about processes in the sub-lithospheric mantle. As satellite-derived gravity models now provide the long to medium-wavelength gravity field at unprecedented accuracy, techniques used to process gravity data need to be updated. We show that when determining these long-wavelengths, the treatment of topographic and isostatic effects is a likely source of error. We constructed a global isostatic model and calculated global topographic and isostatic effect. These calculations were done for ground stations as well as stations at satellite height. We considered both gravity and gravity gradients. Using these results, we determined how much of the gravity signal comes from distant sources. We find that a significant long-wavelength bias is introduced if far-field effects on the topographic effect are neglected. However, due to isostatic compensation far-field effects of the topographic effect are to a large degree compensated by the far-field isostatic effect. This means that far-field effects can be reduced effectively by always considering topographic masses together with their compensating isostatic masses. We show that to correctly represent the ultra-long wavelengths, a global background model should be used. This is demonstrated both globally and for a continental-scale case area in North America. In the case of regional modeling, where the ultra-long wavelengths are not of prime importance, gravity gradients can be used to help minimize correction errors caused by far-field effects.
Angyan, Janos G.; Gerber, Iann C.; Savin, Andreas; Toulouse, Julien
2005-07-15
Long-range exchange and correlation effects, responsible for the failure of currently used approximate density functionals in describing van der Waals forces, are taken into account explicitly after a separation of the electron-electron interaction in the Hamiltonian into short- and long-range components. We propose a 'range-separated hybrid' functional based on a local density approximation for the short-range exchange-correlation energy, combined with a long-range exact exchange energy. Long-range correlation effects are added by a second-order perturbational treatment. The resulting scheme is general and is particularly well adapted to describe van der Waals complexes, such as rare gas dimers.
Corrections to fringe distortion due to flow density gradients in optical interferometry
NASA Technical Reports Server (NTRS)
Cho, Y. C.; Carr, L. W.; Chandrasekhara, M. S.
1993-01-01
An analytical method is formulated to account for distortions of optical interferograms used for studies of flow over airfoils experiencing dynamic stall. It is shown that such distortions are generated primarily due to optical path deflections in the test flow, caused by large density gradients. Such in-flow optical path deflections are neglected in conventional optical techniques for flow studies. The present method employs a ray analysis to determine these in-flow optical path deflections, and accurately predicts the interferogram distortions.
van de Streek, Jacco; Neumann, Marcus A.
2014-01-01
In 2010 we energy-minimized 225 high-quality single-crystal (SX) structures with dispersion-corrected density functional theory (DFT-D) to establish a quantitative benchmark. For the current paper, 215 organic crystal structures determined from X-ray powder diffraction (XRPD) data and published in an IUCr journal were energy-minimized with DFT-D and compared to the SX benchmark. The on average slightly less accurate atomic coordinates of XRPD structures do lead to systematically higher root mean square Cartesian displacement (RMSCD) values upon energy minimization than for SX structures, but the RMSCD value is still a good indicator for the detection of structures that deserve a closer look. The upper RMSCD limit for a correct structure must be increased from 0.25 Å for SX structures to 0.35 Å for XRPD structures; the grey area must be extended from 0.30 to 0.40 Å. Based on the energy minimizations, three structures are re-refined to give more precise atomic coordinates. For six structures our calculations provide the missing positions for the H atoms, for five structures they provide corrected positions for some H atoms. Seven crystal structures showed a minor error for a non-H atom. For five structures the energy minimizations suggest a higher space-group symmetry. For the 225 SX structures, the only deviations observed upon energy minimization were three minor H-atom related issues. Preferred orientation is the most important cause of problems. A preferred-orientation correction is the only correction where the experimental data are modified to fit the model. We conclude that molecular crystal structures determined from powder diffraction data that are published in IUCr journals are of high quality, with less than 4% containing an error in a non-H atom. PMID:25449625
Brandenburg, Jan Gerit; Alessio, Maristella; Civalleri, Bartolomeo; Peintinger, Michael F; Bredow, Thomas; Grimme, Stefan
2013-09-26
We extend the previously developed geometrical correction for the inter- and intramolecular basis set superposition error (gCP) to periodic density functional theory (DFT) calculations. We report gCP results compared to those from the standard Boys-Bernardi counterpoise correction scheme and large basis set calculations. The applicability of the method to molecular crystals as the main target is tested for the benchmark set X23. It consists of 23 noncovalently bound crystals as introduced by Johnson et al. (J. Chem. Phys. 2012, 137, 054103) and refined by Tkatchenko et al. (J. Chem. Phys. 2013, 139, 024705). In order to accurately describe long-range electron correlation effects, we use the standard atom-pairwise dispersion correction scheme DFT-D3. We show that a combination of DFT energies with small atom-centered basis sets, the D3 dispersion correction, and the gCP correction can accurately describe van der Waals and hydrogen-bonded crystals. Mean absolute deviations of the X23 sublimation energies can be reduced by more than 70% and 80% for the standard functionals PBE and B3LYP, respectively, to small residual mean absolute deviations of about 2 kcal/mol (corresponding to 13% of the average sublimation energy). As a further test, we compute the interlayer interaction of graphite for varying distances and obtain a good equilibrium distance and interaction energy of 6.75 Å and -43.0 meV/atom at the PBE-D3-gCP/SVP level. We fit the gCP scheme for a recently developed pob-TZVP solid-state basis set and obtain reasonable results for the X23 benchmark set and the potential energy curve for water adsorption on a nickel (110) surface. PMID:23947824
Band-gap corrected density functional theory calculations for InAs/GaSb type II superlattices
Wang, Jianwei; Zhang, Yong
2014-12-07
We performed pseudopotential based density functional theory (DFT) calculations for GaSb/InAs type II superlattices (T2SLs), with bandgap errors from the local density approximation mitigated by applying an empirical method to correct the bulk bandgaps. Specifically, this work (1) compared the calculated bandgaps with experimental data and non-self-consistent atomistic methods; (2) calculated the T2SL band structures with varying structural parameters; (3) investigated the interfacial effects associated with the no-common-atom heterostructure; and (4) studied the strain effect due to lattice mismatch between the two components. This work demonstrates the feasibility of applying the DFT method to more exotic heterostructures and defect problems related to this material system.
NASA Technical Reports Server (NTRS)
Ingels, F.; Schoggen, W. O.
1981-01-01
Several methods for increasing bit transition densities in a data stream are summarized, discussed in detail, and compared against constraints imposed by the 2 MHz data link of the space shuttle high rate multiplexer unit. These methods include use of alternate pulse code modulation waveforms, data stream modification by insertion, alternate bit inversion, differential encoding, error encoding, and use of bit scramblers. The psuedo-random cover sequence generator was chosen for application to the 2 MHz data link of the space shuttle high rate multiplexer unit. This method is fully analyzed and a design implementation proposed.
NASA Astrophysics Data System (ADS)
Pederson, Mark R.; Heaton, Richard A.; Harrison, Joseph G.
1989-01-01
The uniform-density electron gas is studied within the framework of the Wannier-function (WF) formulation of the self-interaction-corrected local-spin-density approximation (SIC-LSD). While the results of the present work follow rigorously from a variational formulation, they may also be qualitatively understood in terms of the local-bonding-site concept introduced by Mott in his theory of the metal-insulator transition. SIC-LSD admits metallic-state solutions at ordinary electron densities just as in traditional LSD theory. The result of introducing SIC to the metallic state is an overall downward shift of the LSD eigenvalues. This shift is largest for states near k=0 and vanishes for states near the Fermi energy ɛF. As such, the orbital energies at ɛF are found to be in exact agreement with both the exchange-only version of LSD and Hartree-Fock (HF). Beyond metallic-state solutions, this formulation of SIC-LSD also admits insulator solutions at very low electron densities and may thus have important application to the problem of Wigner crystallization.
NASA Astrophysics Data System (ADS)
Mattsson, A. E.; Mattsson, T. R.; Jennison, D. R.
2003-03-01
Based on the correction scheme presented in [1] we have developed a procedure to correct for the surface self-energy error (both exchange and correlation) in density functional calculations on real systems. The method works equally well for all approximations of the exchange-correlation functional, e.g. the local density and general gradient approximations. It has been successfully applied to Al, Pt, Pd, and Mo vacancy formation energies [2,3] and the Pd(111)/α-alumina work of adhesion [4]. We present the current status of our efforts and discuss how to extend the procedure to general systems. Sandia is a multiprogram laboratory operated by Sandia Corporation, a Lockheed Martin Company, for the United States Department of Energy under Contract DE-AC04-94AL85000. [1] A. E. Mattsson, W. Kohn, J. Chem. Phys. 115, 3441 (2001). [2] K. Carling et al, Phys. Rev. Lett. 85, 3862 (2000). [3] T. R. Mattsson, A. E. Mattsson, Phys. Rev. B (Dec. 2002). [4] A. E. Mattsson, D. R. Jennison, Surf. Sci. 58, L611 (2002).
Lill, Sten O Nilsson
2010-09-01
An evaluation of a dispersion-corrected density functional theory method (B3LYP-DCP) [I.D. Mackie, G.A. DiLabio, Interactions in large, polyaromatic hydrocarbon dimers: application of density functional theory with dispersion corrections, J. Phys. Chem. A 112 (2008) 10968-10976] for three systems of biochemical interest is presented. Firstly, structures and energies of isomers of the tripeptide Phe-Gly-Phe have been compared with CCSD(T)/CBS//RI-MP2/cc-pVTZ literature values. In the system aromatic interactions compete with XH-pi (X=C, N) interactions and hydrogen bonds which makes it a reliable model for proteins. The resulting mean absolute deviation between B3LYP-DCP and CCSD(T)/CBS relative energies is found to be 0.50 kcal mol(-1). Secondly, a phenylalanine derivative featuring a CH-pi interaction has been investigated. A comparison between the optimized geometry and X-ray crystal data shows that B3LYP-DCP accurately predicts the interaction between the two aromatic rings. Thirdly, the dipeptide Ac-Phe-Phe-NH(2) which contains an edge-to-face interaction between two aromatic rings has been studied. The study demonstrates the general applicability of the B3LYP-DCP method on systems which features interactions typically present in biochemical compounds. PMID:20609606
A correction for the Hartree-Fock density of states for jellium without screening.
Blair, Alexander I; Kroukis, Aristeidis; Gidopoulos, Nikitas I
2015-02-28
We revisit the Hartree-Fock (HF) calculation for the uniform electron gas, or jellium model, whose predictions--divergent derivative of the energy dispersion relation and vanishing density of states (DOS) at the Fermi level--are in qualitative disagreement with experimental evidence for simple metals. Currently, this qualitative failure is attributed to the lack of screening in the HF equations. Employing Slater's hyper-Hartree-Fock (HHF) equations, derived variationally, to study the ground state and the excited states of jellium, we find that the divergent derivative of the energy dispersion relation and the zero in the DOS are still present, but shifted from the Fermi wavevector and energy of jellium to the boundary between the set of variationally optimised and unoptimised HHF orbitals. The location of this boundary is not fixed, but it can be chosen to lie at arbitrarily high values of wavevector and energy, well clear from the Fermi level of jellium. We conclude that, rather than the lack of screening in the HF equations, the well-known qualitative failure of the ground-state HF approximation is an artifact of its nonlocal exchange operator. Other similar artifacts of the HF nonlocal exchange operator, not associated with the lack of electronic correlation, are known in the literature. PMID:25725721
A correction for the Hartree-Fock density of states for jellium without screening
NASA Astrophysics Data System (ADS)
Blair, Alexander I.; Kroukis, Aristeidis; Gidopoulos, Nikitas I.
2015-02-01
We revisit the Hartree-Fock (HF) calculation for the uniform electron gas, or jellium model, whose predictions—divergent derivative of the energy dispersion relation and vanishing density of states (DOS) at the Fermi level—are in qualitative disagreement with experimental evidence for simple metals. Currently, this qualitative failure is attributed to the lack of screening in the HF equations. Employing Slater's hyper-Hartree-Fock (HHF) equations, derived variationally, to study the ground state and the excited states of jellium, we find that the divergent derivative of the energy dispersion relation and the zero in the DOS are still present, but shifted from the Fermi wavevector and energy of jellium to the boundary between the set of variationally optimised and unoptimised HHF orbitals. The location of this boundary is not fixed, but it can be chosen to lie at arbitrarily high values of wavevector and energy, well clear from the Fermi level of jellium. We conclude that, rather than the lack of screening in the HF equations, the well-known qualitative failure of the ground-state HF approximation is an artifact of its nonlocal exchange operator. Other similar artifacts of the HF nonlocal exchange operator, not associated with the lack of electronic correlation, are known in the literature.
A correction for the Hartree-Fock density of states for jellium without screening
Blair, Alexander I.; Kroukis, Aristeidis; Gidopoulos, Nikitas I.
2015-02-28
We revisit the Hartree-Fock (HF) calculation for the uniform electron gas, or jellium model, whose predictions—divergent derivative of the energy dispersion relation and vanishing density of states (DOS) at the Fermi level—are in qualitative disagreement with experimental evidence for simple metals. Currently, this qualitative failure is attributed to the lack of screening in the HF equations. Employing Slater’s hyper-Hartree-Fock (HHF) equations, derived variationally, to study the ground state and the excited states of jellium, we find that the divergent derivative of the energy dispersion relation and the zero in the DOS are still present, but shifted from the Fermi wavevector and energy of jellium to the boundary between the set of variationally optimised and unoptimised HHF orbitals. The location of this boundary is not fixed, but it can be chosen to lie at arbitrarily high values of wavevector and energy, well clear from the Fermi level of jellium. We conclude that, rather than the lack of screening in the HF equations, the well-known qualitative failure of the ground-state HF approximation is an artifact of its nonlocal exchange operator. Other similar artifacts of the HF nonlocal exchange operator, not associated with the lack of electronic correlation, are known in the literature.
Zurek, Magdalena; Sladen, Louise; Johansson, Edvin; Olsson, Marita; Jackson, Sonya; Zhang, Hui; Mayer, Gaell; Hockings, Paul D.
2016-01-01
Purpose A magnetic resonance imaging method is presented that allows for the simultaneous assessment of oxygen delivery, oxygen uptake, and parenchymal density. The technique is applied to a mouse model of porcine pancreatic elastase (PPE) induced lung emphysema in order to investigate how structural changes affect lung function. Method Nine-week-old female C57BL6 mice were instilled with saline or PPE at days 0 and 7. At day 19, oxygen delivery, oxygen uptake, and lung density were quantified from T1 and proton-density measurements obtained via oxygen-enhanced magnetic resonance imaging (OE-MRI) using an ultrashort echo-time imaging sequence. Subsequently, the lungs were sectioned for histological observation. Blood-gas analyses and pulmonary functional tests via FlexiVent were performed in separate cohorts. Principal Findings PPE-challenged mice had reduced density when assessed via MRI, consistent with the parenchyma loss observed in the histology sections, and an increased lung compliance was detected via FlexiVent. The oxygenation levels, as assessed via the blood-gas analysis, showed no difference between PPE-challenged animals and control. This finding was mirrored in the global MRI assessments of oxygen delivery and uptake, where the changes in relaxation time indices were matched between the groups. The heterogeneity of the same parameters however, were increased in PPE-challenged animals. When the oxygenation status was investigated in regions of varying density, a reduced oxygen-uptake was found in low-density regions of PPE-challenged mice. In high-density regions the uptake was higher than that of regions of corresponding density in control animals. The oxygen delivery was proportional to the oxygen uptake in both groups. Conclusions The proposed method allowed for the regional assessment of the relationship between lung density and two aspects of lung function, the oxygen delivery and uptake. When compared to global indices of lung function, an
Allodji, Rodrigue S; Thiébaut, Anne C M; Leuraud, Klervi; Rage, Estelle; Henry, Stéphane; Laurier, Dominique; Bénichou, Jacques
2012-12-30
A broad variety of methods for measurement error (ME) correction have been developed, but these methods have rarely been applied possibly because their ability to correct ME is poorly understood. We carried out a simulation study to assess the performance of three error-correction methods: two variants of regression calibration (the substitution method and the estimation calibration method) and the simulation extrapolation (SIMEX) method. Features of the simulated cohorts were borrowed from the French Uranium Miners' Cohort in which exposure to radon had been documented from 1946 to 1999. In the absence of ME correction, we observed a severe attenuation of the true effect of radon exposure, with a negative relative bias of the order of 60% on the excess relative risk of lung cancer death. In the main scenario considered, that is, when ME characteristics previously determined as most plausible from the French Uranium Miners' Cohort were used both to generate exposure data and to correct for ME at the analysis stage, all three error-correction methods showed a noticeable but partial reduction of the attenuation bias, with a slight advantage for the SIMEX method. However, the performance of the three correction methods highly depended on the accurate determination of the characteristics of ME. In particular, we encountered severe overestimation in some scenarios with the SIMEX method, and we observed lack of correction with the three methods in some other scenarios. For illustration, we also applied and compared the proposed methods on the real data set from the French Uranium Miners' Cohort study. PMID:22996087
Yamakawa, Youichi; Kontani, Hiroshi
2015-06-26
We present a microscopic derivation of the nematic charge-density wave (CDW) formation in cuprate superconductors based on the three-orbital d-p Hubbard model by introducing the vertex correction (VC) into the charge susceptibility. The CDW instability at q=(Δ(FS),0), (0,Δ(FS)) appears when the spin fluctuations are strong, due to the strong charge-spin interference represented by the VC. Here, Δ(FS) is the wave number between the neighboring hot spots. The obtained spin-fluctuation-driven CDW is expressed as the "intra-unit-cell orbital order" accompanied by the charge transfer between the neighboring atomic orbitals, which is actually observed by the scanning tunneling microscope measurements. We predict that the cuprate CDW and the nematic orbital order in Fe-based superconductors are closely related spin-fluctuation-driven phenomena. PMID:26197139
NASA Astrophysics Data System (ADS)
Hahn, Torsten; Rückerl, Florian; Liebing, Simon; Pederson, Mark
We present our experimental and theoretical results on novel Picene/F4TCNQ and Manganese-Phthalocyanine/F4TCNQ donor / acceptor systems. We apply the recently developed Fermi-orbital based approach for self-interaction corrected density functional theory (FO-SIC DFT) to these materials and compare the results to standard DFT calculations and to experimental data obtained by photoemission spectroscopy. We focus our analysis on the description of the magnitude of the ground state charge transfer and on the details of the formed hybrid orbitals. Further, we show that for weakly bound donor / acceptor systems the FO-SIC approach delivers a more realistic description of the electronic structure compared to standard DFT calculations Support by DFG FOR1154 is greatly acknowledged.
Baer, Marcel D.; Mundy, Christopher J.; McGrath, Matthew J.; Kuo, I-F W.; Siepmann, Joern I.; Tobias, Douglas J.
2011-09-28
First-principles molecular dynamics simulations, in which the forces are computed from electronic structure calculations, have great potential to provide unique insight into structure, dynamics, electronic properties, and chemistry at aqueous interfaces that is not available from empirical force fields. The majority of current first-principles simulations are driven by forces derived from density functional theory with generalized gradient approximations to the exchange-correlation energy, which do not capture dispersion interactions. We have carried out first-principles molecular dynamics simulations of air-water interfaces employing a particular generalized gradient approximation to the exchange-correlation functional (BLYP), with and without empirical dispersion corrections. We assess the utility of the dispersion corrections by comparison of a variety of structural, dynamic, and thermodynamic properties of bulk and interfacial water with experimental data, as well as other first-principles and force field-based simulations. This work was supported by the US Department of Energy, Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences & Biosciences. Pacific Northwest National Laboratory is operated by Battelle for the US Department of Energy.
NASA Astrophysics Data System (ADS)
Zhong, X.; Rungger, I.; Zapol, P.; Heinonen, O.
2015-03-01
Understanding electronic properties of substoichiometric phases of titanium oxide such as Magnéli phase T i4O7 is crucial in designing and modeling resistive switching devices. Here we present our study on Magnéli phase T i4O7 together with rutile Ti O2 and T i2O3 using density functional theory methods with atomic-orbital-based self-interaction correction (ASIC). We predict a new antiferromagnetic (AF) ground state in the low temperature (LT) phase, and we explain energy difference with a competing AF state using a Heisenberg model. The predicted energy ordering of these states in the LT phase is calculated to be robust in a wide range of modeled isotropic strain. We have also investigated the dependence of the electronic structures of the Ti-O phases on stoichiometry. The splitting of titanium t2 g orbitals is enhanced with increasing oxygen deficiency as Ti-O is reduced. The electronic properties of all these phases can be reasonably well described by applying ASIC with a "standard" value for transition metal oxides of the empirical parameter α of 0.5 representing the magnitude of the applied self-interaction correction.
Rotstein, S.; Lax, I.; Svane, G. )
1990-01-01
The relative electron density of lung tissue was measured from computer tomography (CT) slices in 33 breast cancer patients treated by various techniques of adjuvant radiotherapy. The measurements were made before radiotherapy, 3 months and 9 months after completion of radiation therapy. The changes in lung densities at 3 months and 9 months were compared to radiation induced radiological (CT) findings. In addition, subjective symptoms such as cough and dyspnoea were assessed before and after radiotherapy. It was observed that the mean of the relative electron density of lung tissue varied from 0.25 when the whole lung was considered to 0.17 when only the anterior lateral quarter of the lung was taken into account. In patients with positive radiological (CT) findings the mean lung density of the anterior lateral quarter increased 2.1 times 3 months after radiotherapy and was still increased 1.6 times 6 months later. For those patients without findings, in the CT pictures the corresponding values were 1.2 and 1.1, respectively. The standard deviation of the pixel values within the anterior lateral quarter of the lung increased 3.8 times and 3.2 times at 3 months and 9 months, respectively, in the former group, as opposed to 1.2 and 1.1 in the latter group. Thirteen patients had an increase in either cough or dyspnoea as observed 3 months after completion of radiotherapy. In eleven patients these symptoms persisted 6 months later. No significant correlation was found between radiological findings and subjective symptoms. However, when three different treatment techniques were compared among 29 patients the highest rate of radiological findings was observed in patients in which the largest lung volumes received the target dose. A tendency towards an increased rate of subjective symptoms was also found in this group.
NASA Astrophysics Data System (ADS)
Fujino, Takeshi; Maeda, Hiroshi; Moriizumi, Koichi; Kato, Takaaki; Tsubouchi, Natsuro
1994-12-01
Proximity effect correction for electron beam (EB) direct writing is studied in this paper. An iterative calculation of exposure-dose modulation by equalizing the deposited energy of all figures requires an extremely long calculation time, especially in the case of high EB acceleration voltage. Therefore it is not practical for the correction of LSI. The correction method using pattern-area density, however, could realize the high-speed proximity effect correction. In this paper, this method is first investigated from the standpoint of the correction accuracy. Next, the applicability of this method to two critical cases is examined. One is the patterning on a heavy-metal substrate, on which backscattering yield of electrons is high. The other is the application to the cell-projection exposure, in which it is impossible to modulate exposure dose inside the cell. Lastly, the calculation time of the proximity effect correction is evaluated for 64 Mbit dynamic random access memories (MbDRAMs).
Lo, P.; Young, S.; Kim, H. J.; Brown, M. S.
2016-01-01
Purpose: To investigate the effects of dose level and reconstruction method on density and texture based features computed from CT lung nodules. Methods: This study had two major components. In the first component, a uniform water phantom was scanned at three dose levels and images were reconstructed using four conventional filtered backprojection (FBP) and four iterative reconstruction (IR) methods for a total of 24 different combinations of acquisition and reconstruction conditions. In the second component, raw projection (sinogram) data were obtained for 33 lung nodules from patients scanned as a part of their clinical practice, where low dose acquisitions were simulated by adding noise to sinograms acquired at clinical dose levels (a total of four dose levels) and reconstructed using one FBP kernel and two IR kernels for a total of 12 conditions. For the water phantom, spherical regions of interest (ROIs) were created at multiple locations within the water phantom on one reference image obtained at a reference condition. For the lung nodule cases, the ROI of each nodule was contoured semiautomatically (with manual editing) from images obtained at a reference condition. All ROIs were applied to their corresponding images reconstructed at different conditions. For 17 of the nodule cases, repeat contours were performed to assess repeatability. Histogram (eight features) and gray level co-occurrence matrix (GLCM) based texture features (34 features) were computed for all ROIs. For the lung nodule cases, the reference condition was selected to be 100% of clinical dose with FBP reconstruction using the B45f kernel; feature values calculated from other conditions were compared to this reference condition. A measure was introduced, which the authors refer to as Q, to assess the stability of features across different conditions, which is defined as the ratio of reproducibility (across conditions) to repeatability (across repeat contours) of each feature. Results: The
A human embryonic lung fibroblast with a high density of muscarinic acetylcholine receptors.
André, C; Marullo, S; Convents, A; Lü, B Z; Guillet, J G; Hoebeke, J; Strosberg, D A
1988-01-15
Binding studies with the radiolabeled muscarinic antagonists dexetimide, quinuclidinyl benzilate and N-methylscopolamine showed that the human embryonic lung fibroblast CCL137 possesses approximately 2 X 10(5) muscarinic receptors/cell, i.e. 2.1 pmol/mg membrane protein. These receptors showed a marked stereoselectivity towards dexetimide and levetimide and only low affinity for another antagonist, pirenzepine. The muscarinic agonist carbamylcholine inhibited forskolin-stimulated adenylate cyclase and induced phosphatidylinositide turnover in the intact cells. Both effects were inhibited by the muscarinic antagonist atropine. Affinity labeling with tritiated propylbenzylcholine mustard revealed a protein of 72 kDa. Finally, down-regulation of the membrane receptors following prolonged treatment with the agonist carbamylcholine was assessed by means of the hydrophilic antagonist N-methylscopolamine. PMID:2828056
Aoki, Masahiko; Hirose, Katsumi; Sato, Mariko; Akimoto, Hiroyoshi; Kawaguchi, Hideo; Hatayama, Yoshiomi; Fujioka, Ichitaro; Tanaka, Mitsuki; Ono, Shuichi; Takai, Yoshihiro
2016-07-01
The purpose of this study was to investigate the prognostic significance of average iodine density as assessed by dual-energy computed tomography (DE-CT) for lung tumors treated with stereotactic body radiotherapy (SBRT). From March 2011 to August 2014, 93 medically inoperable patients with 74 primary lung cancers and 19 lung metastases underwent DE-CT prior to SBRT of a total dose of 45-60 Gy in 5-10 fractions. Of these 93 patients, nine patients had two lung tumors. Thus, 102 lung tumors were included in this study. DE-CT was performed for pretreatment evaluation. Regions of interest were set for the entire tumor, and average iodine density was obtained using a dedicated imaging software and evaluated with regard to local control. The median follow-up period was 23.4 months (range, 1.5-54.5 months). The median value of the average iodine density was 1.86 mg/cm(3) (range, 0.40-9.27 mg/cm(3)). Two-year local control rates for the high and low average iodine density groups divided by the median value of the average iodine density were 96.9% and 75.7% (P = 0.006), respectively. Tumors with lower average iodine density showed a worse prognosis, possibly reflecting a hypoxic cell population in the tumor. The average iodine density exhibited a significant impact on local control. Our preliminary results indicate that iodine density evaluated using dual-energy spectral CT may be a useful, noninvasive and quantitative assessment of radio-resistance caused by presumably hypoxic cell populations in tumors. PMID:26826198
Aoki, Masahiko; Hirose, Katsumi; Sato, Mariko; Akimoto, Hiroyoshi; Kawaguchi, Hideo; Hatayama, Yoshiomi; Fujioka, Ichitaro; Tanaka, Mitsuki; Ono, Shuichi; Takai, Yoshihiro
2016-01-01
The purpose of this study was to investigate the prognostic significance of average iodine density as assessed by dual-energy computed tomography (DE-CT) for lung tumors treated with stereotactic body radiotherapy (SBRT). From March 2011 to August 2014, 93 medically inoperable patients with 74 primary lung cancers and 19 lung metastases underwent DE-CT prior to SBRT of a total dose of 45–60 Gy in 5–10 fractions. Of these 93 patients, nine patients had two lung tumors. Thus, 102 lung tumors were included in this study. DE-CT was performed for pretreatment evaluation. Regions of interest were set for the entire tumor, and average iodine density was obtained using a dedicated imaging software and evaluated with regard to local control. The median follow-up period was 23.4 months (range, 1.5–54.5 months). The median value of the average iodine density was 1.86 mg/cm3 (range, 0.40–9.27 mg/cm3). Two-year local control rates for the high and low average iodine density groups divided by the median value of the average iodine density were 96.9% and 75.7% (P = 0.006), respectively. Tumors with lower average iodine density showed a worse prognosis, possibly reflecting a hypoxic cell population in the tumor. The average iodine density exhibited a significant impact on local control. Our preliminary results indicate that iodine density evaluated using dual-energy spectral CT may be a useful, noninvasive and quantitative assessment of radio-resistance caused by presumably hypoxic cell populations in tumors. PMID:26826198
Nahak, Prasant; Nag, Kaushik; Hillier, Ashley; Devraj, Ravi; Thompson, David W; Manna, Kausik; Makino, Kimiko; Ohshima, Hiroyuki; Nakahara, Hiromichi; Shibata, Osamu; Panda, Amiya Kumar
2014-01-01
Lung surfactant is a complex mixture of lipid and protein, responsible for alveolar stability, becomes dysfunctional due to alteration of its structure and function by leaked serum materials in disease. Serum proteins, cholesterol and low density lipoprotein (LDL) were studied with bovine lipid extract surfactant (BLES) using Langmuir films, and bilayer dispersions using Raman spectroscopy. While small amount of cholesterol (10 wt%) and LDL did not significantly affect the adsorption and surface tension lowering properties of BLES. However serum lipids, whole serum as well as higher amounts of cholesterol, and LDL dramatically altered the surface properties of BLES films, as well as gel-fluid structures formed in such films observed using atomic force microscopy (AFM). Raman-spectroscopic studies revealed that serum proteins, LDL and excess cholesterol had fluidizing effects on BLES bilayers dispersion, monitored from the changes in hydrocarbon vibrational modes during gel-fluid thermal phase transitions. This study clearly suggests that patho-physiological amounts of serum lipids (and not proteins) significantly alter the molecular arrangement of surfactant in films and bilayers, and can be used to model lung disease. PMID:25409691
Gebhardt, Julian; Viñes, Francesc; Bleiziffer, Patrick; Hieringer, Wolfgang; Görling, Andreas
2014-03-21
We investigate the capability of low-coordinated sites on small model clusters to act as active centers for hydrogen storage. A set of small magic clusters with the formula (XY)6 (X = Mg, Ba, Be, Zn, Cd, Na, Li, B and Y = O, Se, S, F, I, N) and a "drumlike" hexagonal shape showing a low coordination number of three was screened. Oxide clusters turned out to be the most promising candidates for hydrogen storage. For these ionic compounds we explored the suitability of different van der Waals (vdW) corrections to density-functional calculations by comparing the respective H2 physisorption profile to highly accurate CCSD(T) (Coupled Cluster Singles Doubles with perturbative Triples) calculations. The Grimme D3 vdW correction in combination with the Perdew-Burke-Ernzerhof exchange-correlation functional was found to be the best approach compared to CCSD(T) hydrogen physisorption profiles and is, therefore, suited to study these and other light metal oxide systems. H2 adsorption on sites of oxide model clusters is found to meet the adsorption energy criteria for H2 storage, with bond strengths ranging from 0.15 to 0.21 eV. Energy profiles and estimates of kinetic constants for the H2 splitting reaction reveal that H2 is likely to be adsorbed molecularly on sites of (MgO)6, (BaO)6, and (BeO)6 clusters, suggesting a rapid H2 uptake/release at operating temperatures and moderate pressures. The small mass of beryllium and magnesium makes such systems appealing for meeting the gravimetric criterion for H2 storage. PMID:24499810
NASA Astrophysics Data System (ADS)
Greuter, Marcel J. W.; Groen, Jaap M.; Nicolai, Lieuwe J.; Dijkstra, Hildebrand; Oudkerk, Matthijs
2009-02-01
Objective: To quantify the influence of velocity, calcification density and acquisition time on coronary calcium determination using multi-detector CT, dual-source CT and EBT. Materials and Methods: Artificial arteries with four calcifications of increasing density were attached to a robotic arm to which a linear movement was applied between 0 and 120 mm/s (step 10 mm/s). The phantom was scanned five times on 64-slice MDCT, DSCT and EBT using a standard acquisition protocol and the average Agatston score was determined. Results: Increasing motion artifacts were observed at increasing velocities on all scanners, with increasing severity from EBT to DSCT to 64-slice MDCT. The Agatston score showed a linear dependency on velocity from which a correction factor was derived. This correction factor showed a linear dependency on calcification density (0.92<=R2<=0.95). The slope and offset of this correction factor also showed a linear dependency on acquisition time (0.84<=R2<=0.86). Conclusion: The Agatston score is highly dependent on the average density of individual calcifications. The dependency of the Agatston score on velocity shows a linear behaviour on calcification density. A quantitative method could be derived which corrects the measured calcium score for the influence of velocity, calcification density and acquisition time.
Tsukihara, Masayoshi; Noto, Yoshiyuki; Sasamoto, Ryuta; Hayakawa, Takahide; Saito, Masatoshi
2015-03-15
Purpose: To achieve accurate tissue inhomogeneity corrections in radiotherapy treatment planning, the authors had previously proposed a novel conversion of the energy-subtracted computed tomography (CT) number to an electron density (ΔHU–ρ{sub e} conversion), which provides a single linear relationship between ΔHU and ρ{sub e} over a wide range of ρ{sub e}. The purpose of this study is to present an initial implementation of the ΔHU–ρ{sub e} conversion method for a treatment planning system (TPS). In this paper, two example radiotherapy plans are used to evaluate the reliability of dose calculations in the ΔHU–ρ{sub e} conversion method. Methods: CT images were acquired using a clinical dual-source CT (DSCT) scanner operated in the dual-energy mode with two tube potential pairs and an additional tin (Sn) filter for the high-kV tube (80–140 kV/Sn and 100–140 kV/Sn). Single-energy CT using the same DSCT scanner was also performed at 120 kV to compare the ΔHU–ρ{sub e} conversion method with a conventional conversion from a CT number to ρ{sub e} (Hounsfield units, HU–ρ{sub e} conversion). Lookup tables for ρ{sub e} calibration were obtained from the CT image acquisitions for tissue substitutes in an electron density phantom (EDP). To investigate the beam-hardening effect on dosimetric uncertainties, two EDPs with different sizes (a body EDP and a head EDP) were used for the ρ{sub e} calibration. Each acquired lookup table was applied to two radiotherapy plans designed using the XiO TPS with the superposition algorithm for an anthropomorphic phantom. The first radiotherapy plan was for an oral cavity tumor and the second was for a lung tumor. Results: In both treatment plans, the performance of the ΔHU–ρ{sub e} conversion was superior to that of the conventional HU–ρ{sub e} conversion in terms of the reliability of dose calculations. Especially, for the oral tumor plan, which dealt with dentition and bony structures, treatment
Pokhrel, D; Badkul, R; Jiang, H; Estes, C; Kumar, P; Wang, F
2014-06-01
Purpose: Lung-SBRT uses hypo-fractionated dose in small non-IMRT fields with tissue-heterogeneity corrected plans. An independent MU verification is mandatory for safe and effective delivery of the treatment plan. This report compares planned MU obtained from iPlan-XVM-Calgorithm against spreadsheet-based hand-calculation using most commonly used simple TMR-based method. Methods: Treatment plans of 15 patients who underwent for MC-based lung-SBRT to 50Gy in 5 fractions for PTV V100%=95% were studied. ITV was delineated on MIP images based on 4D-CT scans. PTVs(ITV+5mm margins) ranged from 10.1- 106.5cc(average=48.6cc). MC-SBRT plans were generated using a combination of non-coplanar conformal arcs/beams using iPlan XVM-Calgorithm (BrainLAB iPlan ver.4.1.2) for Novalis-TX consisting of micro-MLCs and 6MV-SRS (1000MU/min) beam. These plans were re-computed using heterogeneity-corrected Pencil-Beam (PB-hete) algorithm without changing any beam parameters, such as MLCs/MUs. Dose-ratio: PB-hete/MC gave beam-by-beam inhomogeneity-correction-factors (ICFs):Individual Correction. For independent-2nd-check, MC-MUs were verified using TMR-based hand-calculation and obtained an average ICF:Average Correction, whereas TMR-based hand-calculation systematically underestimated MC-MUs by ∼5%. Also, first 10 MC-plans were verified with an ion-chamber measurement using homogenous phantom. Results: For both beams/arcs, mean PB-hete dose was systematically overestimated by 5.5±2.6% and mean hand-calculated MU systematic underestimated by 5.5±2.5% compared to XVMC. With individual correction, mean hand-calculated MUs matched with XVMC by - 0.3±1.4%/0.4±1.4 for beams/arcs, respectively. After average 5% correction, hand-calculated MUs matched with XVMC by 0.5±2.5%/0.6±2.0% for beams/arcs, respectively. Smaller dependence on tumor volume(TV)/field size(FS) was also observed. Ion-chamber measurement was within ±3.0%. Conclusion: PB-hete overestimates dose to lung tumor relative to
NASA Astrophysics Data System (ADS)
Chwee, T. S.; Sullivan, M. B.
2012-10-01
The adsorption energies and changes in surface work functions for benzene on unreconstructed Cu(111), Ag (111), and Au (111) at low coverages have been studied within the framework of dispersion corrected Kohn-Sham density functional theory. Corrections to account for long range dispersive effects between the adsorbate and metal substrate were incorporated via the exchange-hole dipole moment method of Becke and Johnson [J. Chem. Phys. 123, 154101 (2005), 10.1063/1.2065267]. We show that the dispersion corrected calculations yield significantly improved adsorption energies and work function shifts that are in good agreement with experimental values.
Successes and failures of Hubbard-corrected density functional theory: The case of Mg doped LiCoO2
NASA Astrophysics Data System (ADS)
Santana, Juan A.; Kim, Jeongnim; Kent, P. R. C.; Reboredo, Fernando A.
2014-10-01
We have evaluated the successes and failures of the Hubbard-corrected density functional theory approach to study Mg doping of LiCoO2. We computed the effect of the U parameter on the energetic, geometric, and electronic properties of two possible doping mechanisms: (1) substitution of Mg onto a Co (or Li) site with an associated impurity state and (2) formation of impurity-state-free complexes of substitutional Mg and point defects in LiCoO2. We find that formation of impurity states results in changes on the valency of Co in LiCoO2. Variation of the Co U shifts the energy of the impurity state, resulting in energetic, geometric, and electronic properties that depend significantly on the specific value of U. In contrast, the properties of the impurity-state-free complexes are insensitive to U. These results identify reasons for the strong dependence on the doping properties on the chosen value of U and for the overall difficulty of achieving agreement with the experimentally known energetic and electronic properties of doped transition metal oxides such as LiCoO2.
Successes and failures of Hubbard-corrected density functional theory: the case of Mg doped LiCoO2.
Santana, Juan A; Kim, Jeongnim; Kent, P R C; Reboredo, Fernando A
2014-10-28
We have evaluated the successes and failures of the Hubbard-corrected density functional theory approach to study Mg doping of LiCoO2. We computed the effect of the U parameter on the energetic, geometric, and electronic properties of two possible doping mechanisms: (1) substitution of Mg onto a Co (or Li) site with an associated impurity state and (2) formation of impurity-state-free complexes of substitutional Mg and point defects in LiCoO2. We find that formation of impurity states results in changes on the valency of Co in LiCoO2. Variation of the Co U shifts the energy of the impurity state, resulting in energetic, geometric, and electronic properties that depend significantly on the specific value of U. In contrast, the properties of the impurity-state-free complexes are insensitive to U. These results identify reasons for the strong dependence on the doping properties on the chosen value of U and for the overall difficulty of achieving agreement with the experimentally known energetic and electronic properties of doped transition metal oxides such as LiCoO2. PMID:25362331
Successes and failures of Hubbard-corrected density functional theory. The case of Mg doped LiCoO2
Santana Palacio, Juan A.; Kim, Jeongnim; Kent, Paul R.; Reboredo, Fernando A.
2014-10-28
We have evaluated the successes and failures of the Hubbard-corrected density functional theory approach to study Mg doping of LiCoO2. We computed the effect of the U parameter on the energetic, geometric, and electronic properties of two possible doping mechanisms: (1) substitution of Mg onto a Co (or Li) site with an associated impurity state and (2) formation of impurity-state-free complexes of substitutional Mg and point defects in LiCoO2. We find that formation of impurity states results in changes on the valency of Co in LiCoO2. Variation of the Co U shifts the energy of the impurity state, resulting inmore » energetic, geometric, and electronic properties that depend significantly on the specific value of U. In contrast, the properties of the impurity-state-free complexes are insensitive to U. These results identify reasons for the strong dependence on the doping properties on the chosen value of U and for the overall difficulty of achieving agreement with the experimentally known energetic and electronic properties of doped transition metal oxides such as LiCoO2.« less
Zhang, Igor Ying; Xu, Xin
2013-05-16
An unbiased understanding of complex molecular systems from first-principles critically demands theoretical methods with uniform accuracy for diverse interactions with different natures covering short-, medium-, and long-range correlations. Among the state-of-the-art density functional approximations (DFAs), doubly hybrid (DH) DFAs (e.g., XYG3 in this Letter) provide a remarkable improvement over the conventional DFAs (e.g., B3LYP in this Letter). Even though XYG3 works quite well in many cases of noncovalent bonding interactions (NCIs), it is incomplete in describing the pure long-range dispersive interactions. Here, we address such concerns by adding a scaled long-range contribution from the second-order perturbation theory (PT2). The long-range-corrected XYG3 (lrc-XYG3) is proposed without reparameterizing the three parameters in the original XYG3. Due to its overall excellent performance for all testing sets constructed for various purposes, lrc-XYG3 is the recommended method, which is expected to provide a balanced description of diverse interactions in complex molecular systems. PMID:26282977
Bozkaya, Uğur
2016-04-14
An efficient implementation of the asymmetric triples correction for the coupled-cluster singles and doubles [ΛCCSD(T)] method [S. A. Kucharski and R. J. Bartlett, J. Chem. Phys. 108, 5243 (1998); T. D. Crawford and J. F. Stanton, Int. J. Quantum Chem. 70, 601 (1998)] with the density-fitting [DF-ΛCCSD(T)] approach is presented. The computational time for the DF-ΛCCSD(T) method is compared with that of ΛCCSD(T). Our results demonstrate that the DF-ΛCCSD(T) method provide substantially lower computational costs than ΛCCSD(T). Further application results show that the ΛCCSD(T) and DF-ΛCCSD(T) methods are very beneficial for the study of single bond breaking problems as well as noncovalent interactions and transition states. We conclude that ΛCCSD(T) and DF-ΛCCSD(T) are very promising for the study of challenging chemical systems, where the coupled-cluster singles and doubles with perturbative triples method fails. PMID:27083709
NASA Astrophysics Data System (ADS)
Bozkaya, Uǧur
2016-04-01
An efficient implementation of the asymmetric triples correction for the coupled-cluster singles and doubles [ΛCCSD(T)] method [S. A. Kucharski and R. J. Bartlett, J. Chem. Phys. 108, 5243 (1998); T. D. Crawford and J. F. Stanton, Int. J. Quantum Chem. 70, 601 (1998)] with the density-fitting [DF-ΛCCSD(T)] approach is presented. The computational time for the DF-ΛCCSD(T) method is compared with that of ΛCCSD(T). Our results demonstrate that the DF-ΛCCSD(T) method provide substantially lower computational costs than ΛCCSD(T). Further application results show that the ΛCCSD(T) and DF-ΛCCSD(T) methods are very beneficial for the study of single bond breaking problems as well as noncovalent interactions and transition states. We conclude that ΛCCSD(T) and DF-ΛCCSD(T) are very promising for the study of challenging chemical systems, where the coupled-cluster singles and doubles with perturbative triples method fails.
A Correction for the IRI Topside Electron Density Model Based on Alouette/ISIS Topside Sounder Data
NASA Technical Reports Server (NTRS)
Bilitza, D.
2004-01-01
The topside segment of the International Reference Ionosphere (IRI) electron density model (and also of the Bent model) is based on the limited amount of topside data available at the time (40,OOO Alouette 1 profiles). Being established from such a small database it is therefore not surprising that the models have well-known shortcomings, for example, at high solar activities. Meanwhile a large data base of close to 200,000 topside profiles from Alouette 1,2, and ISIS I, 2 has become available online. A program of automated scaling and inversion of a large volume of digitized ionograms adds continuously to this data pool. We have used the currently available ISIs/Alouette topside profiles to evaluate the IRI topside model and to investigate ways of improving the model. The IRI model performs generally well at middle latitudes and shows discrepancies at low and high latitudes and these discrepancies are largest during high solar activity. In the upper topside IRI consistently overestimates the measurements. Based on averages of the data-model ratios we have established correction factors for the IRI model. These factors vary with altitude, modified dip latitude, and local time.
Gohr, Sebastian; Grimme, Stefan; Söhnel, Tilo; Paulus, Beate; Schwerdtfeger, Peter
2013-11-01
First-principles density functional theory (DFT) is used to study the solid-state modifications of carbon dioxide up to pressures of 60 GPa. All known molecular CO2 structures are investigated in this pressure range, as well as three non-molecular modifications. To account for long-range van der Waals interactions, the dispersion corrected DFT method developed by Grimme and co-workers (DFT-D3) is applied. We find that the DFT-D3 method substantially improves the results compared to the uncorrected DFT methods for the molecular carbon dioxide crystals. Enthalpies at 0 K and cohesive energies support only one possibility of the available experimental solutions for the structure of phase IV: the R3c modification, proposed by Datchi and co-workers [Phys. Rev. Lett. 103, 185701 (2009)]. Furthermore, comparing bulk moduli with experimental values, we cannot reproduce the quite large--rather typical for covalent crystal structures--experimental values for the molecular phases II and III. PMID:24206310
APOM and High-Density Lipoprotein are associated with Lung Function and Percent Emphysema
Burkart, Kristin M; Manichaikul, Ani; Wilk, Jemma B; Ahmed, Firas S; Burke, Gregory L; Enright, Paul; Hansel, Nadia N; Haynes, Demondes; Heckbert, Susan R; Hoffman, Eric A; Kaufman, Joel D; Kurai, Jun; Loehr, Laura; London, Stephanie J; Meng, Yang; O’Connor, George T; Oelsner, Elizabeth; Petrini, Marcy; Pottinger, Tess D; Powell, Charles A; Redline, Susan; Rotter, Jerome I; Smith, Lewis J; Artigas, María Soler; Tobin, Martin D; Tsai, Michael Y; Watson, Karol; White, Wendy; Young, Taylor R; Rich, Stephen S; Barr, R Graham
2014-01-01
Chronic obstructive pulmonary disease (COPD) is linked to cardiovascular disease; however, there are few studies on the associations of cardiovascular genes with COPD. We assessed the association of lung function with 2,100 genes selected for cardiovascular diseases among 20,077 European-Americans and 6,900 African-Americans. We performed replication of significant loci in the other racial group and an independent consortium of Europeans, tested the associations of significant loci with percent emphysema, and examined gene expression in an independent sample. We then tested the association of a related lipid biomarker with FEV1/FVC and percent emphysema. We identified one new polymorphism for FEV1/FVC (rs805301) in European-Americans (p=1.3×10−6) and a second (rs707974) in the combined European-American and African-American analysis (p=1.38×10−7). Both SNPs flank the gene for apolipoprotein M (apoM), a component of HDL. Both replicated in an independent cohort. SNPs in a second gene related to apoM and HDL, PCSK9, were associated with FEV1/FVC among African-Americans. rs707974 was associated with percent emphysema among European-Americans and African-Americans, and APOM expression was related to FEV1/FVC and percent emphysema. Higher HDL levels were associated with lower FEV1/FVC and greater percent emphysema. These findings suggest a novel role for the APOM/HDL pathway in the pathogenesis of COPD and emphysema. PMID:23900982
Hao, Feng Mattsson, Ann E.; Armiento, Rickard
2014-05-14
We have previously proposed that further improved functionals for density functional theory can be constructed based on the Armiento-Mattsson subsystem functional scheme if, in addition to the uniform electron gas and surface models used in the Armiento-Mattsson 2005 functional, a model for the strongly confined electron gas is also added. However, of central importance for this scheme is an index that identifies regions in space where the correction provided by the confined electron gas should be applied. The electron localization function (ELF) is a well-known indicator of strongly localized electrons. We use a model of a confined electron gas based on the harmonic oscillator to show that regions with high ELF directly coincide with regions where common exchange energy functionals have large errors. This suggests that the harmonic oscillator model together with an index based on the ELF provides the crucial ingredients for future improved semi-local functionals. For a practical illustration of how the proposed scheme is intended to work for a physical system we discuss monoclinic cupric oxide, CuO. A thorough discussion of this system leads us to promote the cell geometry of CuO as a useful benchmark for future semi-local functionals. Very high ELF values are found in a shell around the O ions, and take its maximum value along the Cu–O directions. An estimate of the exchange functional error from the effect of electron confinement in these regions suggests a magnitude and sign that could account for the error in cell geometry.
NASA Astrophysics Data System (ADS)
Nooruddin, Hasan A.; Anifowose, Fatai; Abdulraheem, Abdulazeez
2014-03-01
Soft computing techniques are recently becoming very popular in the oil industry. A number of computational intelligence-based predictive methods have been widely applied in the industry with high prediction capabilities. Some of the popular methods include feed-forward neural networks, radial basis function network, generalized regression neural network, functional networks, support vector regression and adaptive network fuzzy inference system. A comparative study among most popular soft computing techniques is presented using a large dataset published in literature describing multimodal pore systems in the Arab D formation. The inputs to the models are air porosity, grain density, and Thomeer parameters obtained using mercury injection capillary pressure profiles. Corrected air permeability is the target variable. Applying developed permeability models in recent reservoir characterization workflow ensures consistency between micro and macro scale information represented mainly by Thomeer parameters and absolute permeability. The dataset was divided into two parts with 80% of data used for training and 20% for testing. The target permeability variable was transformed to the logarithmic scale as a pre-processing step and to show better correlations with the input variables. Statistical and graphical analysis of the results including permeability cross-plots and detailed error measures were created. In general, the comparative study showed very close results among the developed models. The feed-forward neural network permeability model showed the lowest average relative error, average absolute relative error, standard deviations of error and root means squares making it the best model for such problems. Adaptive network fuzzy inference system also showed very good results.
Arora, Priya; Moudgil, R. K.; Bhukal, Nisha
2015-05-15
Static density-density correlation function has been calculated for a spin-polarized two-dimensional quantum electron fluid by including the first-order exchange and self-energy corrections to the random-phase approximation (RPA). This is achieved by determining these corrections to the RPA linear density-density response function, obtained by solving the equation of motion for the single-particle Green’s function. Resulting infinite hierarchy of equations (involving higher-order Green’s functions) is truncated by factorizing the two-particle Green’s function as a product of the single-particle Green’s function and one-particle distribution function. Numerical results of correlation function are compared directly against the quantum Monte Carlo simulation data due to Tanatar and Ceperley for different coupling parameter (r{sub s}) values. We find almost exact agreement for r{sub s}=1, with a noticeable improvement over the RPA. Its quality, however, deteriorates with increasing r{sub s}, but correction to RPA is quite significant.
Mammoto, Tadanori; Chen, Zhao; Jiang, Amanda; Jiang, Elisabeth; Ingber, Donald E; Mammoto, Akiko
2016-01-01
Angiogenesis, the growth of new blood vessels, plays a key role in organ development, homeostasis, and regeneration. The cooperation of multiple angiogenic factors, rather than a single factor, is required for physiological angiogenesis. Recently, we have reported that soluble platelet-rich plasma (PRP) extract, which contains abundant angiopoietin-1 and multiple other angiogenic factors, stimulates angiogenesis and maintains vascular integrity in vitro and in vivo. In this report, we have demonstrated that mouse PRP extract increases phosphorylation levels of the Wnt coreceptor low-density lipoprotein receptor-related protein 5 (LRP5) and thereby activates angiogenic factor receptor Tie2 in endothelial cells (ECs) and accelerates EC sprouting and lung epithelial cell budding in vitro. PRP extract also increases phosphorylation levels of Tie2 in the mouse lungs and accelerates compensatory lung growth and recovery of exercise capacity after unilateral pneumonectomy in mice, whereas soluble Tie2 receptor or Lrp5 knockdown attenuates the effects of PRP extract. Because human PRP extract is generated from autologous peripheral blood and can be stored at -80°C, our findings may lead to the development of novel therapeutic interventions for various angiogenesis-related lung diseases and to the improvement of strategies for lung regeneration. PMID:26091161
Huang, C; Ouyang, J; Reese, T G; Wu, Y; El Fakhri, G; Ackerman, J L
2015-10-21
Due to the lack of signal from solid bone in normal MR sequences for the purpose of MR-based attenuation correction, investigators have proposed using the ultrashort echo time (UTE) pulse sequence, which yields signal from bone. However, the UTE-based segmentation approach might not fully capture the intra- and inter-subject bone density variation, which will inevitably lead to bias in reconstructed PET images. In this work, we investigated using the water- and fat-suppressed proton projection imaging (WASPI) sequence to obtain accurate and continuous attenuation for bones. This approach is capable of accounting for intra- and inter-subject bone attenuation variations. Using data acquired from a phantom, we have found that that attenuation correction based on the WASPI sequence is more accurate and precise when compared to either conventional MR attenuation correction or UTE-based segmentation approaches. PMID:26405761
NASA Astrophysics Data System (ADS)
Zhou, Ri-Feng; Wang, Jue; Chen, Wei-Min
2009-07-01
Due to X-ray attenuation being approximately proportional to material density, it is possible to measure the inner density through Industrial Computed Tomography (ICT) images accurately. In practice, however, a number of factors including the non-linear effects of beam hardening and diffuse scattered radiation complicate the quantitative measurement of density variations in materials. This paper is based on the linearization method of beam hardening correction, and uses polynomial fitting coefficient which is obtained by the curvature of iron polychromatic beam data to fit other materials. Through theoretical deduction, the paper proves that the density measure error is less than 2% if using pre-filters to make the spectrum of linear accelerator range mainly 0.3 MeV to 3 MeV. Experiment had been set up at an ICT system with a 9 MeV electron linear accelerator. The result is satisfactory. This technique makes the beam hardening correction easy and simple, and it is valuable for measuring the ICT density and making use of the CT images to recognize materials.
NASA Astrophysics Data System (ADS)
Blake, Nick P.; Metiu, Horia
1999-04-01
Electrons that are confined to zeolite cavities are modeled using a simplified pseudopotential scheme to represent the interaction of the electrons with both the sodalite framework and the Na+ ions. By comparing theory with recent experimental studies of G centers in Na-doped NaBr-SOD it is demonstrated that restricted forms of density functional theory, where two electrons are forced to pair in the same Kohn-Sham orbital, fail to correctly predict the true nature of the singlet, (spin unpolarized), G center. Electron confinement leads to generalized gradient corrections to the exchange of 0.74 eV and self-interaction corrections (SIC) of 0.7 eV over calculations performed in the local spin density approximation (LSDA). Only the self-interaction corrected generalized gradient approximation and the unrestricted Hartree-Fock approximation are in accord with experiment for the relative stability of the triplet (spin polarized) state. The unrestricted Hartree-Fock method is used to show that G-center absorptions will be blueshifted with respect to absorptions due to the isolated F centers. Constructing a Hubbard Hamiltonian we show that the exchange coupling ranges in values from 2.3 meV(UHF) to 3.6 meV(SIC-LSDA) corresponding to Neel temperatures that range from 27 to 41 K in agreement with experiment.
NASA Astrophysics Data System (ADS)
Tarpanov, D.; Toivanen, J.; Dobaczewski, J.; Carlsson, B. G.
2014-01-01
Background: Models based on using perturbative polarization corrections and mean-field blocking approximation give conflicting results for masses of odd nuclei. Purpose: We systematically investigate the polarization and mean-field models, implemented within self-consistent approaches that use identical interactions and model spaces, to find reasons for the conflicts between them. Methods: For density-dependent interactions and with pairing correlations included, we derive and study links between the mean-field and polarization results obtained for energies of odd nuclei. We also identify and discuss differences between the polarization-correction and full particle-vibration-coupling (PVC) models. Numerical calculations are performed for the mean-field ground-state properties of deformed odd nuclei and then compared to the polarization corrections determined using the approach that conserves spherical symmetry. Results: We have identified and numerically evaluated self-interaction (SI) energies that are at the origin of different results obtained within the mean-field and polarization-correction approaches. Conclusions: Mean-field energies of odd nuclei are polluted by the SI energies, and this makes them different from those obtained using polarization-correction methods. A comparison of both approaches allows for the identification and determination of the SI terms, which then can be calculated and removed from the mean-field results, giving the self-interaction-free energies. The simplest deformed mean-field approach that does not break parity symmetry is unable to reproduce full PVC effects.
Yamashiro, Tsuneo; Moriya, Hiroshi; Tsubakimoto, Maho; Matsuoka, Shin; Murayama, Sadayuki
2016-01-01
Purpose Four-dimensional dynamic-ventilation computed tomography (CT) imaging demonstrates continuous movement of the airways and lungs, which cannot be depicted with conventional CT. We aimed to investigate continuous changes in lung density and airway dimensions and to assess the correlation with spirometric values in smokers. Materials and methods This retrospective study was approved by the Institutional Review Board, and informed consent was waived. Twenty-one smokers including six patients with COPD underwent four-dimensional dynamic-ventilation CT during free breathing (160 mm in length). The mean lung density (MLD) of the scanned lung and luminal areas (Ai) of fixed points in the trachea and the right proximal bronchi (main bronchus, upper bronchus, bronchus intermedius, and lower bronchus) were continuously measured. Concordance between the time curve of the MLD and that of the airway Ai values was expressed by cross-correlation coefficients. The associations between these quantitative measurements and the forced expiratory volume in 1 second/forced vital capacity (FEV1/FVC) values were assessed by Spearman’s rank correlation analysis. Results On the time curve for the MLD, the Δ-MLD1.05 values between the peak inspiratory frame to the later third frame (1.05 seconds later) were strongly correlated with the FEV1/FVC (ρ=0.76, P<0.0001). The cross-correlation coefficients between the airway Ai and MLD values were significantly correlated with the FEV1/FVC (ρ=−0.56 to −0.66, P<0.01), except for the right upper bronchus. This suggested that the synchrony between the airway and lung movement was lost in patients with severe airflow limitation. Conclusion Respiratory changes in the MLD and synchrony between the airway Ai and the MLD measured with dynamic-ventilation CT were correlated with patient’s spirometric values. PMID:27110108
Kocman, Mikuláš; Jurečka, Petr; Dubecký, Matúš; Otyepka, Michal; Cho, Yeonchoo; Kim, Kwang S
2015-03-01
Hydrogen storage in carbonaceous materials and their derivatives is currently a widely investigated topic. The rational design of novel adsorptive materials is often attempted with the help of computational chemistry tools, in particular density functional theory (DFT). However, different exchange-correlation functionals provide a very wide range of hydrogen binding energies. The aim of this article is to offer high level QM reference data based on coupled-cluster singles and doubles calculations with perturbative triple excitations, CCSD(T), and a complete basis set limit estimate that can be used to assess the accuracy of various DFT-based predictions. For one complex, the CCSD(T) result is verified against diffusion quantum Monte Carlo calculations. Reference binding curves are calculated for two model compounds representing weak and strong hydrogen adsorption: coronene (-4.7 kJ mol(-1) per H2), and coronene modified with boron and lithium (-14.3 kJ mol(-1)). The reference data are compared to results obtained with widely used density functionals including pure DFT, M06, DFT-D3, PBE-TS, PBE + MBD, optB88-vdW, vdW-DF, vdW-DF2 and VV10. We find that whereas DFT-D3 shows excellent results for weak hydrogen adsorption on coronene, most of the less empirical density based dispersion functionals except VV10 overestimate this interaction. On the other hand, some of the less empirical density based dispersion functionals better describe stronger binding in the more polar coroB2Li22H2 complex which is one of realistic models for high-capacity hydrogen storage materials. Our results may serve as a guide for choosing suitable DFT methods for quickly evaluating hydrogen binding potential and as a reference for assessing the accuracy of the previously published DFT results. PMID:25655486
Song, Jong-Won; Hirao, Kimihiko
2015-10-14
Since the advent of hybrid functional in 1993, it has become a main quantum chemical tool for the calculation of energies and properties of molecular systems. Following the introduction of long-range corrected hybrid scheme for density functional theory a decade later, the applicability of the hybrid functional has been further amplified due to the resulting increased performance on orbital energy, excitation energy, non-linear optical property, barrier height, and so on. Nevertheless, the high cost associated with the evaluation of Hartree-Fock (HF) exchange integrals remains a bottleneck for the broader and more active applications of hybrid functionals to large molecular and periodic systems. Here, we propose a very simple yet efficient method for the computation of long-range corrected hybrid scheme. It uses a modified two-Gaussian attenuating operator instead of the error function for the long-range HF exchange integral. As a result, the two-Gaussian HF operator, which mimics the shape of the error function operator, reduces computational time dramatically (e.g., about 14 times acceleration in C diamond calculation using periodic boundary condition) and enables lower scaling with system size, while maintaining the improved features of the long-range corrected density functional theory.
Song, Jong-Won; Hirao, Kimihiko
2015-10-14
Since the advent of hybrid functional in 1993, it has become a main quantum chemical tool for the calculation of energies and properties of molecular systems. Following the introduction of long-range corrected hybrid scheme for density functional theory a decade later, the applicability of the hybrid functional has been further amplified due to the resulting increased performance on orbital energy, excitation energy, non-linear optical property, barrier height, and so on. Nevertheless, the high cost associated with the evaluation of Hartree-Fock (HF) exchange integrals remains a bottleneck for the broader and more active applications of hybrid functionals to large molecular and periodic systems. Here, we propose a very simple yet efficient method for the computation of long-range corrected hybrid scheme. It uses a modified two-Gaussian attenuating operator instead of the error function for the long-range HF exchange integral. As a result, the two-Gaussian HF operator, which mimics the shape of the error function operator, reduces computational time dramatically (e.g., about 14 times acceleration in C diamond calculation using periodic boundary condition) and enables lower scaling with system size, while maintaining the improved features of the long-range corrected density functional theory. PMID:26472368
Cheng, Nai-Ming; Fang, Yu-Hua Dean; Tsan, Din-Li
2016-01-01
Purpose We compared attenuation correction of PET images with helical CT (PET/HCT) and respiration-averaged CT (PET/ACT) in patients with non-small-cell lung cancer (NSCLC) with the goal of investigating the impact of respiration-averaged CT on 18F FDG PET texture parameters. Materials and Methods A total of 56 patients were enrolled. Tumors were segmented on pretreatment PET images using the adaptive threshold. Twelve different texture parameters were computed: standard uptake value (SUV) entropy, uniformity, entropy, dissimilarity, homogeneity, coarseness, busyness, contrast, complexity, grey-level nonuniformity, zone-size nonuniformity, and high grey-level large zone emphasis. Comparisons of PET/HCT and PET/ACT were performed using Wilcoxon signed-rank tests, intraclass correlation coefficients, and Bland-Altman analysis. Receiver operating characteristic (ROC) curves as well as univariate and multivariate Cox regression analyses were used to identify the parameters significantly associated with disease-specific survival (DSS). A fixed threshold at 45% of the maximum SUV (T45) was used for validation. Results SUV maximum and total lesion glycolysis (TLG) were significantly higher in PET/ACT. However, texture parameters obtained with PET/ACT and PET/HCT showed a high degree of agreement. The lowest levels of variation between the two modalities were observed for SUV entropy (9.7%) and entropy (9.8%). SUV entropy, entropy, and coarseness from both PET/ACT and PET/HCT were significantly associated with DSS. Validation analyses using T45 confirmed the usefulness of SUV entropy and entropy in both PET/HCT and PET/ACT for the prediction of DSS, but only coarseness from PET/ACT achieved the statistical significance threshold. Conclusions Our results indicate that 1) texture parameters from PET/ACT are clinically useful in the prediction of survival in NSCLC patients and 2) SUV entropy and entropy are robust to attenuation correction methods. PMID:26930211
Cady, John W.
1977-01-01
A computer program is presented which performs, for one or more bodies, along a profile perpendicular to strike, both forward calculations for the magnetic and gravity anomaly fields and independent gravity and magnetic inverse calculations for density and susceptibility or remanent magnetization.
NASA Technical Reports Server (NTRS)
Mertens, C. J.; Xu, X.; Fernandez, J. R.; Bilitza, D.; Russell, J. M., III; Mlynczak, M. G.
2009-01-01
Auroral infrared emission observed from the TIMED/SABER broadband 4.3 micron channel is used to develop an empirical geomagnetic storm correction to the International Reference Ionosphere (IRI) E-region electron densities. The observation-based proxy used to develop the storm model is SABER-derived NO+(v) 4.3 micron volume emission rates (VER). A correction factor is defined as the ratio of storm-time NO+(v) 4.3 micron VER to a quiet-time climatological averaged NO+(v) 4.3 micron VER, which is linearly fit to available geomagnetic activity indices. The initial version of the E-region storm model, called STORM-E, is most applicable within the auroral oval region. The STORM-E predictions of E-region electron densities are compared to incoherent scatter radar electron density measurements during the Halloween 2003 storm events. Future STORM-E updates will extend the model outside the auroral oval.
NASA Astrophysics Data System (ADS)
Güler-Kılıç, Sümeyra; Kılıç, ćetin
2015-06-01
Semilocal and dispersion-corrected density-functional calculations have been performed to study the crystal structure, equation of state, and electronic structure of metal tellurohalides with chemical formula MeTeI where Me=Bi, Au, or Pd. A comparative investigation of the results of these calculations is conducted which reveals the role of van der Waals attraction. It is shown that the prediction of crystal structure of metal tellurohalides is systematically improved thanks to the inclusion of van der Waals dispersion. It is found for BiTeI and AuTeI that the energy versus volume curve is anomalously flat in the vicinity of equilibrium volume and the calculated equation of state has an excessively steep slope in the low-pressure region; these are also fixed in the dispersion-corrected calculations. Analysis based on the computation of the volume and axial compressibilities shows that predicting the anisotropy of BiTeI via the semilocal calculations yields an unrealistic result, whereas the results of dispersion-corrected calculations agree with the experimental compressibility data. Our calculations render that BiTeI (AuTeI) is a narrow band gap semiconductor with Rashba-type spin splitting at the band edges (with an indirect band gap) while PdTeI is a metal with relatively low density of states at the Fermi level. The band gaps of BiTeI and AuTeI obtained via semilocal (dispersion-corrected) calculations are found to be greater (smaller) than the respective experimental values, which is against (in line with) the expected trend. Similarly, the Rashba parameters of BiTeI are bracketed by the respective values obtained via semilocal and dispersion-corrected calculations, e.g., a larger value for the Rashba parameter αR is obtained in association with the reduction of the band gap caused by modification of the crystal structure owing to van der Waals attraction. Excellent agreement with the experimental Rashba parameters is obtained via interpolation of the
Gudmundsdóttir, Hildur; Zhang, Yao; Weber, Peter M; Jónsson, Hannes
2014-12-21
Theoretical calculations of Rydberg excited states of molecular clusters consisting of N,N-dimethylisopropylamine molecules using a Perdew-Zunger self-interaction corrected energy functional are presented and compared with results of resonant multiphoton ionization measurements. The binding energy of the Rydberg electron in the monomer is calculated to be 2.79 eV and 2.27 eV in the 3s and 3p state, respectively, which compares well with measured values of 2.88 eV and 2.21 eV. Three different stable configurations of the dimer in the ground state were found using an energy functional that includes van der Waals interaction. The lowest ground state energy conformation has the two N-atoms widely separated, by 6.2 Å, while the Rydberg state energy is lowest for a configuration where the N-atoms of the two molecules come close together, separated by 3.7 Å. This conformational change is found to lower the Rydberg electron binding energy by 0.2 eV. The self-interaction corrected functional gives a highly localized hole on one of the two molecules, unlike results obtained using the PBE functional or the hybrid B3LYP functional which give a delocalized hole. For the trimer, the self-interaction corrected calculation gives a Rydberg electron binding energy lowered further by 0.13 eV as compared with the dimer. The calculated results compare well with trends observed in experimental measurements. The reduction of the Rydberg electron binding energy with cluster size can be ascribed to an effective delocalization of the positive charge of the hole by the induced and permanent dipole moments of the neighboring molecules. A further decrease observed to occur on a time scale of tens of ps can be ascribed to a structural rearrangement of the clusters in the Rydberg state where molecules rotate to orient their dipoles in response to the formation of the localized hole. PMID:25527936
Gudmundsdóttir, Hildur; Zhang, Yao; Weber, Peter M.; Jónsson, Hannes
2014-12-21
Theoretical calculations of Rydberg excited states of molecular clusters consisting of N,N-dimethylisopropylamine molecules using a Perdew-Zunger self-interaction corrected energy functional are presented and compared with results of resonant multiphoton ionization measurements. The binding energy of the Rydberg electron in the monomer is calculated to be 2.79 eV and 2.27 eV in the 3s and 3p state, respectively, which compares well with measured values of 2.88 eV and 2.21 eV. Three different stable configurations of the dimer in the ground state were found using an energy functional that includes van der Waals interaction. The lowest ground state energy conformation has the two N-atoms widely separated, by 6.2 Å, while the Rydberg state energy is lowest for a configuration where the N-atoms of the two molecules come close together, separated by 3.7 Å. This conformational change is found to lower the Rydberg electron binding energy by 0.2 eV. The self-interaction corrected functional gives a highly localized hole on one of the two molecules, unlike results obtained using the PBE functional or the hybrid B3LYP functional which give a delocalized hole. For the trimer, the self-interaction corrected calculation gives a Rydberg electron binding energy lowered further by 0.13 eV as compared with the dimer. The calculated results compare well with trends observed in experimental measurements. The reduction of the Rydberg electron binding energy with cluster size can be ascribed to an effective delocalization of the positive charge of the hole by the induced and permanent dipole moments of the neighboring molecules. A further decrease observed to occur on a time scale of tens of ps can be ascribed to a structural rearrangement of the clusters in the Rydberg state where molecules rotate to orient their dipoles in response to the formation of the localized hole.
Lao, Ka Un; Herbert, John M.
2014-01-28
The performance of second-order symmetry-adapted perturbation theory (SAPT) calculations using Kohn-Sham (KS) orbitals is evaluated against benchmark results for intermolecular interactions. Unlike previous studies of this “SAPT(KS)” methodology, the present study uses non-empirically tuned long-range corrected (LRC) functionals for the monomers. The proper v{sub xc} (r)→0 asymptotic limit is achieved by tuning the range separation parameter in order to satisfy the condition that the highest occupied KS energy level equals minus the molecule's ionization energy, for each monomer unit. Tests for He{sub 2}, Ne{sub 2}, and the S22 and S66 data sets reveal that this condition is important for accurate prediction of the non-dispersion components of the energy, although errors in SAPT(KS) dispersion energies remain unacceptably large. In conjunction with an empirical dispersion potential, however, the SAPT(KS) method affords good results for S22 and S66, and also accurately predicts the whole potential energy curve for the sandwich isomer of the benzene dimer. Tuned LRC functionals represent an attractive alternative to other asymptotic corrections that have been employed in density-functional-based SAPT calculations, and we recommend the use of tuned LRC functionals in both coupled-perturbed SAPT(DFT) calculations and dispersion-corrected SAPT(KS) calculations.
NASA Astrophysics Data System (ADS)
Bartók, Albert P.; Gillan, Michael J.; Manby, Frederick R.; Csányi, Gábor
2013-08-01
We show how machine learning techniques based on Bayesian inference can be used to enhance the computer simulation of molecular materials, focusing here on water. We train our machine-learning algorithm using accurate, correlated quantum chemistry, and predict energies and forces in molecular aggregates ranging from clusters to solid and liquid phases. The widely used electronic-structure methods based on density functional theory (DFT) by themselves give poor accuracy for molecular materials like water, and we show how our techniques can be used to generate systematically improvable one- and two-body corrections to DFT with modest extra resources. The resulting corrected DFT scheme is considerably more accurate than uncorrected DFT for the relative energies of small water clusters and different ice structures and significantly improves the description of the structure and dynamics of liquid water. However, our results for ice structures and the liquid indicate that beyond-two-body DFT errors cannot be ignored, and we suggest how our machine-learning methods can be further developed to correct these errors.
Silvestrelli, Pier Luigi; Ambrosetti, Alberto
2014-03-28
The Density Functional Theory (DFT)/van der Waals-Quantum Harmonic Oscillator-Wannier function (vdW-QHO-WF) method, recently developed to include the vdW interactions in approximated DFT by combining the quantum harmonic oscillator model with the maximally localized Wannier function technique, is applied to the cases of atoms and small molecules (X=Ar, CO, H{sub 2}, H{sub 2}O) weakly interacting with benzene and with the ideal planar graphene surface. Comparison is also presented with the results obtained by other DFT vdW-corrected schemes, including PBE+D, vdW-DF, vdW-DF2, rVV10, and by the simpler Local Density Approximation (LDA) and semilocal generalized gradient approximation approaches. While for the X-benzene systems all the considered vdW-corrected schemes perform reasonably well, it turns out that an accurate description of the X-graphene interaction requires a proper treatment of many-body contributions and of short-range screening effects, as demonstrated by adopting an improved version of the DFT/vdW-QHO-WF method. We also comment on the widespread attitude of relying on LDA to get a rough description of weakly interacting systems.
Sriwayu, Wa Ode; Haryanto, Freddy; Khotimah, Siti Nurul; Latief, Fourier Dzar Eljabbar
2015-04-16
We have designed and fabricated phantom mimicking breast cancer composition known as a region that has low contrast density. The used compositions are a microcalcifications, fatty tissues and tumor mass by using Al{sub 2}O{sub 3}, C{sub 27}H{sub 46}O, and hard nylon materials. Besides, phantom also has a part to calculate low cost criteria /CNR (Contrast to Noise Ratio). Uniformity will be measured at water distillation medium located in a part of phantom scale contrast. Phantom will be imaged by using micro ct-sky scan 1173 high energy type, and then also can be quantified CT number to examine SkyScan 1173 performance in imaging low contrast density materials. Evaluation of CT number is done at technique configuration parameter using voltage of 30 kV, exposure 0.160 mAs, and camera resolution 560x560 pixel, the effect of image quality to reconstruction process is evaluated by varying image processing parameters in the form of beam hardening corrections with amount of 25%, 66% and100% with each smoothing level S10,S2 and S7. To obtain the better high quality image, the adjustment of beam hardening correction should be 66% and smoothing level reach maximal value at level 10.
NASA Astrophysics Data System (ADS)
Sriwayu, Wa Ode; Haryanto, Freddy; Khotimah, Siti Nurul; Latief, Fourier Dzar Eljabbar
2015-04-01
We have designed and fabricated phantom mimicking breast cancer composition known as a region that has low contrast density. The used compositions are a microcalcifications, fatty tissues and tumor mass by using Al2O3, C27H46O, and hard nylon materials. Besides, phantom also has a part to calculate low cost criteria /CNR (Contrast to Noise Ratio). Uniformity will be measured at water distillation medium located in a part of phantom scale contrast. Phantom will be imaged by using micro ct-sky scan 1173 high energy type, and then also can be quantified CT number to examine SkyScan 1173 performance in imaging low contrast density materials. Evaluation of CT number is done at technique configuration parameter using voltage of 30 kV, exposure 0.160 mAs, and camera resolution 560x560 pixel, the effect of image quality to reconstruction process is evaluated by varying image processing parameters in the form of beam hardening corrections with amount of 25%, 66% and100% with each smoothing level S10,S2 and S7. To obtain the better high quality image, the adjustment of beam hardening correction should be 66% and smoothing level reach maximal value at level 10.
Silvestrelli, Pier Luigi; Ambrosetti, Alberto
2014-03-28
The Density Functional Theory (DFT)/van der Waals-Quantum Harmonic Oscillator-Wannier function (vdW-QHO-WF) method, recently developed to include the vdW interactions in approximated DFT by combining the quantum harmonic oscillator model with the maximally localized Wannier function technique, is applied to the cases of atoms and small molecules (X=Ar, CO, H2, H2O) weakly interacting with benzene and with the ideal planar graphene surface. Comparison is also presented with the results obtained by other DFT vdW-corrected schemes, including PBE+D, vdW-DF, vdW-DF2, rVV10, and by the simpler Local Density Approximation (LDA) and semilocal generalized gradient approximation approaches. While for the X-benzene systems all the considered vdW-corrected schemes perform reasonably well, it turns out that an accurate description of the X-graphene interaction requires a proper treatment of many-body contributions and of short-range screening effects, as demonstrated by adopting an improved version of the DFT/vdW-QHO-WF method. We also comment on the widespread attitude of relying on LDA to get a rough description of weakly interacting systems. PMID:24697424
Bučko, Tomáš; Hafner, Jürgen; Lebègue, Sébastien; Ángyán, János G
2012-04-28
Periodic dispersion corrected DFT calculations have been performed to study the spin-crossover transition of Fe(phen)(2)(NCS)(2) in the molecular and in the crystalline state. We show that London dispersion interactions play a crucial role in the cohesion of the crystals. Based on calculations of vibrational eigenstates of the isolated molecule and of the crystalline phase in both the low- and high-spin states, the transition entropies and enthalpies have been calculated. We demonstrate that, due to the stabilization of the low-spin state by intermolecular dispersion forces, the transition enthalpy at the transition temperature is larger for the crystalline phase in comparison with an isolated molecule. The effective coordination number of the nitrogen atoms of the ligands around the iron atom has been identified as the order parameter driving the quasi-reversible low-spin to high-spin transition in the crystal. Finally, using constrained geometry relaxations at fixed values of the coordination number, we computed the energy barrier of the LS to HS transition and found it to be in a reasonable agreement with the experimental value. PMID:22415338
NASA Astrophysics Data System (ADS)
Gritsenko, Oleg; Baerends, Evert Jan
2004-07-01
Time-dependent density functional theory (TDDFT) calculations of charge-transfer excitation energies ωCT are significantly in error when the adiabatic local density approximation (ALDA) is employed for the exchange-correlation kernel fxc. We relate the error to the physical meaning of the orbital energy of the Kohn-Sham lowest unoccupied molecular orbital (LUMO). The LUMO orbital energy in Kohn-Sham DFT—in contrast to the Hartree-Fock model—approximates an excited electron, which is correct for excitations in compact molecules. In CT transitions the energy of the LUMO of the acceptor molecule should instead describe an added electron, i.e., approximate the electron affinity. To obtain a contribution that compensates for the difference, a specific divergence of fxc is required in rigorous TDDFT, and a suitable asymptotically correct form of the kernel fxcasymp is proposed. The importance of the asymptotic correction of fxc is demonstrated with the calculation of ωCT(R) for the prototype diatomic system HeBe at various separations R(He-Be). The TDDFT-ALDA curve ωCT(R) roughly resembles the benchmark ab initio curve ωCTCISD(R) of a configuration interaction calculation with single and double excitations in the region R=1-1.5 Å, where a sizable He-Be interaction exists, but exhibits the wrong behavior ωCT(R)≪ωCTCISD(R) at large R. The TDDFT curve obtained with fxcasymp however approaches ωCTCISD(R) closely in the region R=3-10 Å. Then, the adequate rigorous TDDFT approach should interpolate between the LDA/GGA ALDA xc kernel for excitations in compact systems and fxcasymp for weakly interacting fragments and suitable interpolation expressions are considered.
NASA Astrophysics Data System (ADS)
Liu, Bin; Cheng, Lei; Curtiss, Larry; Greeley, Jeffrey
2014-04-01
The hydrogenation of furfural to furfuryl alcohol on Pd(111), Cu(111) and Pt(111) is studied with both standard Density Functional Theory (DFT)-GGA functionals and with van der Waals-corrected density functionals. VdW-DF functionals, including optPBE, optB88, optB86b, and Grimme's method, are used to optimize the adsorption configurations of furfural, furfuryl alcohol, and related intermediates resulting from hydrogenation of furfural, and the results are compared to corresponding values determined with GGA functionals, including PW91 and PBE. On Pd(111) and Pt(111), the adsorption geometries of the intermediates are not noticeably different between the two classes of functionals, while on Cu(111), modest changes are seen in both the perpendicular distance and the orientation of the aromatic ring with respect to the planar surface. In general, the binding energies increase substantially in magnitude as a result of van der Waals contributions on all metals. In contrast, however, dispersion effects on the kinetics of hydrogenation are relatively small. It is found that activation barriers are not significantly affected by the inclusion of dispersion effects, and a Brønsted-Evans-Polanyi relationship developed solely from PW91 calculations on Pd(111) is capable of describing corresponding results on Cu(111) and Pt(111), even when the dispersion effects are included. Finally, the reaction energies and barriers derived from the dispersion-corrected and pure GGA calculations are used to plot simple potential energy profiles for furfural hydrogenation to furfuryl alcohol on the three considered metals, and an approximately constant downshift of the energetics due to the dispersion corrections is observed.
Liu, Bin; Cheng, Lei; Curtiss, Larry A.; Greeley, Jeffrey P.
2014-04-01
The hydrogenation of furfural to furfuryl alcohol on Pd(111), Cu(111) and Pt(111) is studied with both standard Density Functional Theory (DFT)-GGA functionals and with van der Waals-corrected density functionals. VdWDF functionals, including optPBE, optB88, optB86b, and Grimme's method, are used to optimize the adsorption configurations of furfural, furfuryl alcohol, and related intermediates resulting from hydrogenation of furfural, and the results are compared to corresponding values determined with GGA functionals, including PW91 and PBE. On Pd(111) and Pt(111), the adsorption geometries of the intermediates are not noticeably different between the two classes of functionals, while on Cu(111), modest changes are seen in both the erpendicular distance and the orientation of the aromatic ringwith respect to the planar surface. In general, the binding energies increase substantially in magnitude as a result of van derWaals contributions on all metals. In contrast, however, dispersion effects on the kinetics of hydrogenation are relatively small. It is found that activation barriers are not significantly affected by the inclusion of dispersion effects, and a Brønsted–Evans–Polanyi relationship developed solely fromPW91 calculations on Pd(111) is capable of describing corresponding results on Cu(111) and Pt(111), even when the dispersion effects are included. Finally, the reaction energies and barriers derived from the dispersion-corrected and pure GGA calculations are used to plot simple potential energy profiles for furfural hydrogenation to furfuryl alcohol on the three considered metals, and an approximately constant downshift of the energetics due to the dispersion corrections is observed.
NASA Astrophysics Data System (ADS)
Ding, Hong; Ray, Keith G.; Ozolins, Vidvuds; Asta, Mark
2012-01-01
Structural and vibrational properties of α-MoO3 are studied employing two recently proposed methodologies for incorporating van der Waals (vdW) contributions in density functional theory (DFT) based calculations. The DFT-D2 [S. Grimme, J. Comput. Chem.JCCHDD0192-865110.1002/jcc.20495 27, 1787 (2006)] and optB88 vdW-DFT [J. Klimeš , J. Phys.: Condens. MatterPRBMDO0953-898410.1088/0953-8984/22/2/022201 22, 022201 (2010)] methods are shown to give rise to increased accuracy in predicted lattice parameters, relative to conventional DFT methods. Calculated vibrational frequencies agree with measurements to within 5% and 10% for modes involving bonded and nonbonded interactions in this compound, respectively.
Mukherjee, Sanchita; Kailasam, Senthilkumar; Bansal, Manju; Bhattacharyya, Dhananjay
2014-01-01
Double helical structures of DNA and RNA are mostly determined by base pair stacking interactions, which give them the base sequence-directed features, such as small roll values for the purine-pyrimidine steps. Earlier attempts to characterize stacking interactions were mostly restricted to calculations on fiber diffraction geometries or optimized structure using ab initio calculations lacking variation in geometry to comment on rather unusual large roll values observed in AU/AU base pair step in crystal structures of RNA double helices. We have generated stacking energy hyperspace by modeling geometries with variations along the important degrees of freedom, roll, and slide, which were chosen via statistical analysis as maximally sequence dependent. Corresponding energy contours were constructed by several quantum chemical methods including dispersion corrections. This analysis established the most suitable methods for stacked base pair systems despite the limitation imparted by number of atom in a base pair step to employ very high level of theory. All the methods predict negative roll value and near-zero slide to be most favorable for the purine-pyrimidine steps, in agreement with Calladine's steric clash based rule. Successive base pairs in RNA are always linked by sugar-phosphate backbone with C3'-endo sugars and this demands C1'-C1' distance of about 5.4 Å along the chains. Consideration of an energy penalty term for deviation of C1'-C1' distance from the mean value, to the recent DFT-D functionals, specifically ωB97X-D appears to predict reliable energy contour for AU/AU step. Such distance-based penalty improves energy contours for the other purine-pyrimidine sequences also. © 2013 Wiley Periodicals, Inc. Biopolymers 101: 107-120, 2014. PMID:23722519
Wang, Jianwei; Zhang, Yong; Wang, Lin-Wang
2015-07-31
We propose a systematic approach that can empirically correct three major errors typically found in a density functional theory (DFT) calculation within the local density approximation (LDA) simultaneously for a set of common cation binary semiconductors, such as III-V compounds, (Ga or In)X with X = N,P,As,Sb, and II-VI compounds, (Zn or Cd)X, with X = O,S,Se,Te. By correcting (1) the binary band gaps at high-symmetry points , L, X, (2) the separation of p-and d-orbital-derived valence bands, and (3) conduction band effective masses to experimental values and doing so simultaneously for common cation binaries, the resulting DFT-LDA-based quasi-first-principles methodmore » can be used to predict the electronic structure of complex materials involving multiple binaries with comparable accuracy but much less computational cost than a GW level theory. This approach provides an efficient way to evaluate the electronic structures and other material properties of complex systems, much needed for material discovery and design.« less
Kruse, Holger; Grimme, Stefan
2012-04-21
A semi-empirical counterpoise-type correction for basis set superposition error (BSSE) in molecular systems is presented. An atom pair-wise potential corrects for the inter- and intra-molecular BSSE in supermolecular Hartree-Fock (HF) or density functional theory (DFT) calculations. This geometrical counterpoise (gCP) denoted scheme depends only on the molecular geometry, i.e., no input from the electronic wave-function is required and hence is applicable to molecules with ten thousands of atoms. The four necessary parameters have been determined by a fit to standard Boys and Bernadi counterpoise corrections for Hobza's S66×8 set of non-covalently bound complexes (528 data points). The method's target are small basis sets (e.g., minimal, split-valence, 6-31G*), but reliable results are also obtained for larger triple-ζ sets. The intermolecular BSSE is calculated by gCP within a typical error of 10%-30% that proves sufficient in many practical applications. The approach is suggested as a quantitative correction in production work and can also be routinely applied to estimate the magnitude of the BSSE beforehand. The applicability for biomolecules as the primary target is tested for the crambin protein, where gCP removes intramolecular BSSE effectively and yields conformational energies comparable to def2-TZVP basis results. Good mutual agreement is also found with Jensen's ACP(4) scheme, estimating the intramolecular BSSE in the phenylalanine-glycine-phenylalanine tripeptide, for which also a relaxed rotational energy profile is presented. A variety of minimal and double-ζ basis sets combined with gCP and the dispersion corrections DFT-D3 and DFT-NL are successfully benchmarked on the S22 and S66 sets of non-covalent interactions. Outstanding performance with a mean absolute deviation (MAD) of 0.51 kcal/mol (0.38 kcal/mol after D3-refit) is obtained at the gCP-corrected HF-D3/(minimal basis) level for the S66 benchmark. The gCP-corrected B3LYP-D3/6-31G* model
NASA Astrophysics Data System (ADS)
Tong, Xiao-Min; Chu, Shih-I.
1998-02-01
We present a self-interaction-free relativistic density-functional theory (DFT). The theory is based on the extension of our recent nonrelativistic DFT treatment with optimized effective potential (OEP) and self-interaction correction (SIC) [Phys. Rev. A 55, 3406 (1997)] to the relativistic domain. Such a relativistic OEP-SIC procedure yields an orbital-independent single-particle local potential with proper long-range Coulombic (-1/r) behavior. The method is applied to the ground-state energy calculations for atoms with Z=2-106. A comparison with the corresponding nonrelativistic OEP-SIC calculations and other relativistic calculations is made. It is shown that the ionization potentials (obtained from the highest occupied orbital energies) and individual orbital binding energies determined by the present relativistic OEP-SIC method agree well with the experimental data across the Periodic Table.
NASA Astrophysics Data System (ADS)
Dinh, Phuong Mai; Gao, Cong Zhang; Klüpfel, Peter; Reinhard, Paul-Gerhard; Suraud, Eric; Vincendon, Marc; Wang, Jing; Zhang, Feng Shou
2014-08-01
We present a detailed analysis of ground state and optical response properties of small metal water complexes. Such complexes represent prototypical systems for analysing chromophore effects in relation to irradiation in a biological environment. The mixing of a metal atom with organic ones leads to the coexistence of covalent and metallic bondings which requires an elaborate treatment of the self interaction correction (SIC) within density functional theory (DFT). This is a particularly key issue in the context of time dependent DFT which represents the natural tool of investigation of irradiation scenarios in such systems. We show that these complexes require a highly elaborate treatment of the SIC which can be attributed to the mixing of bonding types.
NASA Technical Reports Server (NTRS)
Mertens, Christoper J.; Winick, Jeremy R.; Russell, James M., III; Mlynczak, Martin G.; Evans, David S.; Bilitza, Dieter; Xu, Xiaojing
2007-01-01
The response of the ionospheric E-region to solar-geomagnetic storms can be characterized using observations of infrared 4.3 micrometers emission. In particular, we utilize nighttime TIMED/SABER measurements of broadband 4.3 micrometers limb emission and derive a new data product, the NO+(v) volume emission rate, which is our primary observation-based quantity for developing an empirical storm-time correction the IRI E-region electron density. In this paper we describe our E-region proxy and outline our strategy for developing the empirical storm model. In our initial studies, we analyzed a six day storm period during the Halloween 2003 event. The results of this analysis are promising and suggest that the ap-index is a viable candidate to use as a magnetic driver for our model.
NASA Astrophysics Data System (ADS)
Verma, A. K.; Fienga, A.; Laskar, J.; Issautier, K.; Manche, H.; Gastineau, M.
2013-02-01
The Mars Global Surveyor (MGS), Mars Express (MEX), and Venus Express (VEX) experienced several superior solar conjunctions. These conjunctions cause severe degradations of radio signals when the line of sight between the Earth and the spacecraft passes near to the solar corona region. The primary objective of this work is to deduce a solar corona model from the spacecraft navigation data acquired at the time of solar conjunctions and to estimate its average electron density. The corrected or improved data are then used to fit the dynamical modeling of the planet motions, called planetary ephemerides. We analyzed the radio science raw data of the MGS spacecraft using the orbit determination software GINS. The range bias, obtained from GINS and provided by ESA for MEX and VEX, are then used to derive the electron density profile. These profiles are obtained for different intervals of solar distances: from 12 R⊙ to 215 R⊙ for MGS, 6 R⊙ to 152 R⊙ for MEX, and from 12 R⊙ to 154 R⊙ for VEX. They are acquired for each spacecraft individually, for ingress and egress phases separately and both phases together, for different types of solar winds (fast, slow), and for solar activity phases (minimum, maximum). We compared our results with the previous estimations that were based on in situ measurements, and on solar type III radio and radio science studies made at different phases of solar activity and at different solar wind states. Our results are consistent with estimations obtained by these different methods. Moreover, fitting the planetary ephemerides including complementary data that were corrected for the solar corona perturbations, noticeably improves the extrapolation capability of the planetary ephemerides and the estimation of the asteroids masses. Tables 5, 6 and Appendix A are available in electronic form at http://www.aanda.org
NASA Astrophysics Data System (ADS)
Mardirossian, Narbe; Head-Gordon, Martin
2014-05-01
The limit of accuracy for semi-empirical generalized gradient approximation (GGA) density functionals is explored by parameterizing a variety of local, global hybrid, and range-separated hybrid functionals. The training methodology employed differs from conventional approaches in 2 main ways: (1) Instead of uniformly truncating the exchange, same-spin correlation, and opposite-spin correlation functional inhomogeneity correction factors, all possible fits up to fourth order are considered, and (2) Instead of selecting the optimal functionals based solely on their training set performance, the fits are validated on an independent test set and ranked based on their overall performance on the training and test sets. The 3 different methods of accounting for exchange are trained both with and without dispersion corrections (DFT-D2 and VV10), resulting in a total of 491 508 candidate functionals. For each of the 9 functional classes considered, the results illustrate the trade-off between improved training set performance and diminished transferability. Since all 491 508 functionals are uniformly trained and tested, this methodology allows the relative strengths of each type of functional to be consistently compared and contrasted. The range-separated hybrid GGA functional paired with the VV10 nonlocal correlation functional emerges as the most accurate form for the present training and test sets, which span thermochemical energy differences, reaction barriers, and intermolecular interactions involving lighter main group elements.
Mardirossian, Narbe; Head-Gordon, Martin
2014-05-14
The limit of accuracy for semi-empirical generalized gradient approximation (GGA) density functionals is explored by parameterizing a variety of local, global hybrid, and range-separated hybrid functionals. The training methodology employed differs from conventional approaches in 2 main ways: (1) Instead of uniformly truncating the exchange, same-spin correlation, and opposite-spin correlation functional inhomogeneity correction factors, all possible fits up to fourth order are considered, and (2) Instead of selecting the optimal functionals based solely on their training set performance, the fits are validated on an independent test set and ranked based on their overall performance on the training and test sets. The 3 different methods of accounting for exchange are trained both with and without dispersion corrections (DFT-D2 and VV10), resulting in a total of 491 508 candidate functionals. For each of the 9 functional classes considered, the results illustrate the trade-off between improved training set performance and diminished transferability. Since all 491 508 functionals are uniformly trained and tested, this methodology allows the relative strengths of each type of functional to be consistently compared and contrasted. The range-separated hybrid GGA functional paired with the VV10 nonlocal correlation functional emerges as the most accurate form for the present training and test sets, which span thermochemical energy differences, reaction barriers, and intermolecular interactions involving lighter main group elements.
NASA Astrophysics Data System (ADS)
Verma, Prakash; Bartlett, Rodney J.
2016-07-01
Core excitation energies are computed with time-dependent density functional theory (TD-DFT) using the ionization energy corrected exchange and correlation potential QTP(0,0). QTP(0,0) provides C, N, and O K-edge spectra to about an electron volt. A mean absolute error (MAE) of 0.77 and a maximum error of 2.6 eV is observed for QTP(0,0) for many small molecules. TD-DFT based on QTP (0,0) is then used to describe the core-excitation spectra of the 22 amino acids. TD-DFT with conventional functionals greatly underestimates core excitation energies, largely due to the significant error in the Kohn-Sham occupied eigenvalues. To the contrary, the ionization energy corrected potential, QTP(0,0), provides excellent approximations (MAE of 0.53 eV) for core ionization energies as eigenvalues of the Kohn-Sham equations. As a consequence, core excitation energies are accurately described with QTP(0,0), as are the core ionization energies important in X-ray photoionization spectra or electron spectroscopy for chemical analysis.
Verma, Prakash; Bartlett, Rodney J
2016-07-21
Core excitation energies are computed with time-dependent density functional theory (TD-DFT) using the ionization energy corrected exchange and correlation potential QTP(0,0). QTP(0,0) provides C, N, and O K-edge spectra to about an electron volt. A mean absolute error (MAE) of 0.77 and a maximum error of 2.6 eV is observed for QTP(0,0) for many small molecules. TD-DFT based on QTP (0,0) is then used to describe the core-excitation spectra of the 22 amino acids. TD-DFT with conventional functionals greatly underestimates core excitation energies, largely due to the significant error in the Kohn-Sham occupied eigenvalues. To the contrary, the ionization energy corrected potential, QTP(0,0), provides excellent approximations (MAE of 0.53 eV) for core ionization energies as eigenvalues of the Kohn-Sham equations. As a consequence, core excitation energies are accurately described with QTP(0,0), as are the core ionization energies important in X-ray photoionization spectra or electron spectroscopy for chemical analysis. PMID:27448875
NASA Astrophysics Data System (ADS)
Ambrosetti, Alberto; Silvestrelli, Pier Luigi
2016-07-01
The cohesive energy, equilibrium lattice constant, and bulk modulus of Au, Ag, and Cu noble metals are computed by different van der Waals (vdW)-corrected density functional theory (DFT) methods, including vdW-DF, vdW-DF2, vdW-DF-cx, rVV10, and PBE-D. Two specifically designed methods are also developed in order to effectively include dynamical screening effects: the DFT/vdW-WF2p method, based on the generation of maximally localized Wannier functions, and the RPAp scheme (in two variants), based on a single-oscillator model of the localized electron response. Comparison with results obtained without explicit inclusion of van der Waals effects, such as with the local density approximation (LDA), PBE, PBEsol, or the hybrid PBE0 functional, elucidates the importance of a suitable description of screened van der Waals interactions even in the case of strong metal bonding. Many-body effects are also quantitatively evaluated within the RPAp approach.
NASA Astrophysics Data System (ADS)
Nguyen, Chuong V.; Hieu, Nguyen N.; Nguyen, Duong T.
2015-11-01
Strain-dependent structural and electronic properties of MoS2 materials are investigated using first principles calculations. The structural and electronic band structures of the MoS2 with relaxed unit cells are optimized and calculated by the dispersion-corrected density functional theory (DFT-D2). Calculations within the local density approximation (LDA) and GGA using PAW potentials were also performed for specific cases for the purpose of comparison. The effect of strain on the band gap and the dependence of formation energy on strain of MoS2 are also studied and discussed using the DFT-D2 method. In bulk MoS2, the orbitals shift towards the higher/lower energy area when strain is applied along the z/ x direction, respectively. The energy splitting of Mo4 d states is in the range from 0 to 2 eV, which is due to the reduction of the electronic band gap of MoS2.
NASA Astrophysics Data System (ADS)
Silvestrelli, Pier Luigi; Ambrosetti, Alberto
2016-02-01
The DFT/vdW-WF2s1 method based on the generation of localized Wannier functions, recently developed to include the van der Waals interactions in the density functional theory and describe adsorption processes on metal surfaces by taking metal-screening effects into account, is applied to the case of the interaction of Xe with noble-metal surfaces, namely Ag(111), Au(111), and Cu(111). The study is also repeated by adopting the DFT/vdW-QHO-WF variant relying on the quantum harmonic oscillator model which describes well many body effects. Comparison of the computed equilibrium binding energies and distances, and the C_3 coefficients characterizing the adatom-surface van der Waals interactions, with available experimental and theoretical reference data shows that the methods perform well and elucidates the importance of properly including screening effects. The results are also compared with those obtained by other vdW-corrected DFT schemes, including PBE-D, vdW-DF, vdW-DF2, rVV10, and by the simpler local density approximation and semi-local (PBE) generalized gradient approximation approaches.
DiLabio, Gino A.; Koleini, Mohammad
2014-05-14
Dispersion-correcting potentials (DCPs) are atom-centered Gaussian functions that are applied in a manner that is similar to effective core potentials. Previous work on DCPs has focussed on their use as a simple means of improving the ability of conventional density-functional theory methods to predict the binding energies of noncovalently bonded molecular dimers. We show in this work that DCPs developed for use with the LC-ωPBE functional along with 6-31+G(2d,2p) basis sets are capable of simultaneously improving predicted noncovalent binding energies of van der Waals dimer complexes and covalent bond dissociation enthalpies in molecules. Specifically, the DCPs developed herein for the C, H, N, and O atoms provide binding energies for a set of 66 noncovalently bonded molecular dimers (the “S66” set) with a mean absolute error (MAE) of 0.21 kcal/mol, which represents an improvement of more than a factor of 10 over unadorned LC-ωPBE/6-31+G(2d,2p) and almost a factor of two improvement over LC-ωPBE/6-31+G(2d,2p) used in conjunction with the “D3” pairwise dispersion energy corrections. In addition, the DCPs reduce the MAE of calculated X-H and X-Y (X,Y = C, H, N, O) bond dissociation enthalpies for a set of 40 species from 3.2 kcal/mol obtained with unadorned LC-ωPBE/6-31+G(2d,2p) to 1.6 kcal/mol. Our findings demonstrate that broad improvements to the performance of DFT methods may be achievable through the use of DCPs.
Zibordi-Besse, Larissa; Tereshchuk, Polina; Chaves, Anderson S; Da Silva, Juarez L F
2016-06-23
Transition-metal (TM) nanoparticles supported on oxides or carbon black have attracted much attention as potential catalysts for ethanol steam reforming reactions for hydrogen production. To improve the performance of nanocatalysts, a fundamental understanding of the interaction mechanism between water and ethanol with finite TM particles is required. In this article, we employed first-principles density functional theory with van der Waals (vdW) corrections to investigate the interaction of ethanol and water with TM13 clusters, where TM = Ni, Cu, Pd, Ag, Pt, and Au. We found that both water and ethanol bind via the anionic O atom to onefold TM sites, while at higher-energy structures, ethanol binds also via the H atom from the CH2 group to the TM sites, which can play an important role at real catalysts. The putative global minimum TM13 configurations are only slightly affected upon the adsorption of water or ethanol; however, for few systems, the compact higher-energy icosahedron structure changes its configuration upon ethanol or water adsorption. That is, those configurations are only shallow local minimums in the phase space. Except few deviations, we found similar trends for the magnitude of the adsorption energies of water and ethanol, that is, Ni13 > Pt13 > Pd13 and Cu13 > Au13 > Ag13, which is enhanced by the addition of the vdW correction (i.e., from 4% to 62%); however, the trend is the same. We found that the magnitude of the adsorption energy increases by shifting the center of gravity of the d-states toward the highest occupied molecular orbital. On the basis of the Mulliken and Hirshfeld charge analysis, as well as electron density differences, we identified the location of the charge redistribution and a tiny charge transfer (from 0.01 e to 0.19 e) from the molecules to the TM13 clusters. Our vibrational analysis indicates the red shifts in the OH modes upon binding of both water and ethanol molecules to the TM13 clusters, suggesting a weakening of
Redden, Robert A; Doolin, Edward J
2006-01-01
Fibroblasts are important cellular components in wound healing, scar formation, and fibrotic disorders; and the fibroblast-populated collagen-gel (FPCG) model allows examination of fibroblast behavior in an in vitro three-dimensional environment similar to that in vivo. Contraction of free-floating FPCGs depends on an active and dynamic cytoskeleton, and the contraction dynamics are highly influenced by cell density. We investigated mechanistic differences between high- and low-cell density FPCG contraction by evaluating contraction dynamics in detail, using specific cytoskeletal disruptors. Collagen gels were seeded with human lung fibroblasts at either high (HD) or low (LD) density, and incubated with or without cytoskeletal disruptors colchicine (microtubules) or cytochalasin D (microfilaments). Gel area was measured daily. FPCG contraction curves were essentially sigmoidal, featuring an initial period of no contraction (lag phase), followed by a period of rapid contraction (log phase). Contraction curves of HD-FPCGs were distinct from those of LD-FPCGs. For example, HD-FPCGs had a negligible lag phase (compared with 3 d for LD-FPCGs) and exhibited a higher rate of log-phase contraction. Both colchicine and cytochalasin dose-dependently inhibited contraction but specifically affected different phases of contraction in HD- and LD-FPCGs; and colchicine inhibited LD-FPCGs much more than HD-FPCGs. The data indicate that LD- and HD-FPCGs contract through different primary mechanisms. Microtubules and microfilaments are both complementarily and dynamically involved in the contraction of FPCGs, and cell density influences primary cytoskeletal mechanisms. These results provide valuable information about fibroblast behavior in healing and fibrosis, and may suggest novel treatment options. PMID:16759151
Volume correction in computed tomography densitometry for follow-up studies on pulmonary emphysema.
Stoel, Berend C; Putter, Hein; Bakker, M Els; Dirksen, Asger; Stockley, Rob A; Piitulainen, Eeva; Russi, Erich W; Parr, David; Shaker, Saher B; Reiber, Johan H C; Stolk, Jan
2008-12-15
Lung densitometry in drug evaluation trials can be confounded by changes in inspiration levels between computed tomography (CT) scans, limiting its sensitivity to detect changes over time. Therefore our aim was to explore whether the sensitivity of lung densitometry could be improved by correcting the measurements for changes in lung volume, based on the estimated relation between density (as measured with the 15th percentile point) and lung volume. We compared four correction methods, using CT data of 143 patients from five European countries. Patients were scanned, generally twice per visit, at baseline and after 2.5 years. The methods included one physiological model and three linear mixed-effects models using a volume-density relation: (1) estimated over the entire population with one scan per visit (model A) and two scans per visit (model B); and (2) estimated for each patient individually (model C). Both log-transformed and original volume and density values were evaluated and the differences in goodness-of-fit between methods were tested. Model C fitted best (P < 0.0001, P < 0.0001, and P = 0.064), when two scans were available. The most consistent progression estimation was obtained between sites, when both volume and density were log-transformed. Sensitivity was improved using repeated CT scans by applying volume correction to individual patient data. Volume correction reduces the variability in progression estimation by a factor of two, and is therefore recommended. PMID:19056717
Santana Palacio, Juan A.; Kim, Jeongnim; Kent, Paul R.; Reboredo, Fernando A.
2014-10-28
We have evaluated the successes and failures of the Hubbard-corrected density functional theory approach to study Mg doping of LiCoO_{2}. We computed the effect of the U parameter on the energetic, geometric, and electronic properties of two possible doping mechanisms: (1) substitution of Mg onto a Co (or Li) site with an associated impurity state and (2) formation of impurity-state-free complexes of substitutional Mg and point defects in LiCoO_{2}. We find that formation of impurity states results in changes on the valency of Co in LiCoO_{2}. Variation of the Co U shifts the energy of the impurity state, resulting in energetic, geometric, and electronic properties that depend significantly on the specific value of U. In contrast, the properties of the impurity-state-free complexes are insensitive to U. These results identify reasons for the strong dependence on the doping properties on the chosen value of U and for the overall difficulty of achieving agreement with the experimentally known energetic and electronic properties of doped transition metal oxides such as LiCoO_{2}.
Nakata, Ayako; Tsuneda, Takao; Hirao, Kimihiko
2011-12-14
A long-range corrected (LC) time-dependent density functional theory (TDDFT) incorporating relativistic effects with spin-orbit couplings is presented. The relativistic effects are based on the two-component zeroth-order regular approximation Hamiltonian. Before calculating the electronic excitations, we calculated the ionization potentials (IPs) of alkaline metal, alkaline-earth metal, group 12 transition metal, and rare gas atoms as the minus orbital (spinor) energies on the basis of Koopmans' theorem. We found that both long-range exchange and spin-orbit coupling effects are required to obtain Koopmans' IPs, i.e., the orbital (spinor) energies, quantitatively in DFT calculations even for first-row transition metals and systems containing large short-range exchange effects. We then calculated the valence excitations of group 12 transition metal atoms and the Rydberg excitations of rare gas atoms using spin-orbit relativistic LC-TDDFT. We found that the long-range exchange and spin-orbit coupling effects significantly contribute to the electronic spectra of even light atoms if the atoms have low-lying excitations between orbital spinors of quite different electron distributions.
NASA Astrophysics Data System (ADS)
Anatole von Lilienfeld, O.
2013-08-01
Generalised gradient approximated (GGA) density functional theory (DFT) typically overestimates polarisability and bond-lengths, and underestimates force constants of covalent bonds. To overcome this problem we show that one can use empirical force correcting atom centred potentials (FCACPs), parametrised for every nuclear species. Parameters are obtained through minimisation of a penalty functional that explicitly encodes hybrid DFT forces and static polarisabilities of reference molecules. For hydrogen, fluorine, chlorine and carbon the respective reference molecules consist of H2, F2, Cl2 and CH4. The transferability of this approach is assessed for harmonic frequencies in a small set of chlorofluorocarbon molecules. Numerical evidence, gathered for CF4, CCl4, CCl3F, CCl2F2, CClF3, ClF, HF, HCl, CFH3, CF2H2, CF3H, CHCl3, CH2Cl2 and CH3Cl indicates that the GGA+FCACP level of theory yields harmonic frequencies that are significantly more consistent with hybrid DFT values, as well as slightly reduced molecular polarisability.
Rykounov, A. A.
2015-06-07
The influence of pressure on the thermodynamic, structural, and elastic properties of the 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) molecular crystal at T = 0 is systematically studied. Calculations are carried out using density functional theory methods in a plane wave basis set with dispersion corrections for the exchange-correlation part of total energy, and ultrasoft pseudopotentials. The equilibrium unit cell parameters, the cold compression curve in the pressure range of 0–50 GPa and the sound speeds are computed. The effect of finite pressure on the molecular structure of TATB is elucidated from the analysis of relative changes in the intra- and intermolecular geometrical parameters. For the first time, the full set of elastic constants of this crystal at zero and non-zero pressures is determined from ab initio calculations. The resulted structural, elastic, and acoustic properties of TATB are shown to be in a good agreement with available experimental and theoretical data.
Lu, Hongxiu; Zhang, Hongguang; Zhang, Dong; Lu, Hongwang; Ma, Dedong
2015-01-01
Background Early detection of cancer is critical and is expected to contribute significantly to the success of cancer therapy and improvement of patient survival rates. Material/Methods A biocompatible, reconstituted, high-density lipoprotein (rHDL)-based nano-system containing calcium carbonate and near-infrared fluorescence dye (NIRF), methylene blue (MB), was fabricated and characterized by particle size, zeta potential, and morphology observation. The safety profile was confirmed by bovine serum albumin (BSA) challenge assay, hemolysis test, MTT assay, and in vivo long-term toxicity assay. The tumor targetability was assessed by cellular uptake, competitive inhibition experiments, and in vivo imaging assay. Results The self-assembled rHDL/MB/CCPs exhibited desirable and homogenous particle size, neutral surface charges, high bovine serum albumin stability, low hemolytic activity, and negligible cytotoxicity in vitro. The results obtained from confocal scanning laser microscopy and flow cytometry indicated that SR-BI coating exerted tumor-targeting function, which induced high and specific cellular uptake of rHDL/MB/CCPs. In vivo investigation in an A549 tumor xenografts-bearing mouse model revealed that rHDL/MB/CCPs possessed strong tumor targetability. Conclusions rHDL/MB/CCPs could be a safe tumor-targeting probe for cancer detection. PMID:26365043
... the NHLBI on Twitter. What Is a Lung Transplant? A lung transplant is surgery to remove a person's diseased lung ... a healthy lung from a deceased donor. Lung transplants are used for people who are likely to ...
Shi, Hongchang; Wang, Yilei; Hua, Ruimao
2015-12-28
Correction for 'Acid-catalyzed carboxylic acid esterification and ester hydrolysis mechanism: acylium ion as a sharing active intermediate via a spontaneous trimolecular reaction based on density functional theory calculation and supported by electrospray ionization-mass spectrometry' by Hongchang Shi et al., Phys. Chem. Chem. Phys., 2015, DOI: 10.1039/c5cp02914g. PMID:26583937
Xiang, H; Li, B; Behrman, R; Russo, G; Kachnic, L; Lu, H; Fernando, H
2015-06-15
Purpose: To measure the CT density model variations between different CT scanners used for treatment planning and impact on the accuracy of MC dose calculation in lung SBRT. Methods: A Gammex electron density phantom (RMI 465) was scanned on two 64-slice CT scanners (GE LightSpeed VCT64) and a 16-slice CT (Philips Brilliance Big Bore CT). All three scanners had been used to acquire CT for CyberKnife lung SBRT treatment planning. To minimize the influences of beam hardening and scatter for improving reproducibility, three scans were acquired with the phantom rotated 120° between scans. The mean CT HU of each density insert, averaged over the three scans, was used to build the CT density models. For 14 patient plans, repeat MC dose calculations were performed by using the scanner-specific CT density models and compared to a baseline CT density model in the base plans. All dose re-calculations were done using the same plan beam configurations and MUs. Comparisons of dosimetric parameters included PTV volume covered by prescription dose, mean PTV dose, V5 and V20 for lungs, and the maximum dose to the closest critical organ. Results: Up to 50.7 HU variations in CT density models were observed over the baseline CT density model. For 14 patient plans examined, maximum differences in MC dose re-calculations were less than 2% in 71.4% of the cases, less than 5% in 85.7% of the cases, and 5–10% for 14.3% of the cases. As all the base plans well exceeded the clinical objectives of target coverage and OAR sparing, none of the observed differences led to clinically significant concerns. Conclusion: Marked variations of CT density models were observed for three different CT scanners. Though the differences can cause up to 5–10% differences in MC dose calculations, it was found that they caused no clinically significant concerns.
Greuter, M. J. W.; Groen, J. M.; Nicolai, L. J.; Dijkstra, H.; Oudkerk, M.
2009-11-15
Purpose: The objective of this study is to quantify the influence of linear motion, calcification density, and temporal resolution on coronary calcium determination using multidetector computed tomography (MDCT), dual source CT (DSCT), and electron beam tomography (EBT) and to find a quantitative method which corrects for the influences of these parameters using a linear moving cardiac phantom. Methods: On a robotic arm with artificial arteries with four calcifications of increasing density, a linear movement was applied between 0 and 120 mm/s (step of 10 mm/s). The phantom was scanned five times on 64-slice MDCT, DSCT, and EBT using a standard acquisition protocol. The average Agatston, volume, and mass scores were determined for each velocity, calcification, and scanner. Susceptibility to motion was quantified using a cardiac motion susceptibility (CMS) index. Resemblance to EBT and physical volume and mass was quantified using a {Delta} index. Results: Increasing motion artifacts were observed at increasing velocities on all scanners, with increasing severity from EBT to DSCT to 64-slice MDCT. The calcium score showed a linear dependency on motion from which a correction factor could be derived. This correction factor showed a linear dependency on the mean calcification density with a good fit for all three scoring methods and all three scanners (0.73{<=}R{sup 2}{<=}0.95). The slope and offset of this correction factor showed a linear dependency on temporal resolution with a good fit for all three scoring methods and all three scanners (0.83{<=}R{sup 2}{<=}0.98). CMS was minimal for EBT and increasing values were observed for DSCT and highest values for 64-slice MDCT. CMS was minimal for mass score and increasing values were observed for volume score and highest values for Agatston score. For all densities and scoring methods DSCT showed on average the closest resemblance to EBT calcium scores. When using the correction factor, CMS index decreased on average by
NASA Astrophysics Data System (ADS)
Casida, Mark E.; Gutierrez, Fabien; Guan, Jingang; Gadea, Florent-Xavier; Salahub, Dennis; Daudey, Jean-Pierre
2000-11-01
Time-dependent density-functional theory (TDDFT) is an increasingly popular approach for calculating molecular excitation energies. However, the TDDFT lowest triplet excitation energy, ωT, of a closed-shell molecule often falls rapidly to zero and then becomes imaginary at large internuclear distances. We show that this unphysical behavior occurs because ωT2 must become negative wherever symmetry breaking lowers the energy of the ground state solution below that of the symmetry unbroken solution. We use the fact that the ΔSCF method gives a qualitatively correct first triplet excited state to derive a "charge-transfer correction" (CTC) for the time-dependent local density approximation (TDLDA) within the two-level model and the Tamm-Dancoff approximation (TDA). Although this correction would not be needed for the exact exchange-correlation functional, it is evidently important for a correct description of molecular excited state potential energy surfaces in the TDLDA. As a byproduct of our analysis, we show why TDLDA and LDA ΔSCF excitation energies are often very similar near the equilibrium geometries. The reasoning given here is fairly general and it is expected that similar corrections will be needed in the case of generalized gradient approximations and hybrid functionals.
Burns, Lori A; Sherrill, David; Vazquez-Mayagoitia, Alvaro; Sumpter, Bobby G
2011-01-01
A systematic study of techniques for treating non-covalent interactions within the computationally efficient density functional theory (DFT) framework is presented through comparison to benchmark-quality evaluations of binding strength com- piled for molecular complexes of diverse size and nature. In particular, the effi- cacy of functionals deliberately crafted to encompass long-range forces, a posteri- ori DFT+dispersion corrections (DFT-D2 and DFT-D3), and exchange-hole dipole moment (XDM) theory is assessed against a large collection (469 energy points) of reference interaction energies at the CCSD(T) level of theory extrapolated to the estimated complete basis set limit. The established S22 and JSCH test sets of minimum-energy structures, as well as collections of dispersion-bound (NBC10) and hydrogen-bonded (HBC6) dissociation curves and a pairwise decomposition of a protein-ligand reaction site (HSG), comprise the chemical systems for this work. From evaluations of accuracy, consistency, and efficiency for PBE-D, BP86-D, B97-D, PBE0-D, B3LYP-D, B970-D, M05-2X, M06-2X, B97X-D, B2PLYP-D, XYG3, and B3LYP-XDM methodologies, it is concluded that distinct, often contrasting, groups of these elicit the best performance within the accessible double- or robust triple- basis set regimes and among hydrogen-bonded or dispersion-dominated complexes. For overall results, M05-2X, B97-D3, and B970-D2 yield superior values in conjunc- tion with aug-cc-pVDZ, for a mean absolute deviation of 0.41 0.49 kcal/mol, and B3LYP-D3, B97-D3, B97X-D, and B2PLYP-D3 dominate with aug-cc-pVTZ, af- fording, together with XYG3/6-311+G(3df,2p), a mean absolute deviation of 0.33 0.38 kcal/mol.
Zhu, Wei; Germain, Claire; Liu, Zheng; Sebastian, Yinong; Devi, Priyanka; Knockaert, Samantha; Brohawn, Philip; Lehmann, Kim; Damotte, Diane; Validire, Pierre; Yao, Yihong; Valge-Archer, Viia; Hammond, Scott A; Dieu-Nosjean, Marie-Caroline; Higgs, Brandon W
2015-01-01
T and B cell receptor (TCR and BCR, respectively) Vβ or immunoglobulin heavy chain complementarity-determining region 3 sequencing allows monitoring of repertoire changes through recognition, clonal expansion, affinity maturation, and T or B cell activation in response to antigen. TCR and BCR repertoire analysis can advance understanding of antitumor immune responses in the tumor microenvironment. TCR and BCR repertoires of sorted CD4+, CD8+ or CD19+ cells in tumor, non-tumoral distant tissue (NT), and peripheral compartments (blood/draining lymph node [P]) from 47 non-small cell lung cancer (NSCLC) patients (agemedian = 68 y) were sequenced. The clonotype spectra were assessed among different tissues and correlated with clinical and immunological parameters. In all tissues, CD4+ and CD8+ TCR repertoires had greater clonality relative to CD19+ BCR. CD4+ T cells exhibited greater clonality in NT compared to tumor (p = 0.002) and P (p < 0.001), concentrated among older patients (age > 68). Younger patients exhibited greater CD4+ T cell diversity in P compared to older patients (p = 0.05), and greater CD4+ T cell clonality in tumor relative to P (p < 0.001), with fewer shared clonotypes between tumor and P than older patients (p = 0.04). More interestingly, greater CD4+ and CD8+ T cell clonality in tumor and P, respectively (both p = 0.05), correlated with high density of tumor-associated tertiary lymphoid structure (TLS) B cells, a biomarker of higher overall survival in NSCLC. Results indicate distinct adaptive immune responses in NSCLC, where peripheral T cell diversity is modulated by age, and tumor T cell clonal expansion is favored by the presence of TLSs in the tumor microenvironment. PMID:26587322
Glubb, Dylan M; Cerri, Elisa; Giese, Alexandra; Zhang, Wei; Mirza, Osman; Thompson, Emma E.; Chen, Peixian; Das, Soma; Jassem, Jacek; Rzyman, Witold; Lingen, Mark W.; Salgia, Ravi; Hirsch, Fred R.; Dziadziuszko, Rafal; Ballmer-Hofer, Kurt; Innocenti, Federico
2011-01-01
Purpose VEGFR-2 plays a crucial role in mediating angiogenic endothelial cell responses via the VEGF pathway and angiogenesis inhibitors targeting VEGFR-2 are in clinical use. As angiogenesis is a host-driven process, functional heritable variation in KDR, the gene encoding VEGFR-2, may affect VEGFR-2 function, and ultimately, the extent of tumor angiogenesis. Experimental Design We resequenced KDR using 24 DNAs each from healthy Caucasian, African American and Asian groups. Non-synonymous genetic variants were assessed for function using phosphorylation assays. Luciferase reporter gene assays were used to examine effects of variants on gene expression. KDR mRNA and protein expression, and microvessel density (MVD) were measured in non-small cell lung cancer (NSCLC) tumor samples and matching patient DNA samples were genotyped to test for associations with variants of interest. Results KDR resequencing led to the discovery of 120 genetic variants, of which 25 had not been previously reported. Q472H had increased VEGFR-2 protein phosphorylation and associated with increased MVD in NSCLC tumor samples. −2854C and −2455A increased luciferase expression and associated with higher KDR mRNA levels in NSCLC samples. −271A reduced luciferase expression and associated with lower VEGFR-2 levels in NSCLC samples. −906C and 23408G, associated with higher KDR mRNA levels in NSCLC samples. Conclusions This study has defined KDR genetic variation in three populations and identified common variants that impact on tumoral KDR expression and vascularization. These findings may have important implications for understanding the molecular basis of genetic associations between KDR variation and clinical phenotypes related to VEGFR-2 function. PMID:21712447
Müller, Philipp; Rothschild, Sacha I; Arnold, Walter; Hirschmann, Petra; Horvath, Lukas; Bubendorf, Lukas; Savic, Spasenija; Zippelius, Alfred
2016-01-01
Tumor-infiltrating lymphocytes play an important role in cell-mediated immune destruction of cancer cells and tumor growth control. We investigated the heterogeneity of immune cell infiltrates between primary non-small cell lung carcinomas (NSCLC) and corresponding metastases. Formalin-fixed, paraffin-embedded primary tumors and corresponding metastases from 34 NSCLC patients were analyzed by immunohistochemistry for CD4, CD8, CD11c, CD68, CD163 and PD-L1. The percentage of positively stained cells within the stroma and tumor cell clusters was recorded and compared between primary tumors and metastases. We found significantly fewer CD4(+) and CD8(+) T cells within tumor cell clusters as compared with the stromal compartment, both in primary tumors and corresponding metastases. CD8(+) T cell counts were significantly lower in metastatic lesions than in the corresponding primary tumors, both in the stroma and the tumor cell islets. Of note, the CD8/CD4 ratio was significantly reduced in metastatic lesions compared with the corresponding primary tumors in tumor cell islets, but not in the stroma. We noted significantly fewer CD11c(+) cells and CD68(+) as well as CD163(+) macrophages in tumor cell islets compared with the tumor stroma, but no difference between primary and metastatic lesions. Furthermore, the CD8/CD68 ratio was higher in primary tumors than in the corresponding metastases. We demonstrate a differential pattern of immune cell infiltration in matched primary and metastatic NSCLC lesions, with a significantly lower density of CD8(+) T cells in metastatic lesions compared with the primary tumors. The lower CD8/CD4 and CD8/CD68 ratios observed in metastases indicate a rather tolerogenic and tumor-promoting microenvironment at the metastatic site. PMID:26541588
... Emergencies Cardiac Emergencies Eye Emergencies Lung Emergencies Surgeries Lung Emergencies People with Marfan syndrome can be at ... should be considered an emergency. Symptoms of sudden lung collapse (pneumothorax) Symptoms of a sudden lung collapse ...
... version of this page please turn Javascript on. Lung Cancer What is Lung Cancer? How Tumors Form The body is made ... button on your keyboard.) Two Major Types of Lung Cancer There are two major types of lung ...
Metastases to the lung; Metastatic cancer to the lung ... Metastatic tumors in the lungs are cancers that developed at other places in the body (or other parts of the lungs) and spread through the ...
Aisner, J.
1985-01-01
This book contains 13 chapters. Some of the chapter titles are: The Pathology of Lung Cancer; Radiotherapy for Non-Small-Cell Cancer of the Lung; Chemotherapy for Non-Small-Cell Lung Cancer; Immunotherapy in the Management of Lung Cancer; Preoperative Staging and Surgery for Non-Small-Cell Lung Cancer; and Prognostic Factors in Lung Cancer.
Lung B-line artefacts and their use
Mathis, Gebhard; Blaivas, Michael; Volpicelli, Giovanni; Seibel, Armin; Wastl, Daniel; Atkinson, Nathan S. S.; Cui, Xin-Wu; Fan, Mei; Yi, Dong
2016-01-01
Background The analysis of lung artefacts has gained increasing importance as markers of lung pathology. B-line artefact (BLA), caused by a reverberation phenomenon, is the most important lung artefact. In this review, we discuss the current role of BLA in pneumology and explore open questions of the published consensus. Methods We summarized current literature about BLA. Also, we presented observations on healthy subjects and patients with interstitial syndrome (pulmonary fibrosis and edema), to investigate technical factors influencing BLA visualization. Results BLA imaging is influenced by more factors than recently assumed. When multiple BLA is visualized in the lung, they represent a sign of increased density due to the loss of aeration in the lung periphery. This condition may indicate different diseases including cardiogenic pulmonary edema, diffuse or focal interstitial lung diseases (ILD), infections and acute respiratory distress syndrome (ARDS). Correct interpretation of BLA in lung ultrasound is strongly influenced by associated sonographic signs and careful integration of all relevant clinical information. Conclusions BLA is useful to monitor clinical response, and may become crucial in directing the diagnostic process. Further research is warranted to clarify technical adjustments, different probe and machine factors that influence the visualization of BLA. PMID:27293860
Thierfelder, Christian; Schwerdtfeger, Peter; Saue, Trond
2007-09-15
The electric field gradient in late transition metal compounds is incorrectly determined by most density functionals. We show that the coupling of short-range density functional based with long-range wave function based methods using a reparametrization of the Coulomb-attenuated Becke three-parameter Lee-Yang-Parr approximation gives reliable results for the electric field gradients of copper and gold for a series of compounds. This results in nuclear quadrupole moments of -0.208 b for {sup 63}Cu and +0.526 b for {sup 197}Au in good agreement with experimental values of -0.220(15) and +0.547(16)b, respectively.
Technology Transfer Automated Retrieval System (TEKTRAN)
Impaired cell cholesterol trafficking in Niemann-Pick type C (NPC) disease results in the first known instance of impaired regulation of the ATP-binding cassette transporter A1 (ABCA1), a lipid transporter mediating the rate-limiting step in high density lipoprotein (HDL) formation, as a cause of lo...
... many disorders affecting the lungs, such as asthma, COPD, infections like influenza, pneumonia and tuberculosis, lung cancer, and many other breathing problems. Some lung diseases can lead to respiratory failure. Dept. of Health and Human Services Office on Women's Health
A collapsed lung happens when air enters the pleural space, the area between the lung and the chest wall. If it is a ... is called pneumothorax. If only part of the lung is affected, it is called atelectasis. Causes of ...
... balloon-like tissues (blebs) that cause lung collapse ( pneumothorax ) Wedge resection, to remove part of a lobe ... Treat injuries that cause lung tissue to collapse ( pneumothorax or hemothorax ) Treat permanently collapsed lung tissue ( atelectasis ) ...
A collapsed lung happens when air enters the pleural space, the area between the lung and the chest wall. If it is a total collapse, it is called pneumothorax. If only part of the lung is affected, ...
... the lungs to take in oxygen and release carbon dioxide. People with this type of lung disorder often ... the lungs to take up oxygen and release carbon dioxide. These diseases may also affect heart function. An ...
NASA Astrophysics Data System (ADS)
Jones, Andrew Osler
There is an increasing interest in the use of inhomogeneity corrections for lung, air, and bone in radiotherapy treatment planning. Traditionally, corrections based on physical density have been used. Modern algorithms use the electron density derived from CT images. Small fields are used in both conformal radiotherapy and IMRT, however their beam characteristics in inhomogeneous media have not been extensively studied. This work compares traditional and modern treatment planning algorithms to Monte Carlo simulations in and near low-density inhomogeneities. Field sizes ranging from 0.5 cm to 5 cm in diameter are projected onto a phantom containing inhomogeneities and depth dose curves are compared. Comparisons of the Dose Perturbation Factors (DPF) are presented as functions of density and field size. Dose Correction Factors (DCF), which scale the algorithms to the Monte Carlo data, are compared for each algorithm. Physical scaling algorithms such as Batho and Equivalent Pathlength (EPL) predict an increase in dose for small fields passing through lung tissue, where Monte Carlo simulations show a sharp dose drop. The physical model-based collapsed cone convolution (CCC) algorithm correctly predicts the dose drop, but does not accurately predict the magnitude. Because the model-based algorithms do not correctly account for the change in backscatter, the dose drop predicted by CCC occurs further downstream compared to that predicted by the Monte Carlo simulations. Beyond the tissue inhomogeneity all of the algorithms studied predict dose distributions in close agreement with Monte Carlo simulations. Dose-volume relationships are important in understanding the effects of radiation to the lung. Dose within the lung is affected by a complex function of beam energy, lung tissue density, and field size. Dose algorithms vary in their abilities to correctly predict the dose to the lung tissue. A thorough analysis of the effects of density, and field size on dose to the lung
Development and application of a random lung model for dose calculations in radiotherapy
NASA Astrophysics Data System (ADS)
Liang, Liang
Radiotherapy requires accurate dose calculations in the human body, especially in disease sites with large variations of electron density in neighboring tissues, such as the lung. Currently, the lung is modeled by a voxelized geometry interpolated from computed tomography (CT) scans to various resolutions. The simplest such voxelized lung, the atomic mix model, is a homogenized whole lung with a volume-averaged bulk density. However, according traditional transport theory, even the relatively fine CT voxelization of the lung is not valid, due to the extremely small mean free path (MFP) of the electrons. The purpose of this thesis is to study the impact of the lung's heterogeneities on dose calculations in lung treatment planning. We first extend the traditional atomic mix theory for charged particles by approximating the Boltzmann equation for electrons to its Fokker-Planck (FP) limit, and then applying a formal asymptotic analysis to the BFP equation. This analysis raises the length scale for homogenizing a heterogeneous medium from the electron mean free path (MFP) to the much larger electron transport MFP. Then, using the lung's anatomical data and our new atomic mix theory, we build a realistic 2 1/2-D random lung model. The dose distributions for representative realizations of the random lung model are compared to those from the atomic mix approximation of the random lung model, showing that significant perturbations may occur with small field sizes and large lung structures. We also apply our random lung model to a more realistic lung phantom and investigate the effect of CT resolutions on lung treatment planning. We show that, compared to the reference 1 x 1 mm2 CT resolution, a 2 x 2 mm2 CT resolution is sufficient to voxelize the lung, while significant deviations in dose can be observed with a larger 4 x 4 mm 2 CT resolution. We use the Monte Carlo method extensively in this thesis, to avoid systematic errors caused by inaccurate heterogeneity corrections