Photon Transport Research Articles

Luminescent Anthracene-based Elastic Mixed Molecular Crystals for Flexible and FRET-Assisted Optical Waveguides

Abstract Optical waveguides based on elastic molecular crystals are of interest as flexible and compact optical communication materials. Low emission efficiency is often a problem in terms of communication signal strength, and an increase in the loss factor α due to fluorescence reabsorption in the crystal reduces the photon transport efficiency. Here, we report the development and improvement of Förster resonance energy transfer (FRET)-assisted optical waveguides using anthracene-based elastic mixed molecular crystals. 9,10-Dicyanoanthracene was selected as the dopant for solution-grown crystallization of 9,10-dibromoanthracene. 1H NMR measurements of the obtained crystals showed that the acceptor doping to be 1%-5%. Elastic behavior was observed even when doped with a few percent of the acceptor. The quantum efficiency was 0.016, a dramatic improvement over the previous luminescent EMMC (0.004). The α value (92 dB/cm) of this crystal containing 1%-doped 9,10-dicyanoanthracene is much lower than that of the crystal consisting of only 9,10-dibromoanthracene (1258 dB/cm) due to the reduced reabsorption in the FRET system. We have demonstrated a practical approach towards developing improved fluorescent, highly efficient, and flexible optical waveguides by constructing the mixed crystal structure by selecting appropriate acceptor molecules and their low doping ratios.

A framework to model charge sharing and pulse pileup for virtual imaging trials of photon-counting CT

Objective.This study describes the development, validation, and integration of a detector response model that accounts for the combined effects of x-ray crosstalk, charge sharing, and pulse pileup in photon-counting detectors.Approach.The x-ray photon transport was simulated using Geant4, followed by analytical charge sharing simulation in MATLAB. The analytical simulation models charge clouds with Gaussian-distributed charge densities, which are projected on a 3×3 pixel neighborhood of interaction location to compute detected counts. For pulse pileup, a prior analytical method for redistribution of energy-binned counts was implemented for delta pulses. The x-ray photon transport and charge sharing components were validated using experimental data acquired on the CdTe-based Pixirad-1/Pixie-III detector using monoenergetic beams at 26, 33, 37, and 50 keV. The pulse pileup implementation was verified with a comparable Monte Carlo simulation. The model output without pulse pileup was used to generate spatio-energetic response matrices for efficient simulation of scanner-specific photon-counting CT (PCCT) images with DukeSim, with pulse pileup modeled as a post-processing step on simulated projections. For analysis, images for the Gammex multi-energy phantom and the XCAT chest phantom were simulated at 120 kV, both with and without pulse pileup for a range of doses (27-1344 mAs). The XCAT images were evaluated qualitatively at 120 mAs, while images for the Gammex phantom were evaluated quantitatively for all doses using measurements of attenuation coefficients and Calcium concentrations.Main results.Reasonable agreement was observed between simulated and experimental spectra with Mean Absolute Percentage Error Values (MAPE) between 10%and 31%across all incident energies and detector modes. The increased pulse pileup from increased dose affected attenuation coefficients and calcium concentrations, with an effect on calcium quantification as high as MAPE of 28%.Significance.The presented approach demonstrates the viability of the model for enabling VITs to assess and optimize the clinical performance of PCCT.

An implicit unified gas-kinetic particle method with large time steps for gray radiation transport

