A novel method for simulating nuclear explosion with chemical explosion to form an approximate plane wave: Field test and numerical simulation

A nuclear explosion in the rock mass medium can produce strong shock waves,seismic shocks,and other destructive effects,which can cause extreme damage to the underground protection infrastructures.With the increase in nuclear explosion power,underground protection engineering enabled by explosion-proof impact theory and technology ushered in a new challenge.This paper proposes to simulate nuclear explosion tests with on-site chemical explosion tests in the form of multi-hole explosions.First,the mechanism of using multi-hole simultaneous blasting to simulate a nuclear explosion to generate approximate plane waves was analyzed.The plane pressure curve at the vault of the underground protective tunnel under the action of the multi-hole simultaneous blasting was then obtained using the impact test in the rock mass at the site.According to the peak pressure at the vault plane,it was divided into three regions:the stress superposition region,the superposition region after surface reflection,and the approximate plane stress wave zone.A numerical simulation approach was developed using PFC and FLAC to study the peak particle velocity in the surrounding rock of the underground protective cave under the action of multi-hole blasting.The time-history curves of pressure and peak pressure partition obtained by the on-site multi-hole simultaneous blasting test and numerical simulation were compared and analyzed,to verify the correctness and rationality of the formation of an approximate plane wave in the simulated nuclear explosion.This comparison and analysis also provided a theoretical foundation and some research ideas for the ensuing study on the impact of a nuclear explosion.

Genetic evaluations of island populations of Rhesus Macaque (Macaca mulatta) in China: Implications for conservation management

DEAR EDITOR,The macaques belongs to the genus Macaca,consisting of at least 23 species(Roos et al.,2019).Among all congeners,rhesus macaque(M.mulatta)is regarded as the widest distributed non-human primate species in the world.Its native range spans in East Asia,northern part of Southeast Asia and Indian subcontinent(Liu et al.,2018).Listed as“Least Concern”on the IUCN Red List,this species is locally threatened due to habitat loss and degradation in China and Thailand(Lu et al.,2018).Nevertheless,pet release resulting in hybridization with other congeners(e.g.,rhesus macaque×crab-eating macaque(M.fascicularis))was documented in Hong Kong SAR,China(Wong&Ni,2000),threatening genetic integrity of wild populations.

Hydromechanical characterization of gas transport amidst uncertainty for underground nuclear explosion detection

Given the challenge of definitively discriminating between chemical and nuclear explosions using seismic methods alone,surface detection of signature noble gas radioisotopes is considered a positive identification of underground nuclear explosions(UNEs).However,the migration of signature radionuclide gases between the nuclear cavity and surface is not well understood because complex processes are involved,including the generation of complex fracture networks,reactivation of natural fractures and faults,and thermo-hydro-mechanical-chemical(THMC)coupling of radionuclide gas transport in the subsurface.In this study,we provide an experimental investigation of hydro-mechanical(HM)coupling among gas flow,stress states,rock deformation,and rock damage using a unique multi-physics triaxial direct shear rock testing system.The testing system also features redundant gas pressure and flow rate measurements,well suited for parameter uncertainty quantification.Using porous tuff and tight granite samples that are relevant to historic UNE tests,we measured the Biot effective stress coefficient,rock matrix gas permeability,and fracture gas permeability at a range of pore pressure and stress conditions.The Biot effective stress coefficient varies from 0.69 to 1 for the tuff,whose porosity averages 35.3%±0.7%,while this coefficient varies from 0.51 to 0.78 for the tight granite(porosity<1%,perhaps an underestimate).Matrix gas permeability is strongly correlated to effective stress for the granite,but not for the porous tuff.Our experiments reveal the following key engineering implications on transport of radionuclide gases post a UNE event:(1)The porous tuff shows apparent fracture dilation or compression upon stress changes,which does not necessarily change the gas permeability;(2)The granite fracture permeability shows strong stress sensitivity and is positively related to shear displacement;and(3)Hydromechanical coupling among stress states,rock damage,and gas flow appears to be stronger in tight granite than in porous tuff.

Seismic response characteristics of nuclear island structure at generic soil and rock sites

