SEARCH · Engineering Papers

Results for “PIC”

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.

SourceResearch connectionFrom yearThrough year

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

At least 127 records · Page 7

All sources NASA NTRS DOE OSTI

NASA NTRSDate not supplied

Photonic Integrated Circuit TUned for Reconnaissance and Exploration (PICTURE)

The mid-infrared (MIR) spectral range (3-5 μm) is of particular interest for remotely sensing gaseous molecules such as H2O, CO2, CH4, N2O, CO, NH3, and other compounds. The infrared spectra of planets, moons, comets, and asteroids are rich in information, including gas composition and surface mineralogy. Significant advances in technology have emerged for ground-based telescopes, including the higher spectral resolution permitted by cross-dispersed instruments and heterodyne techniques. New technologies offer the opportunity to break this barrier, by using solid-state photonics. In the past three years, under the NASA ROSES PICASSO program, we have been developing key components for a revolutionary MIR spectrometer: Photonic Integrated Circuit TUned for Reconnaissance and Exploration (PICTURE), which is based on integrated photonics technology that offers ultra-small size, weight, and power (SWaP) with non-moving parts and low cost for future planetary missions. In this paper, we will describe our science and technology development progress leading to demonstrating the concept and functionalities of the PICTURE instrument. The photonic integrated circuit spectrometer (PICS) of the PICTURE instrument uses an integrated heterodyne detection scheme to significantly reduce SWaP and improve sensitivity. Our program goals are to advance the building blocks needed for the PICS, which include arrayed waveguide gratings (AWGs), quantum cascade lasers (QCLs) as local oscillators, and quantum cascade detectors (QCDs) as heterodyne detectors. In addition, we are leveraging a NASA SBIR program to develop the toolsets needed to fabricate MIR photonics lanterns (PL). The PL is a critical component that enables PICTURE to bring a signal from the collecting telescope to the PICS, which requires single optical mode inputs. PICTURE focuses on the CO spectral band at 4.6-4.8 μm that is of key importance to cometary science. We also continue to explore future instrument concepts that exploit the broad wavelength potential of the PICS design to perform spectroscopy spanning the full MIR and longwave IR (LWIR) bands. This will have wide applicability in planetary science, for example, to probe the strongest CO2, H2O, and CH4 transitions that are difficult (CH4) or impossible (CO2, H2O) to detect using Earth-based telescopes due to atmospheric opacity. Each integrated-photonics-spectrometer chip will feature a single heterodyne room-temperature laser with a wide wavelength tuning range, or multiple local-oscillator lasers for much broader spectral coverage. The fully developed PICTURE instrument will provide a cost-effective, high-resolution spectrometer for future space applications.

Anthony W Yu↗

NASA NTRSDate not supplied

Photonic Integrated Circuit TUned for Reconnaissance and Exploration (PICTURE)

Anthony W. Yu↗

NASA NTRSDate not supplied

Multi-Gbps Fiber-Optic Wavefront Sensing Coherent Optical Receiver

We report recent progress on technology developments at the NASA Glenn Research Center using photonic lanterns for coherent optical communications applications. In particular, the development of a spatial mode-diversity optical receiver using a photonic lantern combined with a photonic integrated circuit (PIC). The PIC is designed for compatibility with the NASA Laser Communication Relay Demonstration’s (LCRD) differential phase-shift keying signaling format. We report on the PIC design and future characterization and packaging efforts. Additionally, we report on the development of a 19-channel photonic lantern and provide updates on efforts to use the photonic lantern as a wavefront sensing device. The photonic lantern insertion loss, and spatial-mode transfer characteristics are detailed, as well as the development of a predictive framework in which the intensities in the 19 single-mode output channels are used to reconstruct the wavefront of the light entering the multi-moded input side of the lantern. Finally, we discuss concepts for potential future systems using photonic lanterns and PICs that combine both wavefront sensing and coherent communications onto a single integrated platform.

Optical Communications↗

DOE OSTI2025

Coarse Graining Discrete Element Method Information in Particle-in-Cell Length Scales Using a Machine Learning Approach

