Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ghost points”

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.

13 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Multiple interaction points in ghost collisions

The Ghost Collider makes use of unique “ghost” bunches, which are electrically neutral combinations of electrons and positrons within the same RF bucket, eliminates the beam-beam effects typically present at the interaction point (IP) in conventional colliders. This allows for the novel possibility of placing multiple interaction regions in series, achieving additive luminosity without introducing significant disruption. However, to get higher luminosity, the beta functions at the IP reaches millimeter scale, which in turn adds significant chromatic contribution to the collider. Correcting these chromatic effects is essential to maintain beam stability and ensure high luminosity during collider operation. By carefully adjusting the phase advance between two IRs that are placed in series, it becomes possible to cancel chromaticity globally, enabling stable collider operation while preserving high luminosity. In this paper we discuss the design of such IR/IRs to be used in a ghost collider.

Accelerator Physics↗
DOE OSTI2019

Fourth Order Finite Difference Methods for the Wave Equation with Mesh Refinement Interfaces

In this work, we analyze two types of summation-by-parts finite difference operators for approximating the second derivative with variable coefficient. The first type uses ghost points, while the second type does not use any ghost points. A previously unexplored relation between the two types of summation-by-parts operators is investigated. By combining them we develop a new fourth order accurate finite difference discretization with hanging nodes on the mesh refinement interface. We take the model problem as the two-dimensional acoustic wave equation in second order form in terms of acoustic pressure, and we prove energy stability for the proposed method. Compared to previous approaches using ghost points, the proposed method leads to a smaller system of linear equations that needs to be solved for the ghost point values. Another attractive feature of the proposed method is that the explicit time step does not need to be reduced relative to the corresponding periodic problem. Numerical experiments, both for smoothly varying and discontinuous material properties, demonstrate that the proposed method converges to fourth order accuracy. A detailed comparison of the accuracy and the time-step restriction with the simultaneous-approximation-term penalty method is also presented. (An erratum is attached.)

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Elastic Wave Propagation in Curvilinear Coordinates with Mesh Refinement Interfaces by a Fourth Order Finite Difference Method

In this work, we develop a fourth order accurate finite difference method for the three dimensional elastic wave equation in isotropic media with the piecewise smooth material property. In our model, the material property can be discontinuous at curved interfaces. The governing equations are discretized in second order form on curvilinear meshes by using a fourth order finite difference operator satisfying a summation-by-parts property. The method is energy stable and high order accurate. The highlight is that mesh sizes can be chosen according to the velocity structure of the material so that computational efficiency is improved. At the mesh refinement interfaces with hanging nodes, physical interface conditions are imposed by using ghost points and interpolation. With a fourth order predictor-corrector time integrator, the fully discrete scheme is energy conserving. Numerical experiments are presented to verify the fourth order convergence rate and the energy conserving property.

58 GEOSCIENCES↗
DOE OSTI2025

A Block-Structured Adaptive Mesh Framework to Solve Radiation Transfer Equation in Irregular Embedded Geometries

