-
Particle-In-Cell Code Comparison for Ion Acceleration: EPOCH and Smilei
Authors:
Soham Banerjee,
Joseph R. Smith,
Chris Orban
Abstract:
Particle-in-Cell (PIC) codes are a popular tool to model laser-plasma interactions. Many different PIC codes already exist, and many new PIC codes are being developed constantly. It is therefore important to compare different PIC codes to ascertain which code is best suited for a particular kind of physical problem. In a paper by Smith et al. (2021) they compared three different codes on a problem…
▽ More
Particle-in-Cell (PIC) codes are a popular tool to model laser-plasma interactions. Many different PIC codes already exist, and many new PIC codes are being developed constantly. It is therefore important to compare different PIC codes to ascertain which code is best suited for a particular kind of physical problem. In a paper by Smith et al. (2021) they compared three different codes on a problem relating to proton acceleration in the Target Normal Sheath Acceleration regime from a normal incidence ultra-intense laser pulse. Smith et al. (2021) included in their study the widely used EPOCH code. However, they did not include results from the Smilei code, which is another popular PIC code in the plasma community with a variety of features and physics packages. In the present work, we compare the Smilei code to the EPOCH code for the same test case as Smith et al. (2021). Broadly we find the two codes to be highly consistent with agreement in total ion, electron, and field energy at a percent level or better. The electron and ion energy distribution functions agree well at lower energies and the differences at higher energies (e.g. because of the finite number of macroparticles) are similar to what Smith et al. (2021) saw for other codes. We found that Smilei consumed 25% more RAM than EPOCH did but the execution time was 30% less for Smilei on one processor. We include the input files to encourage future comparisons.
△ Less
Submitted 25 December, 2023;
originally announced December 2023.
-
Applying Machine Learning Methods to Laser Acceleration of Protons: Lessons Learned from Synthetic Data
Authors:
Ronak Desai,
Thomas Zhang,
Ricky Oropeza,
John J. Felice,
Joseph R. Smith,
Alona Kryshchenko,
Chris Orban,
Michael L. Dexter,
Anil K. Patnaik
Abstract:
Researchers in the field of ultra-intense laser science are beginning to embrace machine learning methods. In this study we consider three different machine learning methods -- a two-hidden layer neural network, Support Vector Regression and Gaussian Process Regression -- and compare how well they can learn from a synthetic data set for proton acceleration in the Target Normal Sheath Acceleration…
▽ More
Researchers in the field of ultra-intense laser science are beginning to embrace machine learning methods. In this study we consider three different machine learning methods -- a two-hidden layer neural network, Support Vector Regression and Gaussian Process Regression -- and compare how well they can learn from a synthetic data set for proton acceleration in the Target Normal Sheath Acceleration regime. The synthetic data set was generated from a previously published theoretical model by Fuchs et al. 2005 that we modified. Once trained, these machine learning methods can assist with efforts to maximize the peak proton energy, or with the more general problem of configuring the laser system to produce a proton energy spectrum with desired characteristics. In our study we focus on both the accuracy of the machine learning methods and the performance on one GPU including the memory consumption. Although it is arguably the least sophisticated machine learning model we considered, Support Vector Regression performed very well in our tests.
△ Less
Submitted 16 April, 2024; v1 submitted 29 July, 2023;
originally announced July 2023.
-
Methods to Simplify Object Tracking in Video Data
Authors:
Chris Orban,
Scott Zimmerman,
Jessica T. Kulp,
Jennifer Boughton,
Zach Perrico,
Brianna Rapp,
Richelle Teeling-Smith
Abstract:
Recent years have seen an explosion of interest in analyzing the motion of objects in video data as a way for students to connect the concepts of physics to something tangible like a video recording of an experiment. A variety of software exists for students to look at individual frames and click on the object to infer the x,y position. Some of these tools include a capability to automatically ide…
▽ More
Recent years have seen an explosion of interest in analyzing the motion of objects in video data as a way for students to connect the concepts of physics to something tangible like a video recording of an experiment. A variety of software exists for students to look at individual frames and click on the object to infer the x,y position. Some of these tools include a capability to automatically identify the position of the object in the frame. But it is not unusual, especially when inexperienced users are recording the video and configuring the program, for these algorithms to struggle to "lock on" to the moving object. In this paper, we both include some general advice to help object tracking algorithms locate an object and we provide our own algorithms that are simpler and potentially more effective than the sophisticated image processing algorithms that are currently being used. These algorithms focus not on a "template image" of the moving object but instead distinguish between the object and the background by analyzing only the colors of individual pixels. These algorithms are built into a free and open source program called the STEMcoding Object Tracker (http://go.osu.edu/objecttracker) which works in the browser (without any downloads) and is compatible with a variety of operating systems including chromebooks.
△ Less
Submitted 14 October, 2022; v1 submitted 16 June, 2022;
originally announced June 2022.
-
A Particle-In-Cell Code Comparison for Ion Acceleration: EPOCH, LSP, and WarpX
Authors:
Joseph R. Smith,
Chris Orban,
Nashad Rahman,
Brendan McHugh,
Ricky Oropeza,
Enam A. Chowdhury
Abstract:
There are now more Particle-in-Cell (PIC) codes than ever before that researchers use to simulate intense laser-plasma interactions. To date, there have been relatively few direct comparisons of these codes in the literature, especially for relativistic intensity lasers interacting with thin overdense targets. To address this we perform a code comparison of three PIC codes: EPOCH, LSP, and WarpX f…
▽ More
There are now more Particle-in-Cell (PIC) codes than ever before that researchers use to simulate intense laser-plasma interactions. To date, there have been relatively few direct comparisons of these codes in the literature, especially for relativistic intensity lasers interacting with thin overdense targets. To address this we perform a code comparison of three PIC codes: EPOCH, LSP, and WarpX for the problem of laser-driven ion acceleration in a 2D(3v) geometry for a $10^{20}$ W cm$^{-2}$ intensity laser. We examine the plasma density, ion energy spectra, and laser-plasma coupling of the three codes and find strong agreement. We also run the same simulation 20 times with different random seeds to explore statistical fluctuations of the outputs. We then compare the execution times and memory usage of the codes (without "tuning" to improve performance) using between 1 and 48 processors on one node. We provide input files to encourage larger and more frequent code comparisons in this field.
△ Less
Submitted 30 June, 2021; v1 submitted 31 March, 2021;
originally announced March 2021.
-
Particle-in-Cell modeling of a potential demonstration experiment for double pulse enhanced target normal sheath acceleration
Authors:
Nashad Rahman,
Joseph R. Smith,
Gregory Ngirmang,
Chris Orban
Abstract:
Ultra intense lasers are a promising source of energetic ions for various applications. An interesting approach described in Ferri et al. 2019 argues from Particle-in-Cell simulations that using two laser pulses of half energy (half intensity) arriving with close to 45 degrees angle of incidence is significantly more effective at accelerating ions than one pulse at full energy (full intensity). Fo…
▽ More
Ultra intense lasers are a promising source of energetic ions for various applications. An interesting approach described in Ferri et al. 2019 argues from Particle-in-Cell simulations that using two laser pulses of half energy (half intensity) arriving with close to 45 degrees angle of incidence is significantly more effective at accelerating ions than one pulse at full energy (full intensity). For a variety of reasons, at the time of this writing there has not yet been a true experimental confirmation of this enhancement. In this paper we perform 2D Particle-in-Cell simulations to examine if a milliJoule class, 5x10^18 W cm^-2 peak intensity laser system could be used for such a demonstration experiment. Laser systems in this class can operate at a kHz rate which should be helpful for addressing some of the challenges of performing this experiment. Despite investigating a 3.5 times lower intensity than Ferri et al. 2019 did, we find that the double pulse approach enhances the peak proton energy and the energy conversion to protons by a factor of about three compared to a single laser pulse with the same total laser energy. We also comment on the nature of the enhancement and describe simulations that examine how the enhancement may depend on the spatial or temporal alignment of the two pulses.
△ Less
Submitted 8 October, 2021; v1 submitted 12 January, 2021;
originally announced January 2021.
-
Code-to-Code Comparison and Validation of the Radiation-Hydrodynamics Capabilities of the FLASH Code Using a Laboratory Astrophysical Jet
Authors:
Chris Orban,
Milad Fatenejad,
Don Q. Lamb
Abstract:
The potential for laser-produced plasmas to yield fundamental insights into high energy density physics (HEDP) and deliver other useful applications can sometimes be frustrated by uncertainties in modeling the properties and behavior of these plasmas using radiation-hydrodynamics codes. In an effort to overcome this and to corroborate the accuracy of the HEDP capabilities in the publicly available…
▽ More
The potential for laser-produced plasmas to yield fundamental insights into high energy density physics (HEDP) and deliver other useful applications can sometimes be frustrated by uncertainties in modeling the properties and behavior of these plasmas using radiation-hydrodynamics codes. In an effort to overcome this and to corroborate the accuracy of the HEDP capabilities in the publicly available FLASH radiation-hydrodynamics code, we present detailed code-to-code comparisons between FLASH and the HYDRA code developed at Lawrence Livermore National Laboratory using previously published HYDRA simulations from Grava et al. (2008). That study describes a laser experiment that produced a jet-like feature that the authors compare to astrophysical jets. Importantly, the Grava et al. (2008) experiment included detailed x-ray interferometric measurements of electron number densities and a time-integrated measurement of the soft x-ray spectrum. Despite markedly different methods for treating the computational mesh, and different equation of state and opacity models, the FLASH results resemble the results from HYDRA and, most importantly, the experimental measurements of electron density. Having validated the FLASH code in this way, we use the code to further investigate and understand the formation of the jet seen in the Grava et al. (2008) experiment and discuss its relation to the Wan et al. (1997) experiment at the NOVA laser.
△ Less
Submitted 18 April, 2022; v1 submitted 17 June, 2020;
originally announced June 2020.
-
Optimizing Laser-Plasma Interactions for Ion Acceleration using Particle-in-Cell Simulations and Evolutionary Algorithms
Authors:
Joseph R. Smith,
Chris Orban,
John T. Morrison,
Kevin M. George,
Gregory K. Ngirmang,
Enam A. Chowdhury,
W. Mel Roquemore
Abstract:
The development of ultra-intense laser-based sources of high energy ions is an important goal, with a variety of potential applications. One of the barriers to achieving this goal is the need to maximize the conversion efficiency from laser energy to ion energy. We apply a new approach to this problem, in which we use an evolutionary algorithm to optimize conversion efficiency by exploring variati…
▽ More
The development of ultra-intense laser-based sources of high energy ions is an important goal, with a variety of potential applications. One of the barriers to achieving this goal is the need to maximize the conversion efficiency from laser energy to ion energy. We apply a new approach to this problem, in which we use an evolutionary algorithm to optimize conversion efficiency by exploring variations of the target density profile with thousands of one-dimensional particle-in-cell (PIC) simulations. We then compare this "optimal" target identified by the one-dimensional PIC simulations to more conventional choices, such as with an exponential scale length pre-plasma, with fully three-dimensional PIC simulations. The optimal target outperforms the conventional targets in terms of maximum ion energy by 20% and show a significant enhancement of conversion efficiency to high energy ions. This target geometry enhances laser coupling to the electrons, while still allowing the laser to strongly reflect from an effectively thin target. These results underscore the potential for this statistics-driven approach to guide research into optimizing laser-plasma simulations and experiments.
△ Less
Submitted 15 May, 2020;
originally announced May 2020.
-
Using Virtual Reality in Electrostatics Instruction: The Impact of Training
Authors:
Chris D. Porter,
Joseph R. H. Smith,
Eric M. Stagar,
Amber Simmons,
Megan Nieberding,
Jonathan R. Brown,
Abigail Ayers,
Chris M. Orban
Abstract:
Recent years have seen a resurgence of interest in using Virtual Reality (VR) technology to benefit instruction, especially in physics and related subjects. As VR devices improve and become more widely available, there remains a number of unanswered questions regarding the impact of VR on student learning and how best to use this technology in the classroom. On the topic of electrostatics, for exa…
▽ More
Recent years have seen a resurgence of interest in using Virtual Reality (VR) technology to benefit instruction, especially in physics and related subjects. As VR devices improve and become more widely available, there remains a number of unanswered questions regarding the impact of VR on student learning and how best to use this technology in the classroom. On the topic of electrostatics, for example, a large, controlled, randomized study performed by Smith et al. 2017\cite{smith17}, found that VR-based instruction had an overall negligible impact on student learning compared to videos or images. However, they did find a strong trend for students who reported frequent video game play to learn better from VR than other media. One possible interpretation of this result is that extended videogame play provides a kind of "training" that enables a student to learn more comfortably in the virtual environment. In the present work we consider if a VR training activity that is unrelated to electrostatics can help prepare students to learn electrostatics from subsequent VR instruction. We find that preliminary VR training leads to a small but statistically significant improvement in student performance on our electrostatics assessment. We also find that student reported game play is still correlated with higher scores on this metric.
△ Less
Submitted 1 July, 2020; v1 submitted 22 January, 2020;
originally announced January 2020.
-
Experimental observation and computational modeling of radial Weibel instability in high intensity laser-plasma interactions
Authors:
Gregory K. Ngirmang,
John T. Morrison,
Kevin M. George,
Joseph R. Smith,
Kyle D. Frische,
Chris Orban,
Enam A. Chowdhury,
W. Mel Roquemore
Abstract:
When a relativistic intensity laser interacts with the surface of a solid density target, suprathermal electron currents are subject to Weibel instability filamentation when propagating through the thermal population of the bulk target. We present time resolved shadowgraphy of radial ionization front expansion and Weibel instability filamentation within a thin, sub-micron, sheet initiated by irrad…
▽ More
When a relativistic intensity laser interacts with the surface of a solid density target, suprathermal electron currents are subject to Weibel instability filamentation when propagating through the thermal population of the bulk target. We present time resolved shadowgraphy of radial ionization front expansion and Weibel instability filamentation within a thin, sub-micron, sheet initiated by irradiation with a short pulse, high intensity laser. High temporal (100 fs) and spatial (1 $μ$m) resolution shadowgraphy of the interaction reveals a relativistic expansion of the ionization front within a 120 $μ$m diameter region surrounding the laser-target interaction, corroborated by simulations to expand at $0.77c$, where $c$ is the speed of light. Filamentation within the patch persists for several picoseconds and seeds the eventual recombination and heating dynamics on the nanosecond timescale. Particle-in-cell simulations were conducted to elucidate the electron dynamics leading to the radial expansion of the critical surface. Computational results report the ionization expansion is due to field ionization of the expanding hot electron population. Filamentation within the expansion is due to the Weibel instability which is supported by the magnetic fields present.
△ Less
Submitted 28 October, 2019;
originally announced October 2019.
-
Computational Thinking in Introductory Physics
Authors:
Chris Orban,
Richelle Teeling-Smith
Abstract:
Computational Thinking (CT) is still a relatively new term in the lexicon of learning objectives and science standards. There is not yet widespread agreement on the precise definition or implementation of CT, and efforts to assess CT are still maturing, even as more states adopt K-12 computer science standards. In this article we will try to summarize what CT means for a typical introductory (i.e.…
▽ More
Computational Thinking (CT) is still a relatively new term in the lexicon of learning objectives and science standards. There is not yet widespread agreement on the precise definition or implementation of CT, and efforts to assess CT are still maturing, even as more states adopt K-12 computer science standards. In this article we will try to summarize what CT means for a typical introductory (i.e. high school or early college) physics class. This will include a discussion of the ways that instructors may already be incorporating elements of CT in their classes without knowing it. Our intention in writing this article is to provide a helpful, concise and readable introduction to this topic for physics instructors. We also put forward some ideas for what the future of CT in introductory physics may look like.
△ Less
Submitted 20 November, 2019; v1 submitted 18 July, 2019;
originally announced July 2019.
-
A controlled study of virtual reality in first-year magnetostatics
Authors:
Chris D. Porter,
Jonathan Brown,
Joseph R. Smith,
Amber Simmons,
Megan Nieberding,
Abigail E. Ayers,
Chris Orban
Abstract:
Stereoscopic virtual reality (VR) has experienced a resurgence due to flagship products such as the Oculus Rift, HTC Vive and smartphone-based VR solutions like Google Cardboard. This is causing the question to resurface: how can stereoscopic VR be useful in instruction, if at all, and what are the pedagogical best practices for its use? To address this, and to continue our work in this sphere, we…
▽ More
Stereoscopic virtual reality (VR) has experienced a resurgence due to flagship products such as the Oculus Rift, HTC Vive and smartphone-based VR solutions like Google Cardboard. This is causing the question to resurface: how can stereoscopic VR be useful in instruction, if at all, and what are the pedagogical best practices for its use? To address this, and to continue our work in this sphere, we performed a study of 289 introductory physics students who were sorted into three different treatment types: stereoscopic virtual reality, WebGL simulation, and static 2D images, each designed to provide information about magnetic fields and forces. Students were assessed using preliminary items designed to focus on heavily-3D systems. We report on assessment reliability, and on student performance. Overall, we find that students who used VR did not significantly outperform students using other treatment types. There were significant differences between sexes, as other studies have noted. Dependence on students' self-reported 3D videogame play was observed, in keeping with previous studies, but this dependence was not restricted to the VR treatment.
△ Less
Submitted 11 July, 2019;
originally announced July 2019.
-
Particle-in-Cell Simulations of Density Peak Formation and Ion Acceleration from Short Pulse Laser-Driven Ponderomotive Steepening
Authors:
Joseph R. Smith,
Chris Orban,
Gregory K. Ngirmang,
John T. Morrison,
Kevin M. George,
Enam A. Chowdhury,
W. M. Roquemore
Abstract:
We use particle-in-cell (PIC) simulations and simple analytic models to investigate the laser-plasma interaction known as ponderomotive steepening. When normally incident laser light reflects at the critical surface of a plasma, the resulting standing electromagnetic wave modifies the electron density profile via the ponderomotive force, which creates peaks in the electron density separated by app…
▽ More
We use particle-in-cell (PIC) simulations and simple analytic models to investigate the laser-plasma interaction known as ponderomotive steepening. When normally incident laser light reflects at the critical surface of a plasma, the resulting standing electromagnetic wave modifies the electron density profile via the ponderomotive force, which creates peaks in the electron density separated by approximately half of the laser wavelength. What is less well studied is how this charge imbalance accelerates ions towards the electron density peaks, modifying the ion density profile of the plasma. Idealized PIC simulations with an extended underdense plasma shelf are used to isolate the dynamics of ion density peak growth for a 42 fs pulse from an 800 nm laser with an intensity of 10$^{18}$ W cm$^{-2}$. These simulations exhibit sustained longitudinal electric fields of 200 GV m$^{-1}$, which produce counter-steaming populations of ions reaching a few keV in energy. We compare these simulations to theoretical models, and we explore how ion energy depends on factors such as the plasma density and the laser wavelength, pulse duration, and intensity. We also provide relations for the strength of longitudinal electric fields and an approximate timescale for the density peaks to develop. These conclusions may be useful investigating the phenomenon of ponderomotive steepening as advances in laser technology allow shorter and more intense pulses to be produced at various wavelengths. We also discuss the parallels with other work studying the interference from two counter-propagating laser pulses.
△ Less
Submitted 2 May, 2019;
originally announced May 2019.
-
High repetition rate (>= kHz) targets and optics from liquid microjets for the study and application of high intensity laser-plasma interactions
Authors:
K. M. George,
J. T. Morrison,
S. Feister,
G. Ngirmang,
J. R. Smith,
A. J. Klim,
J. Snyder,
D. Austin,
W. Erbsen,
K. D. Frische,
J. Nees,
C. Orban,
E. A. Chowdhury,
W. M. Roquemore
Abstract:
High intensity laser-plasma interactions produce a wide array of energetic particles and beams with promising applications. Unfortunately, high repetition rate and high average power requirements for many applications are not satisfied by the lasers, optics, targets, and diagnostics currently employed. Here, we address the need for high repetition rate targets and optics through the use of liquids…
▽ More
High intensity laser-plasma interactions produce a wide array of energetic particles and beams with promising applications. Unfortunately, high repetition rate and high average power requirements for many applications are not satisfied by the lasers, optics, targets, and diagnostics currently employed. Here, we address the need for high repetition rate targets and optics through the use of liquids. A novel nozzle assembly is used to generate high-velocity, laminar-flowing liquid microjets which are compatible with a low-vacuum environment, generate little to no debris, and exhibit precise positional and dimensional tolerances. Jets, droplets, submicron thick sheets, and other configurations are characterized with pump-probe shadowgraphy to evaluate their use as targets. To demonstrate a high repetition rate, consumable liquid optical element, we present a plasma mirror created by a submicron thick liquid sheet. This plasma mirror provides etalon-like anti-reflection properties in the low-field of 0.1% and high reflectivity as a plasma, 69%, at a repetition rate of 1 kHz. Practical considerations of fluid compatibility, in-vacuum operation, and estimates of maximum repetition rate, in excess of 10 kHz, are addressed. The targets and optics presented here enable the use of relativistically intense lasers at high average power and make possible many long proposed applications.
△ Less
Submitted 20 February, 2019; v1 submitted 12 February, 2019;
originally announced February 2019.
-
A Hybrid Approach for Using Programming Exercises in Introductory Physics
Authors:
Chris Orban,
Richelle M. Teeling-Smith,
Joseph R. H. Smith,
Chris D. Porter
Abstract:
Incorporating computer programming exercises into introductory physics is a delicate task that involves a number of choices that may have an effect on student learning. We present a "hybrid" approach that speaks to a number of common concerns regarding cognitive load which arise when using programming exercises in introductory physics classes where many students are absolute beginner programmers.…
▽ More
Incorporating computer programming exercises into introductory physics is a delicate task that involves a number of choices that may have an effect on student learning. We present a "hybrid" approach that speaks to a number of common concerns regarding cognitive load which arise when using programming exercises in introductory physics classes where many students are absolute beginner programmers. This "hybrid" approach provides the student with a highly interactive web-based visualization, not unlike a PhET or Physlet interactive, but importantly the student is shown only the subset of the code where the initial conditions are set and the system variables are evolved. We highlight results from a coding activity that resembles the classic game Asteroids. The goals of this activity are to show how a simple 1D code can be modified into a 2D code, and to reinforce ideas about the relationship between force, velocity, and acceleration vectors. Survey results from four semesters of introductory physics classes at the Ohio State University's Marion campus, in which a high percentage of the students are weak or absolute beginner programmers, provide evidence that most students can complete coding tasks without severe difficulty. Other survey results are promising for future work where conceptual learning will be assessed in a direct way using metrics like the Animated Force Concept Inventory (Dancy & Beichner 2006). The exercise highlighted here and others from our group are available for general use at http://compadre.org/PICUP
△ Less
Submitted 23 June, 2018;
originally announced June 2018.
-
A controlled study of stereoscopic virtual reality in freshman electrostatics
Authors:
Joseph R. Smith,
Amber Byrum,
Timothy M. McCormick,
Nick Young,
Chris Orban,
Chris D. Porter
Abstract:
Virtual reality (VR) has long promised to revolutionize education, but with little follow-through. Part of the reason for this is the prohibitive cost of immersive VR headsets or caves. This has changed with the advent of smartphone-based VR (along the lines of Google cardboard) which allows students to use smartphones and inexpensive plastic or cardboard viewers to enjoy stereoscopic VR simulatio…
▽ More
Virtual reality (VR) has long promised to revolutionize education, but with little follow-through. Part of the reason for this is the prohibitive cost of immersive VR headsets or caves. This has changed with the advent of smartphone-based VR (along the lines of Google cardboard) which allows students to use smartphones and inexpensive plastic or cardboard viewers to enjoy stereoscopic VR simulations. We have completed the largest-ever such study on 627 students enrolled in calculus-based freshman physics at The Ohio State University. This initial study focused on student understanding of electric fields. Students were split into three treatments groups: VR, video, and static 2D images. Students were asked questions before, during, and after treatment. Here we present a preliminary analysis including overall post-pre improvement among the treatment groups, dependence of improvement on gender, and previous video game experience. Results on select questions are discussed. Several electric field visualizations similar to those used in this study are freely available on Google Play http://go.osu.edu/BuckeyeVR
△ Less
Submitted 5 July, 2017;
originally announced July 2017.
-
A novel approach for using programming exercises in electromagnetism coursework
Authors:
Chris Orban,
Chris D. Porter,
Nash K. Brecht,
Richelle M. Teeling-Smith,
Kathy A. Harper
Abstract:
While there exists a significant number of web interactives for introductory physics, students are almost never shown the computer code that generates these interactives even when the physics parts of these programs are relatively simple. Building off of a set of carefully-designed classical mechanics programming exercises that were constructed with this goal in mind, we present a series of electr…
▽ More
While there exists a significant number of web interactives for introductory physics, students are almost never shown the computer code that generates these interactives even when the physics parts of these programs are relatively simple. Building off of a set of carefully-designed classical mechanics programming exercises that were constructed with this goal in mind, we present a series of electromagnetism programming exercises in a browser-based framework called p5.js. Importantly, this framework can be used to highlight the physics aspects of an interactive simulation code while obscuring other details. This approach allows absolute beginner programmers to gain experience in modifying and running the program without becoming overwhelmed. We plan to probe the impact on student conceptual learning using the Brief Electricity and Magnetism Assessment and other questions. We invite collaborators and teachers to adopt this framework in their high school or early undergraduate classes. All exercises are available at http://compadre.org/PICUP
△ Less
Submitted 1 July, 2017;
originally announced July 2017.
-
Particle-in-Cell Simulations of Electron Beam Production from Infrared Ultra-intense Laser Interactions
Authors:
Gregory K. Ngirmang,
Chris Orban,
Scott Feister,
John T. Morrison,
Enam A. Chowdhury,
W. M. Roquemore
Abstract:
Advances in ultra-intense laser technology are enabling, for the first time, relativistic intensities at mid-infrared (mid-IR) wavelengths. Anticipating further experimental research in this domain, we present high-resolution two dimensional Particle-in-Cell (PIC) simulation results using the Large- Scale Plasma (LSP) code that explore intense mid-IR laser interactions with dense targets. We prese…
▽ More
Advances in ultra-intense laser technology are enabling, for the first time, relativistic intensities at mid-infrared (mid-IR) wavelengths. Anticipating further experimental research in this domain, we present high-resolution two dimensional Particle-in-Cell (PIC) simulation results using the Large- Scale Plasma (LSP) code that explore intense mid-IR laser interactions with dense targets. We present the results of thirty PIC simulations over a wide range of intensities (0.03 < $a_0$ < 39) and wavelengths (λ =780 nm, 3 μm, and 10 μm). Earlier studies, limited to λ =780 nm and $a_0 \sim$ 1 [1, 2], identified super-ponderomotive electron acceleration in the laser specular direction for normal- incidence laser interactions with dense targets. We extend this research to mid-IR wavelengths and find a more general result that normal-incidence super-ponderomotive electron acceleration occurs provided that the laser intensity is not highly relativistic ($a_0 \lesssim 1$) and that the pre-plasma scale length is similar to or longer than the laser wavelength. Under these conditions, ejected electron angular and energy distributions are similar to expectations from an analytic model used in [2]. We also find that, for $a_0 \sim 1$, the mid-IR simulations exhibit a classic ponderomotive steepening pattern with multiple peaks in the ion and electron density distribution. Experimental validation of this basic laser-plasma interaction process will be possible in the near future using mid-IR laser technology and interferometry.
△ Less
Submitted 26 April, 2017;
originally announced April 2017.
-
A Game-Centered, Interactive Approach for Using Programming Exercises in Introductory Physics
Authors:
Chris Orban,
Chris Porter,
Joseph R. H. Smith,
Nash K. Brecht,
Chris A. Britt,
Richelle M. Teeling-Smith,
Kathy A. Harper
Abstract:
Incorporating computer programming exercises in introductory physics is a delicate task that involves a number of choices that may have a strong affect on student learning. We present an approach that speaks to a number of common concerns that arise when using programming exercises in introductory physics classes where most students are absolute beginner programmers. These students need an approac…
▽ More
Incorporating computer programming exercises in introductory physics is a delicate task that involves a number of choices that may have a strong affect on student learning. We present an approach that speaks to a number of common concerns that arise when using programming exercises in introductory physics classes where most students are absolute beginner programmers. These students need an approach that is (1) simple, involving 75 or fewer lines of well-commented code, (2) easy to use, with browser-based coding tools, (3) interactive, with a high frame rate to give a video-game like feel, (4) step-by-step with the ability to interact with intermediate stages of the "correct" program and (5) thoughtfully integrated into the physics curriculum, for example, by illustrating velocity and acceleration vectors throughout. We present a set of hour-long activities for classical mechanics that resemble well-known games such as "asteroids", "lunar lander" and "angry birds". Survey results from the first activity from four semesters of introductory physics classes at OSU in which a high percentage of the students are weak or absolute beginner programmers seems to confirm that the level of difficulty is appropriate for this level and that the students enjoy the activity. These exercises are available for general use at http://compadre.org/PICUP In the future we plan to assess conceptual knowledge using an animated version of the Force Concept Inventory originally developed by M. Dancy.
△ Less
Submitted 31 May, 2017; v1 submitted 7 January, 2017;
originally announced January 2017.
-
3D PIC simulations of electron beams created via reflection of intense laser light from a water target
Authors:
Gregory K. Ngirmang,
Chris Orban,
Scott Feister,
John T. Morrison,
Kyle Frische,
Enam A. Chowdhury,
W. M. Roquemore
Abstract:
We present 3D Particle-in-Cell (PIC) modeling of an ultra-intense laser experiment by the Extreme Light group at the Air Force Research Laboratory (AFRL) using the PIC code LSP. This is the first time PIC simulations have been performed in 3D for this experiment which involves an ultra-intense, short-pulse (30 fs) laser interacting with a water jet target at normal incidence. These 3D PIC simulati…
▽ More
We present 3D Particle-in-Cell (PIC) modeling of an ultra-intense laser experiment by the Extreme Light group at the Air Force Research Laboratory (AFRL) using the PIC code LSP. This is the first time PIC simulations have been performed in 3D for this experiment which involves an ultra-intense, short-pulse (30 fs) laser interacting with a water jet target at normal incidence. These 3D PIC simulation results are compared to results from 2D(3$v$) PIC simulations for both $5.4~\cdot~10^{17}$ W cm$^{-2}$ and $3~\cdot~10^{17}$ W cm$^{-2}$ intensities. Comparing the 2D(3$v$) and 3D simulation results, the laser-energy-to-ejected-electron-energy conversion efficiencies were comparable, but the angular distribution of ejected electrons show interesting differences with qualitative differences at higher intensity. An analytic plane-wave model is discussed which provides some explanation for the angular distribution and energies of ejected electrons in the 2D(3$v$) simulations. We also performed a 3D simulation with circularly polarized light and found a significantly higher conversion efficiency and peak electron energy, which is promising for future experiments.
△ Less
Submitted 5 January, 2016; v1 submitted 16 October, 2015;
originally announced October 2015.
-
Relativistic electron acceleration by mJ-class kHz lasers normally incident on liquid targets
Authors:
Scott Feister,
Drake R. Austin,
John T. Morrison,
Kyle D. Frische,
Chris Orban,
Gregory Ngirmang,
Abraham Handler,
Joseph R. H. Smith,
Mark Schillaci,
Jay A. LaVerne,
Enam A. Chowdhury,
R. R. Freeman,
W. M. Roquemore
Abstract:
We report observation of kHz-pulsed-laser-accelerated electron energies up to 3 MeV in the -$k_\text{laser}$ (backward) direction from a 3 mJ laser interacting at normal incidence with a solid density, flowing-liquid target. The electrons/MeV/s.r. >1 MeV recorded here using a mJ-class laser exceeds or equals that of prior super-ponderomotive electron studies employing lasers at lower repetition-ra…
▽ More
We report observation of kHz-pulsed-laser-accelerated electron energies up to 3 MeV in the -$k_\text{laser}$ (backward) direction from a 3 mJ laser interacting at normal incidence with a solid density, flowing-liquid target. The electrons/MeV/s.r. >1 MeV recorded here using a mJ-class laser exceeds or equals that of prior super-ponderomotive electron studies employing lasers at lower repetition-rates and oblique incidence. Focal intensity of the 40-fs-duration laser is 1.5 $\cdot$ 10$^{18}$ W cm$^{-2}$, corresponding to only ~80 keV electron ponderomotive energy. Varying laser intensity confirms electron energies in the laser-reflection direction well above what might be expected from ponderomotive scaling in normal-incidence laser-target geometry. This direct, normal-incidence energy spectrum measurement is made possible by modifying the final focusing off-axis-paraboloid (OAP) mirror with a central hole that allows electrons to pass, and restoring laser intensity through adaptive optics. A Lanex-based, optics-free high-acquisition rate (>100 Hz) magnetic electron-spectrometer was developed for this study to enable shot-to-shot statistical analysis and real-time feedback, which was leveraged in finding optimal pre-plasma conditions. 3D Particle-in-cell simulations of the interaction show qualitative super-ponderomotive spectral agreement with experiment. The demonstration of a high-repetition-rate, high-flux source containing >MeV electrons from a few-mJ, 40 fs laser and a simple liquid target encourages development of future $\geq$kHz-repetition, fs-duration electron-beam applications.
△ Less
Submitted 20 July, 2017; v1 submitted 28 August, 2015;
originally announced August 2015.
-
Backward-propagating MeV electrons from $10^{18}$ W/cm$^2$ laser interactions with water
Authors:
John T. Morrison,
Enam A. Chowdhury,
Kyle D. Frische,
Scott Feister,
Vladimir M. Ovchinnikov,
John A. Nees,
Chris Orban,
Richard R. Freeman,
W. Melvyn Roquemore
Abstract:
We present an experimental study of the generation of $\sim$MeV electrons opposite to the direction of laser propagation following the relativistic interaction at normal incidence of a $\sim$3 mJ, $10^{18}$ W/cm$^2$ short pulse laser with a flowing 30 $μ$m diameter water column target. Faraday cup measurements record hundreds of pC charge accelerated to energies exceeding 120 keV, and energy-resol…
▽ More
We present an experimental study of the generation of $\sim$MeV electrons opposite to the direction of laser propagation following the relativistic interaction at normal incidence of a $\sim$3 mJ, $10^{18}$ W/cm$^2$ short pulse laser with a flowing 30 $μ$m diameter water column target. Faraday cup measurements record hundreds of pC charge accelerated to energies exceeding 120 keV, and energy-resolved measurements of secondary x-ray emissions reveal an x-ray spectrum peaking above 800 keV, which is significantly higher energy than previous studies with similar experimental conditions and more than five times the $\sim$110 keV ponderomotive energy scale for the laser. We show that the energetic x-rays generated in the experiment result from backward-going, high-energy electrons interacting with the focusing optic and vacuum chamber walls with only a small component of x-ray emission emerging from the target itself. We also demonstrate that the high energy radiation can be suppressed through the attenuation of the nanosecond-scale pre-pulse. These results are supported by 2D Particle-in-Cell (PIC) simulations of the laser-plasma interaction that exhibit beam-like backward-propagating MeV electrons.
△ Less
Submitted 15 March, 2015; v1 submitted 9 January, 2015;
originally announced January 2015.
-
A Novel Femtosecond-Gated, High-Resolution, Frequency-Shifted Shearing Interferometry Technique for Probing Pre-Plasma Expansion in Ultra-Intense Laser Experiments
Authors:
S. Feister,
J. A. Nees,
J. T. Morrison,
K. D. Frische,
C. Orban,
E. A. Chowdhury,
W. M. Roquemore
Abstract:
Ultra-intense laser-matter interaction experiments (>10$^{18}$ W/cm$^{2}$) with dense targets are highly sensitive to the effect of laser "noise" (in the form of pre-pulses) preceding the main ultra-intense pulse. These system-dependent pre-pulses in the nanosecond and/or picosecond regimes are often intense enough to modify the target significantly by ionizing and forming a plasma layer in front…
▽ More
Ultra-intense laser-matter interaction experiments (>10$^{18}$ W/cm$^{2}$) with dense targets are highly sensitive to the effect of laser "noise" (in the form of pre-pulses) preceding the main ultra-intense pulse. These system-dependent pre-pulses in the nanosecond and/or picosecond regimes are often intense enough to modify the target significantly by ionizing and forming a plasma layer in front of the target. Time resolved interferometry offers a robust way to characterize the expanding plasma during this period. We have developed a novel pump-probe interferometry system for an ultra-intense laser experiment that uses two short-pulse amplifiers synchronized by one ultra-fast seed oscillator to achieve 40-femtosecond time resolution over hundreds of nanoseconds, using a variable delay line and other techniques. The first of these amplifiers acts as the pump and delivers maximal energy to the interaction region. The second amplifier is frequency shifted and then frequency doubled to generate the femtosecond probe pulse. After passing through the laser-target interaction region, the probe pulse is split and recombined in a laterally sheared Michelson interferometer. Importantly, the frequency shift in the probe allows strong plasma self-emission at the second harmonic of the pump to be filtered out, allowing plasma expansion near the critical surface and elsewhere to be clearly visible in the interferograms. Three-phase reconstructions of the electron densities are then inferred by Abel inversion. This interferometric system delivers precise measurements of pre-plasma expansion that can identify the condition of the target at the moment that the ultra-intense pulse arrives. Such measurements are indispensable for correlating laser pre-pulse measurements with instantaneous plasma profiles and for enabling realistic Particle-in-Cell simulations of the ultra-intense laser-matter interaction.
△ Less
Submitted 17 July, 2014; v1 submitted 13 June, 2014;
originally announced June 2014.
-
Backward-Propagating MeV Electrons in Ultra-Intense Laser Interactions: Standing Wave Acceleration and Coupling to the Reflected Laser Pulse
Authors:
Chris Orban,
John T. Morrison,
Enam D. Chowdhury,
John A. Nees,
Kyle Frische,
Scott Feister,
W. Melvyn Roquemore
Abstract:
Laser-accelerated electron beams have been created at a kHz repetition rate from the {\it reflection} of intense ($\sim10^{18}$ W/cm$^2$), $\sim$40 fs laser pulses focused on a continuous water-jet in an experiment at the Air Force Research Laboratory. This paper investigates Particle-in-Cell (PIC) simulations of the laser-target interaction to identify the physical mechanisms of electron accelera…
▽ More
Laser-accelerated electron beams have been created at a kHz repetition rate from the {\it reflection} of intense ($\sim10^{18}$ W/cm$^2$), $\sim$40 fs laser pulses focused on a continuous water-jet in an experiment at the Air Force Research Laboratory. This paper investigates Particle-in-Cell (PIC) simulations of the laser-target interaction to identify the physical mechanisms of electron acceleration in this experiment. We find that the standing-wave pattern created by the overlap of the incident and reflected laser is particularly important because this standing wave can "inject" electrons into the reflected laser pulse where the electrons are further accelerated. We identify two regimes of standing wave acceleration: a highly relativistic case ($a_0~\geq~1$), and a moderately relativistic case ($a_0~\sim~0.5$) which operates over a larger fraction of the laser period. In previous studies, other groups have investigated the highly relativistic case for its usefulness in launching electrons in the forward direction. We extend this by investigating electron acceleration in the {\it specular (back reflection) direction} and over a wide range of intensities ($10^{17}-10^{19}$ W cm$^{-2}$).
△ Less
Submitted 22 January, 2015; v1 submitted 24 May, 2014;
originally announced May 2014.
-
Cosmological Perturbation Theory as a Tool for Estimating Box-Scale Effects in N-body Simulations
Authors:
Chris Orban
Abstract:
In performing cosmological N-body simulations, it is widely appreciated that the growth of structure on the largest scales within a simulation box will be inhibited by the finite size of the simulation volume. Following ideas set forth in Seto (1999), this paper shows that standard (a.k.a. 1-loop) cosmological perturbation theory (SPT) can be used to predict, in an approximate way, the deleterious…
▽ More
In performing cosmological N-body simulations, it is widely appreciated that the growth of structure on the largest scales within a simulation box will be inhibited by the finite size of the simulation volume. Following ideas set forth in Seto (1999), this paper shows that standard (a.k.a. 1-loop) cosmological perturbation theory (SPT) can be used to predict, in an approximate way, the deleterious effect of the box scale on the power spectrum of density fluctuations in simulation volumes. Alternatively, this approach can be used to quickly estimate post facto the effect of the box scale on power spectrum results from existing simulations. In this way SPT can help determine whether larger box sizes or other more-sophisticated methods are needed to achieve a particular level of precision for a given application (e.g. simulations to measure the non-linear evolution of baryon acoustic oscillations). I focus on SPT in this note and show that its predictions differ only by about a factor of two or less from the measured suppression inferred from both powerlaw and $Λ$CDM $N$-body simulations. It should be possible to improve the accuracy of these predictions through using more-sophisticated perturbation theory models. An appendix compares power spectrum measurements from the powerlaw simulations at outputs where box-scale effects are minimal to perturbation theory models and previously-published fitting functions. These power spectrum measurements are included with this paper to aid efforts to develop new perturbation theory models.
△ Less
Submitted 16 June, 2014; v1 submitted 17 December, 2013;
originally announced December 2013.
-
A Radiation-Hydrodynamics Code Comparison for Laser-Produced Plasmas: FLASH versus HYDRA and the Results of Validation Experiments
Authors:
Chris Orban,
Milad Fatenejad,
Sugreev Chawla,
Scott C. Wilks,
Donald Q. Lamb
Abstract:
The potential for laser-produced plasmas to yield fundamental insights into high energy density physics (HEDP) and deliver other useful applications can sometimes be frustrated by uncertainties in modeling the properties and expansion of these plasmas using radiation-hydrodynamics codes. In an effort to overcome this and to corroborate the accuracy of the HEDP capabilities recently added to the pu…
▽ More
The potential for laser-produced plasmas to yield fundamental insights into high energy density physics (HEDP) and deliver other useful applications can sometimes be frustrated by uncertainties in modeling the properties and expansion of these plasmas using radiation-hydrodynamics codes. In an effort to overcome this and to corroborate the accuracy of the HEDP capabilities recently added to the publicly available FLASH radiation-hydrodynamics code, we present detailed comparisons of FLASH results to new and previously published results from the HYDRA code used extensively at Lawrence Livermore National Laboratory. We focus on two very different problems of interest: (1) an Aluminum slab irradiated by 15.3 and 76.7 mJ of "pre-pulse" laser energy and (2) a mm-long triangular groove cut in an Aluminum target irradiated by a rectangular laser beam. Because this latter problem bears a resemblance to astrophysical jets, Grava et al., Phys. Rev. E, 78, (2008) performed this experiment and compared detailed x-ray interferometric measurements of electron number densities to HYDRA simulations. Thus, the former problem provides an opportunity for code-to-code comparison, while the latter provides an opportunity for both code-to-code comparison and validation. Despite radically different schemes for determining the computational mesh, and different equation of state and opacity models, the HYDRA and FLASH codes yield results that are in excellent agreement for both problems and with the experimental data for the latter. Having validated the FLASH code in this way, we use the code to further investigate the formation of the jet seen in the Grava et al. (2008) experiment and discuss its relation to the Wan et al. (1997) experiment at the NOVA laser.
△ Less
Submitted 6 June, 2013;
originally announced June 2013.
-
Keeping It Real: Revisiting a Real-Space Approach to Running Ensembles of Cosmological N-body Simulations
Authors:
Chris Orban
Abstract:
In setting up initial conditions for ensembles of cosmological N-body simulations there are, fundamentally, two choices: either maximizing the correspondence of the initial density field to the assumed fourier-space clustering or, instead, matching to real-space statistics and allowing the DC mode (i.e. overdensity) to vary from box to box as it would in the real universe. As a stringent test of b…
▽ More
In setting up initial conditions for ensembles of cosmological N-body simulations there are, fundamentally, two choices: either maximizing the correspondence of the initial density field to the assumed fourier-space clustering or, instead, matching to real-space statistics and allowing the DC mode (i.e. overdensity) to vary from box to box as it would in the real universe. As a stringent test of both approaches, I perform ensembles of simulations using power law and a "powerlaw times a bump" model inspired by baryon acoustic oscillations (BAO), exploiting the self-similarity of these initial conditions to quantify the accuracy of the matter-matter two-point correlation results. The real-space method, which was originally proposed by Pen 1997 and implemented by Sirko 2005, performed well in producing the expected self-similar behavior and corroborated the non-linear evolution of the BAO feature observed in conventional simulations, even in the strongly-clustered regime (sigma8 >= 1). In revisiting the real-space method championed by Sirko 2005, it was also noticed that this earlier study overlooked an important integral constraint correction to the correlation function in results from the conventional approach that can be important in LambdaCDM simulations with Lbox =< 1 Gpc/h and on scales r >= Lbox / 10. Rectifying this shows that the fourier space and real space methods are about equally accurate and efficient for modeling the evolution and growth of the correlation function, contrary to previous claims. An appendix provides a useful independent-of-epoch analytic formula for estimating the importance of the integral constraint bias on correlation function measurements in LambdaCDM simulations.
△ Less
Submitted 17 April, 2013; v1 submitted 10 January, 2012;
originally announced January 2012.
-
Self-Similar Bumps and Wiggles: Isolating the Evolution of the BAO Peak with Power-law Initial Conditions
Authors:
Chris Orban,
David H. Weinberg
Abstract:
Motivated by cosmological surveys that demand accurate theoretical modeling of the baryon acoustic oscillation (BAO) feature in galaxy clustering, we analyze N-body simulations in which a BAO-like gaussian bump modulates the linear theory correlation function ξ_L(r)=(r_0/r)^{n+3} of an underlying self-similar model with initial power spectrum P(k)=A k^n. These simulations test physical and analyti…
▽ More
Motivated by cosmological surveys that demand accurate theoretical modeling of the baryon acoustic oscillation (BAO) feature in galaxy clustering, we analyze N-body simulations in which a BAO-like gaussian bump modulates the linear theory correlation function ξ_L(r)=(r_0/r)^{n+3} of an underlying self-similar model with initial power spectrum P(k)=A k^n. These simulations test physical and analytic descriptions of BAO evolution far beyond the range of most studies, since we consider a range of underlying power spectra (n=-0.5, -1, -1.5) and evolve simulations to large effective correlation amplitudes (equivalent to σ_8=4-12 for r_bao = 100 Mpc/h). In all cases, non-linear evolution flattens and broadens the BAO bump in ξ(r) while approximately preserving its area. This evolution resembles a "diffusion" process in which the bump width σ_bao is the quadrature sum of the linear theory width and a length proportional to the rms relative displacement Σ_pair(r_bao}) of particle pairs separated by r_bao. For n=-0.5 and n=-1, we find no detectable shift of the location of the BAO peak, but the peak in the n=-1.5 model shifts steadily to smaller scales, following r_peak}/r_bao = 1-1.08(r_0/r_bao)^{1.5}. The "SimpleRG" perturbation theory scheme and, to a lesser extent, standard 1-loop perturbation theory are fairly successful at explaining the non-linear evolution of the fourier power spectrum of our models. Analytic models also explain why the ξ(r) peak shifts much more for n=-1.5 than for n >= -1, though no ab initio model we have examined reproduces all of our numerical results. Simulations with L_box = 10 r_bao} and L_box = 20 r_bao yield consistent results for ξ(r) at the BAO scale, provided one corrects for the integral constraint imposed by the uniform density box.
△ Less
Submitted 5 July, 2011; v1 submitted 7 January, 2011;
originally announced January 2011.
-
Delving Deeper into the Tumultuous Lives of Galactic Dwarfs: Modeling Star Formation Histories
Authors:
Chris Orban,
Oleg Y. Gnedin,
Daniel R. Weisz,
Evan D. Skillman,
Andrew E. Dolphin,
John A. Holtzman
Abstract:
The paucity of observed dwarf galaxies in the Local Group relative to the abundance of predicted dark matter halos remains one of the greatest puzzles of the LCDM paradigm. Solving this puzzle now requires not only matching the numbers of objects but also understanding the details of their star formation histories. We present a summary of such histories derived from the HST data using the color-…
▽ More
The paucity of observed dwarf galaxies in the Local Group relative to the abundance of predicted dark matter halos remains one of the greatest puzzles of the LCDM paradigm. Solving this puzzle now requires not only matching the numbers of objects but also understanding the details of their star formation histories. We present a summary of such histories derived from the HST data using the color-magnitude diagram fitting method. To reduce observational uncertainties, we condense the data into five cumulative parameters -- the fractions of stellar mass formed in the last 1, 2, 5, and 10 Gyr, and the mean stellar age. We interpret the new data with a phenomenological model based on the mass assembly histories of dark matter halos and the Schmidt law of star formation. The model correctly predicts the radial distribution of the dwarfs and the fractions of stars formed in the last 5 and 10 Gyr. However, in order to be consistent with the observations, the model requires a significant amount of recent star formation in the last 2 Gyr. Within the framework of our model, this prolonged star formation can be achieved by adding a stochastic variation in the density threshold of the star formation law. The model results are not sensitive to late gas accretion, the slope of the Schmidt law, or the details of cosmic reionization. A few discrepancies still remain: our model typically predicts too large stellar masses, only a modest population of ultra-faint dwarfs, and a small number of dwarfs with anomalously young stellar populations. Nevertheless, the observed star formation histories of Local Group dwarfs are generally consistent the expected star formation in cold dark matter halos.
△ Less
Submitted 2 July, 2008; v1 submitted 7 May, 2008;
originally announced May 2008.
-
Potential Vorticity Evolution of a Protoplanetary Disk with An Embedded Protoplanet
Authors:
H. Li,
S. Li,
J. Koller,
B. Wendroff,
R. Liska,
C. Orban,
E. Liang,
D. Lin
Abstract:
We present two-dimensional inviscid hydrodynamic simulations of a protoplanetary disk with an embedded planet, emphasizing the evolution of potential vorticity (the ratio of vorticity to density) and its dependence on numerical resolutions. By analyzing the structure of spiral shocks made by the planet, we show that progressive changes of the potential vorticity caused by spiral shocks ultimatel…
▽ More
We present two-dimensional inviscid hydrodynamic simulations of a protoplanetary disk with an embedded planet, emphasizing the evolution of potential vorticity (the ratio of vorticity to density) and its dependence on numerical resolutions. By analyzing the structure of spiral shocks made by the planet, we show that progressive changes of the potential vorticity caused by spiral shocks ultimately lead to the excitation of a secondary instability. We also demonstrate that very high numerical resolution is required to both follow the potential vorticity changes and identify the location where the secondary instability is first excited. Low-resolution results are shown to give the wrong location. We establish the robustness of a secondary instability and its impact on the torque onto the planet. After the saturation of the instability, the disk shows large-scale non-axisymmetry, causing the torque on the planet to oscillate with large amplitude. The impact of the oscillating torque on the protoplanet's migration remains to be investigated.
△ Less
Submitted 17 March, 2005;
originally announced March 2005.