Environmental two-level systems (TLS) have been identified as significant decoherence sources in Josephson junction (JJ) based circuits. For such quantum devices to be functional, the removal or control of the TLS is a necessity. Understanding the microscopic origins of the 'strongly coupled' TLS type is one current path of investigation to that end. The delocalized oxygen model suggests the atomic position of an oxygen atom is spatially delocalized in the oxide forming the JJ barrier. In this report we extend this model from its previous 2+1D construction to a complete 3D description using a Wick-rotated time-dependent Schrodinger equation to solve for time-independent solutions in three dimensions. We compute experimentally observable parameters for phase qubits and compare the results to the 2+1D framework. We devise a Voronoi classification scheme to investigate oxygen atoms delocalizing within strained and non-strained crystalline lattices, as well as realistic atomic positions a JJ amorphous tunnel barrier constructed in previous density functional studies.
Generalized collision operator for fast electrons interacting with partially ionized impurities
Linnea Hesslow, Ola Embréus, Mathias Hoppe, Timothy C. DuBois, Gergely Papp, Martin Rahm, Tünde Fülöp
Accurate modelling of the interaction between fast electrons and partially ionized atoms is important for evaluating tokamak disruption mitigation schemes based on material injection. This requires accounting for the effect of screening of the impurity nuclei by the cloud of bound electrons. In this paper, we generalize the Fokker–Planck operator in a fully ionized plasma by accounting for the effect of screening. We detail the derivation of this generalized operator, and calculate the effective ion length scales, needed in the components of the collision operator, for a number of ion species commonly appearing in fusion experiments. We show that for high electric fields, the secondary runaway growth rate can be substantially larger than in a fully ionized plasma with the same effective charge, although the growth rate is significantly reduced at near-critical electric fields. Furthermore, by comparison with the Boltzmann collision operator, we show that the Fokker–Planck formalism is accurate even for large impurity content.
Proton acceleration by a pair of successive ultraintense femtosecond laser pulses
Julien Ferri, Lovisa Senje, Malay Dalui, Kristoffer Svensson, Bastian Aurand, Martin Hansson, Anders Persson, Olle Lundh, Claes-Göran Wahlström, Laurent Gremillet, Evangelos Siminos, Timothy C. DuBois, Longqing Yi, Joana Martins and Tünde Fülöp
We investigate the target normal sheath acceleration of protons in thin aluminum targets irradiated at relativistic intensity by two time-separated ultrashort (35 fs) laser pulses. For identical laser pulses and target thicknesses of 3 and 6 μm, we observe experimentally that the second pulse boosts the maximum energy and charge of the proton beam produced by the first pulse for time delays below ∼0.6−1 ps. By using two-dimensional particle-in-cell simulations we examine the variation of the proton energy spectra with respect to the time-delay between the two pulses. We demonstrate that the expansion of the target front surface caused by the first pulse significantly enhances the hot-electron generation by the second pulse arriving after a few hundreds of fs time delay. This enhancement, however, does not suffice to further accelerate the fastest protons driven by the first pulse once three-dimensional quenching effects have set in. This implies a limit to the maximum time delay that leads to proton energy enhancement, which we theoretically determine.
Origins of plateau formation in ion energy spectra under target normal sheath acceleration
Timothy C. DuBois, Evangelos Siminos, Julien Ferri, Laurent Gremillet and Tünde Fülöp
Target normal sheath acceleration (TNSA) is a method employed in laser–matter interaction experiments to accelerate light ions (usually protons). Laser setups with durations of a few 10 fs and relatively low intensity contrasts observe plateau regions in their ion energy spectra when shooting on thin foil targets with thicknesses of order 10 µm. In this paper we identify a mechanism which explains this phenomenon using one dimensional particle-in-cell simulations. Fast electrons generated from the laser interaction recirculate back and forth through the target, giving rise to time-oscillating charge and current densities at the target backside. Periodic decreases in the electron density lead to transient disruptions of the TNSA sheath field: peaks in the ion spectra form as a result, which are then spread in energy from a modified potential driven by further electron recirculation. The ratio between the laser pulse duration and the recirculation period (dependent on the target thickness, including the portion of the pre-plasma which is denser than the critical density) determines if a plateau forms in the energy spectra.
Effect of partially-screened nuclei on fast-electron dynamics
Linnea Hesslow, Ola Embréus, Adam Stahl, Timothy C. DuBois, Gergely Papp, Sarah L. Newton, Tünde Fülöp
We analyze the dynamics of fast electrons in plasmas containing partially ionized impurity atoms, where the screening effect of bound electrons must be included. We derive analytical expressions for the deflection and slowing-down frequencies, and show that they are increased significantly compared to the results obtained with complete screening, already at subrelativistic electron energies. Furthermore, we show that the modifications to the deflection and slowing down frequencies are of equal importance in describing the runaway current evolution. Our results greatly affect fast-electron dynamics and have important implications, e.g., for the efficacy of mitigation strategies for runaway electrons in tokamak devices, and energy loss during relativistic breakdown in atmospheric discharges.
Relativistic Vlasov-Maxwell modelling using finite volumes and adaptive mesh refinement
Benjamin Svedung Wettervik, Timothy C. DuBois, Evangelos Siminos, Tünde Fülöp
The dynamics of collisionless plasmas can be modelled by the Vlasov-Maxwell system of equations. An Eulerian approach is needed to accurately describe processes that are governed by high energy tails in the distribution function, but is of limited efficiency for high dimensional problems. The use of an adaptive mesh can reduce the scaling of the computational cost with the dimension of the problem. Here, we present a relativistic Eulerian Vlasov-Maxwell solver with block-structured adaptive mesh refinement in one spatial and one momentum dimension. The discretization of the Vlasov equation is based on a high-order finite volume method. A flux corrected transport algorithm is applied to limit spurious oscillations and ensure the physical character of the distribution function. We demonstrate a speed-up by a factor of 7x in a typical scenario involving laser-plasma interaction with an underdense plasma due to the use of an adaptive mesh.
Vlasov modelling of laser-driven collisionless shock acceleration of protons
Benjamin Svedung Wettervik, Timothy C. DuBois, Tünde Fülöp
Ion acceleration due to the interaction between a short high-intensity laser pulse and a moderately overdense plasma target is studied using Eulerian Vlasov–Maxwell simulations. The effects of variations in the plasma density profile and laser pulse parameters are investigated, and the interplay of collisionless shock and target normal sheath acceleration is analyzed. It is shown that the use of a layered-target with a combination of light and heavy ions, on the front and rear side, respectively, yields a strong quasi-static sheath-field on the rear side of the heavy-ion part of the target. This sheath-field increases the energy of the shock-accelerated ions while preserving their mono-energeticity.
Constructing ab initio models of ultra-thin Al-AlOx-Al barriers
Timothy C. DuBois, Martin J. Cyster, George Opletal, Salvy P. Russo, Jared H. Cole
The microscopic structure of ultra-thin oxide barriers often plays a major role in modern nano-electronic devices. In the case of superconducting electronic circuits, their operation depends on the electrical non-linearity provided by one or more such oxide layers in the form of ultra-thin tunnel barriers (also known as Josephson junctions). Currently available fabrication techniques manufacture an amorphous oxide barrier, which is attributed as a major noise source within the device. The nature of this noise is currently an open question and requires both experimental and theoretical investigation. Here, we present a methodology for constructing atomic scale computational models of Josephson junctions using a combination of molecular mechanics, empirical and ab initio methods. These junctions consist of ultra-thin amorphous aluminium-oxide layers sandwiched between crystalline aluminium. The stability and structure of these barriers as a function of density and stoichiometry are investigated, which we compare to experimentally observed parameters.
The dynamics of particle transport under the influence of localized high energy anomalies (explosions) is a complicated phenomenon dependent on many physical parameters of both the particle and the medium it resides in. Here we present a conceptual model that establishes simple scaling laws for particle dispersion in relation to the energy released in a blast, properties of the medium, physical properties of particles, and their initial position away from a blast epicenter. These dependencies are validated against numerical simulations and we discuss predictions of the model which can be validated experimentally. Other applications and extensions to the framework are also considered.
Atomic delocalization as a microscopic origin of two-level defects in Josephson junctions
Timothy C. DuBois, Salvy P. Russo and Jared H. Cole<Paste>
Identifying the microscopic origins of decoherence sources prevalent in Josephson junction (JJ) based circuits is central to their use as functional quantum devices. Focussing on so called 'strongly coupled' two-level defects, we construct a theoretical model using the atomic position of the oxygen which is spatially delocalized in the oxide forming the JJ barrier. Using this model, we investigate which atomic configurations give rise to two-level behaviour of the type seen in experiments. We compute experimentally observable parameters for phase qubits and examine defect response under the effects of applied electric field and strain.
Tracer Dispersion in the Turbulent Convective Layer
Alex Skvortsov, Milan Jamriska and Timothy C. DuBois
Experimental results for passive tracer dispersion in the turbulent surface layer under convective conditions are presented. In this case, the dispersion of tracer particles is determined by the interplay of two mechanisms: buoyancy and advection. In the atmospheric surface layer under unstable stratification the buoyancy mechanism dominates when the distance from the ground is greater than the Monin–Obukhov length, resulting in a different exponent in the scaling law for relative separation of Lagrangian particles (a deviation from the celebrated Richardson’s law). This conclusion is supported by atmospheric observations. Exit-time statistics and a probability density function of concentration increments derived from a previously published experimental dataset demonstrate a noticeable difference between tracer dispersion in the convective and neutrally stratified surface layers.
Delocalized Oxygen as the Origin of Two-Level Defects in Josephson Junctions
Timothy C. DuBois, Manolo C. Per, Salvy P. Russo and Jared H. Cole
One of the key problems facing superconducting qubits and other Josephson junction devices is the decohering effects of bistable material defects. Although a variety of phenomenological models exist, the true microscopic origin of these defects remains elusive. For the first time we show that these defects may arise from delocalization of the atomic position of the oxygen in the oxide forming the Josephson junction barrier. Using a microscopic model, we compute experimentally observable parameters for phase qubits. Such defects are charge neutral but have nonzero response to both applied electric field and strain. This may explain the observed long coherence time of two-level defects in the presence of charge noise, while still coupling to the junction electric field and substrate phonons.
Statistical characterisation of bio-aerosol background in an urban environment
Milan Jamriska, Timothy C. DuBois and Alex Skvortsov
In this paper we statistically characterise the bio-aerosol background in an urban environment. To do this we measure concentration levels of naturally occurring microbiological material in the atmosphere over several days. Naturally occurring bio-aerosols can be considered as noise, as they mask the presence of signals coming from biological material of interest (such as an intentionally released biological agent). Analysis of this 'biobackground' was undertaken in the 1—10 µm size range and a 3—9% contribution was found to be biological in origin e values which are in good agreement with other studies reported in the literature. A model based on the physics of turbulent mixing and dispersion was developed and validated against this analysis. The Gamma distribution (the basis of our model) is shown to comply with the scaling laws of the concentration moments of our data, which enables us to universally characterise both biological and non-biological material in the atmosphere. An application of this model is proposed to build a framework for the development of novel algorithms for bio-aerosol detection and rapid characterisation.
Scaling laws of passive tracer dispersion in the turbulent surface layer
Alex Skvortsov, Milan Jamriska and Timothy C. DuBois
Experimental results for passive tracer dispersion in the turbulent surface layer under stable conditions are presented. In this case, the dispersion of tracer particles is determined by the interplay of three mechanisms: relative dispersion (celebrated Richardson’s mechanism), shear dispersion (particle separation due to variation of the mean velocity field) and specific surface-layer dispersion (induced by the gradient of the energy dissipation rate in the turbulent surface layer). The latter mechanism results in the rather slow (ballistic) law for the mean squared particle separation. Based on a simplified Langevin equation for particle separation we found that the ballistic regime always dominates at large times. This conclusion is supported by our extensive atmospheric observations. Exit-time statistics are derived from the experimental data set and show a reasonable match with the simple dimensional asymptotes for different mechanisms of tracer dispersion, as well as predictions of the multifractal model and experimental data from other sources.
On a Delocalised Oxygen Model of Two-Level System Defects in Josephson Junctions
Timothy C. DuBois
Doctor of Philosophy (Physics), RMIT University (2015)
Constructing circuits which invoke the tunnelling effect of Josephson junctions in the superconducting regime manifest controllable, quantum properties. The development of a quantum computer employing these phenomena is just one of the many advantages to be gained from building such devices.
However, the use of quantum bits for computation is dependent on the ability to operate them in an essentially isolated environment. The phenomenon of decoherence refers to the instability of a quantum state of a system when it interacts with the surrounding environment. Superconducting qubits are sensitive to decoherence mechanisms within the readout leads connecting to the device, and more importantly from the materials from which they are manufactured. Removal or control of these imperfections is required before quantum computers using superconducting circuit architecture can be realised.
One identifiable noise source is the so called 'strongly coupled' two-level system (TLS). Comparable resonance frequencies to the qubit; strong coupling strengths and decoherence times long enough to allow coherent oscillations between the qubit and TLS have been experimentally measured. The premise of this thesis is that positional anharmonicity of oxygen atoms arises within the AlOx barrier of the Josephson junction solely due to its amorphous construction. This ansatz allows the existence of various spatial configurations throughout the layer, causing unique TLS properties based solely on atomic positions and rotation in relation to the external electric field.
To validate this conjecture, Josephson junction models are constructed using a hybrid ab initio and molecular mechanics approach, with various stoichiometry and density properties to reflect experimental observations of the barrier. The resultant atomic positions provide input conditions through a Voronoi classification scheme to a framework describing an oxygen atom that has the capacity to become spatially delocalised as bonds perturb away from a crystalline structure. A direct diagonalisation method is developed for low dimensional descriptions and a Wick-rotated time-dependent Schrödinger equation implementation is used for three dimensional investigations.
Calculated properties are compared to many current experimental strongly coupled TLS measurements in phase qubits, which shows the models' capacity to explain how an oxygen atom can generate a large dipole and appropriate ground to first excited state splitting values expected of a TLS by merely migrating from its preferred lattice position.
On Molecular and DFT Modelling of Titanium Defects in α-Quartz With a Focus on Properties Relevant to Geo-Thermometers
Timothy C. DuBois
Bachelor of Science (Applied Sciences) (Honours), RMIT University (2010)
It's pretty terrible. Contact me directly if you really want me to send you a copy.
Geodynamic processes are extremely difficult to investigate with current technologies. Rheology in the lithosphere and asthenosphere is highly nonlinear and elasticity causes complications near the surface when the ambient temperature is less than approximately 600°C. The current contention between Geologists and Geophysicists is that the Earth’s main geodynamic process is that of Mantle Convection; exerting plate tectonics to organise the system. Considering that plate tectonics is only a kinematic description of observations though, and not a fully dynamic description of plate motion, this and other approximations and assumptions continue to build upon one another yielding incomplete and inconsistent models. Geothermobarometry is the science of measuring temperature and pressure histories of intrusive igneous or metamorphic rock in an attempt to understand these processes in a limited fashion.
The recent work complied on a Titanium-in-Quartz geothermometer by Wark & Watson [Contrib. Mineral. Petrol. 152 (2006)] concluded with an empirical formula generated from the cathodoluminescence of a synthesised Quartz system in the presence of Rutile. The absence of a major pressure effect on the calibration is suggested, but could not be verified.
Generation of an ab initio model of a Ti defect in a SiO2 crystal was undertaken to assist in validating the Ti/Quartz relationship at pressure and temperature over ranges of 0—11 GPa and 0—1500 K respectively. A Molecular Mechanics model was also constructed due to issues pertaining to the construction of the former. An analysis and comparison of the two models conclude that the Titanium-in-Quartz defect investigated showed a pressure dependence of less than 0.2% over the pressure range.
Defect Formation Free Energy of Defects in Anatase
Timothy C. DuBois
Bachelor of Science (Physics), RMIT University (2009)
Using Density Functional Theory, Lattice Dynamics and other modelling methods a theoretical 2x2x1 supercell of Anatase is defected: removing one Titanium atom and adding four Hydrogen atoms to bond to the Oxygen atoms at the defect site. The system is then probed to discover its properties. The Defect Formation Free Energy of the supercell is of particular interest to this experiment, as this property predicts that the defected lattice is thermodynamically stable across the calculated temperature range of 0K to 1000K. Individual atom contributions to the Free Energy raise some questions to the stability of the Hydrogen atoms, which suggests more calculations are required to verify the need for their introduction into the defect. All results await experimental verification.