Competition Between Self-birth and Catalyzed Death in Aggregation Growth with Catalysis Injection

We propose two irreversible aggregation growth models of aggregates of two distinct species （.4 and B） to study the interactions between virus aggregates and medicine efficacy aggregates in the virus-medicine cooperative evolution system. The A-species aggregates evolve driven by self monomer birth and B-species aggregate-catalyzed monomer death in model I and by self birth, catalyzed death, and self monomer exchange reactions in model II, while the catalyst B-species aggregates are assumed to be injected into the system sustainedly or at a periodic time-dependent rate. The kinetic behaviors of the A-species aggregates are investigated by the rate equation approach based on the mean-field theory with the self birth rate kernel IA（k） = Ik, catalyzed death rate kernel JAB（k） = Jk and self exchange rate kernel KA （k, l） = Kkl. The kinetic behaviors of the A-species aggregates are mainly dominated by the competition between the two effects of the self birth （with the effective rate I） and the catalyzed death （with the effective rate JB0）, while the effects of the self exchanges of the A-species aggregates which appear in an effective rate KAo play important roles in the cases of I 〉 JBo and I = JBo. The evolution behaviors of the total mass M1^A（t） and the total aggregate number MA（t） are obtained, and the aggregate size distribution ak（t） of species A is found to approach a generalized scaling form in the case of I ≥ JBo and a special modified scaling form in the case of I 〈 JB0. The periodical evolution of the B-monomers concentration plays an exponential form of the periodic modulation.     

Critical behavior of a dynamical percolation model

The critical behavior of the dynamical percolation model,which realizes the molecular-aggregation conception and describes the crossover between the hadronic phase and the partonic phase,is studied in detail. The critical percolation distance for this model is obtained by using the probability P∞ of the appearance of an infinite cluster. Utilizing the finite-size scaling method the critical exponents γ/ν and τ are extracted from the distribution of the average cluster size and cluster number density. The influences of two model related factors,i.e. the maximum bond number and the definition of the infinite cluster,on the critical behavior are found to be small.     

Evolution of nitrogen structure in N-doped diamond crystal after high pressure and high temperature annealing treatment

In this paper, we have reported an investigation on the evolution of nitrogen structures in diamond crystals which contain nitrogen donor atoms in the range of 1500 ppm-1600 ppm following an annealing treatment at a high pressure of about 6.5 GPa and high temperatures of 1920 K-2120 K. The annealing treatment was found to completely transform nitrogen atoms originally arranged in a single substitutional form （C-center）, into a pair form （A-center）, indicated from infrared （IR） spectra. The photoluminescence （PL） spectra revealed that a small fraction of nitrogen atoms remained in C-center form, while some nitrogen atoms in A-center form were further transformed into N3 and H3 center structures. In addition, PL spectra have revealed the existence of two newly observed nitrogen-related structures with zero phonon lines at 611 nm and 711 nm. All these findings above are very helpful in understanding the formation mechanism of natural diamond stones of the Ia-type, which contains nitrogen atoms in an aggregated form.     

Kinetic Behaviors of a Competitive Population and Fitness System in Exchange-Driven Growth

We proposed an aggregation model of two species aggregates of fitness and population to study the interaction between the two species in their exchange-driven processes of the same species by introducing the monomer birth of fitness catalyzed by the population, where the fitness aggregates perform self-death process and the population aggregates perform self-birth process. The kinetic behaviors of the aggregate size distributions of the fitness and population were analyzed by the rate equation approach with their exchange rate kernel K1（k,j） = K1kj and K2（k,j） = K2kj, the fitness aggregate＇s self-death rate kernel J1 （ k ） = J1 k, population aggregate＇s self-birth rate kernel J2（ k ） = J2k and population-catalyzed fitness birth rate kernel I（k,j） = Ikj＇. The kinetic behavior of the fitness was found depending crucially on the parameter v, which reflects the dependence of the population-catalyzed fitness birth rate on the size of the catalyst （population） aggregate. （i） In the v ≤ 0 case, the effect of catalyzed-birth of fitness is rather weak and the exchange-driven aggregation and self-death of the fitness dominate the process, and the fitness aggregate size distribution αk（t） does not have scale form. （ii） When v ≥0, the effect of the population-catalyzed birth of fitness gets strong enough, and the catalyzed-birth and self-death of the fitness aggregates, together with the self-birth of the population aggregates dominate the evolution process of the fitness aggregates. The aggregate size distribution αk （t） approaches a generalized scaling form.     

