Nuclear Physics Research Articles

Optimized FOPID controller for nuclear research reactor using enhanced planet optimization algorithm

Nuclear reactor control is pivotal for the safe and efficient operation of nuclear power plants, facilitating the regulation of reactor reactivity. This study introduces an optimized fractional-order proportional-integral-derivative (FOPID) controller tailored for maintaining reactivity levels in nuclear power plants, particularly during load-following operations. The controller adjusts the position of control rod to regulate power output effectively. We enhance FOPID controller's performance using a modification of Planet Optimization Algorithm (POA-M), leveraging the strengths of the Arithmetic Optimization Algorithm (AOA) to improve its exploitation capabilities. We evaluate the efficacy of POA-M-FOPID controller against traditional techniques, including POA, AOA, and Particle Swarm Optimization (PSO). We assess its performance using the Egyptian Testing Research Reactor (ETRR-2) as a case study. Our results demonstrate that the POA-M-FOPID controller outperforms alternative algorithms across various control metrics, exhibiting superior resilience and efficiency. Notably, the utilization of the POA-M-FOPID controller yields remarkable improvements in reactor power performance, achieving significantly reduced settling time (25.27 sec) and maximum overshoot (0.67 %) compared to conventional FOPID controllers incorporating POA, AOA, and PSO methods. These findings underscore the effectiveness of POA-M-FOPID in enhancing nuclear reactor control systems, offering potential benefits for broader nuclear power industry in terms of safety, stability, and operational efficiency.

Pseudo-Spin Symmetry and the Hints for Unstable and Superheavy Nuclei

The pseudo-spin symmetry (PSS) provides an important angle to understand nuclear microscopic structure and the novel phenomena found in unstable nuclei. The relativistic Hartree–Fock (RHF) theory, that takes the important degrees of freedom associated with the π-meson and ρ-tensor (ρ-T) couplings into account, provides an appropriate description of the PSS restoration in realistic nuclei, particularly for the pseudo-spin (PS) doublets with high angular momenta (l˜). The investigations of the PSS within the RHF theory are recalled in this paper by focusing on the effects of the Fock terms. Aiming at common artificial shell closures appearing in previous relativistic mean-field calculations, the mechanism responsible for the PSS restoration of high-l˜ orbits is stressed, revealing the manifestation of nuclear in-medium effects on the PSS, and thus, providing qualitative guidance on modeling the in-medium balance between nuclear attractions and repulsions. Moreover, the essential role played by the ρ-T coupling, that contributes mainly via the Fock terms, is introduced as combined with the relations between the PSS and various nuclear phenomena, including the shell structure and the evolution, novel halo and bubble-like phenomena, and the superheavy magicity. As the consequences of the nuclear force in complicated nuclear many-body systems, the PSS itself and the mechanism therein can not only deepen our understanding of nuclear microscopic structure and relevant phenomena, but also provide special insight into the nature of the nuclear force, which can further enrich our knowledge of nuclear physics.

70Zn-induced fusion reactions for synthesis of superheavy element Z=120

In this work, we delved into the intriguing realm of nuclear physics by investigating fusion reactions initiated by [Formula: see text]Zn projectiles colliding with various stable isotopes of thorium nuclei, including [Formula: see text]Th. Our primary objective was the synthesis of the superheavy element [Formula: see text]. To achieve this, we employed an advanced statistical model to evaluate the evaporation residue cross-sections. Remarkably, we observed that the fusion reaction between [Formula: see text]Zn and [Formula: see text]Th exhibited the highest evaporation residue cross-sections, reaching a peak value. Furthermore, the 3n channel within this fusion reaction yielded substantially larger evaporation cross-sections, with a notable value of 0.89 picobarns. The predicted cross-sections fall within the measurable range, making them highly valuable for experimental endeavors aimed at synthesizing the superheavy element [Formula: see text]. This research not only contributes to our understanding of nuclear reactions but also holds promise for expanding our knowledge of the heaviest elements in the periodic table, potentially opening new frontiers in nuclear science.

