Pixel-by-pixel correction of beam hardening artifacts by bowtie filter in fan-beam CT

Objective. Beam hardening (BH) artifacts in computed tomography (CT) images originate from the polychromatic nature of x-ray photons. In a CT system with a bowtie filter, residual BH artifacts remain when polynomial fits are used. These artifacts lead to worse visuals, reduced contrast, and inaccurate CT numbers. This work proposes a pixel-by-pixel correction (PPC) method to reduce the residual BH artifacts caused by a bowtie filter.Approach. The energy spectrum for each pixel at the detector after the photons pass through the bowtie filter was calculated. Then, the spectrum was filtered through a series of water slabs with different thicknesses. The polychromatic projection corresponding to the thickness of the water slab for each detector pixel could be obtained. Next, we carried out a water slab experiment with a mono energyE= 69 keV to get the monochromatic projection. The polychromatic and monochromatic projections were then fitted with a 2nd-order polynomial. The proposed method was evaluated on digital phantoms in a virtual CT system and phantoms in a real CT machine.Main results. In the case of a virtual CT system, the standard deviation of the line profile was reduced by 23.8%, 37.3%, and 14.3%, respectively, in the water phantom with different shapes. The difference of the linear attenuation coefficients (LAC) in the central and peripheral areas of an image was reduced from 0.010 to 0.003cm-1and 0.007cm-1to 0 in the biological tissue phantom and human phantom, respectively. The method was also validated using CT projection data obtained from Activion16 (Canon Medical Systems, Japan). The difference in the LAC in the central and peripheral areas can be reduced by a factor of two.Significance. The proposed PPC method can successfully remove the cupping artifacts in both virtual and authentic CT images. The scanned object's shapes and materials do not affect the technique.

Measurements of All-Particle Energy Spectrum and Mean Logarithmic Mass of Cosmic Rays from 0.3 to 30 PeV with LHAASO-KM2A

We present the measurements of all-particle energy spectrum and mean logarithmic mass of cosmic rays in the energy range of 0.3-30 PeV using data collected from LHAASO-KM2A between September 2021 and December 2022, which is based on a nearly composition-independent energy reconstruction method, achieving unprecedented accuracy. Our analysis reveals the position of the knee at 3.67±0.05±0.15  PeV. Below the knee, the spectral index is found to be -2.7413±0.0004±0.0050, while above the knee, it is -3.128±0.005±0.027, with the sharpness of the transition measured with a statistical error of 2%. The mean logarithmic mass of cosmic rays is almost heavier than helium in the whole measured energy range. It decreases from 1.7 at 0.3 PeV to 1.3 at 3 PeV, representing a 24% decline following a power law with an index of -0.1200±0.0003±0.0341. This is equivalent to an increase in abundance of light components. Above the knee, the mean logarithmic mass exhibits a power law trend towards heavier components, which is reversal to the behavior observed in the all-particle energy spectrum. Additionally, the knee position and the change in power-law index are approximately the same. These findings suggest that the knee observed in the all-particle spectrum corresponds to the knee of the light component, rather than the medium-heavy components.

Mechanistic modelling of relative biological effectiveness of carbon ion beams and comparison with experiments

Objective.Relative biological effectiveness (RBE) plays a vital role in carbon ion radiotherapy, which is a promising treatment method for reducing toxic effects on normal tissues and improving treatment efficacy. It is important to have an effective and precise way of obtaining RBE values to support clinical decisions. A method of calculating RBE from a mechanistic perspective is reported.Approach.Ratio of dose to obtain the same number of double strand breaks (DSBs) between different radiation types was used to evaluate RBE. Package gMicroMC was used to simulate DSB yields. The DSB inductions were then analyzed to calculate RBE. The RBE values were compared with experimental results.Main results.Furusawa's experiment yielded RBE values of 1.27, 2.22, 3.00 and 3.37 for carbon ion beam with dose-averaged LET of 30.3 keVμm-1, 54.5 keVμm-1, 88 keVμm-1and 137 keVμm-1, respectively. RBE values computed from gMicroMC simulations were 1.75, 2.22, 2.87 and 2.97. When it came to a more sophisticated carbon ion beam with 6 cm spread-out Bragg peak, RBE values were 1.61, 1.63, 2.19 and 2.36 for proximal, middle, distal and distal end part, respectively. Values simulated by gMicroMC were 1.50, 1.87, 2.19 and 2.34. The simulated results were in reasonable agreement with the experimental data.Significance.As a mechanistic way for the evaluation of RBE for carbon ion radiotherapy by combining the macroscopic simulation of energy spectrum and microscopic simulation of DNA damages, this work provides a promising tool for RBE calculation supporting clinical applications such as treatment planning.