Aggregation Processes with Catalysis-Driven Decomposition

We propose a three-species aggregation model with catalysis-driven decomposition. Based on the mean-field rate equations, we investigate the evoIution behavior of the system with the size-dependent catalysis-driven decomposition rate J（i; j; k） = Jijk^v and the constant aggregation rates. The results show that the cluster size distribution of the species without decomposition can always obey the conventional scaling law in the case of 0 ≤v ≤ 1, while the kinetic evolution of the decomposed species depends crucially on the index v. Moreover, the total size of the species without decomposition can keep a nonzero value at large times, while the total size of the decomposed species decreases exponentially with time and vanishes finally.     

Experimental studies of THGEM in different Ar/CO2 mixtures

In this paper, the performance of a type of domestic THGEM （THick Gaseous Electron Multiplier） working in Ar/CO2 mixtures is reported in detail. This kind of single THGEM can provide a gain range from 100 to 1000, which is very suitable for application in neutron detection. In order to study its basic characteristics as a reference for the development of a THGEM based neutron detector, the counting rate plateau, the energy resolution and the gain of the THGEM have been measured in different Ar/CO2 mixtures with a variety of electrical fields. For the Ar/CO2（90%/10%） gas mixture, a wide counting rate plateau is achieved from 720 V to 770 V, with a plateau slope of 2.4%/100 V, and an excellent energy resolution of about 22% is obtained at the 5.9 keV full energy peak of the 55Fe X-ray source.     

Decoherence,Orbit, and Classical Linear Frequency Entropy： Quantum-Classical Correspondence

The decoherence process is analyzed for an open quantum system that is classically chaotic, with a classical linear frequency entropy developed to measure the stability of classical motion. Investigation shows that the decoherence measured by the rate of quantum linear entropy production varies significantly with both the underlying classical orbits and the classical linear frequency entropy. Such correspondence is also supported by the further investigation on the Loschmidt Echo.     

Micro-track structure analysis for 100 MeV Si ions in CR-39 by using atomic force microscopy

To analyze the micro-track structure of heavy ions in a polymer material,parameters including bulk etch rate,track etch rate,etch rate ratio,and track core size were measured.The pieces of CR-39 were exposed to 100 MeV Si ions with normal incidence and were etched in 6.25N NaOH solution at 70 C.Bulk etch rate was read out by a profilemeter after several hours of etching.The other parameters were obtained by using an atomic force microscope(AFM)after a short time of etching.We have measured the second etch pits and minute etch pits to obtain the track growth curve and three dimension track structures to track the core size and etch rate measurements.The local dose of the track core was calculated by theδ-ray theory.In our study,we figure out that the bulk etch rate Vb=(1.58±0.022)μm/h,the track etch rate Vt=(2.90±0.529)μm/h,the etch rate ratio V=1.84±0.031,and the track core radii r≈4.65 nm.In the meantime,we find that the micro-track development violates the traditional track-growth model.For this reason,a scenario is carried out to provide an explanation.     

A prototype of a high rating MRPC

Six-gap resistive plate chamber （MRPC） prototypes with semiconductive glass electrodes （bulk resistivity ～ 10^10Ω·cm） were studied for suitability in time-of-flight （TOF） applications at high rates. These studies were performed using a continuous electron beam of 800 MeV at IHEP and an X-ray machine. Time resolutions of about 100 ps and efficiencies larger than 90% were obtained for flux densities up to 28 kHz/cm^2.     

Study on collimation and shielding of the back-streaming neutrons at the CSNS target

The back-streaming neutrons from the spallation target at CSNS are very intense, and can pose serious damage problems for the devices in the accelerator-target interface region. To tackle the problems, a possible scheme for this region was studied, namely a specially designed optics for the proton beam line produces two beam waists, and two collimators are placed at the two waist positions to maximize the collimation effect of the back-streaming neutrons. Detailed Monte Carlo simulations with the beams in the two different CSNS phases show the effectiveness of the collimation system, and the radiation dose rate decreases largely in the interface section. This can ensure the use of epoxy coils for the last magnets and other devices in the beam transport line with reasonable lifetimes, e.g., thirty years. The design philosophy for such an accelerator-target interface region can also be applicable to other high-power proton beam applications.     

Effects of ultra-strong magnetic field on electron capture rates of 55C0 and 56Ni in the magnetar surrounding

