Search NASA⌕ Search

SEARCH · Search NASA

Results for “COLLISION PARAMETER”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 271 records · Page 15

Simulations of classical three-body thermalization in one dimension

One-dimensional systems, such as nanowires or electrons moving along strong magnetic field lines, have peculiar thermalization physics. The binary collision of pointlike particles, typically the dominant process for reaching thermal equilibrium in higher-dimensional systems, cannot thermalize a 1D system. We study how dilute classical 1D gases thermalize through three-body collisions. We consider a system of identical classical point particles with pairwise repulsive inverse power-law potential V ij ∝ 1/|x i –x j | n or the pairwise Lennard-Jones potential. Using Monte Carlo methods, we compute a collision kernel and use it in the Boltzmann equation to evolve a perturbed thermal state with temperature T toward equilibrium. We explain the shape of the kernel and its dependence on the system parameters. Additionally, we implement molecular dynamics simulations of a many-body gas and show agreement with the Boltzmann evolution in the low-density limit. For the inverse power-law potential, the rate of thermalization is proportional to ρ 2 ⁢T$\frac{1}{2}$ – $\frac{1}{n}$, where ρ is the number density. Furthermore, the corresponding proportionality constant decreases with increasing n.

1-dimensional systems↗

Flavour-dependent chemical freeze-out of light nuclei in relativistic heavy-ion collisions

We study the production of light nuclei in Au+Au collisions at $\sqrt{s_{NN}}$ = 7.7–200 GeV and Pb+Pb collisions at $\sqrt{s_{NN}}$ = 2.76 and 5.02 TeV within a flavour-dependent chemical freeze-out scenario, assuming different flavoured hadrons undergo separate chemical freeze-out. Using the Thermal-FIST package, thermal parameters extracted from fits to various sets of hadron yields, including and excluding light nuclei, are used to calculate the ratios of the yields of light nuclei, namely, d/p, $\overline{d}$/$\overline{p}$, t/p, t/d, 4 He/ 3 He, and 3 H Λ / 3 He. A comparison with experimental data from the STAR and ALICE collaborations shows that a sequential freeze-out scenario provides a better description of light nuclei yield ratios than the traditional single freeze-out approach. These results suggest the flavour-dependent chemical freeze-out for final state light nuclei production persists in heavy-ion collisions at both RHIC and LHC energies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Investigating charm quark energy loss in medium with the nuclear modification factor of D 0 -tagged jets

The nuclear modification factor R AA of charm jets, identified by the presence of a D 0 meson among the jet constituents, has been measured for the first time in Pb–Pb collisions at a centre-of-mass energy per nucleon pair $\sqrt{s_{NN}} = 5.02$ TeV with the ALICE detector at the LHC. The D 0 mesons and their charge conjugates are reconstructed from the hadronic decay D 0 → K – π + . Jets are reconstructed from D 0 -meson candidates and charged particles using the anti-k T algorithm with jet resolution parameter R = 0.3, in the jet transverse momentum (p T ) range 5 < $p^{ch}_{T}$$^{jet}$ < 50 GeV/c and pseudorapidity |η ch jet| < 0.6. A hint of reduced suppression in the charm-jet R AA is observed in comparison to inclusive jets in central Pb–Pb collisions with a significance of about 2σ in 20 < $p^{ch}_{T}$$^{jet}$ < 50 GeV/c, suggesting the in-medium energy loss to depend on both the difference between quark and gluon coupling strength (Casimir colour-charge effect) and quark mass (dead-cone effect). The data are compared with model calculations that include mass effects in the in-medium energy loss. Several state-of-the-art models are consistent with the data, with the LIDO model providing the best description of the data in the common kinematic range of inclusive and D 0 -tagged jets, highlighting the role of mass effects in interpreting the results.

Acharya, S. (ORCID:0000000292135329)↗

Identified charged hadron production in Au+Au collisions at $\sqrt{s_{NN}}$ = 54.4 GeV with the STAR detector