The number of traditionally excellent coastal lithologic nuclear power plants is limited.It is a trend to develop nuclear power plants on soil sites in inland areas.Therefore,the seismic safety and adaptability of non-rock nuclear power plant(NPP)sites are the key concerns of nuclear safety researchers.Although the five site categories are clearly defined in the AP1000 design control documents,the effects of nuclear power five site conditions and soil nonlinearity on the seismic response characteristics of nuclear island buildings have not been systematically considered in previous related studies.In this study,targeting the AP1000 nuclear island structure as the research object,three-dimensional finite element models of a nuclear island structure at five types of sites(firm rock site(FR),soft rock site(SR),soft-to-medium soil site(SMS),upper bound soft-to-medium site(SMS-UB),and soft soil site(SS))are established.The partitioned analysis method of soil-structure interaction(PASSI)in the time-domain is used to investigate the effects of site hardness and nonlinearity on the acceleration,displacement,and acceleration response spectrum of the nuclear island structure under seismic excitation.The incremental equilibrium equation and explicit decoupling method are used to analyze the soil nonlinearity described by the Davidenkov model with simplified loading-reloading rules.The results show that,in the linear case,with the increase of site hardness,the peak ground acceleration(PGA)and the peak of acceleration response spectrum of the nuclear island structure increase except for the FR site,while the maximum displacement decreases.In nonlinear analysis,as the site hardness increases,the PGA,maximum displacement,and the peak of acceleration response spectrum of the nuclear island structure increase.The peak value of the acceleration response spectrum in the nonlinear case is greater than that in the linear case for FR,while smaller for SR and soil sites.The site nonlinearity reduces the peak values of the response spectrum for SR and soil sites much more as the site hardness decreases.The results of this study can provide a reference for design of nuclear island structures on soil and rock sites.

Presence of a long nuclear-localization signal sequence in homeodomain transcription factor Nkx 1.2

Homeodomains,a 60-amino acid sequence encoded by 180 nucleotides,are highly conserved DNA-binding motifs that are present in a variety of transcription factors in species ranging from yeast to humans.The NKX proteins belong to the homeodomain(HD)-containing transcription factor family.They play vital roles in the regulation of morphogenesis.NKX1-2 is one member of the NKX subfamily.At present,information about its nuclear localization signal(NLS)sequence is limited.We studied the NLS sequence of zebrafish Nkx1.2 by introducing sequence changes such as deletion,mutation,and truncation,and identified an NLS motif(QNRRTKWKKQ)that is localized at the C-terminus of the homeodomain.Moreover,the deletion of two amino acid residues(RR)in this NLS motif prevents Nkx1.2 from entering the nucleus,indicating that the two amino acids are essential for Nkx1.2 nuclear localization.However,the NLS motif alone is unable to target cytoplasmic protein glutathione S-transferase(GST)to the nucleus.An intact homeodomain is necessary for mediating the complete nuclear transport of cytoplasmic protein.Unlike most nuclear import proteins with short NLS sequences,a long NLS is present in zebrafish Nkx1.2.We also demonstrated that the sequences of homeodomain of NKX1.2 are well conserved among different species.This study is informative to verify the function of the NKX1.2 protein.

Drivers,Trends,and Patterns of Changing Vegetation-greenness in Nansha Islands,China from 2016 to 2022

Changes in vegetation status generally also represents changes in the ecological health of islands and reefs(IRs).However,studies are limited of drivers and trends of vegetation change of Nansha Islands,China and how they relate to climate change and human activities.To resolve this limitation,we studied changes to the Normalized Difference Vegetation Index(NDVI)vegetation-greenness index for 22 IRs of Nansha Islands during normal and extreme conditions.Trends of vegetation greenness were analyzed using Sen's slope and Mann-Kendall test at two spatial scales(pixel and island),and driving factor analyses were performed by time-lagged partial correlation analyses.These were related to impacts from human activities and climatic factors under normal(temperature,precipitation,radiation,and Normalized Difference Built-up Index(NDBI))and extreme conditions(wind speed and latitude of IRs)from 2016 to 2022.Results showed:1)among the 22 IRs,NDVI increased/decreased significantly in 15/4 IRs,respectively.Huayang Reef had the highest NDVI change-rate(0.48%/mon),and Zhongye Island had the lowest(–0.29%/mon).Local spatial patterns were in one of two forms:dotted-form,and degradation in banded-form.2)Under normal conditions,human activities(characterized by NDBI)had higher impacts on vegetation-greenness than other factors.3)Under extreme conditions,wind speed(R^(2)=0.2337,P<0.05)and latitude(R^(2)=0.2769,P<0.05)provided limited explanation for changes from typhoon events.Our results provide scientific support for the sustainable development of Nansha Islands and the United Nations‘Ocean Decade’initiative.

Local Climate Change Induced by Urbanization on a South China Sea Island

