Simulation of electron density maps for two-dimensional crystal structures using Mathematica

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

CT Image Electron Density Quantification in Regions with Metal Implants: Implications for Radiotherapy Treatment Planning

In radiotherapy planning, computed tomography (CT) images are used to quantify the electron density of tissues and provide spatial anatomical information. Treatment planning systems use these data to calculate the expected spatial distribution of absorbed dose in a patient. CT imaging is complicated by the presence of metal implants which cause increased image noise, produce artifacts throughout the image and can exceed the available range of CT number values within the implant, perturbing electron density estimates in the image. Furthermore, current dose calculation algorithms do not accurately model radiation transport at metal-tissue interfaces. Combined, these issues adversely affect the accuracy of dose calculations in the vicinity of metal implants. As the number of patients with orthopedic and dental implants grows, so does the need to deliver safe and effective radiotherapy treatments in the presence of implants. The Medical Physics group at the Cancer Centre of Southeastern Ontario and Queen's University has developed a Cobalt-60 CT system that is relatively insensitive to metal artifacts due to the high energy, nearly monoenergetic Cobalt-60 photon beam. Kilovoltage CT (kVCT) images, including images corrected using a commercial metal artifact reduction tool, were compared to Cobalt-60 CT images throughout the treatment planning process, from initial imaging through to dose calculation. An effective metal artifact reduction algorithm was also implemented for the Cobalt-60 CT system. Electron density maps derived from the same kVCT and Cobalt-60 CT images indicated the impact of image artifacts on estimates of photon attenuation for treatment planning applications. Measurements showed that truncation of CT number data in kVCT images produced significant mischaracterization of the electron density of metals. Dose measurements downstream of metal inserts in a water phantom were compared to dose data calculated using CT images from kVCT and Cobalt-60 systems with and without artifact correction. The superior accuracy of electron density data derived from Cobalt-60 images compared to kVCT images produced calculated dose with far better agreement with measured results. These results indicated that dose calculation errors from metal image artifacts are primarily due to misrepresentation of electron density within metals rather than artifacts surrounding the implants.

Effect of basis set superposition error on the electron density of molecular complexes

The effect of basis set superposition error (BSSE) on molecular complexes is analyzed. The BSSE causes artificial delocalizations which modify the first order electron density. The mechanism of this effect is assessed for the hydrogen fluoride dimer with several basis sets. The BSSE-corrected first-order electron density is obtained using the chemical Hamiltonian approach versions of the Roothaan and Kohn-Sham equations. The corrected densities are compared to uncorrected densities based on the charge density critical points. Contour difference maps between BSSE-corrected and uncorrected densities on the molecular plane are also plotted to gain insight into the effects of BSSE correction on the electron density

Effect of basis set superposition error on the electron density of molecular complexes

The effect of basis set superposition error (BSSE) on molecular complexes is analyzed. The BSSE causes artificial delocalizations which modify the first order electron density. The mechanism of this effect is assessed for the hydrogen fluoride dimer with several basis sets. The BSSE-corrected first-order electron density is obtained using the chemical Hamiltonian approach versions of the Roothaan and Kohn-Sham equations. The corrected densities are compared to uncorrected densities based on the charge density critical points. Contour difference maps between BSSE-corrected and uncorrected densities on the molecular plane are also plotted to gain insight into the effects of BSSE correction on the electron density

Structure, stability, depolarized light scattering, and vibrational spectra of fullerenols from all-electron density-functional-theory calculations

We have investigated the stability, electronic properties, Rayleigh (elastic), and Raman (inelastic) depolarization ratios, infrared and Raman absorption vibrational spectra of fullerenols [C(60)(OH)(n)] with different degrees of hydroxylation by using all-electron density-functional-theory (DFT) methods. Stable arrangements of these molecules were found by means of full geometry optimizations using Becke's three-parameter exchange functional with the Lee, Yang, and Parr correlation functional. This DFT level has been combined with the 6-31G(d,p) Gaussian-type basis set, as a compromise between accuracy and capability to treat highly hydroxylated fullerenes, e.g., C(60)(OH)(36). Thus, the molecular properties of fullerenols were systematically analyzed for structures with n=1, 2, 3, 4, 8, 10, 16, 18, 24, 32, and 36. From the electronic structure analysis of these molecules, we have evidenced an important effect related to the weak chemical reactivity of a possible C(60)(OH)(24) isomer. To investigate Raman scattering and the vibrational spectra of the different fullerenols, frequency calculations are carried out within the harmonic approximation. In this case a systematic study is only performed for n=1-4, 8, 10, 16, 18, and 24. Our results give good agreements with the expected changes in the spectral absorptions due to the hydroxylation of fullerenes.

A quantum molecular similarity analysis of changes in molecular electron density caused by basis set flotation and electric field application

