-
The effect of inter-track coupling on H$_2$O$_2$ productions
Authors:
Ramin Abolfath,
Sedigheh Fardirad,
Abbas Ghasemizad
Abstract:
Background: Lower production of H$_2$O$_2$ in water is a hallmark of ultra-high dose rate (UHDR) compared to the conventional dose rate (CDR). However, the current computational models based on the predicted yield of H$_2$O$_2$ are in opposite of the experimental data. Methods: We construct an analytical model for the rate equation in the production of H$_2$O$_2$ from \ce{^{.}OH}-radicals and use…
▽ More
Background: Lower production of H$_2$O$_2$ in water is a hallmark of ultra-high dose rate (UHDR) compared to the conventional dose rate (CDR). However, the current computational models based on the predicted yield of H$_2$O$_2$ are in opposite of the experimental data. Methods: We construct an analytical model for the rate equation in the production of H$_2$O$_2$ from \ce{^{.}OH}-radicals and use it as a guide to propose a hypothetical geometrical inhomogeneity in the configuration of particles in the FLASH-UHDR beams. We perform a series of Monte Carlo (MC) simulations of the track structures for a system of charged particles impinging the medium in the form of clusters and/or bunches. Results: We demonstrate the interplay of diffusion, reaction rates, and overlaps in track-spacing attribute to a lower yield of H$_2$O$_2$ at FLASH-UHDR vs. CDR. This trend is reversed if spacing among the tracks becomes larger than a critical value, with a length scale that is proportional to the diffusion length of \ce{^{.}OH}-radicals modulated by a rate of decay due to recombination with other species, available within a track, and the space among the tracks. The latter is substantial on the suppressing of the H$_2$O$_2$ population at FLASH-UHDR relative to CDR. Conclusions: Based on our analysis of the present work, at FLASH-UHDR, the lower yield in H$_2$O$_2$ can be interpreted as a signature of bunching the particles in beams of ionizing radiation. The beams enter the medium in closely packed clusters and form inhomogeneities in the track-structure distribution. Thus the MC simulations based on the assumption of uniformly distributed tracks are unable to explain the experimental data.
△ Less
Submitted 25 March, 2024;
originally announced March 2024.
-
A molecular dynamics simulation of the abrupt changes in the thermodynamic properties of water after formation of nano-bubbles / nano-cavities induced by passage of charged particles
Authors:
Ramin Abolfath,
Niayesh Afshordi,
Sohrab Rahvar,
Adri van Duin,
Martin Radler,
Reza Taleei,
Katia Parodi,
Julie Lascaud,
Radhe Mohan
Abstract:
We present a multi-scale formalism that accounts for the formation of nano-scale bubbles/cavities owing to a burst of water molecules after the passage of high energy charged particles that leads to the formation of hot non-ionizing excitations or thermal spikes (TS). We demonstrate the coexistence of a rapidly growing condensed state of water and a hot spot that forms a stable state of diluted wa…
▽ More
We present a multi-scale formalism that accounts for the formation of nano-scale bubbles/cavities owing to a burst of water molecules after the passage of high energy charged particles that leads to the formation of hot non-ionizing excitations or thermal spikes (TS). We demonstrate the coexistence of a rapidly growing condensed state of water and a hot spot that forms a stable state of diluted water at high temperatures and pressures, possibly at a supercritical phase. Depending on the temperature of TS, the thin shell of a highly dense state of water grows by three to five times the speed of sound in water, forming a thin layer of shock wave (SW) buffer, wrapping around the nano-scale cylindrical symmetric cavity. The stability of the cavity, as a result of the incompressibility of water at ambient conditions and the surface tension, allows the transition of supersonic SW to a subsonic contact discontinuity and dissipation to thermo-acoustic sound waves. We further study the mergers of nanobubbles that lead to fountain-like or jet-flow structures at the collision interface. We introduce a time delay in the nucleation of nano-bubbles, a novel mechanism, responsible for the growth and stability of much larger or even micro-bubbles, possibly relevant to FLASH ultra-high dose rate (UHDR). The current study is potentially significant at FLASH-UHDRs. Our analysis predicts the black-body radiation from the transient supercritical state of water localized in nano-cavities wrapping around the track of charged particles can be manifested in the (indirect) water luminescence spectrum.
△ Less
Submitted 9 March, 2024;
originally announced March 2024.
-
Differential tissue sparing of FLASH ultra high dose rates: an {\it in-silico} study
Authors:
Ramin Abolfath,
Alexander Baikalov,
Sohrab Rahvar,
Alberto Fraile,
Stefan Bartzsch,
Emil Schüler,
Radhe Mohan
Abstract:
Purpose: To propose a theory for the differential tissue sparing of FLASH ultra high dose rates (UHDR) through inter-track reaction-diffusion mechanism. Methods: We calculate the time-evolution of particle track-structures using a system of coupled reaction-diffusion equations on a random network designed for molecular transport in porous and disordered media. The network is representative of the…
▽ More
Purpose: To propose a theory for the differential tissue sparing of FLASH ultra high dose rates (UHDR) through inter-track reaction-diffusion mechanism. Methods: We calculate the time-evolution of particle track-structures using a system of coupled reaction-diffusion equations on a random network designed for molecular transport in porous and disordered media. The network is representative of the intra- and inter-cellular diffusion channels in tissues. Spatial cellular heterogeneities over the scale of track spacing have been constructed by incorporating random fluctuations in the connectivity among network sites. Results: We demonstrate the occurrence of phase separation among the tracks as the complexity in intra- and inter-cellular structural increases. At the weak limit of disorder, such as in water and normal tissue, neighboring tracks melt into each other and form a percolated network of nonreactive species. In contrast, at the strong limit of disorder, tracks evolve individually like isolated islands with negligible inter-track overlap. Thus, the spatio-temporal correlation among the chemical domains decreases as the inter-cellular complexity of the tissue increases (e.g. from normal tissue to fractal-type malignant tissue). Conclusions: The differential sparing of FLASH UHDR in normal and tumor tissue may be explained by differences in inter- and intra-cellular structural complexities between the tissue types. The structural complexities of cancerous cells prevent clustering and chemical interaction of tracks, whereas this interaction prevails and thus leads to sparing in normal tissue.
△ Less
Submitted 7 October, 2022;
originally announced October 2022.
-
Modeling interspur interactions as a potential mechanism of the FLASH effect
Authors:
Alexander Baikalov,
Ramin Abolfath,
Radhe Mohan,
Emil Schüler,
Jan J. Wilkens,
Stefan Bartzsch
Abstract:
The mechanism responsible for the FLASH effect, normal tissue sparing by ultra-high dose rate (UHDR) irradiation with isoeffective tumor control compared to conventional dose rate (CDR) irradiation, remains undetermined. Here we investigate the contribution of interspur interactions (interactions between radiolytic species of individual particle tracks) to overall radiochemical interactions as a f…
▽ More
The mechanism responsible for the FLASH effect, normal tissue sparing by ultra-high dose rate (UHDR) irradiation with isoeffective tumor control compared to conventional dose rate (CDR) irradiation, remains undetermined. Here we investigate the contribution of interspur interactions (interactions between radiolytic species of individual particle tracks) to overall radiochemical interactions as a function of irradiation parameters, and suggest an increase in interspur interaction as a potential mechanism for tissue sparing in FLASH radiation therapy. We construct a model that analytically represents the spatiotemporal distribution of spurs in a target volume as a function of irradiation parameters (e.g. dose, dose rate, linear energy transfer), and quantifies the effect of interspur interactions on the ongoing radiochemistry. Spurs evolve under a simplified reaction-diffusion equation with parameters based on Monte Carlo simulations, and interspur interaction is quantified by calculating the expected values of interspur overlap in the target.
The model demonstrates that for any set of irradiation parameters, a minimum critical dose and dose rate are necessary to induce significant interspur interaction, and that interspur interactions correlate negatively with beam linear energy transfer at a fixed dose. The model suggests optimal beam parameters, including dose, dose rate, linear energy transfer, and pulse structure, to maximize interspur interactions. Depending on the rate of radical scavenging in the target, which limits interspur interaction, this model predicts that the irradiation parameters necessary to elicit the FLASH effect may coincide with an onset of significant interspur interactions, suggesting that interspur interaction may be the underlying mechanism of the FLASH effect.
△ Less
Submitted 25 July, 2022;
originally announced July 2022.
-
Towards the Ionizing Radiation Induced Bond Dissociation Mechanism in Oxygen, Water, Guanine and DNA Fragmentation: A Density Functional Theory Simulation
Authors:
Santosh KC,
Ramin Abolfath
Abstract:
The radiation-induced damages in bio-molecules are ubiquitous processes in radiotherapy, radio-biology and critical to space-projects. In this study we present a precise quantification of the fragmentation mechanisms of deoxyribonucleic acid (DNA) and the molecules surrounding DNA such as oxygen and water under non-equilibrium conditions using the first-principle calculations based on density func…
▽ More
The radiation-induced damages in bio-molecules are ubiquitous processes in radiotherapy, radio-biology and critical to space-projects. In this study we present a precise quantification of the fragmentation mechanisms of deoxyribonucleic acid (DNA) and the molecules surrounding DNA such as oxygen and water under non-equilibrium conditions using the first-principle calculations based on density functional theory (DFT). Our results reveal the structural stability of DNA-bases and backbone that withstand up to a combined threshold of charge and hydrogen abstraction owing to simultaneous direct and indirect ionization processes. We show the hydrogen contents of the molecules significantly control the stability in the presence of radiation. This study provides comprehensive information on the impact of the direct and indirect induced bond dissociations and DNA damage, and introduces a systematic methodology in fine-tuning of the input parameters necessary for the large-scale Monte Carlo simulations of radio-biological responses and mitigation of detrimental effects of ionizing radiation.
△ Less
Submitted 27 January, 2022;
originally announced January 2022.
-
The effect of non-ionizing excitations on the diffusion of ion species and inter-track correlations in FLASH ultra-high dose rate radiotherapy
Authors:
Ramin Abolfath,
Alexander Baikalov,
Stefan Bartzsch,
Niayesh Afshordi,
Radhe Mohan
Abstract:
We present a microscopic mechanism that accounts for the outward burst of "cold" ion species (IS) in a high-energy particle track due to coupling with ''hot" non-ion species (NIS). IS refers to radiolysis products of ionized molecules, whereas NIS refers to non-ionized excitations of molecules in a medium. The interaction is mediated by a quantized field of acoustic phonons, a channel that allows…
▽ More
We present a microscopic mechanism that accounts for the outward burst of "cold" ion species (IS) in a high-energy particle track due to coupling with ''hot" non-ion species (NIS). IS refers to radiolysis products of ionized molecules, whereas NIS refers to non-ionized excitations of molecules in a medium. The interaction is mediated by a quantized field of acoustic phonons, a channel that allows conversion of thermal energy of NIS to kinetic energy of IS, a flow of heat from the outer to the inner core of the track structure. We demonstrate the coexistence of "hot" NIS with "cold" IS in the radiation track structures right after their generation. NIS, concentrated within nano-scales volumes wrapping around IS, are the main source of intensive heat-waves and the outward burst of IS due to femto-second time scale IS-NIS coupling. By comparing the transport of IS coupled to NIS with identical configurations of non-interacting IS in thermal equilibrium at room temperature, we demonstrate that the energy gain of IS due to the surrounding hot nanoscopic volumes of NIS significantly increases their effective diffusion constants. The much higher diffusion constants predicted in the present model suggest higher inter-track chemical reaction rates at FLASH-UHDR, as well as lower intra-track reaction rates.
△ Less
Submitted 19 January, 2022;
originally announced January 2022.
-
Oxygen depletion in FLASH ultra-high-dose-rate radiotherapy: A molecular dynamics simulation
Authors:
Ramin Abolfath,
David Grosshans,
Radhe Mohan
Abstract:
We present a first-principles molecular dynamics (MD) simulation and expound upon a mechanism of oxygen depletion hypothesis to explain the mitigation of normal tissue injury observed in ultra-high-dose-rate (UHDR) FLASH radiotherapy. We simulated damage to a segment of DNA (also representing other bio-molecules such as RNA and proteins) in a simulation box filled with H$_2$O and O$_2$ molecules.…
▽ More
We present a first-principles molecular dynamics (MD) simulation and expound upon a mechanism of oxygen depletion hypothesis to explain the mitigation of normal tissue injury observed in ultra-high-dose-rate (UHDR) FLASH radiotherapy. We simulated damage to a segment of DNA (also representing other bio-molecules such as RNA and proteins) in a simulation box filled with H$_2$O and O$_2$ molecules. Attoseconds physical interactions (ionizations, electronic and vibrational excitations) were simulated by using the Monte Carlo track structure code Geant4-DNA. Immediately after ionization, {\it ab initio} Car-Parrinello molecular dynamics (CPMD) simulation was used to identify which H$_2$O and O$_2$ molecules surrounding the DNA-molecule were converted into reactive oxygen species (ROS). Subsequently, the femto- to nano-second reactions of ROS were simulated by using MD with Reactive Force Field (ReaxFF), to illustrate ROS merging into new types of non-reactive oxygen species (NROS) due to strong coupling among ROS. Time-dependent molecular simulations revealed the formation of metastable and transient spaghetti-like complexes among ROS generated at UHDR. At the higher ROS densities produced under UHDR, stranded chains (i.e., NROS) are produced, mediated through attractive electric polarity forces, hydrogen bonds, and magnetic dipole-dipole interactions among hydroxyl (\ce{^{.}OH}) radicals. NROS tend to be less mobile than cellular biomolecules as opposed to the isolated and sparsely dense ROS generated at conventional dose rates (CDR). We attribute this effect to the suppression of bio-molecular damage induced per particle track. At a given oxygen level, as the dose rate increases, the size and number of NROS chains increase, and correspondingly the population of toxic ROS components decreases.
△ Less
Submitted 1 October, 2020;
originally announced October 2020.
-
A New Approach to Modeling the Microdosimetry of Proton Therapy Beams
Authors:
Ramin Abolfath,
Yusuf Helo,
David J. Carlson,
Robert Stewart,
David Grosshans,
Radhe Mohan
Abstract:
{\bf Introduction}: To revisit the formulation of the mean chord length in microdosimetry and replace it by the particle mean free path appropriate for modelings in radio-biology. {\bf Methods}: We perform a collision-by-collision following by event-by-event Geant4 Monte Carlo simulation and calculate double-averaged stepping-length, $\langle\langle l \rangle\rangle$, for a range of target sizes f…
▽ More
{\bf Introduction}: To revisit the formulation of the mean chord length in microdosimetry and replace it by the particle mean free path appropriate for modelings in radio-biology. {\bf Methods}: We perform a collision-by-collision following by event-by-event Geant4 Monte Carlo simulation and calculate double-averaged stepping-length, $\langle\langle l \rangle\rangle$, for a range of target sizes from mm down to $μ$m and depth in water. We consider $\langle\langle l \rangle\rangle$ to represent the particle mean free path. {\bf Results}: We show that $\langle\langle l \rangle\rangle$ continuously drops as a function of depth and asymptotically saturates to a minimum value in low energies, where it exhibits a universal scaling behavior, independent of particle nominal beam energy. We correlate $\langle\langle l \rangle\rangle$ to linear density of DNA damage, complexities of initial lethal lesions and illustrate a relative difference between predictive RBEs in model calculations using mean chord length vs. the proposed mean free path. We demonstrate consistency between rapid increase in RBE within and beyond the Bragg peak and $\langle\langle l \rangle\rangle$, a decreasing function of depth. {\bf Discussion and conclusion}: An interplay between localities in imparted energy at nano-meter scale and subsequent physio-chemical processes, causalities and pathways in DNA damage requires substitution of geometrical chord length of cell nuclei by mean-free path of proton and charged particles to account for a mean distance among sequential collisions in DNA materials.
△ Less
Submitted 27 March, 2020;
originally announced March 2020.
-
A DNA damage multi-scale model for NTCP in proton and hadron therapy
Authors:
Ramin Abolfath,
Chris Peeler,
Dragan Mirkovic,
Radhe Mohan,
David Grosshans
Abstract:
{\bf Purpose}: To develop a first principle and multi-scale model for normal tissue complication probability (NTCP) as a function of dose and LET for proton and in general for particle therapy with a goal of incorporating nano-scale radio-chemical to macro-scale cell biological pathways, spanning from initial DNA damage to tissue late effects.
{\bf Methods}: The method is combination of analytic…
▽ More
{\bf Purpose}: To develop a first principle and multi-scale model for normal tissue complication probability (NTCP) as a function of dose and LET for proton and in general for particle therapy with a goal of incorporating nano-scale radio-chemical to macro-scale cell biological pathways, spanning from initial DNA damage to tissue late effects.
{\bf Methods}: The method is combination of analytical and multi-scale computational steps including (1) derivation of functional dependencies of NTCP on DNA driven cell lethality in nanometer and mapping to dose and LET in millimeter, and (2) 3D-surface fitting to Monte Carlo data set generated based on post radiation image change and gathered for a cohort of 14 pediatric patients treated by scanning beam of protons for ependymoma. We categorize voxel-based dose and LET associated with development of necrosis in NTCP.
{\bf Result}: Our model fits well the clinical data, generated for post radiation tissue toxicity and necrosis. The fitting procedure results in extraction of in-{\it vivo} radio-biological $α$-$β$ indices and their numerical values.
{\bf Discussion and conclusion}: The NTCP model, explored in this work, allows to correlate the tissue toxicities to DNA initial damage, cell lethality and the properties and qualities of radiation, dose and LET.
△ Less
Submitted 7 January, 2020;
originally announced January 2020.
-
Renormalization of radiobiological response functions by energy loss fluctuations and complexities in chromosome aberration induction: deactivation theory for proton therapy from cells to tumor control
Authors:
Ramin Abolfath,
Yusuf Helo,
Lawrence Bronk,
Alejandro Carabe,
David Grosshans,
Radhe Mohan
Abstract:
We employ a multi-scale mechanistic approach to investigate radiation induced cell toxicities and deactivation mechanisms as a function of linear energy transfer in hadron therapy. Our theoretical model consists of a system of Markov chains in microscopic and macroscopic spatio-temporal landscapes, i.e., stochastic birth-death processes of cells in millimeter-scale colonies that incorporates a coa…
▽ More
We employ a multi-scale mechanistic approach to investigate radiation induced cell toxicities and deactivation mechanisms as a function of linear energy transfer in hadron therapy. Our theoretical model consists of a system of Markov chains in microscopic and macroscopic spatio-temporal landscapes, i.e., stochastic birth-death processes of cells in millimeter-scale colonies that incorporates a coarse-grained driving force to account for microscopic radiation induced damage. The coupling, hence the driving force in this process, stems from a nano-meter scale radiation induced DNA damage that incorporates the enzymatic end-joining repair and mis-repair mechanisms. We use this model for global fitting of the high-throughput and high accuracy clonogenic cell-survival data acquired under exposure of the therapeutic scanned proton beams, the experimental design that considers $γ$-H2AX as the biological endpoint and exhibits maximum observed achievable dose and LET, beyond which the majority of the cells undergo collective biological deactivation processes. An estimate to optimal dose and LET calculated from tumor control probability by extension to $~ 10^6$ cells per $mm$-size voxels is presented. We attribute the increase in degree of complexity in chromosome aberration to variabilities in the observed biological responses as the beam linear energy transfer (LET) increases, and verify consistency of the predicted cell death probability with the in-vitro cell survival assay of approximately 100 non-small cell lung cancer (NSCLC) cells.
△ Less
Submitted 23 January, 2019;
originally announced January 2019.
-
A model for relative biological effectiveness of therapeutic proton beams based on a global fit of cell survival data
Authors:
Ramin Abolfath,
Christopher R. Peeler,
Mark Newpower,
Lawrence Bronk,
David Grosshans,
Radhe Mohan
Abstract:
We introduce an approach for global fitting of the recently published high-throughput and high accuracy clonogenic cell-survival data for therapeutic scanned proton beams. Our fitting procedure accounts for the correlation between the cell-survival, the absorbed (physical) dose and the proton linear energy transfer (LET). The fitting polynomials and constraints have been constructed upon generaliz…
▽ More
We introduce an approach for global fitting of the recently published high-throughput and high accuracy clonogenic cell-survival data for therapeutic scanned proton beams. Our fitting procedure accounts for the correlation between the cell-survival, the absorbed (physical) dose and the proton linear energy transfer (LET). The fitting polynomials and constraints have been constructed upon generalization of the microdosimetric kinetic model (gMKM) adapted to account for the low energy and high lineal-energy spectrum of the beam where the current radiobiological models may underestimate the reported relative biological effectiveness (RBE). The parameters (α,\b{eta}) of the linear-quadratic (LQ) model calculated by the presented method reveal a smooth transition from low to high LETs which is an advantage of the current method over methods previously employed to fit the same clonogenic data. Finally, the presented approach provides insight into underlying microscopic mechanisms which, with future study, may help to elucidate radiobiological responses along the Bragg curve and resolve discrepancies between experimental data and current RBE models.
△ Less
Submitted 22 December, 2018;
originally announced December 2018.
-
A molecular dynamics simulation of DNA damage induction by ionizing radiation
Authors:
Ramin M. Abolfath,
David J. Carlson,
Zhe J. Chen,
Ravinder Nath
Abstract:
We present a multi-scale simulation of early stage of DNA damages by the indirect action of hydroxyl ($^\bullet$OH) free radicals generated by electrons and protons. The computational method comprises of interfacing the Geant4-DNA Monte Carlo with the ReaxFF molecular dynamics software. A clustering method was employed to map the coordinates of $^\bullet$OH-radicals extracted from the ionization t…
▽ More
We present a multi-scale simulation of early stage of DNA damages by the indirect action of hydroxyl ($^\bullet$OH) free radicals generated by electrons and protons. The computational method comprises of interfacing the Geant4-DNA Monte Carlo with the ReaxFF molecular dynamics software. A clustering method was employed to map the coordinates of $^\bullet$OH-radicals extracted from the ionization track-structures onto nano-meter simulation voxels filled with DNA and water molecules. The molecular dynamics simulation provides the time evolution and chemical reactions in individual simulation voxels as well as the energy-landscape accounted for the DNA-$^\bullet$OH chemical reaction that is essential for the first principle enumeration of hydrogen abstractions, chemical bond breaks, and DNA-lesions induced by collection of ions in clusters less than the critical dimension which is approximately 2-3 Å. We show that the formation of broken bonds leads to DNA base and backbone damages that collectively propagate to DNA single and double strand breaks. For illustration of the methodology, we focused on particles with initial energy of 1 MeV. Our studies reveal a qualitative difference in DNA damage induced by low energy electrons and protons. Electrons mainly generate small pockets of $^\bullet$OH-radicals, randomly dispersed in the cell volume. In contrast, protons generate larger clusters along a straight-line parallel to the direction of the particle. The ratio of the total DNA double strand breaks induced by a single proton and electron track is determined to be $\approx$ 4 in the linear scaling limit. The tool developed in this work can be used in the future to investigate the relative biological effectiveness of light and heavy ions that are used in radiotherapy.
△ Less
Submitted 2 September, 2013;
originally announced September 2013.
-
Multiscale QM/MM Molecular Dynamics Study on the First Steps of Guanine-Damage by Free Hydroxyl Radicals in Solution
Authors:
Ramin M. Abolfath,
P. K. Biswas,
R. Rajnarayanam,
Thomas Brabec,
Reinhard Kodym,
Lech Papiez
Abstract:
Understanding the damage of DNA bases from hydrogen abstraction by free OH radicals is of particular importance to reveal the effect of hydroxyl radicals produced by the secondary effect of radiation. Previous studies address the problem with truncated DNA bases as ab-initio quantum simulation required to study such electronic spin dependent processes are computationally expensive. Here, for the f…
▽ More
Understanding the damage of DNA bases from hydrogen abstraction by free OH radicals is of particular importance to reveal the effect of hydroxyl radicals produced by the secondary effect of radiation. Previous studies address the problem with truncated DNA bases as ab-initio quantum simulation required to study such electronic spin dependent processes are computationally expensive. Here, for the first time, we employ a multiscale and hybrid Quantum-Mechanical-Molecular-Mechanical simulation to study the interaction of OH radicals with guanine-deoxyribose-phosphate DNA molecular unit in the presence of water where all the water molecules and the deoxyribose-phosphate fragment are treated with the simplistic classical Molecular-Mechanical scheme. Our result illustrates that the presence of water strongly alters the hydrogen-abstraction reaction as the hydrogen bonding of OH radicals with water restricts the relative orientation of the OH-radicals with respective to the the DNA base (here guanine). This results in an angular anisotropy in the chemical pathway and a lower efficiency in the hydrogen abstraction mechanisms than previously anticipated for identical system in vacuum. The method can easily be extended to single and double stranded DNA without any appreciable computational cost as these molecular units can be treated in the classical subsystem as has been demonstrated here.
△ Less
Submitted 24 February, 2012;
originally announced February 2012.
-
Dynamical Magnetic and Nuclear Polarization in Complex Spin Systems: Semi-magnetic II-VI Quantum Dots
Authors:
Ramin M. Abolfath,
Anna Trojnar,
Bahman Roostaei,
Thomas Brabec,
Pawel Hawrylak
Abstract:
Dynamical magnetic and nuclear polarization in complex spin systems is discussed on the example of transfer of spin from exciton to the central spin of magnetic impurity in a quantum dot in the presence of a finite number of nuclear spins. The exciton is described in terms of the electron and heavy hole spins interacting via exchange interaction with magnetic impurity, via hypeprfine interaction w…
▽ More
Dynamical magnetic and nuclear polarization in complex spin systems is discussed on the example of transfer of spin from exciton to the central spin of magnetic impurity in a quantum dot in the presence of a finite number of nuclear spins. The exciton is described in terms of the electron and heavy hole spins interacting via exchange interaction with magnetic impurity, via hypeprfine interaction with a finite number of nuclear spins and via dipole interaction with photons. The time-evolution of the exciton, magnetic impurity and nuclear spins is calculated exactly between quantum jumps corresponding to exciton radiative recombination. The collapse of the wavefunction and the refilling of the quantum dot with new spin polarized exciton is shown to lead to build up of magnetization of the magnetic impurity as well as nuclear spin polarization. The competition between electron spin transfer to magnetic impurity and to nuclear spins simultaneous with the creation of dark excitons is elucidated. The technique presented here opens up the possibility of studying optically induced Dynamical Magnetic and Nuclear Polarization in Complex Spin Systems.
△ Less
Submitted 16 May, 2012; v1 submitted 23 February, 2012;
originally announced February 2012.
-
Computational studies for reduced graphene oxide in hydrogen-rich environment
Authors:
Ramin M. Abolfath,
Kyeongjae Cho
Abstract:
We employ molecular dynamic simulations to study the reduction process of graphene-oxide (GO) in a chemically active environment enriched with hydrogen. We examine the concentration and pressure of hydrogen gas as a function of temperature in which abstraction of oxygen is possible with minimum damage to C-sp$^2$ bonds hence preserving the integrity of the graphene sheet. Through these studies we…
▽ More
We employ molecular dynamic simulations to study the reduction process of graphene-oxide (GO) in a chemically active environment enriched with hydrogen. We examine the concentration and pressure of hydrogen gas as a function of temperature in which abstraction of oxygen is possible with minimum damage to C-sp$^2$ bonds hence preserving the integrity of the graphene sheet. Through these studies we find chemical pathways that demonstrate beneficiary mechanisms for the quality of graphene including formation of water as well as suppression of carbonyl pair holes in favor of hydroxyl and epoxy formation facilitated by hydrogen gas in the environment.
△ Less
Submitted 19 December, 2011;
originally announced December 2011.
-
Spin textures in strongly coupled electron spin and magnetic or nuclear spin systems in quantum dots
Authors:
Ramin M. Abolfath,
Marek Korkusinski,
Thomas Brabec,
Pawel Hawrylak
Abstract:
Controlling electron spins strongly coupled to magnetic and nuclear spins in solid state systems is an important challenege in the field of spintronics and quantum computation. We show here that electron droplets with no net spin in semiconductor quantum dots strongly coupled with magnetic ion/nuclear spin systems break down at low temperature and form a non-trivial anti-ferromagnetic spatially or…
▽ More
Controlling electron spins strongly coupled to magnetic and nuclear spins in solid state systems is an important challenege in the field of spintronics and quantum computation. We show here that electron droplets with no net spin in semiconductor quantum dots strongly coupled with magnetic ion/nuclear spin systems break down at low temperature and form a non-trivial anti-ferromagnetic spatially ordered spin-texture of magneto-polarons. The spatially ordered combined electron-magnetic ion spin-texture, associated with spontaneous symmetry-breaking in the parity of electronic charge density and magnetization of magnetic ions, emerge from both ab-initio density functional approach to the electronic system coupled with mean-field approximation for the magnetic/nuclear spin system and fully mircoscopic exact diagonalization of small systems. The predicted phase diagram determines the critical temperature as a function of coupling strength and identifies possible phases of the strongly coupled spin system. This prediction may arrest fluctuations in spin system and open the way to control, manipulate and prepare magnetic and nuclear spin ensembles in semiconductor nanostructures.
△ Less
Submitted 17 April, 2012; v1 submitted 18 December, 2011;
originally announced December 2011.
-
Spin-blockade effect and coherent control of DNA-damage by free radicals: a proposal on bio-spintronics
Authors:
Ramin M. Abolfath,
Thomas Brabec
Abstract:
Coherent control of OH-free radicals interacting with the spin-triplet state of a DNA molecule is investigated. A model Hamiltonian for molecular spin singlet-triplet resonance is developed. We illustrate that the spin-triplet state in DNA molecules can be efficiently populated, as the spin-injection rate can be tuned to be orders of magnitudes greater than the decay rate due to small spin-orbit c…
▽ More
Coherent control of OH-free radicals interacting with the spin-triplet state of a DNA molecule is investigated. A model Hamiltonian for molecular spin singlet-triplet resonance is developed. We illustrate that the spin-triplet state in DNA molecules can be efficiently populated, as the spin-injection rate can be tuned to be orders of magnitudes greater than the decay rate due to small spin-orbit coupling in organic molecules. Owing to the nano-second life-time of OH free radicals, a non-equilibrium free energy barrier induced by the injected spin triplet state that lasts approximately longer than one-micro second in room temperature can efficiently block the initial Hydrogen abstraction and DNA damage. For a direct demonstration of the spin-blockade effect, a molecular simulation based on an {\em ab-initio} Car-Parrinello molecular dynamics is deployed.
△ Less
Submitted 12 December, 2011;
originally announced December 2011.
-
Reactive Molecular Dynamics study on the first steps of DNA-damage by free hydroxyl radicals
Authors:
Ramin M. Abolfath,
A. C. T. van Duin,
Thomas Brabec
Abstract:
We employ a large scale molecular simulation based on bond-order ReaxFF to simulate the chemical reaction and study the damage to a large fragment of DNA-molecule in the solution by ionizing radiation. We illustrate that the randomly distributed clusters of diatomic OH-radicals that are primary products of megavoltage ionizing radiation in water-based systems are the main source of hydrogen-abstra…
▽ More
We employ a large scale molecular simulation based on bond-order ReaxFF to simulate the chemical reaction and study the damage to a large fragment of DNA-molecule in the solution by ionizing radiation. We illustrate that the randomly distributed clusters of diatomic OH-radicals that are primary products of megavoltage ionizing radiation in water-based systems are the main source of hydrogen-abstraction as well as formation of carbonyl- and hydroxyl-groups in the sugar-moiety that create holes in the sugar-rings. These holes grow up slowly between DNA-bases and DNA-backbone and the damage collectively propagate to DNA single and double strand break.
△ Less
Submitted 12 December, 2011;
originally announced December 2011.
-
Performance engineering of semiconductor spin qubit systems
Authors:
Ramin M. Abolfath,
Thomas Brabec
Abstract:
The performance of a quantum computation system is investigated, with qubits represented by magnetic impurities in coupled quantum dots filled with two electrons. Magnetic impurities are electrically manipulated by electrons. The dominant noise source is the electron mediated indirect coupling between magnetic impurities and host spin bath. As a result of the electron mediated coupling, both noise…
▽ More
The performance of a quantum computation system is investigated, with qubits represented by magnetic impurities in coupled quantum dots filled with two electrons. Magnetic impurities are electrically manipulated by electrons. The dominant noise source is the electron mediated indirect coupling between magnetic impurities and host spin bath. As a result of the electron mediated coupling, both noise properties and the time needed for elementary gate operations, depend on controllable system parameters, such as size and geometry of the quantum dot, and external electric and magnetic fields. We find that the maximum number of quantum operations per coherence time for magnetic impurities increases as electron spin singlet triplet energy gap decreases. The advantage of magnetic impurities over electrons for weak coupling and large magnetic fields will be illustrated.
△ Less
Submitted 28 July, 2010; v1 submitted 14 January, 2010;
originally announced January 2010.
-
Protection Strategy of Sensitive Body Organs in Radiation Therapy
Authors:
Ramin M. Abolfath,
Lech Papiez
Abstract:
In this paper, we investigate protection strategies of sensitive body anatomy against the irradiation to the cancerous moving tumors in intensity modulated radiation therapy. Inspired by optimization techniques developed in statistical physics and dynamical systems, we deploy a method based on variational principles and formulate an efficient genetic algorithm which enable us to search for globa…
▽ More
In this paper, we investigate protection strategies of sensitive body anatomy against the irradiation to the cancerous moving tumors in intensity modulated radiation therapy. Inspired by optimization techniques developed in statistical physics and dynamical systems, we deploy a method based on variational principles and formulate an efficient genetic algorithm which enable us to search for global minima in a complex landscape of irradiation dose delivered to the radiosensitive organs at risk. We take advantage of the internal motion of body anatomy during radiation therapy to reduce the unintentional delivery of the radiation to sensitive organs. We show that the accurate optimization of the control parameters, compare to the conventional IMRT and widely used delivery based on static anatomy assumption, leads to a significant reduction of the dose delivered to the organs at risk.
△ Less
Submitted 28 July, 2009;
originally announced July 2009.
-
Para-ortho transition of artificial H$_2$ molecule in lateral quantum dots doped with magnetic impurities
Authors:
Ramin M. Abolfath
Abstract:
We present the magnetic phase diagram of artificial H$_2$ molecule in lateral quantum dots doped with magnetic impurities as a function of external magnetic field and plunger gate voltage. The onset of Mn-Mn antiferromagnetic-ferromagnetic transition follows the electron spin singlet-triplet transition. We deploy a configuration-interaction method to exactly diagonalize the electron-Mn Hamiltoni…
▽ More
We present the magnetic phase diagram of artificial H$_2$ molecule in lateral quantum dots doped with magnetic impurities as a function of external magnetic field and plunger gate voltage. The onset of Mn-Mn antiferromagnetic-ferromagnetic transition follows the electron spin singlet-triplet transition. We deploy a configuration-interaction method to exactly diagonalize the electron-Mn Hamiltonian and map it to an effective Mn-Mn Heisenberg Hamiltonian. We find that Mn-Mn exchange coupling can be described by RKKY-interaction/magnetic-polaron in weak/strong electron-Mn coupling at low/high magnetic fields.
△ Less
Submitted 27 October, 2009; v1 submitted 27 July, 2009;
originally announced July 2009.
-
Optical control of DNA-base radio-sensitivity
Authors:
Ramin M. Abolfath
Abstract:
{\bf Purpose}: Manipulation of the radio-sensitivity of the nucleotide-base driven by the spin blockade mechanism of diffusive free radicals against ionizing radiation.
{\bf Materials and methods}: We theoretically propose a mechanism which uses the simultaneous application of circularly polarized light and an external magnetic field to control the polarization of the free radicals and create…
▽ More
{\bf Purpose}: Manipulation of the radio-sensitivity of the nucleotide-base driven by the spin blockade mechanism of diffusive free radicals against ionizing radiation.
{\bf Materials and methods}: We theoretically propose a mechanism which uses the simultaneous application of circularly polarized light and an external magnetic field to control the polarization of the free radicals and create S=1 electron-hole spin excitations (excitons) on nucleotide-base. We deploy an ab-initio molecular dynamics model to calculate the characteristic parameters of the light needed for optical transitions.
{\bf Results}: As a specific example, we present the numerical results calculated for a Guanine, in the presence of an OH free radical. To increase the radio-resistivity of this system, a blue light source for the optical pumping and induction of excitons on guanine can be used.
{\bf Conclusions}: The effect of spin-injection on the formation of a free energy barrier in diffusion controlled chemical reaction pathways leads to the control of radiation-induced base damage. The proposed method allows us to manipulate and partially suppress the damage induced by ionizing radiation.
△ Less
Submitted 27 July, 2009;
originally announced July 2009.
-
Electronic States of Magnetic Quantum Dots
Authors:
Ramin M. Abolfath,
Pawel Hawrylak,
Igor Zutic
Abstract:
We study quantum states of electrons in magnetically doped quantum dots as a function of exchange coupling between electron and impurity spins, the strength of Coulomb interaction, confining potential, and the number of electrons. The magnetic phase diagram of quantum dots, doped with a large number of magnetic Mn impurities, can be described by the energy gap in the spectrum of electrons and th…
▽ More
We study quantum states of electrons in magnetically doped quantum dots as a function of exchange coupling between electron and impurity spins, the strength of Coulomb interaction, confining potential, and the number of electrons. The magnetic phase diagram of quantum dots, doped with a large number of magnetic Mn impurities, can be described by the energy gap in the spectrum of electrons and the mean field electron-Mn exchange coupling. A competition between these two parameters leads to a transition between spin-unpolarized and spin-polarized states, in the absence of applied magnetic field. Tuning the energy gap by electrostatic control of nonparabolicity of the confining potential can enable control of magnetization even at the fixed number of electrons. We illustrate our findings by directly comparing Mn-doped quantum dots with parabolic and Gaussian confining potential.
△ Less
Submitted 12 October, 2007;
originally announced October 2007.
-
Piezomagnetic Quantum Dots
Authors:
Ramin M. Abolfath,
Andre Petukhov,
Igor Zutic
Abstract:
We study the influence of deformations on magnetic ordering in quantum dots doped with magnetic impurities. The reduction of symmetry and the associated deformation from circular to elliptical quantum confinement lead to the formation of piezomagnetic quantum dots. The strength of elliptical deformation can be controlled by the gate voltage to change the magnitude of magnetization, at a fixed nu…
▽ More
We study the influence of deformations on magnetic ordering in quantum dots doped with magnetic impurities. The reduction of symmetry and the associated deformation from circular to elliptical quantum confinement lead to the formation of piezomagnetic quantum dots. The strength of elliptical deformation can be controlled by the gate voltage to change the magnitude of magnetization, at a fixed number of carriers and in the absence of applied magnetic field. We reveal a reentrant magnetic ordering with the increase of elliptical deformation and suggest that the piezomagnetic quantum dots can be used as nanoscale magnetic switches.
△ Less
Submitted 18 July, 2007;
originally announced July 2007.
-
Spin transitions induced by a magnetic field in quantum dot molecules
Authors:
Ramin M. Abolfath,
Pawel Hawrylak
Abstract:
We present a theoretical study of magnetic field driven spin transitions of electrons in coupled lateral quantum dot molecules. A detailed numerical study of spin phases of artificial molecules composed of two laterally coupled quantum dots with N=8 electrons is presented as a function of magnetic field, Zeeman energy, and the detuning using real space Hartree-Fock Configuration Interaction (HF-…
▽ More
We present a theoretical study of magnetic field driven spin transitions of electrons in coupled lateral quantum dot molecules. A detailed numerical study of spin phases of artificial molecules composed of two laterally coupled quantum dots with N=8 electrons is presented as a function of magnetic field, Zeeman energy, and the detuning using real space Hartree-Fock Configuration Interaction (HF-CI) technique. A microscopic picture of quantum Hall ferromagnetic phases corresponding to zero and full spin polarization at filling factors $ν=2$ and $ν=1$, and ferrimagnetic phases resulting from coupling of the two dots, is presented.
△ Less
Submitted 30 April, 2008; v1 submitted 20 December, 2006;
originally announced December 2006.
-
Tailoring Magnetism in Quantum Dots
Authors:
Ramin M. Abolfath,
Pawel Hawrylak,
Igor Zutic
Abstract:
We study magnetism in magnetically doped quantum dots as a function of confining potential, particle numbers, temperature, and strength of Coulomb interactions. We explore possibility of tailoring magnetism by controlling the electron-electron Coulomb interaction, without changing the number of particles. The interplay of strong Coulomb interactions and quantum confinement leads to enhanced inho…
▽ More
We study magnetism in magnetically doped quantum dots as a function of confining potential, particle numbers, temperature, and strength of Coulomb interactions. We explore possibility of tailoring magnetism by controlling the electron-electron Coulomb interaction, without changing the number of particles. The interplay of strong Coulomb interactions and quantum confinement leads to enhanced inhomogeneous magnetization which persist at higher temperatures than in the non-interacting case. The temperature of the onset of magnetization can be controlled by changing the number of particles as well as by modifying the quantum confinement and the strength of Coulomb interactions. We predict a series of electronic spin transitions which arise from the competition between the many-body gap and magnetic thermal fluctuations.
△ Less
Submitted 19 December, 2006;
originally announced December 2006.
-
Quantum Hall Ferrimagnetism in lateral quantum dot molecules
Authors:
Ramin M. Abolfath,
Pawel Hawrylak
Abstract:
We demonstrate the existance of ferrimagnetic and ferromagnetic phases in a spin phase diagram of coupled lateral quantum dot molecules in the quantum Hall regime. The spin phase diagram is determined from Hartree-Fock Configuration Interaction method as a function of electron numbers N, magnetic field B, Zeeman energy, and tunneling barrier height. The quantum Hall ferrimagnetic phase correspon…
▽ More
We demonstrate the existance of ferrimagnetic and ferromagnetic phases in a spin phase diagram of coupled lateral quantum dot molecules in the quantum Hall regime. The spin phase diagram is determined from Hartree-Fock Configuration Interaction method as a function of electron numbers N, magnetic field B, Zeeman energy, and tunneling barrier height. The quantum Hall ferrimagnetic phase corresponds to spatially imbalanced spin droplets resulting from strong inter-dot coupling of identical dots. The quantum Hall ferromagnetic phases correspond to ferromagnetic coupling of spin polarization at filling factors between $ν=2$ and $ν=1$.
△ Less
Submitted 25 July, 2006;
originally announced July 2006.
-
Quantum Hall droplets in coupled lateral quantum dots
Authors:
Ramin M. Abolfath,
Pawel Hawrylak,
Michel Pioro-Ladriere,
Andy Sachrajda
Abstract:
We present a theory and Coulomb and Spin Blockade spectroscopy experiments on quantum Hall droplets with controlled electron numbers (N1,N2) in laterally coupled gated quantum dots. The theory is based on the configuration interaction method (CI) coupled with the unrestricted Hartree-Fock (URHF) basis. It allows us to calculate the magnetic field evolution of ground and excited states of coupled…
▽ More
We present a theory and Coulomb and Spin Blockade spectroscopy experiments on quantum Hall droplets with controlled electron numbers (N1,N2) in laterally coupled gated quantum dots. The theory is based on the configuration interaction method (CI) coupled with the unrestricted Hartree-Fock (URHF) basis. It allows us to calculate the magnetic field evolution of ground and excited states of coupled quantum dots with large electron numbers. The method is applied to the spin transitions in the (5,5) droplet. Preliminary experimental results demonstrate the creation of the (5,5) droplet and its Spin Blockade spectra.
△ Less
Submitted 27 March, 2006;
originally announced March 2006.
-
Real Space Hartree-Fock Configuration Interaction Method For Complex Lateral Quantum Dot Molecules
Authors:
Ramin M. Abolfath,
Pawel Hawrylak
Abstract:
We present unrestricted Hartree Fock method coupled with configuration interaction (CI) method (URHF-CI) suitable for the calculation of ground and excited states of large number of electrons localized by complex gate potentials in quasi-two-dimensional quantum dot molecules. The method employs real space finite difference method, incorporating strong magnetic field, for the calculating of singl…
▽ More
We present unrestricted Hartree Fock method coupled with configuration interaction (CI) method (URHF-CI) suitable for the calculation of ground and excited states of large number of electrons localized by complex gate potentials in quasi-two-dimensional quantum dot molecules. The method employs real space finite difference method, incorporating strong magnetic field, for the calculating of single particle states. The Hartree-Fock method is employed for the calculation of direct and exchange interaction contribution to the ground state energy. The effects of correlations are included in energies and directly in the many-particle wavefunctions via configuration interaction (CI) method using a limited set of excitations above the Fermi level. The URHF-CI method and its performance are illustrated on the example of ten electrons confined in a two-dimensional quantum dot molecule.
△ Less
Submitted 27 July, 2006; v1 submitted 30 November, 2005;
originally announced November 2005.
-
Theory of a two-level artificial molecule in laterally coupled quantum Hall droplets
Authors:
Ramin M. Abolfath,
W. Dybalski,
Pawel Hawrylak
Abstract:
We present a theory of laterally coupled quantum Hall droplets with electron numbers (N1,N2) at filling factor $ν=2$. We show that the edge states of each droplet are tunnel coupled and form a two-level artificial molecule. By populating the edge states with one electron each a two electron molecule is formed. We predict the singlet-triplet transitions of the effective two-electron molecule as a…
▽ More
We present a theory of laterally coupled quantum Hall droplets with electron numbers (N1,N2) at filling factor $ν=2$. We show that the edge states of each droplet are tunnel coupled and form a two-level artificial molecule. By populating the edge states with one electron each a two electron molecule is formed. We predict the singlet-triplet transitions of the effective two-electron molecule as a function of the magnetic field, the number of electrons, and confining potential using the configuration interaction method (CI) coupled with the unrestricted Hartree-Fock (URHF) basis. In addition to the singlet-triplet transitions of a 2 electron molecule involving edge states, triplet transitions involving transfer of electrons to the center of individual dots exist for $(N1 \geq 5, N2 \geq 5)$.
△ Less
Submitted 29 March, 2006; v1 submitted 22 September, 2005;
originally announced September 2005.
-
Charge Sensing of an Artificial H2+ Molecule
Authors:
M. Pioro-Ladriere,
M. R. Abolfath,
P. Zawadzki,
J. Lapointe,
S. A. Studenikin,
A. S. Sachrajda,
P. Hawrylak
Abstract:
We report charge detection studies of a lateral double quantum dot with controllable charge states and tunable tunnel coupling. Using an integrated electrometer, we characterize the equilibrium state of a single electron trapped in the doubled-dot (artificial H2+ molecule) by measuring the average occupation of one dot. We present a model where the electrostatic coupling between the molecule and…
▽ More
We report charge detection studies of a lateral double quantum dot with controllable charge states and tunable tunnel coupling. Using an integrated electrometer, we characterize the equilibrium state of a single electron trapped in the doubled-dot (artificial H2+ molecule) by measuring the average occupation of one dot. We present a model where the electrostatic coupling between the molecule and the sensor is taken into account explicitly. From the measurements, we extract the temperature of the isolated electron and the tunnel coupling energy. It is found that this coupling can be tuned between 0 and 60 micro electron-volt in our device.
△ Less
Submitted 28 July, 2005; v1 submitted 1 April, 2005;
originally announced April 2005.