The electron capture rates of 55Co and 56Ni in the ultra-strong magnetic field at four typical temperature- density points have been calculated using the nuclear shell model and Landan energy levels quantized approximate correction. The results show that the electron capture rates of 55Co and 56Ni are increased greatly in the ultra-strong magnetic field, and even exceed two orders of magnitude in the range from 4.414×10^13G to 2.207×10^27 G. The change rate of electron abundance, ye, of 55Co and 56Ni under the condition of B=4.414×10^15G in the magnetar surrounding has been calculated and discussed, the proportions of ye of 55Co and 56Ni in the total Ye have been reduced by 50 percent in all more than the condition without a magnetic field.     

Measurement and analysis of fission rates in a spherical mockup of uranium and polyethylene

Measurements of the reaction rate distribution were carried out using two kinds of Plate Micro Fission Chamber （PMFC）. The first is a depleted uranium chamber and the second an enriched uranium chamber. The material in the depleted uranium chamber is strictly the same as the material in the uranium assembly. With the equation solution to conduct the isotope contribution correction, the fission rate of 238U and 235U were obtained from the fission rate of depleted uranium and enriched uranium. Then, the fission count of 238U and 235U in an individual uranium shell was obtained. In this work, MCNP5 and continuous energy cross sections ENDF/BV.0 were used for the analysis of fission rate distribution and fission count. The calculated results were compared with the experimental ones. The calculation of fission rate of DU and EU were found to agree with the measured ones within 10% except at the positions in polyethylene region and the two positions near the outer surface. Because the fission chamber was not considered in the calculation of the fission counts of 238U and 235U, the calculated results did not agree well with the experimental ones.     

Higher-order differential variational principle and differential equations of motion for mechanical systems in event space

In this paper we study the higher-order differential variational principle and differential equations of motion for mechanical systems in event space. Based on the higher-order d＇Alembert principle of the system, the higher-order velocity energy and the higher-order acceleration energy of the system in event space are defined, the higher-order d＇Alembert- Lagrange principle of the system in event space is established, and the parametric forms of Euler-Lagrange, Nielsen and Appell for this principle are given. Finally, the higher-order differential equations of motion for holonomic systems in event space are obtained.     

Influences of nitrogen flow rate on the structures and properties of Ti and N co-doped diamond-like carbon films deposited by arc ion plating

In this paper, Ti-C-N nanocomposite films are deposited under different nitrogen flow rates by pulsed bias arc ion plating using Ti and graphite targets in the Ar/N2 mixture gas. The surface morphologies, compositions, microstructures, and mechanical properties of the Ti-C-N films are investigated systematically by field emission scanning electron mi- croscopy （FE-SEM）, x-ray photoelectron spectroscopy （XPS）, grazing incident x-ray diffraction （GIXRD）, Raman spectra, and nano-indentation. The results show that the nanocrystalline Ti（C,N） phase precipitates in the film from GIXRD and XPS analysis, and Raman spectra prove the presence of diamond-like carbon, indicating the formation of nanocomposite film with microstructures comprising nanocrystalline Ti（C,N） phase embedded into a diamond-like matrix. The nitrogen flow rate has a significant effect on the composition, structure, and properties of the film. The nano-hardness and elastic modulus first increase and then decrease as nitrogen flow rate increases, reaching a maximum of 34.3 GPa and 383.2 GPa, at a nitrogen flow rate of 90 sccm, respectively.     

Monte Carlo simulation of indoor external exposure due to gamma-emitting radionuclides in building materials

The use of building materials containing naturally occurring radionuclides such as^40K,^238U,^232Th and their progeny, could lead to external exposures to the residents of such buildings. In this paper, a set of models are constructed to calculate the specific effective dose rates（the effective dose rate per Bq/kg of ^40K, the ^238U series,and the ^232Th series） imposed on residents by building materials with the MCNPX code. The effect of chemical composition, position concerned in the room and thickness as well as density of material is analyzed. In order to facilitate more precise assessment of indoor external dose due to gamma-emitting radionuclides in building materials,three regressive expressions are proposed and validated by measured data to calculate specific effective rates for40 K,the238U series and the232 Th series, respectively.     

Extinction Effects of Multiplicative Non-Gaussian L′evy Noise in a Tumor Growth System with Immunization