Here, we present results on the production of 𝜋±, 𝐾±, 𝑝, and $\bar{𝑝}$ in Au+Au collisions at $\sqrt{𝑠_{N⁢N}}$ = 54.4~GeV using the STAR detector at RHIC, at mid-rapidity (|𝑦|< 0.1). Invariant yields of these particles as a function of transverse momentum are shown. We determine bulk properties such as integrated particle yields (𝑑 ⁢𝑁 /𝑑 ⁢𝑦 ), mean transverse momentum (⟨𝑝 𝑇 ⟩), particle ratios, which provide insight into the particle production mechanisms. Additionally, the kinetic freeze-out parameters (𝑇 kin and ⟨𝛽 𝑇 ⟩), which provide information about the dynamics of the system at the time of freeze-out, are obtained. The Bjorken energy density (𝜖 𝐵⁢𝐽 ), which gives an estimate of the energy density in the central rapidity region of the collision zone at the formation time 𝜏, is calculated and presented as a function of multiplicity for various energies. The results are compared with those from the models such as A Multi-Phase Transport (AMPT) and Heavy Ion Jet INteraction Generator (HIJING) for further insights.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Beam energy dependence of net-hyperon yield and its implication on baryon transport mechanism

In the constituent quark model, each quark inside a baryon carries 1/3 unit of the baryon number. An alternative picture exists where the center of a Y-shaped topology of gluon fields, called the baryon junction, carries a unit baryon number. Studying baryon transport over a large rapidity gap (δy) in nuclear collisions provides a possible tool to distinguish these two pictures. A recent analysis of global data on net-proton yield at mid-rapidity in Au+Au collisions showed an exponential dependence on δy and the exponential slope does not vary with event centrality, favoring the baryon junction picture. Since junctions are flavor blind, hyperons – baryons containing valence strange quarks – are expected to exhibit a similar behavior as the proton. This study aims to test this prediction by analyzing hyperon yields in Au+Au collisions at various energies. We observe that net-hyperon yields, after correcting for the strangeness production suppression, adhere to the expected exponential form. The extracted slope parameters for net-Λ, net-$Ξ$ and net-Ω are consistent with each other and with those of net-proton within uncertainties, and exhibit no centrality dependence. Various implementations of the PYTHIA event generator, primarily based on valence quarks for baryon transport, are unable to simultaneously describe the slope parameters for all baryons.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Electric field determination in transient plasmas: in situ & non-invasive methods

Abstract One of the primary basic plasma parameters within transient nonequilibrium plasmas is the reduced electric field strength, roughly understood as the ratio of the electrical energy given to the charged species between two collisions. While physical probes have historically been used for electric field measurements, recent advances in high intensity lasers and sensitive detection methods have allowed for non-invasive optical electric field determination in nearly any discharge configuration with time-resolution up to the sub-nanosecond range and sub-millimeter spatial resolution. This topical review serves to highlight several non-invasive methods for in situ electric field strength determination in transient plasmas ranging from high vacuum environments to atmospheric pressure and above. We will discuss the advantages and proper implementation of (i) laser induced fluorescence dip spectroscopy for measurements in low pressure RF discharges, (ii) optical emission spectroscopy based methods for nitrogen, helium or hydrogen containing discharges, (iii) electric field induced coherent Raman scattering, and (iv) electric field induced second harmonic generation. The physical mechanism for each method will be described as well as basic implementation and highlighting recent results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

First measurement of 𝐷 𝑠⁢1 ⁡(1 + )⁢(2536) + and 𝐷$^{*}_{𝑠⁢2}$⁡(2 + )⁢(2573) + production in proton-proton collisions at $\sqrt{s}$ = 13 TeV at the LHC

The production yields of the orbitally excited charm-strange mesons D s⁢1 ⁡(1 + )⁢(2536)+ and 𝐷$^{*}_{𝑠⁢2}$⁡(2 + )⁢(2573) + were measured for the first time in proton-proton (pp) collisions at a center-of-mass energy of $\sqrt{s}$ =13 TeV with the ALICE experiment at the LHC. The D$^{+}_{s⁢1}$ and D$^{*+}_{s⁢2}$ mesons were measured at midrapidity (|𝑦| < 0.5) in minimum-bias and high-multiplicity pp collisions in the transverse-momentum interval 2< 𝑝 T < 24 GeV/𝑐. Their production yields relative to the D$^{+}_{s}$ ground-state yield were found to be compatible between minimum-bias and high-multiplicity collisions, as well as with previous measurements in e ± ⁢p and e + ⁢e − collisions. The measured D$^{+}_{s⁢1}$/D$^{+}_{s}$ and D$^{*+}_{s⁢2}$/D$^{+}_{s}$ yield ratios are described by statistical hadronization models and can be used to tune the parameters governing the production of excited charm-strange hadrons in Monte Carlo generators, such as pythia 8.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Measurement of substructure-dependent jet suppression in Pb + Pb collisions at 5.02 TeV with the ATLAS detector

