Search NASA⌕ Search

SEARCH · Search NASA

Results for “surface ionization waves”

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 37 records · Page 2

Demonstration of full-scale spatiotemporal diagnostics of solid-density plasmas driven by an ultra-short relativistic laser pulse using an X-ray free-electron laser

Understanding the complex plasma dynamics in ultra-intense relativistic laser–solid interactions is of fundamental importance for applications of laser–plasma-based particle accelerators, the creation of high-energy-density matter, understanding planetary science, and laser-driven fusion energy. However, experimental efforts in this regime have been limited by the lack of accessibility of over-critical densities and the poor spatiotemporal resolution of conventional diagnostics. Over the last decade, the advent of femtosecond brilliant hard X-ray free-electron lasers (XFELs) has opened new horizons to overcome these limitations. Here, for the first time, we present full-scale spatiotemporal measurements of solid-density plasma dynamics, including preplasma generation with tens of nanometer scale length driven by the leading edge of a relativistic laser pulse, ultrafast heating and ionization at the main pulse arrival, the laser-driven blast wave, and transient surface return current-induced compression dynamics up to hundreds of picoseconds after interaction. These observations are enabled by utilizing a novel combination of advanced X-ray diagnostics including small-angle X-ray scattering, resonant X-ray emission spectroscopy, and propagation-based X-ray phase-contrast imaging simultaneously at the European XFEL-HED beamline station.