The South China Sea is a hotspot for regional climate research.Over the past 40 years,considerable improvement has been made in the development and utilization of the islands in the South China Sea,leading to a substantial change in the land-use of the islands.However,research on the impact of human development on the local climate of these islands is lacking.This study analyzed the characteristics of local climate changes on the islands in the South China Sea based on data from the Yongxing Island Observation Station and ERA5 re-analysis.Furthermore,the influence of urbanization on the local climate of the South China Sea islands was explored in this study.The findings revealed that the 10-year average temperature in Yongxing Island increased by approximately 1.11℃from 1961 to 2020,and the contribution of island development and urbanization to the local warming rate over 60 years was approximately 36.2%.The linear increasing trend of the annual hot days from 1961–2020 was approximately 14.84 days per decade.The diurnal temperature range exhibited an increasing trend of 0.05℃per decade,whereas the number of cold days decreased by 1.06days per decade.The rapid increase in construction on Yongxing Island from 2005 to 2021 led to a decrease in observed surface wind speed by 0.32 m s^(-1)per decade.Consequently,the number of days with strong winds decreased,whereas the number of days with weak winds increased.Additionally,relative humidity exhibited a rapid decline from 2001 to 2016 and then rebounded.The study also found substantial differences between the ERA5 re-analysis and observation data,particularly in wind speed and relative humidity,indicating that the use of re-analysis data for climate resource assessment and climate change evaluation on island areas may not be feasible.

An improved analysis method for assessing the nuclear-heating impact on the stability of toroidal field magnets in fusion reactors

The superconducting magnet system of a fusion reactor plays a vital role in plasma confinement,a process that can be dis-rupted by various operational factors.A critical parameter for evaluating the temperature margin of superconducting magnets during normal operation is the nuclear heating caused by D-T neutrons.This study investigates the impact of nuclear heat-ing on a superconducting magnet system by employing an improved analysis method that combines neutronics and thermal hydraulics.In the magnet system,toroidal field(TF)magnets are positioned closest to the plasma and bear the highest nuclear-heat load,making them prime candidates for evaluating the influence of nuclear heating on stability.To enhance the modeling accuracy and facilitate design modifications,a parametric TF model that incorporates heterogeneity is established to expedite the optimization design process and enhance the accuracy of the computations.A comparative analysis with a homogeneous TF model reveals that the heterogeneous model improves accuracy by over 12%.Considering factors such as heat load,magnetic-field strength,and cooling conditions,the cooling circuit facing the most severe conditions is selected to calculate the temperature of the superconductor.This selection streamlines the workload associated with thermal-hydraulic analysis.This approach enables a more efficient and precise evaluation of the temperature margin of TF magnets.Moreover,it offers insights that can guide the optimization of both the structure and cooling strategy of superconducting magnet systems.

Coulomb-assisted nonlocal electron transport between two pairs of Majorana bound states in a superconducting island

We investigate the nonlocal transport modulated by Coulomb interactions in devices comprising two interacting Majorana wires,where both nanowires are in proximity to a mesoscopic superconducting(SC)island.Each Majorana bound state(MBS)is coupled to one lead via a quantum dot with resonant levels.In this device,the nonlocal correlations can be induced in the absence of Majorana energy splitting.We find that the negative differential conductance and giant current noise cross correlation could be induced,due to the interplay between nonlocality of MBSs and dynamical Coulomb blockade effect.This feature may provide a signature for the existence of the MBSs.

U-Pb zircon geochronology of basaltic pyroclastic rocks from the basement beneath the Xisha Islands in the northwestern South China Sea and its geological significance

As one of the micro-blocks dispersed in the South China Sea(SCS), the basement of the Xisha Islands has rarely been drilled because of the thick overlying Cenozoic sediments, which has led to a confused understanding of the pre-Cenozoic basement of the Xisha Islands. Well CK-1, a kilometer-scale major scientific drill in the Xisha Islands in the northwestern SCS, penetrated thick reefal limestone(0–888.4 m) and the underlying basement rocks(888.4–901.4 m). In this study, we present the zircon U-Pb ages of basement basaltic pyroclastic rocks from Well CK-1 in the Xisha Islands of the northwestern SCS to investigate the basement nature of the Xisha microblock. The basement of Well CK-1 consists of basaltic pyroclastic rocks on the seamount. The zircon grains yielded apparent ages ranging from ca. 2 138.9 Ma to ca. 36 Ma. The old group of zircon grains from Well CK-1 was considered to be inherited zircons. Two Cenozoic zircons gave a weighted mean 206Pb/238U age of(36.3 ± 1.1) Ma,Mean Squared Weighted Deviations(MSWD) = 1.2, which may represent the maximum age of the volcano eruption. The Yanshanian inherited zircons(116.9–105.7 Ma and 146.1–130.2 Ma) from Well CK-1 are consistent with the zircons from Well XK-1, indicating that the basement of Chenhang Island may be similar to that of Well XK-1. We propose that the Xisha micro-block may have developed on a uniform Late Jurassic metamorphic crystalline basement, intruded by Cretaceous granitic magma.