For a long time, efficient algorithms for high-dimensional equations, represented by photon radiation transport, have been one important topic in the development of computational methods for particle transport processes. In this paper, we present an implicit unified gas-kinetic particle (IUGKP) method for multiscale gray radiative transfer. Based on the integral solution of the radiative transfer equation, the photon transport processes are categorized into non-equilibrium transport processes with a large photon free path and equilibrium transport processes with a small photon free path. The long-path processes are solved by an implicit Monte Carlo (IMC) method, and the short-path processes are solved by an implicit diffusion system. The closure formulation of photon distribution is derived from the local integral solution of the radiative transfer equation to couple the IMC and diffusion system. The improvement of the proposed IUGKP method over UGKP method is that particles can be tracked continuously instead of just until the first collision, making simulation with large time steps possible. The IUGKP method has the properties of asymptotic-preserving (AP) and regime-adaptive (RA). The AP property states that the IUGKP method converges to the consistent numerical methods for the asymptotic limiting equations of RTE in the limiting regimes. The RA property states that the computational accuracy of the IUGKP method adapts to the regimes. In this paper, the mathematical proof of the AP and RA properties is presented, and the multiscale numerical tests are performed to demonstrate the accuracy and efficiency of the IUGKP method.

Toughening effect and mechanisms of AlN whiskers on a low‐temperature sintered AlNw/AlN ceramics

AbstractConventional aluminum nitride (AlN) ceramics exhibited insufficient mechanical properties expected from their potential applications in the presence of dynamic loads and intensive stresses caused by thermal shock. In this study, the mechanical properties of AlN ceramics were enhanced by toughening them using AlN whiskers (AlNw) via gel casting followed by low‐temperature sintering at 1650°C. The addition of AlNw simultaneously increased the bending strength, toughness and thermal conductivity of the AlN ceramics. The maximum values of the bending strength, toughness, and thermal conductivity of 9.0 wt% AlNw/AlN were 303.92 MPa, 3.92 MPa·m1/2 and 186.53 W·m−1·K−1, respectively, which were higher than those of the AlN ceramics by 39.93%, 61.31% and, 15.61%, respectively. The AlNw in AlNw/AlN ceramics tightly bonded with the ceramic matrix, leading to two characteristic toughening mechanisms in the ceramics: crack pinning and deflection at irregular grain boundaries caused by doped whiskers, as well as bridging and stress relaxation caused by whiskers incorporated with the grains. Moreover, the one‐dimensional morphology of AlNw can provide a channel for quick photon transport, thereby enhancing the thermal conductivity of AlNw/AlN.

Theoretical framework for calibrating the depth-dependent optical scattering in layered human skin using spatially offset measurements.

Spatially offset spectroscopy offers an alternative non-invasive method for enabling deep probing of structures and chemical molecules, which is clinically significant for the characterization of chemical and physical alterations in human skin. However, a more precise depth-resolved quantification using the spatially offset measurements still remains a challenge due to the mixed inhomogeneous scattering. Herein, we report a Monte-Carlo-based quantification modeling platform combined with a novel, to the best of our knowledge, scattering spectrum decomposition method to explore the depth-dependent optical scattering contributions in human skin. In the simplified modeling, human skin was empirically set to be composed of multiple layers, and each layer possessed different photon weights for the spatially offset scattering intensity measurements. The modeling results of photon transportation in-and-out of the layered skin substantially discovered that the layer-dependent scattering contribution was compositely encoded into the spatially offset measurements and varied with the illumination incidence angle. For calibrating the layer-dependent scattering contribution, a modified nonlinear independent component processing algorithm was applied to the spatially offset measurements by decomposing the photon weights of each layer. The calibration results figured out the major scattering contribution of each layer along the offset axis under different incidence angles, which were consistent with previous experimental observations. The proposed theoretical framework establishes a feasible approach for spatially offset optical spectroscopies enabling non-invasive quantitative A-line characterization of the concentrations of skin components.

Nonequilibrium energy transport and nonclassical correlations of photons in a two-mode circuit quantum electrodynamic system

Nonequilibrium energy transport and nonclassical correlations of photons in a two-mode circuit quantum electrodynamic system

Nonreciprocal Photon Transport in a Chiral Optomechanical System