This report details the development of a machine learning (ML)-driven framework to coarse-grain inter-particle collision dynamics from high-fidelity Discrete Element Method (DEM) simulations to Particle-in-Cell (PIC) scales for gas-solid systems. Traditional PIC models, while computationally efficient, rely on empirical granular stress formulations that fail to capture the full complexity of collision physics, particularly the heterogeneity in particle dynamics. This study adopts a bottom-up approach, integrating insights from DEM simulations to improve the physical fidelity and interpretability of PIC-scale models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS2017

Simulating PACE Global Ocean Radiances

The NASA PACE mission is a hyper-spectral radiometer planned for launch in the next decade. It is intended to provide new information on ocean biogeochemical constituents by parsing the details of high resolution spectral absorption and scattering. It is the first of its kind for global applications and as such, poses challenges for design and operation. To support pre-launch mission development and assess on-orbit capabilities, the NASA Global Modeling and Assimilation Office has developed a dynamic simulation of global water-leaving radiances, using an ocean model containing multiple ocean phytoplankton groups, particulate detritus, particulate inorganic carbon (PIC), and chromophoric dissolved organic carbon (CDOC) along with optical absorption and scattering processes at 1 nm spectral resolution. The purpose here is to assess the skill of the dynamic model and derived global radiances. Global bias, uncertainty, and correlation are derived using available modern satellite radiances at moderate spectral resolution. Total chlorophyll, PIC, and the absorption coefficient of CDOC (aCDOC), are simultaneously assimilated to improve the fidelity of the optical constituent fields. A 5-year simulation showed statistically significant (P < 0.05) comparisons of chlorophyll (r = 0.869), PIC (r = 0.868), and a CDOC (r =0.890) with satellite data. Additionally, diatoms (r = 0.890), cyanobacteria (r = 0.732), and coccolithophores (r = 0.716) were significantly correlated with in situ data. Global assimilated distributions of optical constituents were coupled with a radiative transfer model (Ocean-Atmosphere Spectral Irradiance Model, OASIM) to estimate normalized water-leaving radiances at 1 nm for the spectral range 250-800 nm. These unassimilated radiances were within 0.074 mW/sq cm/micron/sr of MODIS-Aqua radiances at 412, 443, 488, 531, 547, and 667 nm. This difference represented a bias of 10.4% (model low). A mean correlation of 0.706 (P < 0.05) was found with global distributions of MODIS radiances. These results suggest skill in the global assimilated model and resulting radiances. The reported error characterization suggests that the global dynamical simulation can support some aspects of mission design and analysis. For example, the high spectral resolution of the simulation supports investigations of band selection. The global nature of the radiance representations supports investigations of satellite observing scenarios. Global radiances at bands not available in current and past missions support investigations of mission capability. PACE, ocean color, water-leaving radiances, biogeochemical model, radiative transfer model

PIC↗

DOE OSTI2024

A Particle-in-Cell Method for Plasmas with a Generalized Momentum Formulation, Part II: Enforcing the Lorenz Gauge Condition