Impact of initial fluctuations and nuclear deformations in isobar collisions

Relativistic isobar^(96)_(44)Ru+^(96)_(44)Ru and^(96)_(40)Zr+^(96)_(40)Zrcollisions have revealed intricate differences in their nuclear size and shape,inspiring unconventional studies of nuclear structure using relativistic heavy ion collisions.In this study,we investigate the relative differences in the mean multiplicityR_(<Nch>)and the secondR_(ε2)and third-order eccentricityR_(ε3)between isobar collisions using initial state Glauber models.It is found that initial fluctuations and nuclear deformations have negligible effects on R_(<Nch>)in most central collisions,while both are important for the R_(ε2)and R_(ε3),the degree of which is sensitive to the underlying nucleonic or sub-nucleonic degree of freedom.These features,compared to real data,may probe the particle production mechanism and the physics underlying nuclear structure.

A fast forward computational method for nuclear measurement using volumetric detection constraints

Owing to the complex lithology of unconventional reservoirs,field interpreters usually need to provide a basis for interpretation using logging simulation models.Among the various detection tools that use nuclear sources,the detector response can reflect various types of information of the medium.The Monte Carlo method is one of the primary methods used to obtain nuclear detection responses in complex environments.However,this requires a computational process with extensive random sampling,consumes considerable resources,and does not provide real-time response results.Therefore,a novel fast forward computational method(FFCM)for nuclear measurement that uses volumetric detection constraints to rapidly calculate the detector response in various complex environments is proposed.First,the data library required for the FFCM is built by collecting the detection volume,detector counts,and flux sensitivity functions through a Monte Carlo simulation.Then,based on perturbation theory and the Rytov approximation,a model for the detector response is derived using the flux sensitivity function method and a one-group diffusion model.The environmental perturbation is constrained to optimize the model according to the tool structure and the impact of the formation and borehole within the effective detection volume.Finally,the method is applied to a neutron porosity tool for verification.In various complex simulation environments,the maximum relative error between the calculated porosity results of Monte Carlo and FFCM was 6.80%,with a rootmean-square error of 0.62 p.u.In field well applications,the formation porosity model obtained using FFCM was in good agreement with the model obtained by interpreters,which demonstrates the validity and accuracy of the proposed method.

Improved nuclear mass formula with an additional term from the Fermi gas model

Nuclear mass is a fundamental property of nuclear physics and a necessary input in nuclear astrophysics.Owing to the complexity of atomic nuclei and nonperturbative strong interactions,conventional physical models cannot completely describe nuclear binding energies.In this study,the mass formula was improved by considering an additional term from the Fermi gas model.All nuclear masses in the Atomic Mass Evaluation Database were reproduced with a root-mean-square deviation(RMSD)of -1.86 MeV(1.92 MeV).The new mass formula exhibits good performance in the neutron-rich nuclear region.The RMSD decreases to 0.393 MeV when the ratio of the neutron number to the proton number is≥1.6.

Reliable calculations of nuclear binding energies by the Gaussian process of machine learning

Reliable calculations of nuclear binding energies are crucial for advancing the research of nuclear physics. Machine learning provides an innovative approach to exploring complex physical problems. In this study, the nuclear binding energies are modeled directly using a machine-learning method called the Gaussian process. First, the binding energies for 2238 nuclei with Z > 20 and N > 20 are calculated using the Gaussian process in a physically motivated feature space, yielding an average deviation of 0.046 MeV and a standard deviation of 0.066 MeV. The results show the good learning ability of the Gaussian process in the studies of binding energies. Then, the predictive power of the Gaussian process is studied by calculating the binding energies for 108 nuclei newly included in AME2020. The theoretical results are in good agreement with the experimental data, reflecting the good predictive power of the Gaussian process. Moreover, the α-decay energies for 1169 nuclei with 50 ≤ Z ≤ 110 are derived from the theoretical binding energies calculated using the Gaussian process. The average deviation and the standard deviation are, respectively, 0.047 MeV and 0.070 MeV. Noticeably, the calculated α-decay energies for the two new isotopes ^ (204 )Ac(Huang et al. Phys Lett B 834, 137484(2022)) and ^ (207) Th(Yang et al. Phys Rev C 105, L051302(2022)) agree well with the latest experimental data. These results demonstrate that the Gaussian process is reliable for the calculations of nuclear binding energies. Finally, the α-decay properties of some unknown actinide nuclei are predicted using the Gaussian process. The predicted results can be useful guides for future research on binding energies and α-decay properties.