AbstractChiral interaction between light and quantum emitters leads to emergence development of chiral quantum optics and stimulates a wide range of practical applications in quantum regime, such as single‐photon isolation and photon unidirectional emission. Cavity optomechanics studying the interaction between optical and mechanical resonators plays an important role in the field of quantum optics. However, how to achieve the chiral interaction between light and mechanical oscillators and explore the applications of the chiral optomechanical systems are still difficult encountered in cavity optomechanics. Here, a method is proposed to achieve chiral optomechanical interaction by exploiting directional squeezed light in a multimode optomechanical system. Based on the chiral interaction between photon and phonon, the nonreciprocal photon transport at a single‐photon level can be realized. An isolation ratio of and a negligible insertion loss for the photonic isolator are obtained. This method paves the way to realize chiral optomechanical interaction for conducting chiral optomechanics and opens up the prospect of exploring and utilizing chiral photon–phonon manipulation in the quantum regime.

Dynamic tunable nonreciprocal single-photon scattering mediated by a giant atom assisted with a time-modulated cavity

Nonreciprocal single-photon scattering in a one-dimensional waveguide coupled to a giant two-level atom assisted with a time-modulated single-mode cavity is investigated. The analytic expressions of the single-photon scattering amplitudes are derived by using an effective Floquet Hamiltonian in real space. The scattering characteristics are discussed detail in both the Markovian and the non-Markovian regimes, and the corresponding conditions for achieving perfect nonreciprocal single-photon transmission are obtained. In the Markovian regime, a frequency-tunable single-photon diode with an ideal transmission contrast ratio can be realized by adjusting the frequency of the cavity mode, the local coupling phase difference, and the accumulated phase between the two coupling points. Furthermore, the influence of the intrinsic energy dissipations on the photon transport is discussed in detail. It is found that the dissipations of the cavity and the giant atom affect discriminatively the nonreciprocal single-photon scattering process. In the non-Markovian regime, the influence of the non-Markovian retarded effect induced by the time delay on the nonreciprocal single-photon scattering is discussed in detail. The results reveal that, although the retarded effect leads to a complex nonreciprocal scattering spectrum, dynamic tunable perfect nonreciprocal transmission with more abundant physical phenomena suitable for photons with different frequencies within a larger range can also be achieved. Such a nonreciprocal single-photon device can be used as an elementary unit for various quantum information processing and may have potential applications in quantum network engineering.

Atmospheric Response for MeV γ Rays Observed with Balloon-borne Detectors

The atmospheric response for MeV γ rays (∼0.1–10 MeV) can be characterized in terms of two observed components. The first component is due to photons that reach the detector without scattering. The second component is due to photons that reach the detector after scattering one or more times. While the former can be determined in a straightforward manner, the latter is much more complex to quantify, as it requires tracking the transport of all source photons that are incident on Earth’s atmosphere. The scattered component can cause a significant energy-dependent distortion in the measured spectrum, which is important to account for when making balloon-borne observations. In this work, we simulate the full response for γ-ray transport in the atmosphere. We find that the scattered component becomes increasingly more significant toward lower energies, and at 0.1 MeV, it may increase the measured flux by as much as a factor of ∼2–4, depending on the photon index and off-axis angle of the source. This is particularly important for diffuse sources, whereas the effect from scattering can be significantly reduced for point sources observed with an imaging telescope.

Chiral Bulk Solitons in Photonic Graphene with Decorated Boundaries

AbstractA chiral bulk soliton in a nonlinear photonic lattice with decorated boundaries, presenting a novel approach to manipulate photonic transport without extensive bulk modifications is proposed. Unlike traditional methods that rely on topological edge and corner modes, this strategy leverages the robust chiral propagation of bulk modes. By introducing nonlinearity into the system, a stable bulk soliton, akin to the topological valley Hall effects is found. The chiral bulk soliton exhibits remarkable stability; the energy does not decay even after a long‐distance propagation; and the corresponding Fourier spectrum confirms the absence of inter‐valley scattering indicating a valley‐locking property. The findings not only contribute to the fundamental understanding of nonlinear photonic systems but also hold significant practical implications for the design and optimization of photonic devices.

Epitaxial Growth of Two-Dimensional Organic Crystals with In-Plane Heterostructured Domain Regulation.