Measurement of Ultra-High-Energy Diffuse Gamma-Ray Emission of the Galactic Plane from 10 TeV to 1 PeV with LHAASO-KM2A

The diffuse Galactic γ-ray emission, mainly produced via interactions between cosmic rays and the interstellar medium and/or radiation field, is a very important probe of the distribution, propagation, and interaction of cosmic rays in the Milky Way. In this Letter, we report the measurements of diffuse γ rays from the Galactic plane between 10 TeV and 1 PeV energies, with the square kilometer array of the Large High Altitude Air Shower Observatory (LHAASO). Diffuse emissions from the inner (15°10  TeV). The energy spectrum in the inner Galaxy regions can be described by a power-law function with an index of -2.99±0.04, which is different from the curved spectrum as expected from hadronic interactions between locally measured cosmic rays and the line-of-sight integrated gas content. Furthermore, the measured flux is higher by a factor of ∼3 than the prediction. A similar spectrum with an index of -2.99±0.07 is found in the outer Galaxy region, and the absolute flux for 10≲E≲60  TeV is again higher than the prediction for hadronic cosmic ray interactions. The latitude distributions of the diffuse emission are consistent with the gas distribution, while the longitude distributions show clear deviation from the gas distribution. The LHAASO measurements imply that either additional emission sources exist or cosmic ray intensities have spatial variations.

A tera-electron volt afterglow from a narrow jet in an extremely bright gamma-ray burst

Some gamma-ray bursts (GRBs) have a tera-electron volt (TeV) afterglow, but the early onset of this has not been observed. We report observations with the Large High Altitude Air Shower Observatory of the bright GRB 221009A, which serendipitously occurred within the instrument field of view. More than 64,000 photons >0.2 TeV were detected within the first 3000 seconds. The TeV flux began several minutes after the GRB trigger, then rose to a peak about 10 seconds later. This was followed by a decay phase, which became more rapid ~650 seconds after the peak. We interpret the emission using a model of a relativistic jet with half-opening angle ~0.8°. This is consistent with the core of a structured jet and could explain the high isotropic energy of this GRB.

Cross-Channel Constraints on Resonant Antikaon-Nucleon Scattering

Chiral perturbation theory and its unitarized versions have played an important role in our understanding of the low-energy strong interaction. Yet, so far, such studies typically deal exclusively with perturbative or nonperturbative channels. In this Letter, we report on the first global study of meson-baryon scattering up to one-loop order. It is shown that covariant baryon chiral perturbation theory, including its unitarization for the negative strangeness sector, can describe meson-baryon scattering data remarkably well. This provides a highly nontrivial check on the validity of this important low-energy effective field theory of QCD. We show that the K[over ¯]N related quantities can be better described in comparison with those of lower-order studies, and with reduced uncertainties due to the stringent constraints from the πN and KN phase shifts. In particular, we find that the two-pole structure of Λ(1405) persists up to one-loop order reinforcing the existence of two-pole structures in dynamically generated states.

New physics searches at kaon and hyperon factories