Quantum molecular similarity (QMS) techniques are used to assess the response of the electron density of various small molecules to application of a static, uniform electric field. Likewise, QMS is used to analyze the changes in electron density generated by the process of floating a basis set. The results obtained show an interrelation between the floating process, the optimum geometry, and the presence of an external field. Cases involving the Le Chatelier principle are discussed, and an insight on the changes of bond critical point properties, self-similarity values and density differences is performed

Barrier for the reaction X20+ + X20+ -> X402+ in alkali-metal clusters related to electron density at the bond midpoint of the supermolecule (X20+)2

Using the extended Thomas-Fermi version of density-functional theory (DFT), calculations are presented for the barrier for the reaction Na20++Na20+¿Na402+. The deviation from the simple Coulomb barrier is shown to be proportional to the electron density at the bond midpoint of the supermolecule (Na20+)2. An extension of conventional quantum-chemical studies of homonuclear diatomic molecular ions is then effected to apply to the supermolecular ions of the alkali metals. This then allows the Na results to be utilized to make semiquantitative predictions of position and height of the maximum of the fusion barrier for other alkali clusters. These predictions are confirmed by means of similar DFT calculations for the K clusters.

Inverting ionospheric D-region electron density pro-file from riometer measurement

The aim of this work is to invert the ionospheric electron density profile from Riometer (Relative Ionospheric opacity meter) measurement. The newly Riometer instrument KAIRA (Kilpisjärvi Atmospheric Imaging Receiver Array) is used to measure the cosmic HF radio noise absorption that taking place in the D-region ionosphere between 50 to 90 km. In order to invert the electron density profile synthetic data is used to feed the unknown parameter Neq using spline height method, which works by taking electron density profile at different altitude. Moreover, smoothing prior method also used to sample from the posterior distribution by truncating the prior covariance matrix. The smoothing profile approach makes the problem easier to find the posterior using MCMC (Markov Chain Monte Carlo) method.

Temporal and Spatial Behavior of Electron Density and Temperature in a Laser-Produced Plasma from YBa2Cu3O7

Spectroscopic studies of laser -induced plasma from a high-temperature superconducting material, viz., YBa2Cu3O7 (YBCO), have been carried out. Electron temperature and electron density measurements were made from spectral data. The Stark broad ening of emission lines was used to determine the electron density, and the ratio of line in tensities was exploited for the determination of electron temperature. An initial electron temperature of 2.35 eV and electron density of 2.5 3 1017 cm2 3 were observed. The dependence on electron temperature and density on different experimental parameters such as distance from the target, delay time after the in itiation of the plasm a, and laser irradiance is also discussed in detail. Index Headings: Laser -plasma spectroscopy; Plasma diagnostics; Emission spectroscop y; YBa2Cu3O7.

Time evolution of the electron density and temperature in laser-produced plasmas from YBa2Cu3O7

Laser radiation at 1.06 µm from a pulsed Nd:YAG laser was focused onto a multielement YBa2Cu3O7 target in vacuum and the plasma thus generated was studied using time-resolved spectroscopic techniques. Line broadening of the Ba I emission line at 553.5 nm was monitored as a function of time elapsed after the incidence of a laser pulse on the target. Measured line profiles of barium species were used to infer the electron density and temperature, and the time evolution of these important plasma parameters has been worked out.

Electron density and temperature measurements in a laser produced carbon plasma

Plasma generated by fundamental radiation from a Nd:YAG laser focused onto a graphite target is studied spectroscopically. Measured line profiles of several ionic species were used to infer electron temperature and density at several sections located in front of the target surface. Line intensities of successive ionization states of carbon were used for electron temperature calculations. Stark broadened profiles of singly ionized species have been utilized for electron density measurements. Electron density as well as electron temperature were studied as functions of laser irradiance and time elapsed after the incidence of laser pulse. The validity of the assumption of local thermodynamic equilibrium is discussed in light of the results obtained.

Electron Density Determination of Laser Induced Plasma from Polymethyl Methacrylate Using Phaseshift Detection Technique

Irradiation of a Polymethyl methacrylate target using a pulsed Nd-YAG laser causes plasma formation in the vicinity of the target. The refractive index gradient due to the presence of the plasma is probed using phase-shift detection technique. The phase-shift technique is a simple but sensitive technique for the determination of laser ablation threshold of solids. The number density of laser generated plasma above the ablation threshold from Polymethyl methacrylate is calculated as a function of laser fluence. The number density varies from 2×1016 cm-3 to 2×1017 cm-3 in the fluence interval 2.8-13 J · cm-2.

A quantum molecular similarity analysis of changes in molecular electron density caused by basis set flotation and electric field application

Quantum molecular similarity (QMS) techniques are used to assess the response of the electron density of various small molecules to application of a static, uniform electric field. Likewise, QMS is used to analyze the changes in electron density generated by the process of floating a basis set. The results obtained show an interrelation between the floating process, the optimum geometry, and the presence of an external field. Cases involving the Le Chatelier principle are discussed, and an insight on the changes of bond critical point properties, self-similarity values and density differences is performed