Complex organic lateral heterostructures (OLHs) with spatial distribution of two or more chemical components are crucial for designing and realizing unique structure-dependent optoelectronic applications. However, the precise design of well-defined OLHs with flexible domain regulation remains a considerable challenge. Herein, we present a stepwise solution self-assembly method to synthesize two-dimensional (2D) OLHs with a central rhombus domain and a lateral region featuring tunable blue and green emission based on the sequential nucleation and growth of 2D crystals. By controlling the initial crystallization time of 2,6-diphenylanthracene, the rhombic length ratio (α) of the multicolor-emissive part of the 2D OLHs is precisely modified. Furthermore, a third lateral layer is constructed on the resulting OLHs, demonstrating scalable lateral regulation. Significantly, these prepared 2D OLHs exhibit great excitation position-dependent waveguide characteristics and enable a 0.06 dB/μm low-loss waveguiding, which are conducive to photon transport and conversion for photonic integrated circuits. This work provides a stepwise strategy for the accurate fabrication of 2D OLHs, fabricating the developments of next-generation optoelectronics devices.

Axion-photon mixing in 3D: classical equations and geometric optics

Light particle-photon mixing in magnetised plasmas plays a vital role in constraining the existence of new physics, especially axions, dark photons, and ultra-high-frequency gravitational waves. Recently, we derived an expression for the resonant conversion of axions to photons in inhomogeneous media using kinetic theory to derive photon transport equations. In this work, we show how the same expression for the conversion probability can be obtained from the classical wave equations of axion-electrodynamics by deriving an equivalent transport equation along the photon worldline. This result provides further corroboration of this expression for the resonant production of photons from light particles, which has also recently been supported by independent numerical simulations of full axion-electrodynamics. In addition, this new approach provides a more general expression that accounts for mixing away from resonance, which is integrated along the whole worldline of the photon in a way that naturally incorporates a curved photon trajectory relevant to refractive media where the photon and light-particle worldlines differ.

Evaluation of TOPAS MC tool performance in optical photon transport and radioluminescence-based dosimetry

Radiation therapy plays a pivotal role in modern cancer treatment, demanding precise and accurate dose delivery to tumor sites while minimizing harm to surrounding healthy tissues. Monte Carlo simulations have emerged as indispensable tools for achieving this precision, offering detailed insights into radiation transport and interaction at the subatomic level. As the use of scintillation and luminescence dosimetry becomes increasingly prevalent in radiation therapy, there arises a need for validated Monte Carlo tools tailored to optical photon transport applications. In this paper, an evaluation process of the TOPAS (TOol for PArticle Simulation) Monte Carlo tool for Cerenkov light generation, optical photon transport and radioluminescence based dosimetry is presented. Three distinct sources of validation data are utilized: one from a published set of experimental results and two others from simulations performed with the Geant4 code. The methodology employed for evaluation includes the selection of benchmark experiments, making use of opt3 and opt4 Geant4 physics models and simulation setup, with observed slight discrepancies within the calculation uncertainties. Additionally, the complexities and challenges associated with modeling optical photons generation through luminescence or Cerenkov radiation and their transport are discussed. The results of our evaluation suggests that TOPAS can be used to reliably predict Cerenkov generation, luminescence phenomenon and the behavior of optical photons in common dosimetry scenarios.

Impacts of vertical variations in soil properties on H*(10) simulations for 137Cs deposition