Rare meson decays are among the most sensitive probes of both heavy and light new physics. Among them, new physics searches using kaons benefit from their small total decay widths and the availability of very large datasets. On the other hand, useful complementary information is provided by hyperon decay measurements. We summarize the relevant phenomenological models and the status of the searches in a comprehensive list of kaon and hyperon decay channels. We identify new search strategies for under-explored signatures, and demonstrate that the improved sensitivities from current and next-generation experiments could lead to a qualitative leap in the exploration of light dark sectors.

Constraints on Heavy Decaying Dark Matter from 570 Days of LHAASO Observations

The kilometer square array (KM2A) of the large high altitude air shower observatory (LHAASO) aims at surveying the northern γ-ray sky at energies above 10 TeV with unprecedented sensitivity. γ-ray observations have long been one of the most powerful tools for dark matter searches, as, e.g., high-energy γ rays could be produced by the decays of heavy dark matter particles. In this Letter, we present the first dark matter analysis with LHAASO-KM2A, using the first 340 days of data from 1/2-KM2A and 230 days of data from 3/4-KM2A. Several regions of interest are used to search for a signal and account for the residual cosmic-ray background after γ/hadron separation. We find no excess of dark matter signals, and thus place some of the strongest γ-ray constraints on the lifetime of heavy dark matter particles with mass between 10^{5} and 10^{9}  GeV. Our results with LHAASO are robust, and have important implications for dark matter interpretations of the diffuse astrophysical high-energy neutrino emission.

Virtual computed-tomography system for deep-learning-based material decomposition

Objective. Material decomposition (MD) evaluates the elemental composition of human tissues and organs via computed tomography (CT) and is indispensable in correlating anatomical images with functional ones. A major issue in MD is inaccurate elemental information about the real human body. To overcome this problem, we developed a virtual CT system model, by which various reconstructed images can be generated based on ICRP110 human phantoms with information about six major elements (H, C, N, O, P, and Ca). Approach. We generated CT datasets labelled with accurate elemental information using the proposed generative CT model and trained a deep learning (DL)-based model to estimate the material distribution with the ICRP110 based human phantom as well as the digital Shepp–Logan phantom. The accuracy in quad-, dual-, and single-energy CT cases was investigated. The influence of beam-hardening artefacts, noise, and spectrum variations were analysed with testing datasets including elemental density and anatomical shape variations. Main results. The results indicated that this DL approach can realise precise MD, even with single-energy CT images. Moreover, noise, beam-hardening artefacts, and spectrum variations were shown to have minimal impact on the MD. Significance. Present results suggest that the difficulty to prepare a large CT database can be solved by introducing the virtual CT system and the proposed technique can be applied to clinical radiodiagnosis and radiotherapy.

Accurate Relativistic Chiral Nucleon-Nucleon Interaction up to Next-to-Next-to-Leading Order

We construct a relativistic chiral nucleon-nucleon interaction up to the next-to-next-to-leading order in covariant baryon chiral perturbation theory. We show that a good description of the np phase shifts up to T_{lab}=200  MeV and even higher can be achieved with a χ[over ˜]^{2}/d.o.f. less than 1. Both the next-to-leading-order results and the next-to-next-to-leading-order results describe the phase shifts equally well up to T_{lab}=200  MeV, but for higher energies, the latter behaves better, showing satisfactory convergence. The relativistic chiral potential provides the most essential inputs for relativistic ab initio studies of nuclear structure and reactions, which has been in need for almost two decades.

Exploring Lorentz Invariance Violation from Ultrahigh-Energy γ Rays Observed by LHAASO

Recently, the LHAASO Collaboration published the detection of 12 ultrahigh-energy γ-ray sources above 100 TeV, with the highest energy photon reaching 1.4 PeV. The first detection of PeV γ rays from astrophysical sources may provide a very sensitive probe of the effect of the Lorentz invariance violation (LIV), which results in decay of high-energy γ rays in the superluminal scenario and hence a sharp cutoff of the energy spectrum. Two highest energy sources are studied in this work. No signature of the existence of the LIV is found in their energy spectra, and the lower limits on the LIV energy scale are derived. Our results show that the first-order LIV energy scale should be higher than about 10^{5} times the Planck scale M_{Pl} and that the second-order LIV scale is >10^{-3}M_{Pl}. Both limits improve by at least one order of magnitude the previous results.