Huang, Lingen [Helmholtz-Zentrum Dresden-Rossendor↗

Strong-field nonsequential double photoionization using virtual-detector theory with path summation

We present an ab initio study of the nonsequential strong-field ionization dynamics of a model two-electron atom with helium character. Single- and double-ionization events are characterized and displayed using detector signals extracted at different points in the two-electron two-dimensional space. The double photoelectron momentum distribution is calculated via coherent path-summation over virtual-particle trajectories. A comparison is made between the momentum distributions obtained with the virtual-detector method, the Schrödinger wave function, and the time-dependent surface flux method developed by Tao and Scrinzi [New J. Phys. 14, 013021 (2012)]. Insights into different ionization and electron recollision pathways are gained from detailed virtual-particle tracking and energy-time readouts. Furthermore, this study demonstrates the extension of virtual-detector theory to strong-field multielectron quantum dynamics and highlights the importance of the evolving quantum phase in quasiclassical electron propagation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Excited state dynamics of cis , cis -1,3-cyclooctadiene: UV pump VUV probe time-resolved photoelectron spectroscopy

We present UV pump, vacuum ultraviolet probe time-resolved photoelectron spectroscopy measurements of the excited state dynamics of cis,cis-1,3-cyclooctadiene. A 4.75 eV deep UV pump pulse launches a vibrational wave packet on the first electronically excited state, and the ensuing dynamics are probed via ionization using a 7.92 eV probe pulse. The experimental results indicate that the wave packet undergoes rapid internal conversion to the ground state in under 100 fs. Comparing the measurements with electronic structure and trajectory surface hopping calculations, we are able to interpret the features in the measured photoelectron spectra in terms of ionization to several states of the molecular cation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Guided plasma jets directed onto wet surfaces: angular dependence and control

The optimal use of atmospheric pressure plasma jets (APPJs) for treatment of surfaces—inorganic, organic and liquid—depends on being able to control the flow of plasma-generated reactive species onto the surface. The typical APPJ is a rare gas mixture (RGM) flowed through a tube to which voltage is applied, producing an RGM plasma plume that extends into the ambient air. The RGM plasma plume is guided by a surrounding shroud of air due to the higher electric field required for an ionization wave (IW) to propagate into the air. Additionally, the mixing of the ambient air with the RGM plasma plume then determines the production of reactive oxygen and nitrogen species (RONS). The APPJ is usually oriented perpendicular to the surface being treated. However, the angle of the APPJ with respect to the surface may be a method to control the production of reactive species to the surface due to the change in APPJ propagation properties and the resulting gas dynamics. In this paper, we discuss results from computational and experimental investigations addressing two points—propagation of IWs in APPJs with and without a guiding gas shroud as a function of angle of the APPJ with respect to the surface; and the use of this angle to control plasma activation of thin water layers. We found that APPJs propagating out of the plasma tube into a same-gas environment lack any of the directional properties of shroud-guided jets, and largely follow electric field lines as the angle of the plasma tube is changed. Guided APPJs propagate coaxially with the tube as the angle is changed, and turn perpendicularly towards the surface only a few mm above the surface. The angle of the APPJ produces different gas dynamic distributions, which enable some degree of control over the content of RONS transferred to thin water layers.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Generation of gamma radiation by a subterawatt ultrashort laser pulse: optimisation of preplasma and pulse duration

We report an experimental and numerical study of the acceleration of electrons in a plasma interacting with a subterawatt laser pulse (intensity of ∼3 × 10{sup 18} W cm{sup −2} at a pulse duration of 50 fs). A preplasma layer on the surface of a molybdenum target is formed by an additional laser pulse with a duration of 8 ns and an intensity of ∼2 × 10{sup 12} W cm{sup −2}. It is shown that an increase in the laser pulse duration to 1700 fs at a constant energy (and a proportional decrease in intensity) leads to an increase in the yield of bremsstrahlung γ-radiation by more than an order of magnitude when the nanosecond pulse is ahead of the femtosecond one by 15 − 25 ns. Interferometry data and results of diagnostics of optical and γ-radiation of a plasma demonstrate that the collisional ionisation of atoms by electrons oscillating in the field of such a laser pulse plays an essential role in the formation of electron density profile. The sensitivity of the described effect to the level of amplified spontaneous emission is determined, despite the nanosecond pulse impact. Numerical simulations show that at a large pulse duration, the acceleration of electrons is stipulated by the beaking of plasma waves excited during stimulated Raman scattering of laser radiation. (paper)

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Molecular surface programming of rectifying junctions between InAs colloidal quantum dot solids

Heavy-metal-free III–V colloidal quantum dots (CQDs) show promise in optoelectronics: Recent advancements in the synthesis of large-diameter indium arsenide (InAs) CQDs provide access to short-wave infrared (IR) wavelengths for three-dimensional ranging and imaging. In early studies, however, we were unable to achieve a rectifying photodiode using CQDs and molybdenum oxide/polymer hole transport layers, as the shallow valence bandedge (5.0 eV) was misaligned with the ionization potentials of the widely used transport layers. This occurred when increasing CQD diameter to decrease the bandgap below 1.1 eV. Here, we develop a rectifying junction among InAs CQD layers, where we use molecular surface modifiers to tune the energy levels of InAs CQDs electrostatically. Previously developed bifunctional dithiol ligands, established for II-VI and IV-VI CQDs, exhibit slow reaction kinetics with III-V surfaces, causing the exchange to fail. Here, we study carboxylate and thiolate binding groups, united with electron-donating free end groups, that shift upward the valence bandedge of InAs CQDs, producing valence band energies as shallow as 4.8 eV. Photophysical studies combined with density functional theory show that carboxylate-based passivants participate in strong bidentate bridging with both In and As on the CQD surface. The tuned CQD layer incorporated into a photodiode structure achieves improved performance with EQE (external quantum efficiency) of 35% (>1 μm) and dark current density < 400 nA cm -2 , a >25% increase in EQE and >90% reduced dark current density compared to the reference device. This work represents an advance over previous III-V CQD short-wavelength IR photodetectors (EQE < 5%, dark current > 10,000 nA cm -2 ).

36 MATERIALS SCIENCE↗

Observation of breakdown wave mechanism in avalanche ionization produced atmospheric plasma generated by a picosecond CO 2 laser

Understanding the formation and long-timescale evolution of atmospheric plasmas produced by ultrashort, long-wavelength infrared (LWIR) pulses is an important but partially understood problem. Of particular interest are plasmas produced in air with a peak laser intensity ~10 12 W/cm 2 , the so-called clamping intensity observed in LWIR atmospheric guiding experiments where tunneling and multi-photon ionization operative at near-IR or shorter wavelengths are inoperative. We find that avalanche breakdown on the surface of aerosol (dust) particles can act to seed the breakdown of air observed above the 200 GW/cm 2 threshold when a train of 3 ps 10.6 μm laser pulses separated by 18 ps is used. The breakdown first appears at the best focus but propagates backward toward the focusing optic as the plasma density approaches critical density and makes forward propagation impossible. The velocity of the backward propagating breakdown can be as high as 10 9 cm/s, an order of magnitude greater than measured with ns pulse-produced breakdown, and can be explained rather well by the so-called breakdown wave mechanism. Transverse plasma expansion with a similar velocity is assisted by UV photoionization and is observed as a secondary longitudinal breakdown mechanism in roughly 10% of the shots. When a cm-size, TW power beam is propagated, interception of aerosol particles is guaranteed and several (40 cm -3 ) breakdown sites appear, each initially producing a near-critical density plasma. On a 10 ns–1 μs timescale, shockwaves from each site expand radially and coalesce to produce a large hot gas channel. The radial velocity of the expansion agrees well with the prediction of the blast wave theory developed for ultrafast atmospheric detonations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Electron beam injection from a hollow cathode plasma into a downstream reactive environment: Characterization of secondary plasma production and Si 3 N 4 and Si etching

A material etching system was developed by combining beam electron injection from a direct current (DC) hollow-cathode (HC) electron source with the downstream reactive environment of a remote CF 4 /O 2 low temperature plasma. The energy of the injected beam electrons is controlled using an acceleration electrode biased positively relative to the HC Argon discharge. For an acceleration voltage greater than the ionization potential of Ar, the extracted primary electrons can produce a secondary plasma in the process chamber. We characterized the properties of the secondary plasma by performing Langmuir probe measurements of the electron energy probability function (EEPF) 2.5 cm below the extraction ring. The data indicate the existence of two major groups of electrons, including electrons with a primary beam electron energy which varies as the acceleration voltage is varied, along with low energy electrons produced by ionization of the Ar gas atoms in the process chamber by the injected beam electrons. When combining the HC Ar beam electron with a remote CF 4 /O 2 electron cyclotron wave resonance (ECWR) plasma, the EEPF of both the low energy plasma electron and beam electron components decreases. Additionally, we studied surface etching of Si3N4 and poly-crystalline Si (poly-Si) thin films as a function of process parameters, including the acceleration voltage (0-70 V), discharge current of the HC discharge (1-2 A), pressure (2- 100 mTorr), source to substrate distance (2.5 to 5 cm), and feed gas composition (with or without CF 4 /O 2 ). The direction of the incident beam electrons was perpendicular to the surface. Si 3 N 4 and poly-crystalline silicon etching are seen and indicate an electronneutral synergy effect. No/little remote plasma spontaneous etching was observed for the conditions used in this study, and the etching is confined to the substrate area irradiated by the injected beam electrons. The electron etched Si 3 N 4 surface etching rate profile distribution is confined within a ~30 mm diameter circle, which is slightly broader than the area for which poly-Si etching is seen, and coincides closely with the spatial profile of beam electrons as determined by the Langmuir probe measurements. The magnitude of poly-Si etching rate is by a factor of two times smaller than the Si 3 N 4 etching rate. In this paper, we discuss possible explanations of the data and the role of surface charging.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The Radar Echo Telescope for Cosmic Rays: Pathfinder experiment for a next-generation neutrino observatory

The Radar Echo Telescope for Cosmic Rays (RET-CR) is a recently funded experiment designed to detect the englacial cascade of a cosmic ray-initiated air shower via in-ice radar, toward the goal of a full-scale, next-generation experiment to detect ultrahigh energy neutrinos in polar ice. For cosmic rays with a primary energy greater than 10 PeV, roughly 10% of an air shower’s energy reaches the surface of a high elevation ice sheet (≳2 ˙km) concentrated into a radius of roughly 10 cm. This penetrating shower core creates an in-ice cascade orders of magnitude more dense than the preceding in-air cascade. Additionally, this dense cascade can be detected via the radar echo technique, where transmitted radio waves are reflected from the ionization deposit left in the wake of the cascade. RET-CR will test the radar echo method in nature, with the in-ice cascade of a cosmic ray-initiated air shower serving as a test beam. In this work, we present the projected event rate and sensitivity based upon a three part simulation using corsika, geant4, and radioscatter. RET-CR expects ~1 radar echo event per day.

79 ASTRONOMY AND ASTROPHYSICS↗

Multi-energy ion irradiation effects and distinctive features of color center generation in yttria-stabilized zirconia driven by structural discrepancy

Regarding a requirement for inert matrix fuel (IMF) to burn minor actinides and plutonium and high-level waste immobilization, Yttria-stabilized zirconia (ZrO 2 -6.5 wt% Y 2 O 3 ) (YSZ) crystals with cubic phase are selected to comprehensively evaluate the irradiation resistance by irradiation-induced structural discrepancies across multiple energy conditions. Electronic excitation and nuclear collisions reveal distinct color centers, influenced by varying component ratios that reflect their formation mechanisms. Key factors driving ultra-fast structural transitions include collision cascades, pressure waves, and energy dissipation from electronic excitation. In this work, we find that in the electronic energy loss (E ele )-dominant region, internal latent tracks with unique surface nanostructures emerge, distinct from the defect clusters seen in the nuclear energy loss (E nuc )-dominant region, while the threshold for latent track formation and melting are identified. In conclusion, structural discrepancies driven by different energy-loss mechanisms promote the generation of singly ionized and nearest-neighbor doubly ionized oxygen vacancies, which respectively dominate the formation of F⁺ and T-centers, resulting in bandgap narrowing (4.80 eV → 4.73 eV), and furthermore triggering visible light emission (2.17 eV → 2.20 eV) in irradiated YSZ crystals, thus, enabling the design of novel irradiation-tailored material functionalities.

36 MATERIALS SCIENCE↗

The achievement of the T e,div feedback control by CD 4 seeding on EAST

A multi-function divertor feedback control system has been built on Experimental Advanced Superconducting Tokamak (EAST) to treat the divertor heat load issue. With the real-time data of the Langmuir probes and the impurity seeding, the divertor electron temperature (T e,div ) is well controlled to achieve the partial detachment phase. The first trial of the CD 4 seeding for the T e,div reduction has been achieved on EAST long-pulse discharge. In the seeding phase, the T e,div was maintained close to 5 eV, and the surface temperature of the target plate (T surface,div ) had a reduction of about 150° C. The plasma stored energy had a reduction in the control phase, so it is necessary to find a way to keep the good plasma confinement in the next step. The CD 4 injection also mitigated the low hybrid wave coupling rate in some degree. The big volume of the CD 4 injection lifted the Greenwald density fraction from ~0.4 to ~0.7, which made the SOL into high recycling state. As a result, most of the injected carbon particles were in the high ionized state, and with the lower of the T e,div , the tungsten line emission was suppressed obviously.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A self-consistent hybrid model of kinetic striations in low-current argon discharges

A self-consistent hybrid model of standing and moving striations was developed for low-current DC discharges in noble gases. We introduced the concept of surface diffusion in phase space ( r , u ) (where u denotes the electron kinetic energy) described by a tensor diffusion in the nonlocal Fokker–Planck kinetic equation for electrons in the collisional plasma. Electrons diffuse along surfaces of constant total energy ε = u - eφ ( r ) between energy jumps in inelastic collisions with atoms. Numerical solutions of the 1d1u kinetic equation for electrons were obtained by two methods and coupled to ion transport and Poisson solver. We studied the dynamics of striation formation in Townsend and glow discharges in argon gas at low discharge currents using a two-level excitation-ionization model and a ‘full-chemistry’ model, which includes stepwise and Penning ionization. Standing striations appeared in Townsend and glow discharges at low currents, and moving striations were obtained for the discharge currents exceeding a critical value. These waves originate at the anode and propagate towards the cathode. We have seen two types of moving striations with the two-level and full-chemistry models, which resemble the s and p striations previously observed in the experiments. Simulations indicate that processes in the anode region could control moving striations in the positive column plasma. The developed model helps clarify the nature of standing and moving striations in DC discharges of noble gases at low discharge currents and low gas pressures.

Physics↗

A NICER View of the Massive Pulsar PSR J0740+6620 Informed by Radio Timing and XMM-Newton Spectroscopy

We report on Bayesian estimation of the radius, mass, and hot surface regions of the massive millisecond pulsar PSR J0740+6620, conditional on pulse-profile modeling of Neutron Star Interior Composition Explorer X-ray Timing Instrument event data. We condition on informative pulsar mass, distance, and orbital inclination priors derived from the joint North American Nanohertz Observatory for Gravitational Waves and Canadian Hydrogen Intensity Mapping Experiment/Pulsar wideband radio timing measurements of Fonseca et al. We use XMM-Newton European Photon Imaging Camera spectroscopic event data to inform our X-ray likelihood function. The prior support of the pulsar radius is truncated at 16 km to ensure coverage of current dense matter models. We assume conservative priors on instrument calibration uncertainty. We constrain the equatorial radius and mass of PSR J0740+6620 to be ${12.39}_{-0.98}^{+1.30}$ km and ${2.072}_{-0.066}^{+0.067}$ M ⊙ respectively, each reported as the posterior credible interval bounded by the 16% and 84% quantiles, conditional on surface hot regions that are non-overlapping spherical caps of fully ionized hydrogen atmosphere with uniform effective temperature; a posteriori, the temperature is ${\mathrm{log}}_{10}(T\,[{\rm{K}}])={5.99}_{-0.06}^{+0.05}$ for each hot region. All software for the X-ray modeling framework is open-source and all data, model, and sample information is publicly available, including analysis notebooks and model modules in the Python language. Furthermore, our marginal likelihood function of mass and equatorial radius is proportional to the marginal joint posterior density of those parameters (within the prior support) and can thus be computed from the posterior samples.

79 ASTRONOMY AND ASTROPHYSICS↗

Few-femtosecond resolved imaging of laser-driven nanoplasma expansion

Abstract The free expansion of a planar plasma surface is a fundamental non-equilibrium process relevant for various fields but as-yet experimentally still difficult to capture. The significance of the associated spatiotemporal plasma motion ranges from astrophysics and controlled fusion to laser machining, surface high-harmonic generation, plasma mirrors, and laser-driven particle acceleration. Here, we show that x-ray coherent diffractive imaging can surpass existing approaches and enables the quantitative real-time analysis of the sudden free expansion of laser-heated nanoplasmas. For laser-ionized SiO 2 nanospheres, we resolve the formation of the emerging nearly self-similar plasma profile evolution and expose the so far inaccessible shell-wise expansion dynamics including the associated startup delay and rarefaction front velocity. Our results establish time-resolved diffractive imaging as an accurate quantitative diagnostic platform for tracing and characterizing plasma expansion and indicate the possibility to resolve various laser-driven processes including shock formation and wave-breaking phenomena with unprecedented resolution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Shake loss intensities in x-ray photoelectron spectroscopy: Theory, experiment, and atomic composition accuracy for MgO and related compounds

The relative intensities of XPS core levels, scaled by their photoionization cross sections, are regularly used to determine sample atomic composition. Cross sections, however, give the intensity to all possible final states for the core ionizations, not just to the main peak. This includes all intrinsic satellite structure (shake states and, for open shell systems, the different ionic multiplets). In practice, for solids, this is usually experimentally impossible to determine accurately because such a satellite structure sits on the inelastically scattered electron background and cannot be easily separated. Therefore, usually, only the intensity of the main peak is used. This limits the ultimate possible accuracy of XPS composition determination. The purpose of the present paper is to examine the contributions that a theoretical analysis of losses of intensity can make to improve quantitation. For an MgO single crystal, we show that the correct stoichiometry of 1:1 can be recovered using the theoretical analysis of the experimental MgO peak ratio intensities. For materials with a sufficient bandgap for the XPS main peaks to be separated from the scattered background, the intensity of main peaks can often be accurately determined. Thus, if one uses theory to calculate that fraction of the total intensity lost from a main peak into all its satellite structure, the intensity of just main peaks could then be used to more accurately determine relative atom % composition. This work tests this approach using a single crystal MgO (50% Mg, 50% O) standard. Ab initio electronic structure theory of representative MgO clusters is used to determine Hartree–Fock wave functions for the ground state and final ionized states corresponding to the main Mg 2p and O1s XPS peaks of the oxide. The sudden approximation, SA, is used to determine the fractional losses from these main peaks to shake satellites, which is found to be greater for O1s than Mg2p. This results in predicted “apparent composition” for stoichiometric MgO of 55.2% Mg, 44.8% O instead of the true 50% Mg, 50% O. Equivalent theory for CaO results in a predicted apparent Ca value of 53.4%. Experimentally, using Mg2s or 2p intensity ratio to O1s, we find values between 52.2% and 56.0% Mg using two crystals and four different instrument electron pass energies. The average value of the measurements is 54.5% Mg when corrected for the presence of an adventitious carbon overlayer and slight surface hydroxide. Though this agreement with theory may be somewhat fortuitous, given the potential experimental errors, which are fully discussed, it is similar to that in our earlier study on LiF. We also present preliminary experimental data on Mg(OH) 2 and MgSO 4 , which show a similar trend of apparently higher than 50% Mg, but we have no theory values. We are not yet able to experimentally test for validation of the difference between apparent composition for MgO (55.2% Mg) and CaO (53.4% Ca), owing to significant carbonate formation at the surface of the single crystal CaO. Here, an important conclusion is that the theoretical determination of shake losses, obtained with ab initio wavefunctions and the SA, is likely to be a useful way to calibrate the accuracy and reliability of compositions obtained from XPS intensities and merits further study.

47 OTHER INSTRUMENTATION↗

Deflagration to Detonation Transition Update: XDDT Code Modularization

A legacy FORTRAN 77 implementation of the Baer–Nunziato two-phase mixture theory for deflagration-to-detonation transition (DDT) in reactive granular materials—hereafter the XDDT (eXplosive DDT) code—has been modularized to Fortran 90 with modular structure, external input files, and adaptive mesh capability. During validation, two code defects were identified and corrected: an inconsistency in the nodal solid pressure evaluation and a nonphysical burn-front tracking criterion. The ignition criterion was also corrected to use the granular surface temperature from the interface heat transfer model, matching the original Baer implementation. An initial attempt to validate against Figure 3 of the original Baer and Nunziato (1986) paper revealed that the code’s detonation velocity on a 201-node mesh (5.5 km/s) was approximately 21% below the expected Chapman–Jouguet value for 70% TMD HMX (∼7 km/s). Validation was redirected to the piston-driven DDT experiments of McAfee et al. (1989), Shot B-9036, for which well-characterized ionization-pin data are available. With the compaction-burn coefficient calibrated to 𝐶 𝛼 = 75, the XDDT code reproduces the DDT transition time to within 0.4% and produces a steady-state detonation velocity within 4% of the McAfee experimental value of 6.36 km/s. The burn model was generalized to support pressure-dependent exponents, enabling application to nitrocellulose-based ball propellants (TS3659) with a cube-root pressure dependence. Validation against the Sandusky/Baer PDC82 piston-impact experiment yielded a reactive wave velocity of 2.3–2.8 km/s, in good agreement with the experimental value of ∼2.2 km/s, and wave coalescence within 5% of the experimental timing. The mathematical model, input parameter requirements, and a roadmap for extending XDDT to PETN with an autocatalytic burn model are presented.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

A thermodynamic analysis of hypervelocity impacts on metals

Hypervelocity impacts on solid surfaces generate and propagate shock waves and create extreme density, pressure, and temperature conditions. Understanding hypervelocity impacts requires a multi-physics approach at the intersection of fluid dynamics, solid mechanics, and quantum theory, which has the capability of spanning multiple different dominating physics regimes. These strong shock conditions found in hypervelocity impacts are studied in the context of high-speed collisions between spacecraft and meteoroids and debris. We study the problem of hypervelocity impacts on metals with the lens of understanding the processes that lead to the creation of an impact generated plasma. In order to characterize the impacting particle as well as its threat to spacecraft from the subsequent radiating plasma an analysis was performed that highlights how the kinetic energy of the impactor is deposited into a target, how the shock wave and shocked material in the target evolve, and how these effects can lead towards the creation of an impact plasma through heating. The analysis suggests that at low velocities the majority of the energy from a hypervelocity impact is stored in reversible elastic energy and does not contribute to strong heating. At higher impact velocities, nucleus and electron thermal components dominate, increasing the energy available for ionization. Overall, this work provides the complete set of state variable initial conditions for any study of hypervelocity impact plasma expansion and its effects.

42 ENGINEERING↗

Cohesion and excitations of diamond-structure silicon by quantum Monte Carlo: Benchmarks and control of systematic biases

We have carried out quantum Monte Carlo (QMC) calculations of silicon crystal focusing on the accuracy and systematic biases that affect the electronic structure characteristics. The results show that 64 and 216 atom supercells provide an excellent consistency for extrapolated energies per atom in the thermodynamic limit for ground, excited, and ionized states. We have calculated the ground state cohesion energy with both systematic and statistical errors below ≈ 0.05 eV. The ground state exhibits a fixed-node error of only 1.3(2)% of the correlation energy, suggesting an unusually high accuracy of the corresponding single-reference trial wave function. We obtain a very good agreement between optical and quasiparticle gaps that affirms the marginal impact of excitonic effects. Our most accurate results for band gaps differ from the experiments by about 0.2 eV. This difference is assigned to a combination of residual finite-size and fixed-node errors. In conclusion, we have estimated the crystal Fermi level referenced to vacuum that enabled us to calculate the edges of valence and conduction bands in agreement with experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