The ATLAS detector at the Large Hadron Collider has been used to measure jet substructure modification and suppression in Pb + Pb collisions at a nucleon–nucleon center-of-mass energy $\sqrt{{^S}_{NN}}$ = 5.02 TeV in comparison with proton–proton (pp) collisions at $\sqrt{{^S}}$ = 5.02 TeV. The Pb + Pb data, collected in 2018, have an integrated luminosity of 1.72 nb -1 , while the pp data, collected in 2017, have an integrated luminosity of 260 pb -1 . Jets used in this analysis are clustered using the anti- k t algorithm with a radius parameter R = 0.4. The jet constituents, defined by both tracking and calorimeter information, are used to determine the angular scale r g of the first hard splitting inside the jet by reclustering them using the Cambridge–Aachen algorithm and employing the soft-drop grooming technique. The nuclear modification factor, R AA , used to characterize jet suppression in Pb + Pb collisions, is presented differentially in r g , jet transverse momentum, and in intervals of collision centrality. The R AA value is observed to depend significantly on jet r g . Jets produced with the largest measured r g are found to be twice as suppressed as those with the smallest r g in central Pb + Pb collisions. The R AA values do not exhibit a strong variation with jet p T in any of the r g intervals. The r g and p T dependence of jet R AA is qualitatively consistent with a picture of jet quenching arising from coherence and provides the most direct evidence in support of this approach.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Silica nanocluster binding rate coefficients from molecular dynamics trajectory calculations

Oxide nanoparticle growth from vapor phase precursors occurs in high temperature aerosol reactors first via the formation of nanoclusters, which are nanometer-scale condensed-phase species composed of 10 1 -10 2 atoms. The binding rate for nanoclusters, defined as the rate at which two nanoclusters collide with and stick to one another, is hence of interest in predicting nanoparticle size distribution functions in an aerosol. We have utilized molecular dynamics (MD) trajectory calculations to determine the homogeneous (equal-sized) and heterogeneous (disparate-sized) binding rate coefficients of SiO2 (silica) nanoclusters composed of 18, 144, and 333 atoms. MD trajectory calculations incorporated all-atom models of nanoclusters using a combined Born–Huggins–Mayer-Lennard-Jones potential model, which accounts for both short range interactions and electrostatic interactions. MD calculations were utilized to determine the binding probability as a function of both initial relative velocity and impact parameter; integration of the binding probability across impact parameter and relative velocity space yields the binding rate coefficient. MD trajectory calculations reveal that the most common type of encounters between nanoclusters leading to binding are grazing collisions, i.e. instances where collision would not occur without induced dipole potential influences. The resulting binding rate coefficients are found to be extremely weakly dependent on system temperature, which is in contrast to the use of rate coefficient models which are the product of a hard-sphere collision rate coefficient and a constant enhancement factor (leading to a rate coefficient proportional to the square root of temperature). Enhancement factors defined from MD trajectory calculations fall in the range of 3–9 as system temperature decreases from 1500 K to 300 K. Furthermore, such large values suggest that potential interactions need to be considered when modeling oxide nanocluster growth in the gas phase.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Morphology of shocked lateral outflows in colliding hydrodynamic flows

Supersonic interacting flows occurring in phenomena, such as protostellar jets, give rise to strong shocks and have been demonstrated in several laboratory experiments. To study such colliding flows, we use the AstroBEAR AMR code to conduct hydrodynamic simulations in three dimensions. Here, we introduce variations in the flow parameters of density, velocity, and cross-sectional radius of the colliding flows in order to study the propagation and conical shape of the bow shock formed by collisions between two, not necessarily symmetric, hypersonic flows. We find that the motion of the interaction region is driven by imbalances in ram pressure between the two flows, while the conical structure of the bow shock is a result of shocked lateral outflows being deflected from the horizontal when the flows are of differing cross sections.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A Rate Theory Model of Radiation-Induced Swelling in an Austenitic Stainless Steel

