G.G. Marmitt , I. Alencar , H. Trombini , F.F. Selau , B. Konrad , P.L. Grande
{"title":"PowerMEIS 3: A versatile tool for simulating ion and electron scattering","authors":"G.G. Marmitt , I. Alencar , H. Trombini , F.F. Selau , B. Konrad , P.L. Grande","doi":"10.1016/j.cpc.2025.109639","DOIUrl":"10.1016/j.cpc.2025.109639","url":null,"abstract":"<div><div>The aggressive roadmap for nanotechnology is driving the development of characterization techniques capable of providing nanometric resolution while preserving structural and chemical information for increasingly complex samples. Ion and electron scattering have emerged as powerful methodologies to meet these demands. However, due to the sophistication of modern samples, data interpretation heavily relies on advanced simulations. In this context, we have developed the <span>PowerMEIS<!--> <!-->3</span> computer program, a versatile Monte Carlo tool for simulating the scattering spectra of ions and electrons. This program has been rewritten from its previous versions and incorporates several new features. A detailed description of its implementation is provided after introducing the necessary physical principles. The program's wide range of applications is illustrated through several examples, including Medium Energy Ion Scattering, Rutherford Backscattering Spectrometry, molecular ion scattering, Nuclear Reaction Profiling, and Reflection Electron Energy Loss Spectroscopy. In particular, we demonstrate three distinct strategies for calculating the path integral: the Single Scattering, <em>Connected Trajectory</em>, and <em>Direct Trajectory</em> approaches, all based on a voxel representation of the target sample. Additionally, we compare the performance of <span>PowerMEIS<!--> <!-->3</span> program to other established programs, such as <span>TRBS</span> and <span>SIMNRA</span>. The <em>Connected Trajectory</em> approach is a novel feature in scattering simulations and significantly reduces the simulation time for Multiple Scattering calculations. Moreover, it enables simulations of nanostructures at any incidence angle, a capability not supported by other programs. The program also offers the option to run simulations remotely on a server hosted at <em>Universidade Federal do Rio Grande do Sul</em> (UFRGS). Finally, we discuss the limitations of the <em>Connected Trajectory</em> approach at lower energies, primarily due to the time reversal approximation employed, and highlight possibilities for further development.</div></div>","PeriodicalId":285,"journal":{"name":"Computer Physics Communications","volume":"313 ","pages":"Article 109639"},"PeriodicalIF":7.2,"publicationDate":"2025-04-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143904582","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Multi-expression programming for enhancing MHD heat transfer in a nanofluid-filled enclosure with heat generation and viscous dissipation","authors":"Naeem Ullah , Aneela Bibi , Dianchen Lu","doi":"10.1016/j.cpc.2025.109649","DOIUrl":"10.1016/j.cpc.2025.109649","url":null,"abstract":"<div><div>Efficient thermal management is a critical challenge in various engineering configuration where overheating affects performance, such as electronics, industrial cooling, and HVAC applications. Traditional cooling methods often struggle with confined enclosures, leading to inefficiencies. Nanofluids and optimized heating mechanisms offer a promising solution, but their complex thermal behavior requires precise predictive modeling. This study addresses this challenge by conducting a numerical analysis of heat transfer in nanofluid-filled enclosures with sinusoidal heating. This study employs multi-expression programming technique to improve thermal performance by analyzing heating design and electromagnetic interactions. In this exploration a square enclosure filled with water-based copper oxide nanofluid is evaluated, featuring a centrally located sinusoidal heated element. The enclosure is also partially heated from below, cooled along the sidewalls, while the upper and remaining lower portions are insulated. The numerical simulation explores flow-controlling variables, including nanoparticles volume fraction, heating element amplitude, magnetic field strength and its orientation, viscous dissipation, and heat generation, to assess their impact on flow dynamics and thermal performance. The findings indicate that the Nusselt number increases by <span><math><mrow><mn>26.68</mn><mo>%</mo></mrow></math></span> when nanoparticle concentration reaches <span><math><mrow><mn>4</mn><mo>%</mo></mrow></math></span>, while a rise in Rayleigh number from <span><math><msup><mrow><mn>10</mn></mrow><mn>3</mn></msup></math></span> to <span><math><msup><mrow><mn>10</mn></mrow><mn>6</mn></msup></math></span> results in an approximate <span><math><mrow><mn>75.40</mn><mo>%</mo></mrow></math></span> increase. Moreover, the average percentage decrease in Nusselt number against <span><math><msub><mi>Q</mi><mi>g</mi></msub></math></span> from 0 to 30 is 20.71% while for <span><math><mrow><mi>H</mi><mi>a</mi></mrow></math></span> (10 to 100) it is 42.61%.The multi-expression programming model accurately predicts convective heat transfer trends, achieving a high correlation coefficient (<span><math><mrow><msub><mi>C</mi><mi>R</mi></msub><mo>=</mo><mn>0.99</mn></mrow></math></span> for training, <span><math><mrow><msub><mi>C</mi><mi>R</mi></msub><mo>=</mo><mn>0.94</mn></mrow></math></span> for testing) and low error metrics (<span><math><mrow><mi>R</mi><mi>M</mi><mi>S</mi><mi>E</mi><mspace></mspace><mo>=</mo><mspace></mspace><mn>0.02</mn><mo>,</mo><mspace></mspace><mi>M</mi><mi>A</mi><mi>E</mi><mspace></mspace><mo>=</mo><mspace></mspace><mn>0.03</mn><mo>,</mo><mspace></mspace><mi>P</mi><mi>I</mi><mspace></mspace><mo>=</mo><mspace></mspace><mn>0.06</mn></mrow></math></span> for training), ensuring strong agreement with numerical results.</div></div>","PeriodicalId":285,"journal":{"name":"Computer Physics Communications","volume":"313 ","pages":"Article 109649"},"PeriodicalIF":7.2,"publicationDate":"2025-04-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143908098","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Minimal autocorrelation in hybrid Monte Carlo simulations using exact Fourier acceleration","authors":"Johann Ostmeyer , Pavel Buividovich","doi":"10.1016/j.cpc.2025.109624","DOIUrl":"10.1016/j.cpc.2025.109624","url":null,"abstract":"<div><div>The hybrid Monte Carlo (HMC) algorithm is a ubiquitous method in computational physics with applications ranging from condensed matter to lattice QCD and beyond. However, HMC simulations often suffer from long autocorrelation times, severely reducing their efficiency. In this work two of the main sources of autocorrelations are identified and eliminated. The first source is the sampling of the canonical momenta from a sub-optimal normal distribution, the second is a badly chosen trajectory length. Analytic solutions to both problems are presented and implemented in the exact Fourier acceleration (EFA) method. It completely removes autocorrelations for near-harmonic potentials and consistently yields (close-to-) optimal results for numerical simulations of the Su-Schrieffer-Heeger and the Ising models as well as in lattice gauge theory, in some cases reducing the autocorrelation by multiple orders of magnitude. EFA is advantageous for and easily applicable to any HMC simulation of an action that includes a quadratic part.</div></div>","PeriodicalId":285,"journal":{"name":"Computer Physics Communications","volume":"313 ","pages":"Article 109624"},"PeriodicalIF":7.2,"publicationDate":"2025-04-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143868630","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Elizabeth Suehr , Manuel Gale , Ramon Lopez , Raymond L. Fontenot , Peter Liever , Jennifer S Curtis
{"title":"Multi-sphere rigid-body particles in a parallelized LEBC with LIGGGHTS","authors":"Elizabeth Suehr , Manuel Gale , Ramon Lopez , Raymond L. Fontenot , Peter Liever , Jennifer S Curtis","doi":"10.1016/j.cpc.2025.109636","DOIUrl":"10.1016/j.cpc.2025.109636","url":null,"abstract":"<div><div>A method for the Message Passing Interface (MPI) parallelization of the Lees-Edwards boundary condition (LEBC) within the LIGGGHTS framework for multi-sphere rigid particles was created, allowing for the simulation of very detailed complex shapes. Double-send and double-receive communication was added to LIGGGHTS to allow for shared information across disjointed processor domains along the shearing boundary of the LEBC. The verification of this method is performed via 3D shearing simulations of single spheres and sphere clumps and rods with aspect ratios 2, 4, and 6. The predicted shear stress employing the new parallelized LEBC method matches stress values from granular kinetic theory and previously published simulation results. No LEBC simulations for DEM or multi-sphere rigid particles are known to be parallelized, allowing for computationally difficult LEBC multi-sphere simulations to be performed for the first time.</div></div>","PeriodicalId":285,"journal":{"name":"Computer Physics Communications","volume":"313 ","pages":"Article 109636"},"PeriodicalIF":7.2,"publicationDate":"2025-04-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143908096","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Julio Almansa , Francesc Salvat-Pujol , Gloria Díaz-Londoño , Artur Carnicer , Antonio M. Lallena , Francesc Salvat
{"title":"PENGEOM – A general-purpose geometry package for Monte Carlo simulation of radiation transport in complex material structures (New Version Announcement)","authors":"Julio Almansa , Francesc Salvat-Pujol , Gloria Díaz-Londoño , Artur Carnicer , Antonio M. Lallena , Francesc Salvat","doi":"10.1016/j.cpc.2025.109634","DOIUrl":"10.1016/j.cpc.2025.109634","url":null,"abstract":"<div><div>A new version of the code system <span>pengeom</span>, which provides a complete set of tools to handle different geometries in Monte Carlo simulations of radiation transport, is presented. The distribution package consists of a set of Fortran subroutines and a Java graphical user interface that allows building and debugging the geometry-definition file, and producing images of the geometry in two- and three dimensions. A detailed description of these tools is given in the original paper [<em>Comput. Phys. Commun.</em> <strong>199</strong> (2016) 102–113] and in the code manual included in the distribution package. The present new version differs from the previous one in that 1) it implements a more systematic handling of round-off errors, 2) the set of examples has been updated, and 3) it allows including a single voxelized box as a geometry module. With the last optional feature, a Monte Carlo code can readily be used for describing irradiation processes with complex material structures, such as medical treatments.</div></div>","PeriodicalId":285,"journal":{"name":"Computer Physics Communications","volume":"313 ","pages":"Article 109634"},"PeriodicalIF":7.2,"publicationDate":"2025-04-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143868631","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Grover's search meets Ising models: A quantum algorithm for finding low-energy states","authors":"A.A. Zhukov , A.V. Lebedev , W.V. Pogosov","doi":"10.1016/j.cpc.2025.109627","DOIUrl":"10.1016/j.cpc.2025.109627","url":null,"abstract":"<div><div>We propose a methodology for implementing Grover's algorithm in the digital quantum simulation of disordered Ising models. The core concept revolves around using the evolution operator for the Ising model as the quantum oracle within Grover's search. This operator induces phase shifts for the eigenstates of the Ising Hamiltonian, with the most pronounced shifts occurring for the lowest and highest energy states. Determining these states for a disordered Ising Hamiltonian using classical methods presents an exponentially complex challenge with respect to the number of spins (or qubits) involved. Within our proposed approach, we determine the optimal evolution time by ensuring a phase flip for the target states. This method yields a quadratic speedup compared to classical computation methods and enables the identification of the lowest and highest energy states (or neighboring states) with a high probability ≲1.</div></div>","PeriodicalId":285,"journal":{"name":"Computer Physics Communications","volume":"313 ","pages":"Article 109627"},"PeriodicalIF":7.2,"publicationDate":"2025-04-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143868629","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Dennis Herb , Marco Trenti , Marilena Mantela , Constantinos Simserides , Joachim Ankerhold , Mirko Rossini
{"title":"QuantumDNA: A python package for analyzing quantum charge dynamics in DNA and exploring its biological relevance","authors":"Dennis Herb , Marco Trenti , Marilena Mantela , Constantinos Simserides , Joachim Ankerhold , Mirko Rossini","doi":"10.1016/j.cpc.2025.109626","DOIUrl":"10.1016/j.cpc.2025.109626","url":null,"abstract":"<div><div>The study of DNA charge dynamics is a highly interdisciplinary field that bridges physics, chemistry, biology, and medicine, and plays a critical role in processes such as DNA damage detection, protein-DNA interactions, and DNA-based nanotechnology. However, despite significant progress in each of these areas, knowledge often remains inaccessible to researchers in other scientific communities, limiting the broader impact of advances across disciplines. To bridge this gap, we present QuantumDNA, an open-source Python package for simulating DNA charge transfer and excited state dynamics using quantum physical methods. QuantumDNA combines an efficient Linear Combination of Atomic Orbitals (LCAO) approach combined with tight-binding models and incorporates open quantum systems techniques to account for environmental effects. This approach allows for a rapid yet sufficiently accurate analysis of large DNA ensembles, enabling statistical studies of genetic and epigenetic phenomena. To ensure accessibility, the package features a graphical user interface, making it suitable for researchers across disciplines.</div></div><div><h3>Program summary</h3><div><em>Program Title:</em> QuantumDNA</div><div><em>CPC Library link to program files:</em> <span><span>https://doi.org/10.17632/5mw48c7gbb.1</span><svg><path></path></svg></span></div><div><em>Developer's repository link:</em> <span><span>https://github.com/dehe1011/QuantumDNA</span><svg><path></path></svg></span></div><div><em>Licensing provisions:</em> BSD 3-clause</div><div><em>Programming language:</em> Python</div><div><em>Nature of the Problem:</em> Over the past 60 years, a variety of advanced simulation methods have been employed to explore charge dynamics of DNA. However, most of these approaches are computationally too expensive for the large-scale statistical screening required in genetics and epigenetics. Therefore, theoretical and computational results are often restricted to specialized fields, limiting their accessibility and reproducibility to researchers across disciplines. <em>Solution Method:</em> QuantumDNA combines computational methods from quantum physics and theoretical chemistry to facilitate high-throughput analysis of DNA charge dynamics with sufficient accuracy. Unlike computationally expensive <em>ab initio</em> methods, it utilizes the efficiency of LCAO and TB models to simulate charge dynamics while considering environmental effects, making simulations more accessible and encouraging interdisciplinary research. <em>Additional comments:</em> QuantumDNA is an open-source package featuring a graphical user interface, tutorial Jupyter notebooks, and a dedicated documentation website. The package also supports CPU parallel computing.</div></div>","PeriodicalId":285,"journal":{"name":"Computer Physics Communications","volume":"313 ","pages":"Article 109626"},"PeriodicalIF":7.2,"publicationDate":"2025-04-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143886151","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Bloch-state optimal basis sets: An efficient approach for electronic structure interpolation","authors":"Sasawat Jamnuch , John Vinson","doi":"10.1016/j.cpc.2025.109635","DOIUrl":"10.1016/j.cpc.2025.109635","url":null,"abstract":"<div><div>We present an efficient implementation of the <em>k</em>-space interpolation for electronic structure based on the optimal basis method originally proposed by Shirley [Phys. Rev. B <strong>54</strong>, 16,464 (1996)] The method allows interpolation onto any <em>k</em>-point from a minimal set of input density functional theory (DFT) wavefunctions. Numerically interpolated eigenvalues have an accuracy within 0.01 eV with very small computational cost. The interpolated wavefunctions were used in the Bethe-Salpeter equation to simulate x-ray absorption spectra of different systems, ranging from small bulk crystals to large intermetallic supercells. The method is extensively tested in terms of verification and accuracy for best practices. The approach is shown to be robust and will greatly help accelerate high-throughput DFT studies.</div></div>","PeriodicalId":285,"journal":{"name":"Computer Physics Communications","volume":"313 ","pages":"Article 109635"},"PeriodicalIF":7.2,"publicationDate":"2025-04-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143863532","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
MohammadAli Daeian , Punya Cheema , W. Spencer Smith , Zahra Keshavarz-Motamed
{"title":"A multi-domain lattice Boltzmann method for mesh refinement with curved boundary interfaces","authors":"MohammadAli Daeian , Punya Cheema , W. Spencer Smith , Zahra Keshavarz-Motamed","doi":"10.1016/j.cpc.2025.109637","DOIUrl":"10.1016/j.cpc.2025.109637","url":null,"abstract":"<div><div>Multi-domain grid refinement is a well-known method for mesh refinement in Lattice Boltzmann Methods (LBM). However, the method in three-dimensional cases is currently limited to problems in which the interface between domains can only be surfaces with straight boundaries, and no 3D multi-domain LBM method is specifically tailored for cases with domain interface on a complex curved boundary. Complex geometries like this are frequently observed in blood flow in cardiovascular systems. In this paper, an LBM multi-domain method was developed for grid refinement with curved boundary interfaces. The proposed method is based on using an interpolative second-order wall boundary condition in conjunction with a new image-based ghost node method for near-wall treatment at the interface. The method was verified to show second-order accuracy in space at the domain interface in a circular Poiseuille flow. The methodology was further employed in three different cases: steady idealized stenosis flow, pulsatile flow in the carotid bifurcation, and pulsatile flow in an intracranial aneurysm. The results were compared to single-resolution simulations for each case. For the resolutions used in these cases, the relative L<sup>2</sup> norm of the difference between the multi-domain and fine single-resolution simulation had the following values: 0.005 for velocity magnitude at the stenosis center, 0.002 for mass flowrate splitting in the bifurcation, and 0.008 for wall shear stress in peak systole in the aneurysm dome. For these examples, the method demonstrated up to 65 % speedup for the bifurcation simulation and 39 % speedup for the aneurysm simulation compared to single-resolution simulations.</div></div>","PeriodicalId":285,"journal":{"name":"Computer Physics Communications","volume":"313 ","pages":"Article 109637"},"PeriodicalIF":7.2,"publicationDate":"2025-04-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143894785","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Hao-Yan Liu , Guang-Yu Sun , Yue-Lin Liu , Shu Zhang , Chang-Chun Qi , Sheng Zhou , Wen-Rui Li , Guan-Jun Zhang
{"title":"FlashSim: A novel particle-in-cell numerical model for vacuum surface flashover simulation based on finite element method","authors":"Hao-Yan Liu , Guang-Yu Sun , Yue-Lin Liu , Shu Zhang , Chang-Chun Qi , Sheng Zhou , Wen-Rui Li , Guan-Jun Zhang","doi":"10.1016/j.cpc.2025.109625","DOIUrl":"10.1016/j.cpc.2025.109625","url":null,"abstract":"<div><div>In vacuum-dielectric insulation systems, surface flashover is commonly considered as an interfacial breakdown induced by high electric field, posing a significant threat to the stable and safe operation of numerous vacuum equipment. To clarify its development mechanism and propose relevant suppression strategies, a novel Particle-in-Cell (PIC) vacuum surface flashover simulation model, FlashSim, based on finite element method (FEM), is introduced. The model computes the real-time electric field at cathode triple junction (CTJ) where seed electrons are generated. Electron collisions with the dielectric surface, including elastic backscattering, inelastic backscattering, and true secondary electron emission, are considered. Firstly, the flashover simulation with flat dielectric surface and parallel-plate electrodes shows that electric field at CTJ is enhanced due to surface charge accumulation, leading to an increase in field electron emission (FEE). Low-energy (<53 eV) electrons concentrate near the dielectric surface due to positive charging in the dielectric layer, while high-energy (53eV-1.9 keV) electrons can escape into higher vertical positions. Furthermore, the model employs adaptive mesh, enabling high simulation accuracy without compromising the computational time. Notably, the adopted FEM mesh generator can handle complex boundary geometries, which is validated by simulations with surface grooving for promoting the flashover capability. The performance of grooves is evaluated through electron distribution, anode current density, average surface charge density and electron flux, demonstrating higher flashover strength. The proposed FlashSim simulation model is expected to assist in further revealing the flashover mechanism and developing novel strategies for flashover suppression, and will be further upgraded to include outgassing and plasma discharge.</div></div>","PeriodicalId":285,"journal":{"name":"Computer Physics Communications","volume":"313 ","pages":"Article 109625"},"PeriodicalIF":7.2,"publicationDate":"2025-04-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143868632","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}