The extinction phenomenon induced by multiplicative non-Gaussian Levy noise in a tumor growth model with immune response is discussed. Under the influence of the stochastic immune rate, the model is analyzed in terms of a stochastic differential equation with multiplicative noise. By means of the theory of the infinitesimal generator of Hunt processes, the escape probability, which is used to measure the noise-induced extinction probability of tumor cells, is explicitly expressed as a function of initial tumor cell density, stability index and noise intensity. Based on the numerical calculations, it is found that for different initial densities of tumor cells, noise parameters play opposite roles on the escape probability. The optimally selected values of the multiplicative noise intensity and the stability index are found to maximize the escape probability.     

Influence of thyroid volume reduction on absorbed dose in 131I therapy studied by using Geant4 Monte Carlo simulation

A simulation study has been performed to quantify the effect of volume reduction on the thyroid absorbed dose per decay and to investigate the variation of energy deposition per decay due to β- and γ-activity of 131I with volume/mass of thyroid, for water, ICRP- and ICRU-soft tissue taken as thyroid material. A Monte Carlo model of the thyroid, in the Geant4 radiation transport simulation toolkit was constructed to compute the β- and γ-absorbed dose in the simulated thyroid phantom for various values of its volume. The effect of the size and shape of the thyroid on energy deposition per decay has also been studied by using spherical, ellipsoidal and cylindrical models for the thyroid and varying its volume in 1-25 cm3 range. The relative differences of Geant4 results for different models with each other and MCNP results lie well below 1.870%. The maximum relative difference among the Geant4 estimated results for water with ICRP and ICRU soft tissues is not more than 0.225%. S-values for ellipsoidal, spherical and cylindrical thyroid models were estimated and the relative difference with published results lies within 3.095%. The absorbed fraction values for beta particles show a good agreement with published values within 2.105% deviation. The Geant4 based simulation results of absorbed fractions for gammas again show a good agreement with the corresponding MCNP and EGS4 results （±6.667%） but have 29.032% higher values than that of MIRD calculated values. Consistent with previous studies, the reduction of the thyroid volume is found to have a substantial effect on the absorbed dose. Geant4 simulations confirm dose dependence on the volume/mass of thyroid in agreement with MCNP and EGS4 computed values but are substantially different from MIRD8 data. Therefore, inclusion of size/mass dependence is indicated for 131I radiotherapy of the thyroid.     

First-Principles Calculations of Atomic and Electronic Properties of T1 and In on Si（111）

The atomic and electronic structures of T1 and In on Si（111） surfaces are investigated using the firstprinciples total energy calculations. Total energy optimizations show that the energetically favored structure is 1/3 ML T1 adsorbed at the T4 sites on Si（111） surfaces. The adsorption energy difference of one T1 adatom between （√3 × √3） and （1 × 1） is less than that of each In adatom. The DOS indicates that TI 6p and Si 3p electrons play a very important role in the formation of the surface states. It is concluded that the bonding of TI adatoms on Si（111） surfaces is mainly polar covalent, which is weaker than that of In on Si（111）. So T1 atom is more easy to be migrated than In atom in the same external electric field and the structures of T1 on Si（111） is prone to switch between （√3 × √3） and （1 × 1）.     

Study on the pre-chopper in CSNS LEBT

Physical designing of the pre-chopper in CSNS LEBT is carried out, which includes the deflecting voltage, the length and the width of the deflecting plates, and the gap between the deflecting plates. The most outstanding feature of the design is that both the gap and the width vary with the beam envelope size. So both the requried deflecting voltage and the loaded capacitance are lowered. In order to avoid destruction of the space charge neutralization by the pre-chopper in the whole LEBT, an electron-trapping electrode is arranged to confine the electrostatic field of the pre-chopper to the local area. To examine the reliability of the pre-chopping design in CSNS LEBT, a similar pre-chopping design in ADS RFQ LEBT is set up and an experiment on the pre-chopper is prepared. 3-dimensional simulations are carried out to determine the loaded capacitance and the applied voltage of the electron-trapping electrode.     

Synchronization in a Mutualism Ecosystem Induced by Noise Correlation

Understanding the cause of the synchronization of population evolution is an important issue for ecological improvement. Here we present a Lotka-Volterra-type model driven by two correlated environmental noises and show, via theoretical analysis and direct simulation, that noise correlation can induce a synchronization of the mutualists. The time series of mutual species exhibit a chaotic-like fluctuation, which is independent of the noise correlation, however, the chaotic fluctuation of mutual species ratio decreases with the noise correlation. A quantitative parameter defined for characterizing chaotic fluctuation provides a good approach to measure when the complete synchronization happens.