In a previous paper Christlieb et al. (A particle-in-cell method for plasmas with a generalized momentum formulation, part I: Model formulation, 2024), we developed a new particle-in-cell (PIC) method for the relativistic Vlasov–Maxwell system in which the electromagnetic fields and the equations of motion for the particles were cast in terms of scalar and vector potentials through a Hamiltonian formulation. This new method evolved the potentials under the Lorenz gauge using integral equation methods. New methods to construct spatial derivatives of the potentials that converge at the same rates as the fields were also presented. The new particle method was compared against standard explicit discretizations, including the well-known FDTD-PIC method, for a range of applications involving sheaths and particle beams. Here, this paper extends this new class of methods by focusing on the enforcement the Lorenz gauge condition in both exact and approximate forms using co-located meshes. A time-consistency property of the proposed field solver for the vector potential form of Maxwell’s equations is established, which is shown to preserve the equivalence between the semi-discrete Lorenz gauge condition and the analogous semi-discrete continuity equation. Using this property, we present three methods to enforce a semi-discrete gauge condition. The first method introduces an update for the continuity equation that is consistent with the discretization of the Lorenz gauge condition. Both the finite difference and spectral implementations satisfy this discrete gauge condition to machine precision. The second approach we propose enforces a semi-discrete continuity equation using the boundary integral solution to the field equations. The potential benefit of this approach is that it eliminates spatial derivatives that appear on the particle data, namely the current density, which is often calculated by linear combinations of low-order spline basis functions. This method is ideally suited to boundary integral equation methods that invert multi-dimensional operators without dimensional splitting techniques and will be the subject of future work. The third approach introduces a gauge correcting method that makes direct use of the gauge condition to modify the scalar potential and uses local maps for both the charge and current densities. This results in a gauge error, as the maps do not enforce the continuity equation. The vector potential coming from the current density is taken to be exact, and using the Lorenz gauge, we compute a correction to the scalar potential that makes the two potentials satisfy the gauge condition. This method also enforces the gauge condition to machine precision. We demonstrate two of the proposed methods in the context of periodic domains. Problems defined on bounded domains, including those with complex geometric features remain an ongoing effort. However, this work shows that it is possible to design computationally efficient methods that can effectively enforce the Lorenz gauge condition in a non-staggered PIC formulation.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2026

Machine learning surrogates for ion energy–angle distributions in thermal and RF plasma sheaths

Ion energy–angle distributions (IEADs) at material surfaces are a critical input for plasma–material interaction (PMI) studies in fusion devices, yet they are computationally expensive to obtain using particle-in-cell (PIC) simulations. In this work, we develop a machine learning surrogate based on a deep deconvolutional neural network (DDeCNN) trained on large databases generated with the hPIC2 code. The surrogate is capable of reconstructing IEADs from sheath parameters for both thermal and radio-frequency (RF) plasmas, including cases with multiple ion species. Across thousands of test cases, the model achieves high accuracy, with over 97 % of predictions classified as good or average based on standard error metrics (MAE, MSE, L2). Even in the more challenging RF and multi-species regimes, the surrogate reliably captures the multi-peak structure of PIC results. Once trained, the surrogate produces IEADs in milliseconds on a common workstation, yielding speedups of six to seven orders of magnitude compared with running a full PIC simulation. This computational gain enables dense parameter scans and direct coupling of IEAD predictions with PMI and erosion models on whole-device scales in fusion-relevant conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Calculating the space-charge-limited current density for nonplanar geometries by simulating the charge-free electric field

Calculating the space-charge-limited-current density (SCLCD) for a complicated diode geometry often requires computationally expensive particle-in-cell (PIC) codes. Here, this paper addresses this issue by using the charge-free electric field $E_0$ calculated using COMSOL Multiphysics to determine local and global SCLCD. The SCLCD obtained by using the surface average of $|E_0|^2$ on the cathode recovers theoretical results for one-dimensional (1D) planar, cylindrical, and tip-to-tip geometries in appropriate limits. We next compared tip-to-tip calculations with the SCLCD obtained using the PIC code Empire. The SCLCD calculated using COMSOL agreed well with Empire for flatter 1D tip-to-tip geometries and diverged with increasing sharpness. Physically, Empire predicts lower SCLCD than COMSOL because the electrons spread due to concentrated space-charge at the tip, whereas theory assumes that the electrons follow the charge-free electric field lines. We further assess the behavior of the SCLCD for tips protruding from the centers of flat, circular plates of various areas. Larger plate areas with constant tip size recover the 1D planar SCLCD globally and 1D tip-to-tip SCLCD locally, while reducing the difference between Empire and COMSOL calculations since larger plates capture more of the emitted electrons, reducing SCLCD suppression due to beam spreading. These results show that charge-free electric field simulations can be used to determine the SCLCD without needing to simulate particle dynamics in PIC.