Factors Influencing the Spatial Variability of Air Temperature Urban Heat Island Intensity in Chinese Cities

Few studies have investigated the spatial patterns of the air temperature urban heat island(AUHI)and its controlling factors.In this study,the data generated by an urban climate model were used to investigate the spatial variations of the AUHI across China and the underlying climate and ecological drivers.A total of 355 urban clusters were used.We performed an attribution analysis of the AUHI to elucidate the mechanisms underlying its formation.The results show that the midday AUHI is negatively correlated with climate wetness(humid:0.34 K;semi-humid:0.50 K;semi-arid:0.73 K).The annual mean midnight AUHI does not show discernible spatial patterns,but is generally stronger than the midday AUHI.The urban–rural difference in convection efficiency is the largest contributor to the midday AUHI in the humid(0.32±0.09 K)and the semi-arid(0.36±0.11 K)climate zones.The release of anthropogenic heat from urban land is the dominant contributor to the midnight AUHI in all three climate zones.The rural vegetation density is the most important driver of the daytime and nighttime AUHI spatial variations.A spatial covariance analysis revealed that this vegetation influence is manifested mainly through its regulation of heat storage in rural land.

Gyrokinetic simulation of magnetic-island-induced electric potential vortex mode

Ion temperature gradient(ITG)-driven turbulence with embedded static magnetic islands is simulated by utilizing a gyrokinetic theory-based global turbulence transport code(GKNET)in this work.Different from the traditional equilibrium circular magnetic-surface average(EMSA)method,an advanced algorithm that calculates the perturbed magnetic-surface average(PMSA)of the electric potential has been developed to precisely deal with the zonal flow component in a non-circular magnetic surface perturbed by magnetic islands.Simulations show that the electric potential vortex structure inside islands induced by the magnetic islands is usually of odd parity when using the EMSA method.It is found that the odd symmetry vortex can transfer into an even one after a steep zonal flow gradient,i.e.the flow shear has been built in the vicinity of the magnetic islands by adopting the PMSA algorithm.The phase of the potential vortex in the poloidal cross section is coupled with the zonal flow shear.Such an electric potential vortex mode may be of essential importance in wide topics,such as the turbulence spreading across magnetic islands,neoclassical tearing mode physics,and also the interaction dynamics between the micro-turbulence and MHD activities.

Effect of Moringa oleifera on semen quality profiles in endangered Teressa goat under tropical humid island ecosystem

Objective:To investigate the effect of Moringa oleifera leaf extract on semen quality parameters in Teressa goat buck.Methods:A total of 25 semen samples from six bucks were selected for the study.Spermatozoa of 150×10^(6) were incubated in 0,300,500 and 700μg of leaf extract as groupⅠ,Ⅱ,Ⅲ and Ⅳ,respectively.Liquid stored semen samples were analysed for motility,viability,total sperm abnormalities,plasma membrane,acrosomal and nuclear integrities,seminal plasma intracellular enzymes[aspartate aminotransferase(AST),alanine aminotransferase(ALT)and lactate dehydrogenase(LDH)]and total antioxidant capacity(TAC)and spermatozoa malondialdehyde(MDA)and cholesterol efflux in comparison with the control group(groupⅠ)for up to 72 h.Results:Moringa oleifera leaf extract(500μg/150×10^(6) spermatozoa)-treated semen had significantly higher motility,viability,plasma membrane,acrosomal and nuclear integrities,and TAC(P<0.05),and had significantly lower total sperm abnormalities,AST,ALT,LDH,MDA production and cholesterol efflux compared to those in other Moringa oleifera leaf extract treated(300 or 700μg/150×10^(6) spermatozoa)and control groups at different hours of liquid storage(P<0.05).Semen quality parameters and TAC showed an increasing trend and total sperm abnormalities,MDA production,leakage of intra-cellular enzymes and cholesterol efflux showed a decreasing trend from group Ⅰto group Ⅲand then an opposite trend from group Ⅲ to group Ⅳ at different hours of liquid storage.Conclusions:Moringa oleifera leaf extract 500μg/150×10^(6) spermatozoa can be used for liquid semen preservation in Teressa goat at Andaman and Nicobar Islands.