Physical density estimations of single- and dual-energy CT using material-based forward projection algorithm: a simulation study

OBJECTIVES

This study aims to evaluate the accuracy of physical density prediction in single-energy CT (SECT) and dual-energy CT (DECT) by adapting a fully simulation-based method using a material-based forward projection algorithm (MBFPA).

METHODS

We used biological tissues referenced in ICRU Report 44 and tissue substitutes to prepare three different types of phantoms for calibrating the Hounsfield unit (HU)-to-density curves. Sinograms were first virtually generated by the MBFPA with four representative energy spectra (i.e. 80 kVp, 100 kVp, 120 kVp, and 6 MVp) and then reconstructed to form realistic CT images by adding statistical noise. The HU-to-density curves in each spectrum and their pairwise combinations were derived from the CT images. The accuracy of these curves was validated using the ICRP110 human phantoms.

RESULTS

The relative mean square errors (RMSEs) of the physical density by the HU-to-density curves calibrated with kV SECT nearly presented no phantom size dependence. The kV-kV DECT calibrated curves were also comparable with those from the kV SECT. The phantom size effect became notable when the MV X-ray beams were employed for both SECT and DECT due to beam-hardening effects. The RMSEs were decreased using the biological tissue phantom.

CONCLUSION

Simulation-based density prediction can be useful in the theoretical analysis of SECT and DECT calibrations. The results of this study indicated that the accuracy of SECT calibration is comparable with that of DECT using biological tissues. The size and shape of the calibration phantom could affect the accuracy, especially for MV CT calibrations.

ADVANCES IN KNOWLEDGE

The present study is based on a full simulation environment, which accommodates various situations such as SECT, kV-kV DECT, and even kV-MV DECT. In this paper, we presented the advances pertaining to the accuracy of the physical density prediction when applied to SECT and DECT in the MV X-ray energy range. To the best of our knowledge, this study is the first to validate the physical density estimation both in SECT and DECT using human-type phantoms.

Peta-electron volt gamma-ray emission from the Crab Nebula

The Crab Nebula is a bright source of gamma rays powered by the Crab Pulsar's rotational energy through the formation and termination of a relativistic electron-positron wind. We report the detection of gamma rays from this source with energies from 5 × 10^-4 to 1.1 peta-electron volts with a spectrum showing gradual steepening over three energy decades. The ultrahigh-energy photons imply the presence of a peta-electron volt electron accelerator (a pevatron) in the nebula, with an acceleration rate exceeding 15% of the theoretical limit. We constrain the pevatron's size between 0.025 and 0.1 parsecs and the magnetic field to ≈110 microgauss. The production rate of peta-electron volt electrons, 2.5 × 10³⁶ ergs per second, constitutes 0.5% of the pulsar spin-down luminosity, although we cannot exclude a contribution of peta-electron volt protons to the production of the highest-energy gamma rays.

Review of Deep Learning Based Automatic Segmentation for Lung Cancer Radiotherapy

Lung cancer is the leading cause of cancer-related mortality for males and females. Radiation therapy (RT) is one of the primary treatment modalities for lung cancer. While delivering the prescribed dose to tumor targets, it is essential to spare the tissues near the targets-the so-called organs-at-risk (OARs). An optimal RT planning benefits from the accurate segmentation of the gross tumor volume and surrounding OARs. Manual segmentation is a time-consuming and tedious task for radiation oncologists. Therefore, it is crucial to develop automatic image segmentation to relieve radiation oncologists of the tedious contouring work. Currently, the atlas-based automatic segmentation technique is commonly used in clinical routines. However, this technique depends heavily on the similarity between the atlas and the image segmented. With significant advances made in computer vision, deep learning as a part of artificial intelligence attracts increasing attention in medical image automatic segmentation. In this article, we reviewed deep learning based automatic segmentation techniques related to lung cancer and compared them with the atlas-based automatic segmentation technique. At present, the auto-segmentation of OARs with relatively large volume such as lung and heart etc. outperforms the organs with small volume such as esophagus. The average Dice similarity coefficient (DSC) of lung, heart and liver are over 0.9, and the best DSC of spinal cord reaches 0.9. However, the DSC of esophagus ranges between 0.71 and 0.87 with a ragged performance. In terms of the gross tumor volume, the average DSC is below 0.8. Although deep learning based automatic segmentation techniques indicate significant superiority in many aspects compared to manual segmentation, various issues still need to be solved. We discussed the potential issues in deep learning based automatic segmentation including low contrast, dataset size, consensus guidelines, and network design. Clinical limitations and future research directions of deep learning based automatic segmentation were discussed as well.