Predicting quadrupole deformation via anisotropic flow and transverse momentum spectra in isotopic 54128-135Xe\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\mathbf {{}^{128-135}_{\\qquad \\,54}\ extrm{Xe}}$$\\end{document} collisions at LHC

In the hydrodynamical description of heavy-ion collisions, the elliptic flow v2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${\ extrm{v}}_{2}$$\\end{document} and triangular flow v3\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${\ extrm{v}}_{3}$$\\end{document} are sensitive to the quadrupole deformation β2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${\\beta _{2}}$$\\end{document} of the colliding nuclei. We produce v2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${\ extrm{v}}_{2}$$\\end{document} and v3\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${\ extrm{v}}_{3}$$\\end{document} ratios qualitatively and quantitatively in most-central Xe–Xe collisions at 5.44 TeV. By employing HYDJET++ model, we study the sensitivity of anisotropic flow coefficients and mean transverse momentum to the quadrupole deformation and system-size in isotopic Xe–Xe collisions. Flow observables strongly depend on the strength of nucleon–nucleon scattering occurring in even-A and odd-A nuclei. Flow for odd-A nuclei is suppressed in comparison to flow in even-A collisions. There exists a linear inter-dependence between pT\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${\ extrm{p}}_{\ extrm{T}}$$\\end{document} integrated anisotropic flow and nuclear deformation. Mean transverse momentum signifies the fireball temperature in body–body and tip–tip collisions. There exists a negative linear correlation of ⟨pT⟩\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\langle {\ extrm{p}}_{\ extrm{T}}\\rangle $$\\end{document} with collision system-size and a positive correlation with nuclear deformation. Flow measurements in high-energy, heavy-ion collisions using isotopic collision systems, offer a new precision tool to study nuclear structure physics. Observation of nuclear structure properties like nuclear deformation in a heavy-ion collision such as this would be very interesting.

Measurement of the neutron total cross section of 9Be at the Back-n white neutron source of CSNS

Abstract The neutron total cross section data of 9Be play an important role in nuclear structure model research of light nuclei and nuclear power installations. A measurement of the neutron total cross section of 9Be in the 0.3 eV-120 MeV energy region has been carried out by using time-of-flight method and transmission method with the Neutron Total Cross Sectional Spectrometer (NTOX) based on the multi-cell fast fission chamber at the China Spallation Neutron Source (CSNS)- Back-n white neutron source (Back-n). The fission count-neutron energy distributions of 235U and 238U without sample and with Be samples in three thicknesses were measured in the double-bunch operation mode during beam power of 100 kW. The Bayesian method is used to eliminate the influence of the double-bunch problem on neutron measurement in the energy region above 10 keV. The neutron total cross section of 9Be results is well consistent with ENDF/B-VIII.0 evaluation library data in 0.3 eV-20 MeV energy region. In the energy range of 0.3 eV to 10 keV, the deviation between our results and the evaluation results of ENDF/B-VIII.0 is within 2.5%, which is within 15% from 0.01-20 MeV. In the resonance energy region, the measured resonance energies in our experiment are 0.63 MeV, 0.82 MeV and 2.8 MeV respectively. The results show that the total cross section uncertainties of three thicknesses Be samples are within 2.2% in the energy region below 1 MeV. The total cross section uncertainty of 30 mm Be from 235U is the smallest and less than 5% in the energy region 0.3 eV-120 MeV. The results of this experiment can provide technical support for further data analysis and related nuclear data evaluation.

A new spectrally-invariant approach to the remote sensing of inhomogeneous clouds