Geochemistry of island arc assemblage in the Eastern Desert of Egypt and the role of Pan-African magmatism in crustal growth of the Arabian–Nubian Shield:A review

Neoproterozoic island arc assemblage of the Arabian–Nubian Shield(ANS)in the Eastern Desert(ED)of Egypt comprises juvenile suites of metavolcanics(MV),large amounts of meta-sedimentary rocks(MS),and voluminous metagabbros-diorites(MGD)and syn-tectonic intrusions of older granitoids(OG).We report here the updates of these four rock units in terms of classification,distribution,chemical characteristics,geodynamic evolution,metamorphism,and ages.In addition,we discuss these integrated data to elucidate a reasonable and reliable model for crustal evolution in the ANS.The main features of these rock units indicate their relation to each other and the geodynamic environment dominated by early immature oceanic island arcs to primitive continental arcs.Integrated information of the island arc metavolcanic and plutonic rocks(gabbros,diorites,tonalites,and granodiorites)furnish evidence of the genetic relationships.These include proximity and a coeval nature in the field;all protolith magmas are subalkaline in nature following calc-alkaline series with minor tholeiitic affinities;common geochemical signature of the arc rocks and subduction-related magmatism;their similar enrichment in LREEs;and similar major element compositions with mafic melts derived from metasomatized mantle wedge.The volcano-sedimentary and the OG rocks underwent multiphase deformation events whereas the MGD complexes deformed slightly.Based on the magmatic,sedimentological,and metamorphic evolutions constrained by geochronological data as well as the progressive evolutionary trend from extensional to compressional regimes,a possible gradual decrease in the subducted slab dip angle is the most infl uential in any geodynamic model for arc assemblage in the ED of Egypt.

Benchmark experiment on slab^(238)U with D-T neutrons for validation of evaluated nuclear data

A benchmark experiment on^(238)U slab samples was conducted using a deuterium-tritium neutron source at the China Institute of Atomic Energy.The leakage neutron spectra within energy levels of 0.8-16 MeV at 60°and 120°were measured using the time-of-flight method.The samples were prepared as rectangular slabs with a 30 cm square base and thicknesses of 3,6,and 9 cm.The leakage neutron spectra were also calculated using the MCNP-4C program based on the latest evaluated files of^(238)U evaluated neutron data from CENDL-3.2,ENDF/B-Ⅷ.0,JENDL-5.0,and JEFF-3.3.Based on the comparison,the deficiencies and improvements in^(238)U evaluated nuclear data were analyzed.The results showed the following.(1)The calculated results for CENDL-3.2 significantly overestimated the measurements in the energy interval of elastic scattering at 60°and 120°.(2)The calculated results of CENDL-3.2 overestimated the measurements in the energy interval of inelastic scattering at 120°.(3)The calculated results for CENDL-3.2 significantly overestimated the measurements in the 3-8.5 MeV energy interval at 60°and 120°.(4)The calculated results with JENDL-5.0 were generally consistent with the measurement results.

Nuclear charge radius predictions by kernel ridge regression with odd-even effects

The extended kernel ridge regression(EKRR)method with odd-even effects was adopted to improve the description of the nuclear charge radius using five commonly used nuclear models.These are:(i)the isospin-dependent A^(1∕3) formula,(ii)relativistic continuum Hartree-Bogoliubov(RCHB)theory,(iii)Hartree-Fock-Bogoliubov(HFB)model HFB25,(iv)the Weizsacker-Skyrme(WS)model WS*,and(v)HFB25*model.In the last two models,the charge radii were calculated using a five-parameter formula with the nuclear shell corrections and deformations obtained from the WS and HFB25 models,respectively.For each model,the resultant root-mean-square deviation for the 1014 nuclei with proton number Z≥8 can be significantly reduced to 0.009-0.013 fm after considering the modification with the EKRR method.The best among them was the RCHB model,with a root-mean-square deviation of 0.0092 fm.The extrapolation abilities of the KRR and EKRR methods for the neutron-rich region were examined,and it was found that after considering the odd-even effects,the extrapolation power was improved compared with that of the original KRR method.The strong odd-even staggering of nuclear charge radii of Ca and Cu isotopes and the abrupt kinks across the neutron N=126 and 82 shell closures were also calculated and could be reproduced quite well by calculations using the EKRR method.