Extended Very-High-Energy Gamma-Ray Emission Surrounding PSR J0622+3749 Observed by LHAASO-KM2A

We report the discovery of an extended very-high-energy (VHE) gamma-ray source around the location of the middle-aged (207.8 kyr) pulsar PSR J0622+3749 with the Large High-Altitude Air Shower Observatory (LHAASO). The source is detected with a significance of 8.2σ for E>25  TeV assuming a Gaussian template. The best-fit location is (right ascension, declination) =(95.47°±0.11°,37.92°±0.09°), and the extension is 0.40°±0.07°. The energy spectrum can be described by a power-law spectrum with an index of -2.92±0.17_{stat}±0.02_{sys}. No clear extended multiwavelength counterpart of the LHAASO source has been found from the radio to sub-TeV bands. The LHAASO observations are consistent with the scenario that VHE electrons escaped from the pulsar, diffused in the interstellar medium, and scattered the interstellar radiation field. If interpreted as the pulsar halo scenario, the diffusion coefficient, inferred for electrons with median energies of ∼160  TeV, is consistent with those obtained from the extended halos around Geminga and Monogem and much smaller than that derived from cosmic ray secondaries. The LHAASO discovery of this source thus likely enriches the class of so-called pulsar halos and confirms that high-energy particles generally diffuse very slowly in the disturbed medium around pulsars.

Ultrahigh-energy photons up to 1.4 petaelectronvolts from 12 γ-ray Galactic sources

The extension of the cosmic-ray spectrum beyond 1 petaelectronvolt (PeV; 10¹⁵ electronvolts) indicates the existence of the so-called PeVatrons-cosmic-ray factories that accelerate particles to PeV energies. We need to locate and identify such objects to find the origin of Galactic cosmic rays¹. The principal signature of both electron and proton PeVatrons is ultrahigh-energy (exceeding 100 TeV) γ radiation. Evidence of the presence of a proton PeVatron has been found in the Galactic Centre, according to the detection of a hard-spectrum radiation extending to 0.04 PeV (ref. ²). Although γ-rays with energies slightly higher than 0.1 PeV have been reported from a few objects in the Galactic plane^3-6, unbiased identification and in-depth exploration of PeVatrons requires detection of γ-rays with energies well above 0.1 PeV. Here we report the detection of more than 530 photons at energies above 100 teraelectronvolts and up to 1.4 PeV from 12 ultrahigh-energy γ-ray sources with a statistical significance greater than seven standard deviations. Despite having several potential counterparts in their proximity, including pulsar wind nebulae, supernova remnants and star-forming regions, the PeVatrons responsible for the ultrahigh-energy γ-rays have not yet been firmly localized and identified (except for the Crab Nebula), leaving open the origin of these extreme accelerators.

Role of N*(1535) in the Λ+c → K0ηp decay and the possible фp state in the Λ+c → π0фp decay