Photon transport simulations based on the Monte Carlo method have played a crucial role in assessing and estimating the ambient dose equivalent rates H*(10), resulting from the deposition of 137Cs in soil following the nuclear power plant accident in Fukushima. However, a comprehensive examination of the effect of vertical variations in soil properties on the simulation outcomes has not yet been performed. Disregarding the vertical distribution of soil properties not only leads to potential inaccuracies in the shielding responses of soil layers but also in the determination of the radioactive source inventory, particularly when using the concentration data in Bq/kg. These oversights diminish the reliability of the simulation results. This study addresses several soil property factors that could potentially influence the simulation results, including variations in chemical composition induced by water content, bulk density profile, and estimated inventory profile, all evaluated through an examined simulation model. The results show that inappropriate assignment of the soil density profile can cause considerable errors in the H*(10) simulation outcomes. Furthermore, the sensitivity of H*(10) to variations in soil vertical density is analyzed, with the results indicating that H*(10) can be highly sensitive to changes in the bulk density of the top 0–5 cm soil layers. These results should facilitate the establishment of appropriate simulation strategies and support the reassessment of past simulation results.

A GPU-accelerated Monte Carlo code, RT2 for coupled transport of photon, electron/positron, and neutron

Objective. This work aims to develop a graphics processing unit (GPU)-accelerated Monte Carlo code for the coupled transport of photon, electron/positron and neutron over a broad range of energies for medical applications. Approach. By separating the MC evolution of radiation into source, transport, and interaction kernels, the branch divergence was alleviated. The memory coalescence was achieved by vectorizing the access pattern in which the secondary particles were archived. To accelerate further particle tracking, ray-tracing hardware acceleration in the Nvidia OptiXTM framework was applied. For photon and electron/positron, the EGSnrc interaction modules were ported as a GPU-optimized configuration. For neutron, a group-wised transport based on NJOY21 preprocessed data was implemented. The developed code was validated against CPU-based FLUKA. Neutron, x-ray and electron beams incident on water and ICRP phantoms were simulated. The neutron energy group and the transport parameters of photon and electron were set to be the same in both codes. A single Nvidia RTX 4090 card was used in this code while all 20 threads of a single Intel Core i9-10900K node were used in FLUKA. Main results. The number of histories was set to ensure that statistical uncertainties lower than 2% for all voxels whose doses were larger than 20% of the maximum. In all cases, the dose differences in the voxels between the codes were within 2.5%. For photons and electrons, the developed code was 150–300 times faster than FLUKA in both geometries. For neutrons, the code was respectively 80 and 135 times faster in the water and ICRP phantoms than FLUKA. Significance. This study offers an appropriate solution for uncoalesced memory access and branch divergence commonly encountered in coupled MC transport on the GPU architecture. The formidable acceleration in computing times and accuracy shown in this study can promise a routine clinical use of MC simulations.

Guiding charged particles in vacuum via Lagrange points

We propose a method for guiding charged particles such as electrons and protons, in vacuum, by employing the exotic properties of Lagrange points. This leap is made possible by the dynamics unfolding around these equilibrium points, which stably capture such particles, akin to the way Trojan asteroids are held in Jupiter’s orbit. Unlike traditional methodologies that allow for either focusing or three-dimensional storage of charged particles, the proposed scheme can guide both non-relativistic and relativistic electrons and protons in small cross-sectional areas in an invariant fashion over long distances without any appreciable loss in energy – in a manner analogous to photon transport in optical fibers. Here, particle guiding is achieved by employing twisted electrostatic potentials that in turn induce stable Lagrange points in vacuum. In principle, guidance can be realized within the fundamental mode of the resulting waveguide, thereby presenting a prospect for manipulating these particles in the quantum domain. Our findings may be useful in a wide range of applications in both scientific and technological pursuits. These applications could encompass electron microscopies and lithographies, particle accelerators, quantum and classical communication/sensing systems, as well as methods for shuttling entangled qubits between nodes within a quantum network.

Unconventional Thermophotonic Charge Density Wave.

Charge-order states of broken symmetry, such as charge density wave (CDW), are able to induce exceptional physical properties, however, the precise understanding of the underlying physics is still elusive. Here, we combine fluctuational electrodynamics and density functional theory to reveal an unconventional thermophotonic effect in CDW-bearing TiSe_{2}, referred to as thermophotonic-CDW (tp-CDW). The interplay of plasmon polariton and CDW electron excitations give rise to an anomalous negative temperature dependency in thermal photons transport, offering an intuitive fingerprint for a transformation of the electron order. Additionally, the demonstrated nontrivial features of tp-CDW transition hold promise for a controllable manipulation of heat flow, which could be extensively utilized in various fields such as thermal science and electron dynamics, as well as in next-generation energy devices.