Current operational satellite retrievals of cloud optical and microphysical properties go back to the Nakajima-King technique developed at the end of the 1980 s. This technique is based on library calculations for plane-parallel homogeneous clouds. It often works well for overcast skies but leads to substantial errors for inhomogeneous and broken cloud fields where the plane-parallel-geometry assumption is no longer valid. The basic concept of a new technique was first introduced in nuclear physics to quantify the critical conditions of reactors. Based on the eigenvalues of the radiative transfer equation, their approach provides a powerful means to parameterize the structure of 3D media. This parameterization was later successfully applied to relate surface reflectance spectra to 3D canopy structure and is known as the spectrally-invariant approximation. The proposed approach adapts this technique to the remote sensing of cloud properties such as droplet single scattering albedo and the average number of scatterings (which are the fundamental parameters in radiative transfer theory), with an emphasis on quantifying the associated errors and uncertainties. This retrieval is free from the plane-parallel homogeneous cloud assumption.

Overview of neutral beam injectors for plasma heating and diagnostics developed at Budker INP

An overview of the neutral beam injectors developed at the Budker Institute of Nuclear Physics in Novosibirsk during the last 10 years is presented. These neutral injectors are used for plasma diagnostics, heating and current drive in modern fusion devices with magnetic confinement. An arc or a radio-frequency (RF) discharge generates a plasma in the ion sources of the injectors, and a positive hydrogen or deuterium ion beam is extracted and accelerated by a multiaperture ion-optical system (IOS). The accelerated ion beam is converted into a neutral one in a gas target. The precision multiaperture IOS with spherically concave electrodes provides ballistic focusing of the neutral beam. The high-energy, high-power beam injector based on negative ions, which is currently under development, is described as well. It comprises a RF negative ion source and a wide-aperture electrostatic accelerator separated from the source by a low-energy beam transport line, thereby improving the injector reliability.

Statistical investigation of the angular momentum dependence of the nuclear level density parameter

The nuclear level density (NLD) parameter is crucial for calculating cross-sections in nuclear physics, astrophysics, reactor design, and medical physics. Spin and parity, along with excitation energy, are fundamental properties of an excited nuclear level. Previous investigations into the NLD’s dependence on spin and parity have primarily used approximate methods like parity equidistribution and Gaussian distribution of spins. However, the specific impact of spin and parity on the NLD parameter, a key component in NLD formulation, has not been extensively studied. We examined the spin and parity dependence of the NLD parameter. Our findings demonstrate that the NLD parameter’s dependency on both excitation energy and angular momentum can be accurately characterized by a Laplace distribution, highlighting the complex interplay of these factors in nuclear physics.

Inverse problem in energy-dependent potentials using semiclassical methods

Wave equations with energy-dependent potentials appear in many areas of physics, ranging from nuclear physics to black hole perturbation theory. In this work, we use the semiclassical Wentzel-Kramers-Brillouin (WKB) method to first revisit the computation of bound states of potential wells and reflection/transmission coefficients in terms of the Bohr-Sommerfeld rule and the Gamow formula. We then discuss the inverse problem, in which the latter observables are used as a starting point to reconstruct the properties of the potentials. By extending known inversion techniques to energy-dependent potentials, we demonstrate that so-called width-equivalent or WKB-equivalent potentials are not isospectral anymore. Instead, we explicitly demonstrate that constructing quasi-isospectral potentials with the inverse techniques is still possible. Those reconstructed, energy-independent potentials share key properties with the width-equivalent potentials. We report that including energy-dependent terms allows for a rich phenomenology, particularly for the energy-independent equivalent potentials. Published by the American Physical Society 2024

Survey Analysis of Potential Nuclear Safety Research of Thailand for International Research Collaborative Reinforcement in the 2020s