Many rate theory models of cavity (void) swelling have been published over the past 50 years, all having the same, or similar, structures. A rigorous validation of the models has not been possible because of the dearth of information concerning the microstructures that correspond with the swelling data. Whereas the lack of microstructure information is still an issue for historical swelling data, in the past 10–20 years data have been published on the evolution of the microstructure (point defect yields from collision cascades, cavity number densities, and dislocation densities/yield strengths) allowing certain gaps in information to be filled when considering historic swelling data. With reasonable estimates of key microstructure parameters, a standard rate theory model can be applied, and the model parameter space explored, in connection with historical swelling data. By using published data on: (i) yield strength as a function of dose and temperature (to establish an empirical expression for dislocation density evolution); (ii) cavity number densities as a function of temperature; and (iii) freely migrating defect (FMD) production as a function of primary knock-on atom (PKA) spectrum, the necessary parameter and microstructure inputs that were previously unknown can be used in model development. This paper describes a rate-theory model for void swelling of 316 stainless steel irradiated in the EBR-2 reactor as a function of irradiation temperature and neutron dose.

36 MATERIALS SCIENCE↗

Transition characteristics and electron kinetics in microhollow cathode discharges

We demonstrate the transition characteristics and electron kinetics of microdischarges in a microgap with a cathode having microstructures using two-dimensional particle-in-cell/Monte Carlo collision (2d3v) simulations. It is found that the microdischarge is inside the hollow cavity at lower pressures, forming microhollow cathode discharges (MHCDs), while the dominant discharge moves outside the cavity at higher pressures, locating above the cathode rectangular protrusion tip. The spatial distributions of the microdischarge parameters (e.g., electric potential and electron density) are presented by capturing the transition characteristics. The electron kinetics of the microdischarges are examined based on the moment analysis of the electron Boltzmann equation. The collisional and non-collisional components of the electron power absorption are decomposed, which were seldomly investigated for microdischarges previously. The results provide fundamental insights into MHCD formation with a structured electrode at varying pressure conditions, which could be beneficial for the design and fabrication of microplasma devices in practical applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Probability Density Function for the spatial and intensity distribution of neutron-induced defects in Silicon

The ability to model semiconductor device degradation under neutron irradiation depends upon having a robust modeling capability for the neutron-induced collision cascades as well as a means to analytically fit the resulting probability distributions of defect production and ionizing energy deposition for purposes of extrapolation to low-probability, high-consequence scenarios. In this paper, the widely-utilized binary collision approximation codes MARLOWE and SRIM are deployed in conjunction with a critical examination of their parameterizations as benchmarked against higher-fidelity molecular dynamics simulations. A simple 3-parameter form described by the Generalized Logistic Distribution is shown to be a good fit to Frenkel pair and ionization intensity distributions in bulk silicon. The BCA codes are then applied to simulate cascades in 5 nm layers of a representative gate-all-around nanosheet transistor, where joint probability distributions of threshold levels of damage to multiple layers are evaluated.

36 MATERIALS SCIENCE↗

Turbulent Regimes in Collisions of 3D Alfvén-wave Packets

Using three-dimensional gyrofluid simulations, we revisit the problem of Alfvén-wave (AW) collisions as building blocks of the Alfvénic turbulent cascade and their interplay with magnetic reconnection at magnetohydrodynamic (MHD) scales. Depending on the large-scale value of the nonlinearity parameter χ 0 (the ratio between the AW linear propagation time and nonlinear turnover time), different regimes are observed. For strong nonlinearities (χ 0 ~ 1), turbulence is consistent with a dynamically aligned, critically balanced cascade—fluctuations exhibit a scale-dependent alignment $\sin {\theta }_{{k}_{\perp }}\propto {k}_{\perp }^{-1/4}$, resulting in a ${k}_{\perp }^{-3/2}$ spectrum and ${k}_{\parallel }\propto {k}_{\perp }^{1/2}$ spectral anisotropy. At weaker nonlinearities (small χ 0 ), a spectral break marking the transition between a large-scale weak regime and a small-scale ${k}_{\perp }^{-11/5}$ tearing-mediated range emerges, implying that dynamic alignment occurs also for weak nonlinearities. At χ 0 < 1 the alignment angle ${\theta }_{{k}_{\perp }}$ shows a stronger scale dependence than in the χ 0 ~ 1 regime, namely $\sin {\theta }_{{k}_{\perp }}\propto {k}_{\perp }^{-1/2}$ at χ 0 ~ 0.5, and $\sin {\theta }_{{k}_{\perp }}\propto {k}_{\perp }^{-1}$ at χ 0 ~ 0.1. Dynamic alignment in the weak regime also modifies the large-scale spectrum, scaling approximately as ${k}_{\perp }^{-3/2}$ for χ 0 ~ 0.5 and as ${k}_{\perp }^{-1}$ for χ 0 ~ 0.1. A phenomenological theory of dynamically aligned turbulence at weak nonlinearities that can explain these spectra and the transition to the tearing-mediated regime is provided; at small χ 0 , the strong scale dependence of the alignment angle combines with the increased lifetime of turbulent eddies to allow tearing to onset and mediate the cascade at scales that can be larger than those predicted for a critically balanced cascade by several orders of magnitude. Such a transition to tearing-mediated turbulence may even supplant the usual weak-to-strong transition.