The nonleptonic weak decays of Λ+c → K0ηp and Λ+c → π0фp are investigated from the viewpoint of probing the N*(1535) resonance and the possible фp state. For the Λ+c → K0ηp decay, we study the invariant mass distribution of ηp with both the chiral unitary approach and an effective Lagrangian model. Within the chiral unitary approach, the N*(1535) resonance is dynamically generated from the final state interaction of mesons and baryons in coupled channels. While for the effective Lagrangian model, we take a Breit-Wigner formula for the N*(1535) resonance. We found that the behavior of the N*(1535) resonance in the Λ+c → K0N*(1535) → K0ηp decay within the two approaches is different. For the Λ+c → π0фp decay, we consider a triangle singularity mechanism, where the Λ+c decays into the K*Σ*(1385), the Σ*(1385) decays into the π0Σ/Λ, and then the K*Σ/Λ merge to produce the фp in the final state. This mechanism produces a peak structure around 2020 MeV. In addition, the possibility that there is a hidden-strange pentaquark-like state is also considered by taking into account the final state interactions of K*Λ, K*Σ, and фp. We conclude that it is difficult to search for the hidden-strange state in this decay. However, we do expect nontrivial behavior in the фp invariant mass distribution. The proposed Λ+c decay mechanism here can provide valuable information on the properties of these nuclear resonances and can in principle be tested by experiments such as BESIII, LHCb and Belle-II.

Emergence of a Complete Heavy-Quark Spin Symmetry Multiplet: Seven Molecular Pentaquarks in Light of the Latest LHCb Analysis

A recent analysis by the LHCb Collaboration suggests the existence of three narrow pentaquarklike states-the P_{c}(4312), P_{c}(4440), and P_{c}(4457)-instead of just one in the previous analysis [the P_{c}(4450)]. The closeness of the P_{c}(4312) to the D[over ¯]Σ_{c} threshold and the P_{c}(4440)  and  P_{c}(4457) to the D[over ¯]^{*}Σ_{c} threshold suggests a molecular interpretation of these resonances. We show that these three pentaquarklike resonances can be naturally accommodated in a contact-range effective field theory description that incorporates heavy-quark spin symmetry. This description leads to the prediction of all the seven possible S-wave heavy antimeson-baryon molecules [that is, there should be four additional molecular pentaquarks in addition to the P_{c}(4312), P_{c}(4440), and P_{c}(4457)], providing the first example of a heavy-quark spin symmetry molecular multiplet that is complete. If this is confirmed, it will not only give us an impressive example of the application of heavy-quark symmetries and effective field theories in hadron physics, it will also uncover a clear and powerful ordering principle for the molecular spectrum, reminiscent of the SU(3)-flavor multiplets to which the light hadron spectrum conforms.

Evidence for prevalent Z = 6 magic number in neutron-rich carbon isotopes

The nuclear shell structure, which originates in the nearly independent motion of nucleons in an average potential, provides an important guide for our understanding of nuclear structure and the underlying nuclear forces. Its most remarkable fingerprint is the existence of the so-called magic numbers of protons and neutrons associated with extra stability. Although the introduction of a phenomenological spin-orbit (SO) coupling force in 1949 helped in explaining the magic numbers, its origins are still open questions. Here, we present experimental evidence for the smallest SO-originated magic number (subshell closure) at the proton number six in 13-20C obtained from systematic analysis of point-proton distribution radii, electromagnetic transition rates and atomic masses of light nuclei. Performing ab initio calculations on 14,15C, we show that the observed proton distribution radii and subshell closure can be explained by the state-of-the-art nuclear theory with chiral nucleon-nucleon and three-nucleon forces, which are rooted in the quantum chromodynamics.

Leading SU(3)-breaking corrections to the baryon magnetic moments in chiral perturbation theory

We calculate the baryon magnetic moments using covariant chiral perturbation theory (chiPT) within the extended-on-mass-shell renormalization scheme. By fitting the two available low-energy constants, we improve the Coleman-Glashow description of the data when we include the leading SU(3)-breaking effects coming from the lowest-order loops. This success is in dramatic contrast with previous attempts at the same order using heavy-baryon chiPT and covariant infrared chiPT. We also analyze the source of this improvement with particular attention to the comparison between the covariant results.