Abstract To achieve the long-term challenge of nuclear energy public acceptance in Thailand, nuclear safety research needed to be properly determined in both domestic and international directions, especially in the 2020s which was a period passing the Fukushima disaster over 10 years. Thailand Institute of Nuclear Technology (TINT) has studied nuclear safety research after the Fukushima accident to answer technical and social issues of nuclear power. An update of nuclear safety research from domestic experts and international surveys was needed in order to identify potential collaborative research to serve the goal of public acceptance reinforcement. The objective of this study was to survey, assess and rank the importance and knowledge level of nuclear safety research in Thailand among domestic experts in various fields. The survey was extended to collect the opinion of international participants of the ASEAN Network on Nuclear Power Safety Research (ASEAN NPSR) to analyze the similarity of the nuclear research interest for reinforcing the future collaborative project. As a result, the importance and knowledge level showed diverse important research topics with the priority of research scopes on human factor novel reactor technologies, and risk assessment. According to the ASEAN NPSR survey, the nuclear safety research of severe accidents, risk assessment, and novel reactor technologies were listed as potential collaborative projects. Also, the domestic and ASEAN NPSR survey results helped support the new collaborative research extension session in the annual ASEAN NPSR meeting to together discuss the potential nuclear safety research between members for the 2020s.

Accessing new physics with an undoped, cryogenic CsI CEvNS detector for COHERENT at the SNS

We consider the potential for a 10 kg undoped cryogenic CsI detector operating at the Spallation Neutron Source to measure coherent elastic neutrino-nucleus scattering and its sensitivity to discover new physics beyond the standard model (BSM). Through a combination of increased event rate, lower threshold, and good timing resolution, such a detector would significantly improve on past measurements. We considered tests of several BSM scenarios such as neutrino nonstandard interactions and accelerator-produced dark matter. This detector’s performance was also studied for relevant questions in nuclear physics and neutrino astronomy, namely the weak charge distribution of Cs and I nuclei and detection of neutrinos from a core-collapse supernova. Published by the American Physical Society 2024

The spatial resolution measurements on the small prototype of the Super Charm-Tau Factory drift chamber

Budker Institute of Nuclear Physics is actively evolving the Super Charm-Tau Factory project. A drift chamber with a small hexagonal cell was proposed for this collider. For experimental measurements of a spatial resolution, a drift chamber small prototype consisting of seven cells has been developed. The experiments have been carried out in He/C3H8 (60/40) and He/C2H6 (50/50) gas mixtures at different gas gains using cosmic particles. For this purpose, a calibration procedure for R(t, α) has been developed and fine-tuned, taking into account small prototype features and an angular dependence of isochrones at the edge of the cells. As a result, the average spatial resolutions are less than 100μm for gas gains greater than 5⋅104.

Pulse shape discrimination technique for diffuse supernova neutrino background search with JUNO

Pulse shape discrimination (PSD) is widely used in particle and nuclear physics. Specifically in liquid scintillator detectors, PSD facilitates the classification of different particle types based on their energy deposition patterns. This technique is particularly valuable for studies of the diffuse supernova neutrino background (DSNB), nucleon decay, and dark matter searches. This paper presents a detailed investigation of the PSD technique, applied in the DSNB search performed with the Jiangmen Underground Neutrino Observatory (JUNO). Instead of using conventional cut-and-count methods, we employ methods based on boosted decision trees and neural networks and compare their capability to distinguish the DSNB signals from the atmospheric neutrino neutral-current background events. The two methods demonstrate comparable performance, resulting in a 50–80% improvement in signal efficiency compared to a previous study performed for JUNO (An et al. [JUNO] in J Phys G 43(3):030401, 2016). Moreover, we study the dependence of the PSD performance on the visible energy and final state composition of the events and find a significant dependence on the presence/absence of 11\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$^{11}$$\\end{document}C. Finally, we evaluate the impact of the detector effects (photon propagation, PMT dark noise, and waveform reconstruction) on the PSD performance.

Development of a fully implicit ODE-solver for containment analysis code