79 ASTRONOMY AND ASTROPHYSICS↗

Measurement of the angle between jet axes in Pb-Pb collisions at $\sqrt{s_{NN}}$ =5.02 TeV

This paper presents the first measurement of the angle between different jet axes (denoted as Δ⁢𝑅 axis ) in Pb-Pb collisions. The measurement is carried out in the 0–10 % most-central events at $\sqrt{s_{NN}}$ = 5.02 TeV. Jets are assembled by clustering charged particles at midrapidity using the anti-𝑘 T algorithm with resolution parameters 𝑅 = 0.2 and 0.4 and transverse momenta in the intervals 40 < 𝑝$^{ch jet}_{T}$ < 140 GeV/𝑐 and 80 < 𝑝$^{ch jet}_{T}$ < 140 GeV/𝑐, respectively. Measurements at these low transverse momenta enhance the sensitivity to quark-gluon plasma (QGP) effects. A comparison to models implementing various mechanisms of jet energy loss in the QGP suggests that the observed narrowing of the Pb-Pb distribution relative to 𝑝𝑝 can be explained if quark-initiated jets are more likely to emerge from the medium than gluon-initiated jets. These new measurements disfavor intrajet 𝑝 T broadening as implemented in a simple model calculation with the Baier-Dokshitzer-Mueller-Peigne-Schiff formalism for energy loss in the QGP. The comparison of Pb-Pb and 𝑝𝑝 collisions shows sensitivity to the angular scale at which the QGP can resolve two independent splittings, favoring mechanisms that incorporate incoherent energy loss.

Acharya, S. [LPC] (ORCID:0000000292135329)↗

Measurement of the CKM angle γ in B ± → DK ± and B ± → Dπ ± decays with D → $ {K}_{\mathrm{S}}^0 $h + h -

A measurement of CP -violating observables is performed using the decays B ± → DK ± and B ± → Dπ ± , where the D meson is reconstructed in one of the self-conjugate three-body final states \( {K}_{\mathrm{S}}^0 \) π + π - and \( {K}_{\mathrm{S}}^0 \) K + K - (commonly denoted \( {K}_{\mathrm{S}}^0 \) h + h - ). The decays are analysed in bins of the D -decay phase space, leading to a measurement that is independent of the modelling of the D -decay amplitude. The observables are inter- preted in terms of the CKM angle γ . Using a data sample corresponding to an integrated luminosity of 9 fb - 1 collected in proton-proton collisions at centre-of mass energies of 7, 8, and 13 TeV with the LHCb experiment, γ is measured to be \( \left({68.7}_{-5.1}^{+5.2}\right){}^{\circ} \) . The hadronic parameters \( {r}_B^{D K},{r}_B^{D\pi},{\delta}_B^{D K},\kern0.5em \mathrm{and}\kern0.5em {\delta}_B^{D\pi} \) , which are the ratios and strong-phase differences of the suppressed and favoured B ± decays, are also reported.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Accessing the strong interaction between Λ baryons and charged kaons with the femtoscopy technique at the LHC