Wright, Jack K. [Purdue Univ., West Lafayette, IN ↗

DOE OSTI2025

Modeling betatron radiation using particle-in-cell codes for plasma wakefield accelerator diagnostics

The analysis of plasma wakefield acceleration experimental measurements, particularly in the characterization of photons emitted through the betatron radiation mechanism, requires the development of accurate numerical models. These computational models are crucial for supporting modern instrumentation designed to measure the single-shot, double-differential angular-energy radiation spectra resulting from interactions between beams and plasmas. Motivated by the needs of such applications, this paper presents detailed numerical models of betatron radiation generated in beam-plasma acceleration experiments. These models are based on the integration of the Liénard-Wiechert (LW) potentials, applied to computed particle trajectories. The particle trajectories are generated using three distinct methods: first, by tracking particles through idealized fields in the blowout regime of PWFA; second, by obtaining trajectories using the fast quasistatic particle-in-cell (PIC) code quickpic; and third, obtaining trajectories from the fully self-consistent PIC code osiris. To ensure the accuracy and reliability of these models, the paper includes various benchmark tests using analytical expressions, as well as employing the PIC code epoch, which takes an alternative approach by using a Monte Carlo quantum electrodynamics (QED)-based radiation model. Additionally, the paper presents simulations of the expected experimental betatron radiation spectra, taking into account parameters relevant to PWFA and plasma photocathode experiments at the SLAC FACET-II facility.

Yadav, M. [University of California, Los Angeles, ↗

DOE OSTI2025

Wafer-scale waveguide sidewall roughness scattering loss characterization by image processing

Photonic integrated circuits (PICs) are vital for developing affordable, high-performance optoelectronic devices that can be manufactured at an industrial scale, driving innovation and efficiency in various applications. Optical loss of modes in thin film waveguides and devices is a critical measure of their performance. Thin film growth, lithography, masking, and etching processes are imperfect processes that introduce significant sidewall and top-surface roughness and cause dominating optical losses in waveguides and photonic structures. This roughness, as perturbations couple light from guided to far-field radiation modes, leads to scattering losses that can be estimated from theoretical models. Typically, with UV-based lithography, sidewall roughness is significantly larger than wafer-top surface roughness. Atomic force microscopy (AFM) imaging measurement gives a 3D and high-resolution roughness profile, but the measurement is inconvenient, costly, and unscalable for large-scale PICs and at wafer-scale. Here, we evaluate the sidewall roughness profile based on 2D high-resolution scanning electron microscope (SEM) imaging. We characterized the loss on two homemade nitride and oxide films on 3-inch silicon wafers with 12 waveguide devices on each and correlated the scattering loss estimated from a 2D image-based sidewall profile and theoretical Payne model. The lowest loss of guided fundamental transverse electric (TE 0 ) mode is found at 0.075 dB/cm at 633 nm across 24 devices, a record at visible wavelength. Our work shows 100% success (edge continuity span exceeding 95% of image width/height) in edge detection in image processing of all images to estimate autocorrelation function and optical mode loss. These demonstrations offer valuable insights into waveguide sidewall roughness and a comparison of experimental and 2D SEM image processing based loss estimations with applications in loss characterization at wafer-scale PICs.

Khurana, Mohit (ORCID:0000000207898672)↗

DOE OSTI2025

Scalable low-loss cryogenic packaging of quantum memories in CMOS-foundry processed photonic chips

Optically linked solid-state quantum memories such as color centers in diamond are a promising platform for distributed quantum information processing and networking. Photonic integrated circuits (PICs) have emerged as a crucial enabling technology for these systems, integrating quantum memories with efficient electrical and optical interfaces in a compact and scalable platform. Packaging these hybrid chips into deployable modules while maintaining low optical loss and resiliency to temperature cycling is a central challenge to their practical use. We demonstrate a packaging method for PICs using surface grating couplers and angle-polished fiber arrays that is robust to temperature cycling, offers scalable channel count, applies to a wide variety of PIC platforms and wavelengths, and offers pathways to automated high-throughput packaging. Using this method, we show optically and electrically packaged quantum memory modules integrating all required qubit controls on chip, operating at millikelvin temperatures with <3 dB losses achievable from fiber to quantum memory for the TE 0 mode at a wavelength of 737 nm.

Bernson, Robert [Tyndall National Institute, Cork ↗

DOE OSTI2025

Evaluating Variable-Impedance Magnetically-Insulated Transmission Lines as a Risk-Mitigation Measure for Next-Generation Pulsed Power

This project has produced the first detailed characterizations of power flow resulting from applying the “variable-impedance MITL” concept to real-life systems in Sandia’s pulsed power program (Z and next-generation pulsed power (NGPP)). We present simulation results and analyses for constant-impedance versions of both Z and NGPP and survey the operational viability of several variable-impedance re-designs in the parameter space of linear tapers. Circuit modeling (SCREAMER/Bertha) was used to pinpoint promising candidate designs, and EM-PIC (Empire) simulations were used to evaluate these candidates more rigorously. This approach was particularly successful in the Z regime which resulted in the identification of several viable variable-impedance MITL designs for each level. The approach was more challenged in the operating space NGPP occupies, producing data points that speak to a more restrictive design space due to anode plasma turn-on. In the end, we were able to converge on one viable variable-impedance design for the highest inductance line (level “F”) and one for the highest current line (level “A”). Altogether, the body of simulation evidence presented in this report suggest there does exist flexibility in operating space for magnetically-insulated transmission lines (MITLs) having variable geometric impedance to be a potential enabling technology for safely increasing current delivery (and potentially lowering stack voltage) in pulsed-power drivers by manipulating electron losses; however, operating points for a particular design must be carefully screened. Circuit and EM-PIC modeling provided consistent verdicts in safe operating regimes for operational viability, but additional physics such as anode plasma turn-on which is included in Empire but not in SCREAMER/Bertha was found to be a critical factor affecting power flow that lead to different assessments between the codes. It is not always the case that the occurrence of anode plasma caused a design to fail (some designs turned on anode plasma yet still delivered load currents meeting design targets); the details matter such as how early in the pulse anode surfaces break down (and how large a region). However, in every case that it did fail it was found that the feedback from anode plasma was the cause (i.e., turning off the anode plasma model in Empire restored agreement with the circuit model prediction). As circuit simulations represent an efficient and practical means of surveying design space compared to more computationally-expensive approaches such as EM-PIC, it could be prudent to invest in the research and development of models to include the effects of anode plasma such as ion emission in circuit codes. The variable-impedance MITL design is a new concept that enables controlled manipulation of the initial electron losses in the outer MITL and can be tested on Z today. We encourage follow-on work to explore further optimization (including alternative variable-impedance profiles, e.g., having constant dZ/dR), and to confirm the major findings presented in this report by fielding test hardware on actual Z shots.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2026

Maximum Energy of Particles Accelerated in Gamma-Ray Burst Afterglow Shocks

Particle acceleration in relativistic collisionless shocks remains an open problem in high-energy astrophysics. Particle-in-cell (PIC) simulations predict that electron acceleration in weakly magnetized shocks proceeds via small-angle scattering, leading to a maximum electron energy significantly below the Bohm limit. This upper bound on electron energy manifests observationally as a characteristic synchrotron cutoff, providing a direct probe of the underlying acceleration physics. Gamma-ray burst (GRB) afterglows offer an exceptional laboratory for testing these predictions. Here, we model the spectral evolution of GRB afterglows during the relativistic deceleration phase, incorporating PIC-motivated acceleration prescriptions and self-consistently computing synchrotron and synchrotron self-Compton emission. We find that low-energy bursts in low-density environments, typical of short GRBs, exhibit a pronounced synchrotron cutoff in the GeV band within minutes to hours after the trigger. Applying our framework to GRB 190114C and GRB 130427A, we find that current observations are insufficient to discriminate between PIC-motivated acceleration and the Bohm limit, primarily due to poor photon statistics in the Fermi-Large Area Telescope band. Nevertheless, future MeV–TeV afterglow observations can break model degeneracies and place substantially tighter constraints on the mechanisms responsible for particle acceleration in relativistic shocks. To this end, we simulate a fiducial nearby short GRB as a promising probe of the cutoff location, for which the two acceleration scenarios are cleanly distinguishable and the detection of such an event in the near future remains feasible.

Wu, Zhao-Feng [Purdue University, West Lafayette, ↗

DOE OSTI2026

A Structure-Preserving Decorated Particle Method for the Vlasov-Poisson System

We revisit the Scovel-Weinstein framework (Scovel & Weinstein, CPAM 1994) for reducing the Vlasov-Poisson system while preserving its Hamiltonian structure. Standard particle-in-cell (PIC) algorithms approximate the distribution function by macro-particles with position and velocity. In contrast, Scovel-Weinstein decorated particles involve additional shape degrees of freedom, while maintaining a finite-dimensional reduction with Hamiltonian structure inherited from the continuum model. Although the original work established this structure three decades ago, its computational potential has remained largely unexplored. We present a practical implementation of the Scovel-Weinstein model and compare it with a standard PIC algorithm. Numerical experiments demonstrate that macro-particles in standard PIC can be replaced by far fewer decorated particles while retaining comparable accuracy. This decorated particle approach offers a new structure-preserving paradigm for kinetic plasma simulation.

65M75, 70H05, 70G65↗

NASA NTRS1980

Comparison of the weight loss and adherence of nine different polyimide films thermally aged at 315 C and 350 C in air

Thermal exposure experiments at 315 and 350 C were performed in air on nine different types of polyimides applied to thin 304 stainless steel foils. The tests were conducted to determine which polyimide was the most thermally stable and adherent when subjected to long exposure times at elevated temperatures. One polyimide designated PIC-7 was found to be more thermally stable than the others; however, it did not possess the adherent properties of PIC-2 and PIC-5. It was concluded that as far as thermal stability and adherence are concerned, five of the polyimides are more suitable for high temperature applications than the other four.

Fusaro, R. L.↗

NASA NTRS1986

Prospecting for planets in circumstellar dust - Sifting the evidence from Beta Pictoris

IRAS and near-IR coronographic data for the A5V star Beta Pic are analyzed for evidence of planetary formation. The light scattered from the central star in the system is integrated along the distance to the star to determine the disk magnitude at various distances from the star. A modified gamma distribution is applied to evaluate the scattering coefficient to test a hypothesis that the inner 30 AU of the disk has been swept out by planets. The scattering from the region around Beta Pic is compared with scattering around other A5V stars. Finally, the IR data at 0.89 micron is compared with scattering at the coronograph wavelengths. No evidence is found to support the hypothesis of clearing in the inner disk, although large particle densities can be found very close to Beta Pic. The study illustrates the effectiveness of using scattered light and IR data to discover and characterize matter distributions circumstellar disks.

Diner, D. J.↗

NASA NTRS1990

Simulation of the collapse and dissipation of Langmuir wave packets

Particle-in-cell (PIC) simulations and Zakharov's partial differential equations (PDEs) are used to investigate the collapse of isolated Langmuir wave packets in two dimensions. Collapse thresholds are determined numerically, and the roles of enhanced Langmuir wave damping and nonlinearities not included in the standard Zakharov equations are discussed. The Langmuir wave and ion dynamics in PIC simulations are compared with the predictions of PDE simulations incorporating enhanced Langmuir damping. Electron heating and coherent acceleration in the PIC simulations are discussed and compared with predictions of the transit-time theory.

Newman, D. L.↗

NASA NTRS1993

The circumstellar gas surrounding 51 Ophiuchi - A candidate proto-planetary system similar to Beta Pictoris

Combined archival and recent International Ultraviolet Explorer (IUE) observations of the star, 51 Oph, reveal the presence of variable, accreting gas with velocities as large as + 100 km/s relative to the system. The electron number density of the circumstellar gas is comparable to that observed around the candidate proto-planetary system, Beta Pic. In addition to the cooler gas, absorption from Al III, Si IV, and C IV is present over the velocity range of the accreting gas. The presence of Si IV and C IV in the spectrum of a B9.5 star provides evidence for collisional ionization of the circumstellar gas like that observed in Beta Pic. The combination of H-alpha profiles with double emission peaks to comparable strength, together with detection of transient mass ejection events similar to those observed in other Be stars, suggests that the inclination of the 51 Oph system is within 10-15 deg of the equatorial plane. Collectively these data imply that the 51 Oph system is similar to Beta Pic in both system constituents and orientation, and may be in a similar evolutionary state.

Grady, C. A.↗