The thermal–hydraulic dynamics in containment are governed by a system of stiff ordinary differential equations (ODEs). A fully implicit discretization scheme is adopted to discretize these ODEs in order to mitigate the effects of stiffness. In comparison with explicit or semi-implicit discretization schemes that are subject to Courant limits on time steps, the fully implicit discretization scheme is more suitable for a containment analysis code that focuses on predicting both short-term and long-term thermal–hydraulic parameters after an accident. This study introduces a general-purpose ODE solver for the containment analysis code. The outline of the solver is as follows: The fully implicit discrete equations lead to a large set of nonlinear equations that need to be solved using Newton’s iterative method. The partial derivative components in the Jacobi matrix are calculated by the perturbation method using finite difference approximation, which avoids the complicated derivation of partial derivatives. The scaling modification technique is incorporated into this ODE solver to deal with significant differences in unknown variable magnitudes, and the line search method is introduced to address the difficulty of obtaining an accurate root estimate with Newton’s method when the initial guess is far from the actual root. This proposed ODE solver was applied to two typical stiff ODE problems to test its stiffness-suppressed ability and to demonstrate that this proposed solver can perform calculations with a very large time step. Then, the CASSIA code, a containment analysis code developed by China Nuclear Power Technology Research Institute Co., Ltd (CNPRI), equipped with this ODE solver, was applied to the CSNI (Committee on the Safety of Nuclear Installations) benchmark problem and the Carolinas Virginia Tube Reactor (CVTR) test 3 problem to preliminarily demonstrate that the proposed ODE solver can perform containment thermal–hydraulic analysis correctly. This study could provide references for the development of a home-made containment analysis code.

An analytic description of electron thermalization in kilonovae ejecta

Abstract A simple analytic description is provided of the rate of energy deposition by β-decay electrons in the homologously expanding radioactive plasma ejected in neutron star mergers, valid for a wide range of ejecta parameters- initial entropy, electron fraction {s0, Ye} and scaled (time-independent) density ρt3. The formulae are derived using detailed numerical calculations following the time-dependent composition and β-decay emission spectra (including the effect of delayed deposition). The deposition efficiency depends mainly on ρt3 and only weakly on {s0, Ye}. The time te at which the ratio between the rates of electron energy deposition and energy production drops to 1 − e−1, is given by $t_e=t_{0e}\Big (\frac{\rho t^3}{0.5(\rho t^3)_0}\Big )^a$, where $(\rho t^3)_0=\frac{0.05\, {\rm M}_{\odot }}{4\pi (0.2c)^3}$, t0e(s0, Ye) ≈ 17 days and 0.4 ≤ a(s0, Ye) ≤ 0.5. The fractional uncertainty in te due to nuclear physics uncertainties is $\approx 10~{{\%}}$. The result a ≤ 0.5 reflects the fact that the characteristic β-decay electron energies do not decrease with time (largely due to ‘inverted decay chains’ in which a slowly-decaying isotope decays to a rapidly-decaying isotope with higher end-point energy). We provide an analytic approximation for the time-dependent electron energy deposition rate, reproducing the numerical results to better than 50% (typically $<30~{{\%}}$, well within the energy production rate uncertainty due to nuclear physics uncertainties) over a 3-4 orders-of-magnitude deposition rate decrease with time. Our results may be easily incorporated in calculations of kilonovae light curves (with general density and composition structures), eliminating the need to numerically follow the time-dependent electron spectra. Identifying te, e.g. in the bolometric light curve, will constrain the ejecta density distribution.

The Second-Order Features Adjoint Sensitivity Analysis Methodology for Response-Coupled Forward/Adjoint Linear Systems (2nd-FASAM-L): Mathematical Framework and Illustrative Application to an Energy System

The Second-Order Features Adjoint Sensitivity Analysis Methodology for Response-Coupled Forward/Adjoint Linear Systems (abbreviated as “2nd-FASAM-L”), presented in this work, enables the most efficient computation of exactly obtained mathematical expressions of first- and second-order sensitivities of a generic system response with respect to the functions (“features”) of model parameters. Subsequently, the first- and second-order sensitivities with respect to the model’s uncertain parameters, boundaries, and internal interfaces are obtained analytically and exactly, without needing large-scale computations. Within the 2nd-FASAM-L methodology, the number of large-scale computations is proportional to the number of model features (defined as functions of model parameters), as opposed to being proportional to the number of model parameters. This characteristic enables the 2nd-FASAM-L methodology to maximize the efficiency and accuracy of any other method for computing exact expressions of first- and second-order response sensitivities with respect to the model’s features and/or primary uncertain parameters. The application of the 2nd-FASAM-L methodology is illustrated using a simplified energy-dependent neutron transport model of fundamental significance in nuclear reactor physics.

