An official website of the United States government
Above Threshold Ionization
A fortran code for calculating the energy spectrum for above threshold ionization and tunneling ionization using the window operator method of Schafer PRA 42, 5794 and Schafer CPC 63, 427. This code uses input files from A. L. Harris doi:10.6084/m9.figshare.21084826 (2022)
more »
« less
- Award ID(s):
- 2207209
- PAR ID:
- 10648963
- Publisher / Repository:
- figshare
- Date Published:
- Subject(s) / Keyword(s):
- Atomic and molecular physics
- Format(s):
- Medium: X Size: 20075 Bytes
- Size(s):
- 20075 Bytes
- Right(s):
- CC BY 4.0
- Sponsoring Org:
- National Science Foundation
More Like this
-
Using a code that employs a self-consistent method for computing the effects of photo-ionization on circumstellar gas dynamics, we model the formation of wind-driven nebulae around massive stars. We take into account changes in stellar properties and mass-loss over the star’s evolution. Our simulations show how various properties, such as the density and ionization fraction, change throughout the evolution of the star. The multi-dimensional simulations reveal the presence of strong ionization front instabilities in the main-sequence phase, similar to those seen in galactic ionization fronts. Hydrodynamic instabilities at the interfaces lead to the formation of filaments and clumps that are continually being stripped off and mixed with the low density interior. Even though the winds start out as completely radial, the spherical symmetry is quickly destroyed, and the shocked wind region is manifestly asymmetrical. The simulations demonstrate that it is important to include the effects of the photoionizing photons from the star, and simulations that do not include this may fail to reproduce the observed density profile and ionization structure of wind-blown bubbles around massive stars.more » « less
-
The ponderomotive force has been suggested to be the main mechanism to produce the so-called first ionization potential (FIP) effect—the enrichment of low-FIP elements observed in the outer solar atmosphere, in the solar wind, and in solar energetic events. It is well known that the ionization of these elements occurs within the chromosphere. Therefore, this phenomenon is intimately tied to the plasma state in the chromosphere and the corona. For this study, we combine IRIS observations, a single-fluid 2.5D radiative magnetohydrodynamics (MHD) model of the solar atmosphere, including ion–neutral interaction effects and nonequilibrium (NEQ) ionization effects, and a novel multifluid multispecies numerical model (based on the Ebysus code). Nonthermal velocities of Si IV measured from IRIS spectra can provide an upper limit for the strength of any high-frequency Alfvén waves. With the single-fluid model, we investigate the possible impact of NEQ ionization within the region where the FIP may occur, as well as the plasma properties in those regions. These models suggest that regions with strongly enhanced network and type II spicules are possible sites of large ponderomotive forces. We use the plasma properties of the single-fluid MHD model and the IRIS observations to initialize our multifluid models to investigate the multifluid effects on the ponderomotive force associated with Alfvén waves. Our multifluid analysis reveals that collisions and NEQ ionization effects dramatically impact the behavior of the ponderomotive force in the chromosphere, and existing theories may need to be revisited.more » « less
-
We present the first-of-its-kind coupling of a continuum full- f gyrokinetic turbulence model with a 6D continuum model for kinetic neutrals, carried out using the Gkeyll code. Our objective is to improve the first-principle understanding of the role of neutrals in plasma fueling, detachment, and their interaction with edge plasma profiles and turbulence statistics. Our model includes only atomic hydrogen and incorporates electron-impact ionization, charge exchange, and wall recycling. These features have been successfully verified with analytical predictions and benchmarked with the DEGAS2 Monte Carlo neutral code. We carry out simulations for a scrape-off layer (SOL) with simplified geometry and National Spherical Torus Experiment parameters. We compare these results to a baseline simulation without neutrals and find that neutral interactions reduce the normalized density fluctuation levels and associated skewness and kurtosis, while increasing auto-correlation times. A flatter density profile is also observed, similar to the SOL density shoulder formation in experimental scenarios with high fueling.more » « less
-
null (Ed.)ABSTRACT The extreme ultraviolet region (EUV) provides most of the ionization that creates the high equivalent width (EW) broad and narrow emission lines (BELs and NELs) of quasars. Spectra of hypermassive Schwarzschild black holes (HMBHs; MBH ≥ 1010 M⊙) with α-discs, decline rapidly in the EUV suggesting much lower EWs. Model spectra for BHs of mass 106–1012 M⊙ and accretion rates 0.03 ≤ Lbol/LEdd ≤ 1.0 were input to the cloudy photoionization code. BELs become ∼100 times weaker in EW from MBH ∼ 108 M⊙ to MBH ∼ 1010 M⊙. The high-ionization BELs (O vi 1034 Å, C iv 1549 Å, and He ii 1640 Å) decline in EW from MBH ≥ 106 M⊙, reproducing the Baldwin effect, but regain EW for MBH ≥ 1010 M⊙. The low-ionization lines (Mg ii 2798 Å, H β 4861 Å, and H α 6563 Å) remain weak. Lines for maximally spinning HMBHs behave similarly. Line ratio diagrams for the BELs show that high O vi/H β and low C iv/H α may pick out HMBH, although O vi is often hard to observe. In NEL BPT diagrams, HMBHs lie among star-forming regions, except for highly spinning, high accretion rate HMBHs. In summary, the BELs expected from HMBHs would be hard to detect using the current optical facilities. From 100 to 1012 M⊙, the emission lines used to detect active galactic nuclei (AGNs) only have high EW in the 106–109 M⊙ window, where most AGNs are found. This selection effect may be distorting reported distributions of MBH.more » « less