The interaction between baryons and kaons/antikaons is a crucial ingredient for the strangeness $S$ = 0 and $S$ = -2 sector of the meson–baryon interaction at low energies. In particular, the $Λ\overline{K}$ might help in understanding the origin of states such as the $Ξ$(1620), whose nature and properties are still under debate. Experimental data on $Λ–K$ and $Λ–\overline{K}$ systems are scarce, leading to large uncertainties and tension between the available theoretical predictions constrained by such data. In this Letter we present the measurements of $Λ–K$ + ⊕ $\overline{Λ}–K$ - and $Λ–K$ - ⊕ $\overline{Λ}–K$ + correlations obtained in the high multiplicity triggered data sample in pp collisions at $\sqrt{s}$ = 13 TeV recorded by ALICE at the LHC. The correlation function for both pairs is modeled using the Lednický–Lyuboshits analytical formula and the corresponding scattering parameters are extracted. The $Λ–K$ - ⊕ $\overline{Λ}–K$ + correlations show the presence of several structures at relative momenta k above 200 MeV/c, compatible with the $Ω$ baryon, the $Ξ$(1690), and $Ξ$(1820) resonances decaying into $Λ–K$ - pairs. The low $k^*$ region in the $Λ–K$ - ⊕ $\overline{Λ}–K$ + also exhibits the presence of the $Ξ$(1620) state, expected to strongly couple to the measured pair. The presented data allow to access the $ΛK$ + and $ΛK$ - strong interaction with an unprecedented precision and deliver the first experimental observation of the $Ξ$(1620) decaying into $ΛK$ - .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Collaborative Research: Unravelling the Physics Associated with the Production of Extremely Dense Plasma States of Microscale (Final Report) Nanosecond-pulsed Discharges

The aim of this project is to study and establish the physical mechanisms that contribute to the formation of anomalously dense plasmas in high-pressure nanosecond-pulsed discharges. These discharges have a broad range of applications such as plasma-assisted combustion, plasma flow actuators, biomedical sterilization and exotic materials synthesis. The structure and formation of these discharges, producing high plasma densities of ~10 14 -10 15 cm -3 , are well-studied and understood. Fast-pulsed microscale high-pressure discharges can be driven to even higher densities of > 10 19 cm -3 , approaching warm dense matter conditions. The mechanisms that generate these plasmas have not been understood. Analysis of the warm dense matter state under laboratory conditions is an expensive and non-trivial endeavor. For instance, dense plasmas can be generated by electrical explosions of metal foils and wires. Plasmas generated after the explosion have a short lifetime and often present difficult conditions for diagnostics. Generation of dense plasmas was also achieved during high-voltage nanosecond pulsed discharges when the so-called explosive electron emission is obtained. Unfortunately, this process is very difficult to control for the studies of warm dense matter. In our recent study, we have shown that additional heating of plasma by lasers can further increase the density of plasma and even lead to the fully ionized state. This method, potentially, allows better control of the plasma parameters. In this work, we studied a second stage laser-heated micro-discharge using a self-consistent one-dimensional particle-in-cell Monte Carlo-collision (1D PIC-MCC) model coupled with Maxwell’s equations. We predicted the generation of a fully ionized plasma on the picosecond time scale. However, this model considered the plasma as an ideal gas despite the high pressure and the nearly fully ionized state. The ideal plasma model assumes that the dilute gas approximation is valid, where the inter-particle interactions are negligible. For charged particles this assumption holds as long as the shielded Coulomb potential assumption is valid. For very high plasma densities, this concept breaks down since the Debye sphere surrounding each charged particle no longer contains enough electrons to statistically provide the shielding of the single particle Coulomb interaction potential. At such densities, the plasma can no longer be described as ideal and non-ideal coupling effects need to be considered. In this report, we elucidate our recent work of developing a PIC-MCC model with improvements for non-ideal plasma conditions due to Coulomb coupling at high densities. In particular, we study the interaction of green light radiation and a dense microplasma, and explore the non-ideal plasma effects in this interaction. In this computational model, we implement the two most important non-ideal effects: ionization potential depression (IPD) and enhanced collision cross sections. Our primary goal is to study the physics associated with electromagnetic (EM) wave heating, also called the second-stage wave-heating, and establish the role of plasma non-ideality in this phenomenon. Our secondary goal is to improve the chemistry mechanism of the 1D PIC-MCC model by including a more detailed excited species collision treatment. At high pressures, stepwise ionization from excited species might play an important role in the ionization process. Previously, this ionization mechanism was neglected due to the excitation collision cross section of xenon being smaller than that of ionization. However, a preliminary study showed that the excited species density in the initial microplasma was an order of magnitude higher than the electron density. Therefore, my aim is to determine the significance of this additional ionization pathway to the plasma generation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