Assessment of radiological consequences to a hypothetical accident of the 3-MW TRIGA Mark-II nuclear research reactor of Bangladesh

The radiological consequences to a hypothetical accident are assessed for the 3-MW TRIGA Mark-II Nuclear Research Reactor of Bangladesh for both wet and dry seasons. Accident simulations are crucial for assessing the potentialrisks in the event of a nuclear accident. ORIGEN code was used to carry out the facility's radionuclide inventory. In the event of a hypothetical accident, the Total Effective Dose Equivalent (TEDE) is assessed using HotSpots. The findings indicate that for stability class A during the rainy season, the selected radionuclides arrived at a downwind distance of 0.03 km with an arrival time of less than one minute, and the highest TEDE, respiratory time-integrated air concentration, and ground deposition activity, respectively, were marked at 2.8E-12 Sv, 9.4E-01 Bq-sec/m3, and 2.4E-05 kBq/m2. In the case of I-131 at 0.03 km from the source, the thyroid gland is the most sensitive human organ, followed by the lung, esophagus, thymus, surface bone, skin, and breast. The estimated Committed Effective Dose Equivalent (CEDE) from I-131 is 9.1E-11 Sv. In contrast, the adrenal was the most sensitive organ with a maximum CEDE of 1.3E-13 Sv in the case of Cs-137 at a downwind distance of 0.03 km. The highest TEDE values for the general public and the vicinity around the TRIGA Mark II reactor facility are lower than the annual dose limits.

Commentary on the current status of nuclear medicine medical physics expert support in the UK.

Commentary on the current status of nuclear medicine medical physics expert support in the UK.

Automated thickness measurement system for nuclear fuel plates: A metrological study

Thickness measurements of nuclear fuel plates manufactured at the Nuclear Fuel Center of the Nuclear and Energy Research Institute (IPEN-CNEN/SP) are currently performed using manual external U-shaped frame micrometers with non-rotating spindles and beveled measurement anvils. Results are then manually recorded on dedicated forms. To improve the efficiency and the metrological reliability of the process, a novel thickness measurement system for nuclear fuel plates was developed. This work presents a comprehensive metrological study of the newly designed measurement system, achieved through the construction of a new apparatus with automated mechanisms for synchronized movement, data capture, transcription, and processing. The study results demonstrate the metrological adequacy of the developed system and underscore the importance of the design and structural quality of the measurement equipment. Consequently, the metrological reliability of the results in the dimensional control of nuclear fuel plates is confirmed.

Alanine dosimeters for LET measurement in proton radiotherapy

Proton radiotherapy enables highly conformal dose delivery to the tumor while minimizing exposure to healthy tissues surrounding the target. However, to achieve the most optimal treatment plan, it is necessary to consider various parameters, including variable relative biological effectiveness (RBE) and linear energy transfer (LET) within the target and in the normal tissues. The proton treatment planning systems (TPS) currently under development incorporate both of these parameters. However, before introducing biologically guided treatment planning into the clinic, commissioning and end-to-end tests must be conducted. These procedures require specific tools to verify the planned LET distribution. The same tools will also be necessary for subsequent patient-specific quality assurance (QA).The objective of this work was to evaluate the alanine/EPR dosimetry system developed at the Cyclotron Center Bronowice at the Institute of Nuclear Physics (CCB IFJ PAN) for verifying LET values in therapeutic proton fields. The change in the shape of the EPR spectra of alanine irradiated in proton fields of different LETs, described by the ratio of the intensities of the satellite and central lines (x/y ratio), was used to discriminate LET values. The method does not provide absolute LET values. However, as demonstrated, it can verify the LET values calculated in the TPS after calibration using Monte Carlo (MC) simulation.