Radiation transport arises in various scientific, industrial, and medical fields, and understanding its effect in applications is needed to make accurate predictions, safety assessments and performance optimizations. Solving the Radiation Transport Equation (RTE) is challenging due to its integro-differential nature, which involves both differential and integral terms. The differential term describes the change in radiation intensity due to absorption and emission, while the integral term accounts for scattering. The accurate modeling of radiation is further complicated in many applications due to the complex, irregular geometries. Various methods exist for solving the RTE, including the zonal, Monte Carlo, spherical harmonics, discrete ordinates, and finite volume methods. Traditional mesh-based approaches, which rely on structured or unstructured meshes, struggle with irregular geometries due to: a) the difficulty of conforming structured grids to irregular domains, b) challenges in enforcing boundary conditions correctly, and c) the additional computational cost of unstructured mesh methods. This work presents a second-order accurate method for solving the RTE in irregular geometries. The radiation intensity is discretized using the finite-volume method in both spatial and angular directions on regular Cartesian grid blocks. Leveraging the block-structured adaptive mesh refinement (AMR) framework provided by AMReX, our method refines the grid locally to reduce spatial discretization error, ensuring a converged numerical solution while minimizing computational costs elsewhere. A two-stage deferred correction approach is employed: First, a first-order discretization on grid blocks is solved using an algebraic multigrid method in HYPRE. Second, a correction term is applied explicitly to achieve second-order accuracy. The correction term is calculated by approximating the radiation flux on cell faces using a Total Variation Diminishing (TVD) scheme. This approach ensures quick convergence of the multigrid method while preserving higher-order accuracy of the numerical solution. Irregular geometries are resolved as embedded boundaries (EB), resulting in both cut cells and regular cells. In cut cells, we modify the fluxes using face fractions and incorporate additional contributions from EB boundary conditions. To ensure higher-order convergence near the EB interface, the correction term is modified by interpolating the radiation intensity to fictitious ghost points. The implementation takes advantage of modern supercomputers by leveraging AMReX’sMPI/X parallelization strategy where X can be MPI or a GPU accelerator including CUDA, HIP and DPC++. We validate our solver using classical test cases, both with and without EB, demonstrating accuracy and efficiency. Additionally, we analyze the impact of adaptive mesh refinement on solution accuracy and computational cost, highlighting the advantages of our approach for high-resolution radiation transport simulations.

computational fluid dynamics (CFD)↗
DOE OSTI2026

The Ghost Collider: An Innovative Higgs Factory

The Ghost Collider is an innovative proposal for a 550 GeV centre-of-mass, 275 GeV per beam linear collid-er with four interaction regions, each with the design luminosity. The primary innovation is the use of “ghost bunches” containing equal numbers of electrons and posi-trons so they are electrically neutral. In the linacs, energy is transferred between electrons and positrons in the same bunch, decelerating one type of particle and using the energy to accelerate the other: a new class of Energy Re-covery Linacs. At the interaction points, collisions be-tween two neutral ghost bunches occurs with weakened electromagnetic interaction such as the beam-beam effect or disruption, ensuring that the particles and their energy can be recycled with minimal loss. Four “serial IPs” are incorporated, where chromatic errors produced in one IP are cancelled in the following IP. The neutral beam insta-bility at the IP is also addressed. All interaction points have the nominal luminosity per IP of 2.8 x 1034 cm-2s-1 for a facility luminosity of 1.1 x1035 cm-2s-1 @ 100 MW electrical power for the RF systems. The result is a totally original concept for an electron-positron collider.