Reduced-order model to approximate response matrices for filter stack spectrometers.

We present a reduced-order model to calculate response matrices rapidly for filter stack spectrometers (FSSs). The reduced-order model allows response matrices to be built modularly from a set of pre-computed photon and electron transport and scattering calculations through various filter and detector materials. While these modular response matrices are not appropriate for high-fidelity analysis of experimental data, they encode sufficient physics to be used as a forward model in design optimization studies of FSSs, particularly for machine learning approaches that require sampling and testing a large number of FSS designs.

Technical note: A GPU-based shared Monte Carlo method for fast photon transport in multi-energy x-ray exposures.

The Monte Carlo (MC) method is an accurate technique for particle transport calculation due to the precise modeling of physical interactions. Nevertheless, the MC method still suffers from the problem of expensive computational cost, even with graphics processing unit (GPU) acceleration. Our previous works have investigated the acceleration strategies of photon transport simulation for single-energy CT. But for multi-energy CT, conventional individual simulation leads to unnecessary redundant calculation, consuming more time. This work proposes a novel GPU-based shared MC scheme (gSMC) to reduce unnecessary repeated simulations of similar photons between different spectra, thereby enhancing the efficiency of scatter estimation in multi-energy x-ray exposures. The shared MC method selects shared photons between different spectra using two strategies. Specifically, we introduce spectral region classification strategy to select photons with the same initial energy from different spectra, thus generating energy-shared photon groups. Subsequently, the multi-directional sampling strategy is utilized to select energy-and-direction-shared photons, which have the same initial direction, from energy-shared photon groups. Energy-and-direction-shared photons perform shared simulations, while others are simulated individually. Finally, all results are integrated to obtain scatter distribution estimations for different spectral cases. The efficiency and accuracy of the proposed gSMC are evaluated on the digital phantom and clinical case. The experimental results demonstrate that gSMC can speed up the simulation in the digital case by ∼37.8% and the one in the clinical case by ∼20.6%, while keeping the differences in total scatter results within 0.09%, compared to the conventional MC package, which performs an individual simulation. The proposed GPU-based shared MC simulation method can achieve fast photon transport calculation for multi-energy x-ray exposures.

A Generic Implementation of the D1S Methodology for Dose Rate Assessment by Monte Carlo Codes in Fusion Applications

Dose rate assessment is of utmost importance in fusion reactors because of the maintenance operations that these facilities will require. The standard way to perform this assessment in ITER is use of the Direct 1-Step (D1S) methodology, which consists of performing neutron transport simulation, material activation, and decay photon transport simulation in one step and thus in one simulation only. Usually, implementation of the D1S methodology requires changing the source files in a reference Monte Carlo code. The purpose of the present work is to develop an easy-to-implement method for Monte Carlo codes to help calculate the dose rate in fusion reactors. To do so, the proposal is to act on nuclear data only and not on source files of the simulation codes. This is done by replacing prompt photon production in evaluation files by suitable decay photon production, taking into account radioactive decays of radionuclides and irradiation history. This study was to be applied first on the TRIPOLI-4® Monte Carlo code for the sake of simplicity. TRIPOLI-4 is the reference code for particle transport at CEA. The verification and validation process relies first on a comparison to the reference Rigorous 2-Step (R2S) methodology and then on an experiment, the so-called Frascati Neutron Generator (FNG) dose experiment. The nuclear database to be changed is of the ENDF-6 format, a recurrent format in neutronics studies. The analysis studied both neutron and photon responses to check if the simulation was performing normally in a physical way and to compare the results with references provided by simulations based on the R2S methodology or by the FNG dose experiment. The simulations have proven to be in good agreement with the experimental results.