Hutton, A. [Thomas Jefferson National Accelerator ↗
DOE OSTI2021

Lorentz symmetry in ghost-free massive gravity

The role of Lorentz symmetry in ghost-free massive gravity is studied, emphasizing features emerging in approximately Minkowski spacetime. The static extrema and saddle points of the potential are determined and their Lorentz properties identified. Solutions preserving Lorentz invariance and ones breaking four of the six Lorentz generators are constructed. Locally, globally, and absolutely stable Lorentz-invariant extrema are found to exist for certain parameter ranges of the potential. Gravitational waves in the linearized theory are investigated. Deviations of the fiducial metric from the Minkowski metric are shown to lead to pentarefringence of the five wave polarizations, which can include superluminal modes and subluminal modes with negative energies in certain observer frames. The Newton limit of ghost-free massive gravity is explored. The propagator is constructed and used to obtain the gravitational potential energy between two point masses. Here, the result extends the Fierz-Pauli limit to include corrections generically breaking both rotation and boost invariance.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Chromaticity compensation of a ghost collider

The GHOST collider final focus system is a pure quadrupole-drift beamline targeting ¿* = 2 mm at four serial interaction points, with peak ¿-functions reaching ~65 km in the final triplet. Beam tracking simulations reveal that a nominal bunch develops a pronounced C-shape in longitu- dinal phase space at IP1, with ¿¿ growing from 0.15 mm to ~2.3 mm — a factor of ~15× increase that directly re- duces luminosity. The mechanism is identified as chromat- ically amplified betatron path length: off-momentum par- ticles acquire enlarged betatron amplitudes in the high-¿ final triplet, generating excess path length via the geometric ¿¿ = - 1/2 ¿ (¿'2 + ¿'2) ¿¿ integral. Within the monoenergetic Balandin framework, ¿2¿2 with ¿ = - 1/2 ¿ ¿ ¿¿ dominates ¿2¿2 by over four million times; the full beam ¿¿ is a fur- ther factor of ~27 larger, driven by ¿¿ -induced chromatic amplitude growth. Phase-advance scans confirm that ¿¿ is insensitive to the apochromatic (¿ ˜ 0) condition. When combined with the geometric hourglass effect, this distor- tion poses a challenge to maximizing luminosity within the current lattice design.

Gamage, B. [Thomas Jefferson National Accelerator ↗
DOE OSTI2020

The convolutional double copy: a case study with a point

The double copy relates scattering amplitudes in gauge and gravity theories. It has also been extended to classical solutions, and a number of approaches have been developed for doing so. One of these involves expressing fields in a variety of (super-)gravity theories in terms of convolutions of gauge fields, including also BRST ghost degrees of freedom that map neatly to their corresponding counterparts in gravity. In this paper, we spell out how to use the convolutional double copy to map gauge and gravity solutions in the manifest Lorenz and de Donder gauges respectively. We then apply this to a particular example, namely the point charge in pure gauge theory. As well as clarifying how to use the convolutional approach, our results provide an alternative point of view on a recent discussion concerning whether point charges map to the Schwarzschild solution, or the more general two-parameter JNW solution, which includes a dilaton field. We confirm the latter.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Time-reversed Young’s experiment: Deterministic, diffractionless second-order interference effect

The classic Young’s double-slit experiment exhibits first-order interference, producing alternating bright and dark fringes shaped by the diffraction effect of the slits. In contrast, here we demonstrate that its time-reversed configuration generates an ideal, deterministic second-order ‘ghost’ interference pattern free from diffraction and first-order effects with the use of only a position-fixed detector. The pattern’s size is governed by the dimensions of the ‘effectively extended light source’ formed by point emitter(s). Beyond highlighting the nonreciprocal physics between the two configurations, this system unlocks a range of novel phenomena inaccessible in traditional double-slit experiments. These include fully programmable, digitized interference fringe formations and the ability to align the pattern plane with the source plane on the same side of the setup. Remarkably, our proposed experiment achieves these outcomes without relying on nonclassical correlations or quantum entanglement. Here, by restoring time-reversal symmetry and eliminating diffraction, our theoretical analysis reveals that this approach offers exciting potential for advancing optical imaging and sensing technologies with improved transverse and longitudinal phase-shift sensitivity and resolution limit beyond current limitations.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Algebraic Structures in the Coupling of Gravity to Gauge Theories

Highlights: • Geometric foundation for Quantum General Relativity coupled to Quantum Electrodynamics. • Rigorous introduction to Hopf algebraic renormalization. • Appropriate modifications in case renormalization Hopf algebra is ill-defined. • Generalization of Furry’s Theorem to include gravitons and ghosts. • Discussion of obstructions to multiplicative renormalization. This article is an extension of the author’s second master thesis, Prinz (2017). It aims to introduce to the theory of perturbatively quantized General Relativity coupled to Spinor Electrodynamics, provide the results thereof and set the notation to serve as a starting point for further research in this direction. It includes the differential geometric and Hopf algebraic background, as well as the corresponding Lagrange density and some renormalization theory. Then, a particular problem in the renormalization of Quantum General Relativity coupled to Quantum Electrodynamics is addressed and solved by a generalization of Furry’s Theorem. Next, the restricted combinatorial Green’s functions for all two-loop propagators and all one-loop divergent subgraphs thereof are presented. Finally, relations between these one-loop restricted combinatorial Green’s functions necessary for multiplicative renormalization are discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Chasing ghosts: characterization of artifact generation in coded aperture decoding due to experimental implementation

Coded aperture imaging is a form of lensless aperture imaging that projects multiple overlapping images of the source onto the detector, enhancing signal strength, which is advantageous for low-flux sources or high-resolution imaging. This technique requires decoding of the detector signal to reconstruct the original source, which involves convolving the detector data with the aperture pattern. When the signal is from a centered point source, the reconstructed source image is known as the point spread function (PSF). A clean PSF without artifacts is a Dirac delta function [Appl. Opt. 20, 1858 (1981)]. This paper examines the robustness of the decoding process against variations in experimental tolerances by analyzing artifact growth in the reconstructed PSF. We illustrate the effects of incorrect magnification, rotation, and detector size and find that aperture–detector rotational misalignment about the imaging axis is the most sensitive parameter, with significant artifact generation occurring with angular offsets of less than one degree. We discuss compensation methods for imperfect aperture placement, finding that small detector sizes produce uncompensatable artifact generation, and compare theoretical predictions with experimental PSF measurements of a rank , 6.8 mm thick (less than one mean free path) coded aperture with a 3.5 mm cell size, conducted at the MegaJOuLe Neutron Imaging Radiography dense plasma focus [IEEE Trans. Plasma Sci. 49, 3299 (2021)] using a 2.45 MeV neutron source. Based on our findings, we recommend using magnified coded apertures in the under-sampled regime, which allows for the inclusion of fiducial markers to characterize aperture–detector rotational offsets and the addition of mechanical coupling, where possible, to constrain rotational and magnification offsets.

Selwood, M. P. [Lawrence Livermore National Labora↗
DOE OSTI2026

CHESS 2025: Spectrometer orthorectified at-sensor radiance from NEON AOP imaging spectroscopy surveys

This dataset provides Level 1 (L1) orthorectified at-sensor radiance derived from measurements collected by the Imaging Spectrometer-1 (NIS-1) onboard the NEON (National Ecological Observatory Network) Airborne Observation Platform (AOP) for the 2025 Colorado Headwaters Ecological Spectroscopy Study (CHESS). NIS-1 captures light reflected from the Earth’s surface in 426 discrete wavelength bands as raw digital numbers (DNs; Level 0). These data are then calibrated to physical units (uW/cm²·sr·nm) following the processing steps described in the NEON Imaging Spectrometer Level 1B Calibrated Radiance Algorithm Theoretical Basis Document (ATBD; Gallery 2022). The data delivered here are the primary inputs for the surface reflectance product in “Custom surface reflectance, shade masks, and equivalent water thickness maps for the Colorado Headwaters Ecological Spectroscopy Study” (Carroll et al. 2026). For intertemporal comparison, the radiance data here are most directly relatable to the v2 radiance data in “NEON AOP Imaging Spectroscopy Survey of Upper East River Colorado Watersheds: Raw-Space Radiance and Observational Variable Dataset” (Goulden et al. 2018), to which the same processing methodology was applied. Together, the radiance and reflectance data enable users to exploit the unique reflection signatures of different surface objects for land cover classification, foliar trait mapping, plant vigor assessment, water content estimation, trace-element identification, and other scientific applications. The data were acquired over three study domains in the Upper Gunnison river basin: the upper East River watershed (CRBU); Almont Triangle and Taylor Canyon (ALMO); and Upper Taylor River watershed (UPTA) between 2025-06-13 and 2025-07-15. Within each domain, data are delivered by flightline as orthorectified and calibrated hyperspectral rasters in Hierarchical Data Format version 5 (HDF5) format, with radiance values provided in uW/cm²·sr·nm on a fixed, uniform Universal Transverse Mercator (UTM) grid at 1 meter spatial resolution. The radiance rasters include all 426 NIS-1 spectral bands, along with associated quality-assurance (QA) and diagnostic and ancillary layers needed for atmospheric correction workflows. Orthorectified radiance is produced from pushbroom spectrometer observations by applying NEON’s radiometric calibration (including bad pixel masking, dark subtract, dark pedestal shift correction, electronic panel ghost correction, grating ghost correction, deblur correction and flat-fielding) and spectral calibration (using spectral response function band centers and full-width at half-maximum intensity), followed by geolocation and regridding to the fixed grid. CHESS Project Description: The Colorado Headwaters Ecological Spectroscopy Study (CHESS) comprised a multi-week airborne remote sensing and field observation campaign in the Upper Gunnison Basin, Colorado, conducted in June and July of 2025. Airborne remote sensing was conducted by the National Ecological Observatory Network Airborne Observation Platform (NEON AOP), concurrent with a field campaign run by the Rocky Mountain Biological Laboratory (RMBL), the Lawrence Berkeley National Laboratory (LBNL) and SLAC National Accelerator Laboratory Watershed Function Science Focus Area (SFA), and NASA-JPL (Jet Propulsion Laboratory) Earth Surface Mineral Dust Source Investigation (EMIT) program. Between June 10 and July 18, 2025, the NEON AOP flight team collected high-resolution aerial imaging spectroscopy and Light Detection and Ranging (LiDAR) data over three domains: the Upper East River (CRBU), Almont Triangle (ALMO), and the Upper Taylor Basin (UPTA). In coordination with the flights, a field campaign acquired ground-truth observations, including observations of vegetation composition, foliar traits, forest demography, and subsurface properties in 18 core sampling areas within the domains. Additional surface water observations were taken at over 380 point locations. All CHESS campaign datasets can be found within the CHESS ESS-DIVE data portal: https://data.ess-dive.lbl.gov/portals/chess. Funding Acknowledgement: Field and remote-sensing data acquisition was performed under a grant from the National Aeronautics and Space Administration (80NSSC24K1005). This work was also supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

2018 NEON and 2025 CHESS Campaigns↗
DOE OSTI2025

Neutron Identification Capabilities in MicroBooNE Through the Application of Machine Learning with Blips

Neutrinos (ν) are subatomic particles first observed in 1956 by LANL physicists Clyde Cowan and Frederick Reines but first theorized by Wolfgang Pauli in 1930. Neutrinos are the least massive known particle and are classified as leptons with 3 flavors corresponding to their leptonic counterparts (electron, muon, and tau). We know that they are abundant, 65 billion neutrinos travel through your fingertip every second, and elusive, a single neutrino could fly through a lightyear of lead without interacting at all. Though, there is much still unknown and a better characterization of these “ghostly” particles can give us clues as to the matter/anti-matter asymmetry in the early universe and possible glimpses into new physics. To measure a particle that is extremely light and rarely interacting, physicists have developed an extremely sensitive detector known as a Liquid Argon Time Projection Chamber (LArTPC). The fiducial volume (or TPC) is bombarded with neutrinos, some of which interact with argon (Ar) atoms to produce particles that in turn excite and ionize the Ar. The products of these are free electrons which then drift through the TPC’s applied magnetic field towards a multi-plane wire readout system. The electrons’ charge is collected at this anode and the light from the initial interactions is collected by photomultiplier tubes (PMTs). In conjunction, these mechanisms allow LArTPCs to achieve millimeter spatial resolution and sub-MeV energy thresholds. The detector of interest in this study is the MicroBooNE Experiment at Fermilab. MicroBooNE is an above ground LArTPC with dimensions of approximately 10m × 2.5m × 2.3m, about the size of a school bus. Its purpose is to study neutrinos, so to improve rates of measured ν interactions, the detector is squarely in the path of the Booster Neutrino Beam (BNB) at Fermilab. A major challenge in neutrino studies is energy reconstruction, much of the neutrino’s original energy is lost in interactions that the detector is not sensitive to, often due to low-energy products. The initial goal of this analysis was to better identify neutrons, the main source of poor energy reconstruction in neutrino events. Because neutrons are neutral particles, like neutrinos, we can only directly measure the products of their interactions in LArTPCs. Most of these products are low-energy signals and while each one contributes a negligible amount of energy, collectively these signals make up most of the lost energy in each neutrino event. We define these signals as blips; point-like, isolated depositions of charge in the detector. Blips have MeV-scale energies and are the size of a single hit (charge deposition) or a cluster of a few hits on at least two wire planes. Blips are the principal detector features used to study low-energy physics; thus, they are the key to unlocking information not only about neutrons but gamma photons, supernova and solar neutrinos as well as helping us better identify certain particles. Therefore, this analysis strives to use blips for improved neutron identification (ID) and characterization.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