World Journal of Mechanics
Vol.4 No.6(2014), Article ID:47455,27 pages DOI:10.4236/wjm.2014.46019

Asymptotic Models for Studying Kinetics of Formation of Compact Objects with Strong Internal Bonds

Emil Edipovich Lin

Russian Federal Nuclear Center—All-Russia Research Institute of Experimental Physics, Sarov, Russia

Email: postmaster@ifv.vniief.ru, e.e.lin@ifv.vniief.ru

Copyright © 2014 by author and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY).

http://creativecommons.org/licenses/by/4.0/

Received 30 March 2014; revised 24 April 2014; accepted 23 May 2014

ABSTRACT

An asymptotic method has been developed for investigation of kinetics of formation of compact objects with strong internal bonds. The method is based on the uncertainty relation for a coordinate and a momentum in space of sizes of objects (clusters) with strongly pronounced collective quantum properties resulted from exchange interactions of various physical nature determined by spatial scales of the processes under consideration. The proposed phenomenological approach has been developed by analogy with the all-known ideas about coherent states of quantum mechanical oscillator systems for which a product of coordinate and momentum uncertainties (dispersions) accepts the value, which is minimally possible within uncertainty relations. With such an approach the leading processes are oscillations of components that make up objects, mainly: collective nucleon oscillations in a nucleus and phonon excitations in a mesostructure crystal lattice. This allows us to consider formation and growth of subatomic and mesoscopic objects in the context of a single formalism. The proposed models adequately describe characteristics of formation processes of nuclear matter clusters as well as mesoscopic crystals having covalent and quasi-covalent bonds between atoms.

Keywords:Compact Objects, Formation Kinetics, Clusters, Quantum Properties, Spatial Scales, Nuclear Matter, Mesoscopic Crystals

1. Introduction

In works [1] -[6] various physical objects and their environment are described. They are of interest not only in terms of their individual properties determined by their structure and sizes but also in terms of determining a spatial boundary between macroand microcosm (macrophysics and microphysics). Object sizes “intermediate” between atoms and solid bodies in mesoscopic crystal structures (in which behaviors conform to the rules of both quantum mechanics and classical physics) can reach considerable values, at which the lower complexity limit of a macroscopic crystal is reached. Fundamental considerations and assessments of quantum representation feasibility show that the “coordinate-momentum” uncertainty relation can fit for bodies with the sizes about 10–7 m, i.e., at a level of nanometer scales [2] . Currently more and more papers deal with investigations of formation, evolution and properties of nanostructures and nanostructure-based materials. These investigations cover processes of various physical and chemical nature proceeding in a mesoscopic region of spatial scales up to 10–10 - 10–4 m: from formation of atomic clusters on a molecular level and nanocrystal synthesis to global changes of a macroscopic continuum state [7] -[13] . Mechanical behavior of nanostructured objects is described using molecular dynamics methods and using statistic models allowing for both thermodynamic and kinetic aspects of problems under consideration. These problems cover incompleteness of classical description of process running in mesoscopic scales, necessity of developing new models including phenomenological ones that take into account formation of quantum systems and determine growth regularities and nanostructure consolidation, and determination of phase transformation size dependencies, in particular, ascertainment of the effect of polymorphic transformations on synergy of crystal structures.

In this connection it is important in principle to determine general features of dynamics of object growth both in microcosm and in “mesocosm”. With this purpose this paper investigates kinetics of formation of various in physical nature compact objects having strong internal bonds and considerably different spatial scales: nuclear matter clusters [14] and nanostructured mesoscopic structures with covalent and quasi-covalent bonds between atoms in a crystal lattice [15] . A typical feature of these objects is the existence of purely quantum effect of exchange interaction. Virtual meson exchange occurs between nucleons; electron exchange occurs in crystal structures having covalent bonds. The effect of exchange interaction allows these objects to be viewed as compact clusters with strongly pronounced collective quantum properties (“quantum” clusters). In case of nuclei these quantum properties are connected with strong interactions and manifest themselves in the existence of oscillatory and rotator shells. In mesocrystals quasi-particles (phonons) are excited.

The above creates the necessary precondition for investigating kinetics of object formation from a common point of view. We will consider closed stochastic systems of compact quantum clusters randomly interacting with each other. In accordance with standard ideas of physical kinetics the irreversible aggregation of objects is described using the concept of distribution density wave in the space of cluster sizes a, propagating in time t in direction of cluster size increase. It is noted in Ref. [16] that the behavior of a wave packet of any physical nature is subject to a universal ratio (resulting from Fourier theorem) for half-width of a wave packet and half-width of a spectral line ∆a∙∆k ≥ 1∕4π, where k is a wave number. This results in the uncertainty relation for coordinate and momentum in the space of cluster sizes. For quantum-mechanical systems being in a coherent state (harmonic oscillations), the above product takes the minimum value, the uncertainty relation is fulfilled as an accurate equality, and a process trajectory in a phase space is classical [17] . This means that it is possible to consider relatively large objects, the sizes of which exceed values corresponding to the lower limit of complexity for a macroscopic solid. Thus, the proposed phenomenological approach describes a limiting asymptotic stage of cluster growth process. Under this approach “the leading” processes appear to be oscillations of components that compose objects, for example collective oscillations of nucleons in a nucleus or phonon excitations of mesostructure crystal lattice.

It should be noted that the proposed asymptotic method in no way can replace the developed fundamental methods of investigations and physical models of microcosm and mesocosm dynamics described in the above Refs. [1] -[10] and represented in [18] -[21] . On the other hand it is well known that “when interacting with outer space, physical objects never discover their complete internal complexity potential. Hence incomplete or even phenomenological description of physical phenomena or physical objects sometimes responses better to the main point, and to the understanding of what is going on” [2] .

For example, when solving kinetic equations, the fitting parameters determined using experiments are introduced, as a rule. When planning experiments with new or weakly studied objects, there appears necessity to perform preliminary evaluations of final dimensional and timing characteristics of the investigated processes, i.e., to determine dependence between typical object sizes and the time of their growth. In this respect, the proposed asymptotic models can be useful when first meeting this or that problem, evaluating it or analyzing primary experimental data. For example, when studying cluster radioactivity described in [22] , this method made it possible to determine the spectrum of the most probable mass numbers of relatively light nucleon clusters formed inside heavy nuclei [14] [23] . This method also makes it possible to evaluate both average and maximum mass numbers of super-heavy nuclei and formed inside supernovas, and to determine mass numbers corresponding to the “peaks” of abundance of elements in Galaxy, numerically computed in [24] . In addition, it turned out to be possible to get estimates of time formation and typical sizes of neutron stars at supernova outbursts. In the field of mesoscopics this approach is valid and reasonably effective for materials with profound strong interatomic bonds: carbon nanostructures with covalent bonds as well as for nanostructured particles of light actinides with strong approach of atoms [15] . The expediency of pursuing development of asymptotic models [15] in the field of mesoscopics is connected with both the creation of ultra-hard carbon materials and the problem of phase stability and plutonium polymorphic transformations.

Below are arguments for substantiating and more detailed presenting of formalism of a general kinetic approach to description of formation and growth of objects. Also, it is reasonable to present the expanded fields of application for the proposed asymptotic models.

2. General Kinetic Approach

Let us consider formation and growth of compact clusters in conservative stochastic systems defined as assemblies of total-mass-limited quantum objects interacting with each other randomly. We will describe the process of the irreversible aggregation of objects using the concept of distribution density wave in the space of cluster sizes a. The wave propagates with the time t toward an increase in the cluster size. Such one-dimensional approach allows one not to take into account deviation of a geometric shape of the object from the ideal one. On the assumption of the above the universal relation for a wave packet half-width and a spectral line half-width for coherent processes one can write down the following uncertainty relation for a coordinate and a momentum in the space of cluster sizes [25] :

. (1)

Here, is momentum uncertainty, is cluster mass, is reduced Planck constant. Momentum uncertainty is equal to the momentum itself in the order of magnitude, i.e., interaction of objects either occurs or does not occur. Physical meaning of relation (1) is in the fact that during the time interval of elementary (single) act of interaction between objects the exact cluster size cannot be determined until this interaction is finished either through capture of one object by another or their partial or complete disruption, or elastic scattering. It is connected with the fact that unless the elementary act is finished, it is impossible to determine to which object each of the interacting surface elements (nucleons in case of nuclear matter and atoms in case of crystal mesostructures) is assigned. The sign of approximate equality reflects quasi-coherent behavior of considered systems specified by possible fluctuations connected with the specific character of high intensive processes (high-velocity collisions breaking the initial structure of objects; inharmoniousness of atomic oscillations at environment high temperature).

Formally, it follows from the relation (1) that, i.e., an instantaneous velocity of growth of object size is infinitely large:

Because of random character of aggregation process of investigated objects the analytical derivative cannot be determined, and the index of the fractional derivative is equal to 1/2. To describe such a stochastic process using any analytical function F(a) one should use a differential of the type [26] .

In the closed system of clusters the distribution density function satisfies the condition of total mass conservation:

, (2)

where are current total masses of seeds and clusters, respectively, amax is the maximum cluster size in the system; and the mass of volumetric close-packed cluster can be defined either as, where are size and mass of a cubic shape seed, or as, where is cluster material density, α is geometric factor (for a cube α = 1, for a sphere α = π/6). At the end of aggregation process the total mass of seeds is infinitely small, and hence, the first term in conservation Equation (2) can be neglected. Total cluster mass can be considered as integral of motion, i.e., one can consider that the mass increment at t®tf tends to zero (tf is the time instant of aggregation completion). Revealing increment of a certain integral in expression (2), we obtain

(3)

Since at the end of the process cluster sizes and quantities are large comparing to those of seeds, the “simultaneous” failure of a significant number of clusters at their interaction with seeds is highly unlikely. Frozen distribution of clusters in sizes is established in the system, and an integrand in expression (3) does not change its sign, at least, within the time interval (tf – Dt, tf). Hence from (3) it follows that

Introducing the dimensionless variable a/a0 and revealing the last expression according to rules [26] , for the above differential dF we obtain:

With high the second term in braces can be neglected. It is finally found that the density function of size distribution of large clusters at the end of the aggregation process can be approximated by the following relationship:

(4)

This asymptotic kind of distribution obtained using uncertainty relation (1) is not connected with the initial distribution, and does not depend on the type of a kinetic equation for. Distribution (4) can be regarded as “invariant” of any kinetic equation describing behavior of clusters with quantum properties.

In accordance with the above kind of analytic function differential, evolution of distribution density function during a stochastic process of aggregation can be described in diffusion approximation using FokkerPlanck equation written for the space of cluster sizes a:

(5)

Here, is the average rate of kinematic transport of and diffusion coefficient in the space a, respectively. To linearize this equation and obtain approximate solutions we will accept that the average rate of kinematic transfer of matches in order of magnitude the growth rate of an average cluster size:

In accordance with the ratio (1) the diffusion coefficient can be written as:

We introduce dimensionless variables, where is a typical time scale of object interactions determined by their nature and aggregation mechanism. Then, taking into account the accepted assumptions, Equation (5) takes the following form:

(6)

Here, is dimensionless growth rate of average size of clusters, and. To estimate the value we will draw on ratio (1) and the mass conservation law in an elementary process of interaction between a cluster and a seed in the assumption that seed mass gets equally “smeared” on a cluster surface:

Here, is the density of cluster matter. From the conditions we will find that. Then, from ratio (1) we obtain that. Hence is a small parameter.

Revealing the derivative in the third term of Equation (6), we obtain the following equation:

(7)

The study into asymptotic properties of Equation (7) performed in [25] has shown the following. With high , the law of average size growth becomes:

(8)

The constants are determined by an expression coming from uncertainty relation (1) transcribed for an average size of clusters:

(1')

At, we obtain the following expression:

(9)

Value is given from the condition, i.e., a right-hand portion of expression (9) should not depend on. Specific values of constants in law (8) are related with the value, i.e., they are determined by a mode and mechanism of object interactions in a closed system.

Partial solution of Equation (7) at looks like:

(10)

Here, C(l) is a function of the arbitrary real number l. General solution fitting the law of mass conservation and asymptotic distribution (4) (with C(l) = Bl5, B º const), can be written as:

(11)

Comparing this expression with distribution (4) one can conclude that at t®tf the condition is to be satisfied. At the end of the process of object irreversible aggregation as a result of interaction between large clusters and small seeds in a stochastic conservative system the distribution density of large clusters is in inverse proportion to their masses, and average size of clusters is significantly less than the final maximum size.

Approximate view of -function close to maximum of distribution density in size curve can be determined by omitting in Equation (7) terms with the first variable derivative, since in small neighborhood of local maximum the analytical function changes poorly:. Having written the function as a product of space and time factors and having solved the equation for the spatial factor P by Ventzel-Kramers-Brillouin method [16] , we obtain the following expression:

(12)

where,. The presence of the oscillating factors in {…} in the expression (12) means that function P, in principle, can have variety of local maximums with coordinates. What this means is existence of variety of the most probable (stable) sizes of clusters. The constant in the expression (12) is found from the condition that the first maximum (n = 0) corresponds to seeds. Since, the parameter. Then the condition of maximums looks like:

. (13)

The condition of -function “zeros” with consideration for small and values looks as:

(14)

The presence of many “zeros” in -function means the existence of unstable sizes of clusters.

We consider size distribution of “small” clusters for which kinematic transport of distribution density is expressed weaker than diffusion. Having omitted in Equation (7) the second term responsible for transport and having solved the obtained equation, we will get the following expression for a space factor:

. (15)

Function maximum condition (15) takes the form:

(16)

“Joining” of maximum condition (16) with the first maximum condition in expression (12) allows one to determine the arbitrary constant:

(17)

With the knowledge of, it is possible to determine the coordinate of -function first “zero”:

(18)

“Joining” of expressions (14) and (18) with k = 0 makes it possible to find a constant value in expression (14) and write the following condition for coordinates of function P “zeros”:

. (19)

What this means is existence of unstable clusters with corresponding to expression (19) sizes, the first “unstable” size in this case is equal to.

Expression (15) shows that with the distribution density of seed failure products is approximately proportional to their masses. This is incompliance with the concept that disintegration of seeds into reasonably large fragments is more likely (advantageous in energy), than complete failure into individual components.

The linearization of Fokker-Planck kinetic Equation (5) allows one to establish basic properties of cluster size distribution density function. Expressions (8), (9) can be used for finding laws of time-dependent average size growth of large clusters. We consider two limiting growth modes: 1) seed small flux, when each of the seeds has time to occupy its energy advantageous place on cluster surface before interaction with the following seed begins; 2) seed large flux, when seeds affect the cluster almost simultaneously which is equivalent to cluster collision. In the first case ti parameter does not depend on a cluster size and the law of growth can be written as:

(20)

Here t1 is the time unit in the mode of seed small flux. In case of seed large flux we accept the ti parameter equal to the perturbation back-and-forth travel time in a cluster:, where c0 is velocity of perturbation spread in cluster substance (sonic velocity). Then the law of growth can be written as:

(21)

Here t2 is the time unit in the seed large flux mode (cluster collision).

It should be noted that laws of growth for average size objects can be deduced directly from uncertainty relation (1'), if to substitute in it, determine, and then to solve the differential equation with zero initial condition. With the seed small flux we obtain the following chain:

Solving the last differential relation in quadratures, we will obtain the following expression:

(22)

This exactly agrees with law (20), if to reveal the t1 parameter there. Similarly, with we obtain:

(23)

This agrees with law (21). In the first approximation one can assume that time-dependent laws of growth (22), (23) for average size of objects are invariants of any kinetic equation describing aggregation of clusters with quantum properties.

In number of problems connected with investigations of high-intensive processes with high energy (for example, at high temperatures of crystal objects) it is reasonable to determine the value of time interval of elementary interaction act based on the “energy-time” uncertainty relation:

(24)

Here, is energy level width of an isolated excited state of a quantum mechanical system. Substituting expression (24) in relation (1'), one can write the following chain of differential relations:

(25)

The value is determined by nature of objects and a process mode. For example, in case of phonon excitation of mesoscopic crystals at high temperatures (is Boltzmann constant, is characteristic frequency), the value is directly proportional to number of atoms, covered by phonon excitations: [27] . The value can depend on the average size of an excitation region. Having solved differential relation (25) in quadratures one can find approximate laws of average size growth with time in high-intensive processes of object aggregation followed by multiple production and elimination of phonons at high temperatures.

We consider an issue of upper limits of object sizes. From uncertainty relation (1) and mass conservation condition in elementary interaction between a large cluster and a seed, when, we obtain the following expression for the maximum object size

(26)

Here, is a minimum time interval of an elementary act of object interaction determined by physical nature of the process.

3. Subatomic and Astrophysical Objects

In this section various nuclear phenomena are viewed from unified point of view as processes of formation and growth of compact clusters with pronounced collective quantum properties (strong interaction) in a closed stochastic nucleon system in an excited state. As the excited closed system one can consider a mother nucleus, or two contacting nuclei after their collision (double system). At the initial time such system consists of nuclear matter seeds: alpha-particles and (or) their fragments—tritons and deuterons. As a result of random interaction of these seeds their coalescence takes place. Compact accumulations of nuclear continuum (clusters) are formed. Nucleon oscillations on the surface of adjoining “nuclear” clusters having temporarily (virtually) broken external bonds can result in mutual compensation of these unused bonds. As a result the adjoining clusters consolidate, and larger nuclear objects with various masses (mass numbers A) are formed. Due to Coulomb repulsion these objects can escape from the region of interaction. The proposed model can be applied to the following problems: 1) calculations of the most probable mass numbers A of clusters within the entire range of the existing nuclide systematization [28] [29] ; 2) computed estimation of a discrete spectrum of the most probable numbers A, typical for cluster radioactivity investigated in [30] [31] ; 3) discovering the most probable average mass numbers of clusters formed in processes of spontaneous nuclear fission [32] [33] ; 4) determination of the most probable values A and evaluation of the average value of superheavy clusters formed in processes of nucleosynthesis in stars [24] [34] . Only fast intranuclear processes (r-processes) that go during lifetime of excited states are considered. These states are given by the energy of oscillatory and rotational levels determined in shell model of nucleus. The time of processes of approaching equilibrium is within the interval between the time of direct nuclear reactions (with typical nuclear time of ~10–22 s) and the lifetime of a compound nucleus (10–16 - 10–15 s [29] [35] -[39] ). In addition, there are some attempts to determine separately typical formation times and nuclear matter cluster sizes in astrophysical scales, in particular, for neutron stars.

The basis for analytical description of the processes under consideration is a phenomenological concept (given in Section 2) of the object distribution density wave in space of their sizes a. Such one-dimensional approach allows one not to take into account deviation of a real geometric shape of nucleus from spherical one. Evolution of function during a stochastic process of object aggregation can be described in diffusion approximation using Fokker-Planck Equation (5) written for the space of cluster sizes. To determine the most probable mass numbers we will use expression (12) with an arbitrary real number found out in Section 2 and known expressions [29] [39] connecting the nucleus radius R = a/2 and the mass number:

(27)

where r0 = 1.3 fm is a typical parameter of strong interaction. If one accepts that this expression is approximately fulfilled for seeds as well, then from Equation (12) it is possible to obtain the following expression for the most probable mass numbers An:

(28)

Here, A0 is seed mass number, and the parameter is determined by a seed mass and size and the typical time scale ti of interaction between objects. It is reasonable to accept as such scale a period of high-frequency nucleon oscillations in a nucleus. This parameter can be determined as, where is an average thermal velocity of nucleons in degenerate Fermi gas determined by the known value of average kinetic energy of ideal Fermi gas per a nucleon, equal to EF  = 22 MeV [39] .

To find typical average sizes of cluster nuclides we will use expressions (22) and (23) obtained in Section 2 and transcribe them as follows:

, (29)

. (30)

Here, r is density of nuclear matter. In order to eliminate the parameter a0 from equations, it was accepted that. The absence of a seed size in (29) and (30) reflects an asymptotic stage of the process, when the system “forgets” about initial conditions and its evolution goes on according to internal dynamics of nucleus. This dynamics is connected with excitation of oscillatory and rotational energy levels.

General time of the process is determined by shell energy level G using Heisenberg rule [35] [36] :

. (31)

It is correct for isolated excitation of level in quantum system.

Verification of the proposed formalism can be implemented by the example of “rigid” (deep inelastic) processes involving hadrons described in [40] [41] . We consider these processes as the following chain of events: head-on collision of “nucleon-nucleon” or “lepton-nucleon”®quark knocking-out from gluon cloud®random interactions in intra-hadron medium®new particle formation (hadron jets). In [23] [25] the new particles are considered as clusters consisting of intra-hadron medium components (partons): quarks, gluons, quark-antiquark couples. A condition of generation of new particles can be written as:

Here, is total energy of partons, is total energy of new particles, is total energy of two original colliding particles. At the end of this process all the partons transformed into bound states, , therefore relation (3) is fulfilled, and the proposed in Section 2 stochastic approach is applicable. Then law (21) for describing the evolution of cluster average size provides the following relation between length aunit, time tunit and mass munit units under deep inelastic interaction of fundamental particles:

. (32)

Here c is the maximum velocity of interaction propagation in hadron medium: the speed of light in vacuum. Based on relation (32) and the data from [19] for “conventional” quark masses (uand dquarks) deep inside hadrons and for typical distances, one can obtain the following estimation of the time unit for current quarks (cq) in the state of asymptotic freedom (af):

The quark transition time ttrans from asymptotic freedom to a bound state (confinement) inside a hadron is evaluated according to Formula (21) as . The obtained value corresponds to the time scale of strong interaction.

In processes considered in [40] the minimum value of length unit is. It can be considered as the top estimation of current quark size. Then it goes from relation (32) that, and the time for formation of hadrons having an average size of, in accordance with Formula (21), is. This value is significantly higher than the life time of unstable hadrons (“resonances”):. This reflects formation of stable hadron jets in described processes of catastrophic collisions.

If in accordance with the results of [42] , we accept that perturbation propagation rate in quark-gluon medium is equal to sonic speeds cs = 0.3 c and 0.57 c, then in the above estimates of the time units (the time of quark transformation into a bound state and the time of hadron formation) factors 1.8 and 1.3 will appear. These corrections do not change the estimate of strong interaction time scale and increase the inequality reflecting the formation of stable hadron jets.

The asymptotic distribution j µ m–1 means that: 1) probability of nucleon fragmentation in deep inelastic scattering is higher than probability of their conservation; 2) the number of pions in hadron jet is 1.5 times higher that the number of nucleons. These conclusions are in compliance with concepts of [40] about nature of rigid processes. In particular, at proton frontal collision the following reaction takes place: pp®p0p+p+nn. One can see that there are 3 pions and 2 neutrons in the formed hadron jet.

Consequently, the proposed formalism allows adequate description of an asymptotic form of well-known accidents in the world of elementary particles. This gives grounds to try to apply it in solving the above described problems of intranuclear processes description.

Values An calculated in [14] [25] with the help of Formula (28) with three types of seeds correspond to the mass numbers of the nuclides and their isotopes over the currently known range [28] [29] , including the transfermium elements. In addition, calculated mass numbers of cluster nuclides are also in compliance with standard concepts [29] [43] on isotope “stability valley”. It should be kept in mind that the calculated clusters are formed deep in the sea of nuclear matter and not in particular reactions. In escaping from either the parent nucleus or the double system clusters can emit or capture individual nucleons and seeds of nuclear matter.

The established correspondence allows the characteristic values for various processes of approaching equilibrium (cluster radioactivity, spontaneous fission, synthesis of heavy and superheavy elements) to be extracted from the above values.

Calculated mass numbers of clusters formed inside heavy nuclei in processes of cluster radioactivity make up the following discrete spectrum: A = 6, 7, 9, 11-16, 18, 20, 24, 26, 28, 30, 32, 34, 37-40. In experiments on cluster radioactivity the values A = 14, 24, 28, 32 were recorded. Under discussion is the possibility of recording the outgoing clusters with mass numbers close to stated predicted spectrum (e.g., see [22] [30] [31] ).

In the field of spontaneous fission it is possible to single out two sets of the most probable cluster mass numbers, each being repeated for two or three seed types: A1 = 72, 79, 84, 88, 103, 108 и A2 = 126, 134, 139, 142, 147, 155, 158. The arithmetic means for the two sets of mass numbers are equal to и. The obtained values correspond to the coordinates of the maxima of the total nuclear fission yields [32] [33] . It is possible to describe spontaneous nuclear fission as a result of excitation of the first rotational level with the energy Grot = 100 keV using the Formulae (27) and (30). The total time of an intranuclear process corresponding to the given value Grot, is equal to trot = 6.287 × 10–21 s, and the Formulae (27) and (30) determine an average mass number of light fragments as. An average mass number of heavy fragments < A>heavy is determined from the condition of mass conservation; Amat is mass number of parent nucleus. In case of light actinides (Th, U) we get = 135-138. In diagrams [32] [33] of fragment mass distribution in spontaneous heavy nuclear fission (A = 235, 238) an average mass number of lighter fragments is 90-100, an average mass number of heavier fragments is about 140.

If one assumes that fast nucleosynthesis in stars occurs as a result of transition of nucleon system from high-frequency vibration level with the period ti = 5 × 10–23 s to the first rotational level with the typical energy 100 keV during the lifetime of this level trot = 6.287 × 10–21 s, then Formulae (27) and (29) provide the following estimate of an average mass number of superheavy elements: = 330. The obtained value agrees with the prediction from [24] about possibility of nucleus formations with A > 300 in neutron stars. As for the induced fission as the process competing with nucleosynthesis, in this case one can accept as a typical scale the sonic wave travel time in a nucleus formed as a result of nucleosynthesis: t0 = a/c0. With A = 330, c0 = 5 × 107 m×s–1, we obtain ti = 3.6 × 10–22 s. Then from Formulae (27) and (29) we obtain that with the above lifetime of the first rotational level 6.287 × 10–21 s the average mass number of light fragments of superheavy elements is. In compliance with the mass conservation law the average mass number of heavy fragments is equal to. The obtained values correspond approximately to the first and the last peaks of final abundance of elements in Galaxy [34] . This compliance gives grounds to expect that the above estimate of the average mass number of superheavy elements that can be formed as a result of nucleosynthesis in supernovas and neutron stars is reliable.

One can try to estimate a mass number of the final nuclide Aend formed in r-processes of nucleosynthesis in stars, if he uses expression (26) for the maximum cluster size in a closed system and the known expression (27) for the mass-number-dependence of nucleus radius. If as in (26) we assume the travel time of an electromagnetic wave over a nucleus equal to a/c (c is the velocity of light in vacuum), then from espressions (26) and (27) we will obtain the following equation for a mass number of the final nuclide:

This, with alpha particle (A0 = 4) taken as a seed, yields to Aend = 470. According to the expression (28) mass number 470 corresponds to triton as the seed. This value is by approximately 15% higher than value 408 responsible for hypothetical continuation of 2β-stability line [43] . Hence the obtained estimate is in compliance with the known systematization of nuclei.

The above time trot ~ 10–20 s of intranuclear processes of nucleosynthesis and fission is negligible as compared to the given in [24] [34] times tyield = 0.3 - 2 s, when performing numerical simulation of nucleosynthesis and fission yield in astrophysical phenomena. In accordance with concepts about supernovas and neutron stars published in encyclopedia [44] [45] , the times τyield are connected with hydrodynamic processes of neutron star formations in supernova outbursts as a result of gravitational compression of the star central region.

It is possible to determine times of neutron star formations by an analytical way. For this purpose we will use the analogy with the concept [19] on introduction of phenomenological constants. From considerations of dimensionality and by analogy with the notion of Planck mass, determined through the action constant in microcosm, the gravitation constant GN and the light speed c, we will write the following expression for the action constant KG in the “world” of collapsing stars:

(33)

Here is the Chandrasekhar critical mass [45] [46] above which the star gravitational compression begins, is the mass of the Sun. When postulating the uncertainty principle in gravitational processes by analogy with Heisenberg “energy-time” uncertainty relation for microcosm [46] , one can write the following expression for the lifetime of excited state in gravity system:

. (34)

Here DE is energy level width of an excited system determined by the energy removed by radiation of every type (neutrino, gamma, etc.). With Type I supernova outburst the radiated energy is equal to DE = 1042 - 1043 J [44] [45] . Then from relations (33) and (34) with, GN = 6.67 × 10–11 kg–1∙m3∙s–2, we obtain that the excited state lifetimes of collapsing stars are within the range. The obtained values correspond to generally known astrophysical concepts about times of “cold” neutron star formations as a result of supernova outburst. With Type II supernova outburst the radiated energy is equal to 1041 J [45] . According to relations (33) and (34) this provides the lifetime of the star excited state. This value corresponds to astrophysical concepts about the time of cooling-down of a “hot” neutron star formed as a result of slow escape of supernova shell [45] which is equal in order of magnitude to 10 s. It should be noted that the discovered time range of fast and slow neutron star formation 0.17 - 17 s covers the range (set in [24] [34] ) of times of nucleosynthesis and forced fission yields and times of formation of final element abundance in Galaxy. Hence the time trot ~ 10–20 s corresponds to formation of heavy and superheavy clusters inside composite nuclei, and the times correspond to formation of recorded distribution of elements bulk formed during processes of gravity compression of supernova central regions.

We try to determine a size of the gravitating nuclear matter cluster (neutron star) on the basis of “momentumcoordinate” uncertainty relation, written for astrophysical objects: 

. (35)

Here, M is object mass, V is typical velocity of compression process, as is object localization, i.e., size of star. From here we get the following expression for object localization:

. (36)

In case of relativistic collapse (V ~ c) the object size is equal to. Since M ~ MCh , as is lower than the gravitational radius rg = 2MGN/c2, which means the formation of black hole. If a typical process velocity corresponds in order of magnitude velocity of supernova shell escape (about 20,000 km∙s1 [45] ), then object localization is 1.56 × 104 m. The obtained value corresponds to the minimum size of a typical neutron star equal to 16 km [44] .

4. Mesoscopic Objects

Since the proposed kinetic approach provides an adequate description of processes in microcosm, it is possible to try to use it for considering object growth processes in mesoscopic scales. This section provides a data review are in view on development of the proposed asymptotic method for investigating “mesokinetics” of growth of various nanostructured objects with strong interatomic bonds [15] . Among these are carbon nanostructures, for example, single-wall nanotubes and detonation nanodiamonds [21] [47] (covalent bonds) as well as light actinides with strong approach of atoms [9] (“quasi-covalent” bonds). It is assumed that the growth of nanostructures is caused by their phonon excitations and vibrational interactions, leading to mutual compensation of free electronic bonds of “boundary” (i.e., surface) atoms in adjoining objects and to consolidation of these objects into more large patterns, which are compact clusters.

At the initial stage of crystal objects growth, as a result of random interactions between seeds carrying crystal structure elements, small crystalline particles having the short-range-order are formed. These small particles continue to interact with seeds and between themselves, causing their growth. In the process of this growth the quasi-long range crystal ordering arises in the particle, and the lower limit of complexity for a macroscopic solid is reached. Hence the growth process of mesoobjects with nanocrystalline structure can be considered as formation of clusters in continuum initially consisting of seeds with the set spatial orientation of atoms [15] [25] .

With phenomenological description of single-wall carbon nanotube growth as a result of capture of hexagon seeds based on uncertainty relation (1) the following formula was obtained for estimating a typical nanotube diameter [48] :

where qv is characteristic oscillation temperature of covalent bond tension in C-C carbon atoms, lcov is carbon covalent bond length (is atomic mass unit, is the Boltzmann constant). It follows here from that a typical diameter of a single-wall carbon tube is d = 2.2 nm (qv = 2154 K, lcov = 0.15 nm). This value is in the middle of experimental range d = 1 - 4 nm given in [21] . The expression obtained in [48] for the growth of an average length of a single-wall carbon nanotube looks like:

(37)

The parameter, by analogy with quantum mechanics [49] , can be determined as the upper limit of a diffusion coefficient in the space of cluster sizes. The expression obtained for represents a quantum analog of a classical diffusion solution [50] of a problem about growth of particle of new phase and specifies the least possible time of growth in the system initially consisted of seeds (hexagons). From expression (37) we obtain that the least growth time for a single-wall carbon nanotube with the given in [21] length of 40 nm is 4 µs. It should be noted that in this case we are dealing with the time of growth for an individual nanotube, not with hardware time of synthesis of a bulk of these objects. Of some certain interest is the condition determination for performance of experiments, when single-wall nanotubes with the given average size specified according to expression (37) will grow as a result of phonon effects.

When considering formation and growth of volume-packed objects at average temperatures compared to Debye temperature, we will distinguish between “small” clusters in which atom behavior corresponds to oscillation spectrum of seeds not having a crystalline structure, and “large” clusters in which behavior of atoms corresponds to dynamics of crystal lattice oscillations. A crystal lattice can have either the short-range-order or quasilong-range ordering depended on a cluster size. At the initial stage in processes of collisions between seeds the “first” clusters are formed. As a result of interactions between these clusters and seeds the larger objects are formed in the system. A typical time scale for them can be evaluated as a product of a number of bonded seeds in an object equal by order of magnitude to and the time of excitation of a single seed:

The parameter can be determined as a product of number N of chemical bonds in a seed and the inverse typical frequency of atom oscillations in its structure, namely,. Characteristic oscillation frequency can be expressed through the relevant characteristic temperature based on the relation. Then in compliance with the expressions (8) and (9) the approximate law of growth for an average size in processes of interaction between small clusters and seeds and between them will be written the following way:

. (38)

Here, is a time unit in interaction processes between small clusters. By classical analogy with notions [50] the exponent 1/2 in the Formula (38) corresponds to the diffusion of seeds at the cluster.

Clusters of size progressively increasing in time are formed in the system. They can interact with both seeds and between themselves as well. The minimum object size at which an oscillating spectrum of a small cluster switches over an oscillating spectrum of a large cluster “by a jump”, corresponds to the occurrence of the crystal short-range ordering, and represents the lower boundary of transition region from mesoscopic structures to macroscopic bodies. As crystal structures grow, they tend to the quasi-long-range ordering, and the dynamic system itself “forgets” about initial conditions. It would appear natural that at a certain time some critical average cluster size is reached in this system. This size does not depend on a seed size and is determined only by dynamics of atom oscillations in investigated objects. This corresponds to the quasi-long-range crystal ordering emerged in the system, i.e., to the lower limit of complexity reached in a macroscopic crystal body. The given size corresponds to the upper boundary of a transition region from mesoscopic structures to macroscopic bodies.

The characteristic size of objects, beginning with which the cluster has got the crystal lattice with short-range order, and the characteristic size, such that quasi-long-range crystal ordering arises in the cluster, and the appropriate times of achievement of these sizes were obtained in [15] [25] on the basis of the laws (38), (22, 23) as follows:

(39)

(40)

(41)

(42)

Here, is Debye parameter, n is concentration of atoms in a crystal lattice, , E0 is seed energy in the environment determined by equations for inner energy of nonlinear polyatomic molecules given in [51] ; is seed energy inside a large cluster determined by an oscillating spectrum of a small particle [52] :

are typical temperatures of tension vibrations and flexural vibrations of atom bonds in a molecule, respectively; is number of atoms in a molecule, is number of bonds:

Here, is energy of the basic lattice oscillating state (energy of “zero” oscillations), is internal (dynamic) fractal dimensionality of cluster, is oscillation frequency of atoms in the lattice. Formulae (39), (40) determine the area of average cluster sizes which is a transition region from mesoscopic structures to macroscopic particles.

The growth laws for crystal particles from seeds (polyatomic molecules) look like [15] [53] :

(43)

(44)

Here, M is number of atoms in a seed. Formula (43) describes growth of large clusters under small flux of seeds, Formula (44) corresponds to both large flux of seeds and interaction between large clusters, sound speed being recorded in Debye approximation as.

In case of detonation synthesis of nanodiamonds, a “Lagrangian” layer of shock compressed confined matter (i.e., a layer with a constant mass) is viewed as a closed system [53] . It is assumed that deformation and fragmentation of benzene rings are taking place in a shock wave front that results in production of structural seeds of diamond phase with cubic modification—molecules of cyclohexane carbon skeleton. An effective size of such a seed is determined by the length of a covalent bond of carbon atoms and is equal to a0 = 0.3 nm. The estimated time of seed formation tnf ~ 10–12 s was obtained in [53] from kinetic considerations. As shock compressed layer expands, crystal diamond-like clusters are formed in it as a result of seed coalescence and oscillating interactions between evolving nanostructures and seeds and between themselves.

The minimum average size of diamond-like nanoclusters starting from which they possess a crystal lattice, is estimated as. The cluster of this size has seeds, i.e. about 160 carbon atoms. In accordance with Formula (41) the formation time of such cluster can be estimated as. Using Formula (40) and the time of seed capture determined in [54] as one can estimate the cluster minimum size starting from which it can be regarded as a macroscopic particle:. This cluster includes seeds, i.e., approximately 2000 carbon atoms. From Formula (42) we obtain that the time of synthesizing of this particle is equal to. The further growth of diamond-like particles is described by Formulae (43) and (44). Estimates provided in [53] [54] show that an average size of detonation nanodiamonds about 4 nm can be reached in a time of about (2 − 5) × 10–9 s. These estimates as well as asymptotic mode (4) of cluster size distribution and the idea about existence of local maximum variety within this distribution agree with key experimental data obtained by different authors (see [25] ).

Shock-induced coalescence of detonation nanodiamonds of “water” synthesis under dynamic loading in recoverable ampoule is described in [53] [55] -[57] . Growth of particles is described by the following formulae obtained using the above laws (22), (23):

(45)

(46)

Formula (45) corresponds to a small flux of seeds (nanodiamonds), when an average size growth occurs as a result of “particle-seed”—type interaction. Formula (46) corresponds to a large flux of seeds, when growth occurs as a result of “particle-particle”—type interaction. Because the compression pulse duration is short (t ~ 10–5 s), it is natural to assume that the crystal particles predominantly grow under metastable conditions for the diamond in the phase diagram of carbon: at zero pressure and increased residual temperatures. Cooling-down time for shock heated porous carbon sample is determined using the solution of problem about a heat wave starting from an instant plane source (a thin disk) formed under shock compression of nanodiamond powder and propagating of this wave in material of the ampoule. Estimates performed according to Formulae (45), (46) show that during a cooling-down time of about t ~ 102 s two subsystems of crystal diamond-like particles with average sizes about 1.5 × 10–5 m and 1.5 × 10–4 m are formed in porous sample of nanodiamonds. The mechanism of crystal growth resulted from consecutive connection of “small” seeds to a “large” particle described by Formula (45) provides by an order of magnitude larger sizes than the mechanism of large particle coalescence described by Formula (46). The calculated results agree with the given in [57] experimental bimodal distributions of polycrystalline diamond particles under consideration.

Behavior of light actinides in conditions of internal self-irradiation and external effects in many cases is determined by their phase stability and a character of crystal phase polymorphic transformations [9] . The problem of phase stability and polymorphic transformations of light actinides is considered in terms of phenomenological concepts about evolution of nanoclusters as a result of their phonon excitation and oscillating interactions [15] . The basis for such concepts is, for example, the existence of strong (quasi-covalent) bonds in multiphase plutonium resulted from f-electrons effect on atoms approaching each other [9] .

First, we will single out the following process: a) spontaneous formation of plutonium alpha-phase (a-Pu) from plutonium delta-phase (d-Pu) as a result of sample aging; b) martensitic transition of plutonium delta-phase into alpha-prime (a¢) phase on cooling of samples and reverse transformation on their subsequent heating. We will suppose that at the initial time there is a fairly large number of starting seeds of a certain crystalline phase which could be generated, as a result of temperature fluctuations, on defects of the initial structure either affected by intrinsic radioactive radiation, or as a result of local shear deformations on lattice defects under external action. Since the a¢ phase of Pu differs from the a-phase by the presence of a relatively small number of gallium (Ga) impurity atoms (~ 1% - 2%) and crystal lattice parameters of these phases are close to each other (see, for example, [9] ), we are going to consider the «pure» transitions d ®a ®d  for simplicity.

Therefore, the influence of impurity atoms on kinetics of phase transformations is ignored in a proposed model. As starting seeds of a-Pu we will consider nanoclusters consisting of eight atoms with the only crystallographically possible mutual orientation described in [9] . A characteristic size of an a-phase seed, accepted to be equal to the most linear size of an elementary crystal cell, is a0a = 1 nm; the time of its formation from dimmers of d-phase is evaluated as:

Here, ea is internal energy of a-phase seed, and ed is internal energy of the dimer. Due to connection of two such clusters a unit cell is generated which consists of 16 atoms, and which forms a-Pu monoclinic lattice. Dimmers, tetrahedrons, as well as plane formations of six atoms generating glide planes of face-centered cubic lattice (fcc) can be considered as structural seeds of d-Pu. Irreversible aggregation of seeds results in formation of small (nano-) particles with quasi-crystal (short-range) ordering of atoms in a lattice. Estimates show that the minimum size and time for forming a particle with a short-range ordering in a-Pu crystal lattice are equal to ~ 2 × 10–9 m, tmin ~ 1.9 × 10–10 s, respectively. Estimated parameters determine premartensitic state described in [9] . The following coalescence of such nanocrystal objects brings to formation of larger particles with quasi-long range ordering.

Based on uncertainty relation (1') and the resultant expression (46) in Debye approximation for description of nanoparticle coalescence in [58] , it was stated that at modest temperatures the spontaneous generation of a-Pu polycrystals with sizes of about 10–7 m occurs for the time t = 6 × 10–3 s. Formally it was also obtained that a-Pu “macroscopic” particles with sizes 0.1 - 1 mm can be generated for times of 2 × 10–1 - 2 × 102 years. It should be noted that these estimates are in compliance with the generally known concepts about slow kinetics of spontaneous transformation of alloyed d-Pu in normal storing conditions.

The main mechanism of martensitic transition in the proposed model is addition of starting seeds (polyatomic molecules) to large clusters in the regime of small flux of seeds, when each of them has time to occupy its energy advantageous place on a surface of a large cluster before the interaction with the following seed begins. In this case growth of a martensite particle is described by Formula (22) in which the process finish time, according to concepts [59] about ordering kinetics at phase transformations, is determined as:

(47)

Here, l is a characteristic size of an initial matrix from d-Pu, e.g., a grain average size in a homogenized sample, cd is velocity of propagation of minor elastic perturbations in the initial matrix. A limited average size of particles formed as a result of capturing the seeds with oscillating spectrum of polyatomic molecule [60] automatically follows from (22), (43), (47):

(48)

Here, qDa. is Debye parameter of a-phase. Temperature T0 of the beginning of martensitic transformation is determined from the condition that an average a-phase particle size becomes equal to a characteristic longitudinal size llong of a twin. If this value does not depend on T, then it follows from Formula (48) that with the decreased l, the value T0 has to decrease as well, since in oscillating spectrum Ev of solid particles there is a temperature-independent term connected with zero oscillations of a crystal lattice [15] [52] . Finish temperature T* of martensitic transition is determined from the condition that average particle size of a-phase becomes equal to a critical size at which a quasi-long-range ordering appears in the system, that is from the condition of intersection of the curve (geometric place of points in a plane (, T) determined according to Formula (40)) with the curve in the plane (, T) by Formula (48):. In this case the condition is fulfilled automatically. With lower temperatures T < T*, when, , the quasi-long-range ordering does not have time to be formed, and the system does not reveal explicit macroscopic properties. From the condition we obtain the following equation to determine a “critical” temperature of martensitic transition finish:

From here, it follows that with the decreased l, temperature T* should increase through the presence (in the expression for Ev) of the term responsible for zero oscillations and independent on temperature.

The calculated curves T(l) given in [25] with characteristic temperatures of martensitic transition show that in a diffusion model [60] , with the decreased matrix grain size, transformation “blurring” determined as the difference T0 – T*, decreases also, which is caused by the effect of crystal lattice zero oscillations under significant sample.

Now, we consider the reverse transformation a®d during the heating of the sample preliminary cooled to the “critical” temperature T*. In this case, the vibration spectrum of the small a-Pu particle surrounded by the d-Pu matrix transforms into the vibration spectrum of the environment due to the interaction of atoms at the interface. This process propagates deep into the particle and leads to the transformation of its vibration spectrum in the bulk; i.e., the polymorphic transformation occurs. According to the differential relation (25) when, we can obtain the following relationship for the growth of the reverse transformation region during the heating over the entire volume:

(49)

Here, D(x) is the tabulated Debye function with x = θD/T [61] . Formula (49) determines the lower time boundary of the inverse transformation which is evaluated in [15] [60] as tmin ~ 10–8 s.

As the size of inverse transformation region increases, which is described by Formula (49), this region can capture some individual groups of atoms (atomic clusters) from the initial matrix surface from the d-phase in the region of “contact” with a-phase and recovered d-phase. In addition clusters can be captured from a-phase martensitic particles. The inner energy of clusters can correspond to either intermediate state between two crystal phases (diatomic and polyatomic molecules), or atom oscillation spectrum in a solid. The law of growth of reverse transformation region in the mode of a small flux of seeds described by Formula (43) looks like:

(50)

Here, E0i is the energy of the captured object determined by a number of atoms in a cluster, index i is phase, Evd is energy of the captured object inside the region of inverse transformation, determined by an oscillating spectrum of a solid particle. It is taken into consideration here that as a result of elasticity here, we took into account that the elasticity softening leads to a decrease in the Debye parameter according to the linear law, q0i is Debye parameter at T = 0 [62] . Formula (50) determines the upper time boundary for the inverse transformation within the proposed model estimated as 2.4 s [15] [60] .

In case of capturing seeds directly from a-phase, when their vibration spectrum Eva corresponds to the solid body vibration spectrum, expression (43) for an average size of d-phase particles looks like:

(51)

Within the temperature range of the reverse transformation a®d the linear dependence [62] for Debye parameter can be approximately written as. In case of dimmers as d-phase seeds, the internal cluster energy is written as:

where the first two expansion terms (using Bernoulli numbers Bn [63] ) are taken into consideration:, B1 = 1/6, B2 = 1/30, ···. If one neglects terms of the second (or higher) order of smallness, the temperature factor in Formula (50) in the first approximation can be written as:

.

From the condition of extremeness:

we obtain the following expression for determining the “relevant” temperature:

At this point the second derivative is negative:. The values Ti determine temperatures at which particle sizes growing from dimers of this or that phase reach their maximum in a plane (, T). In case of d-phase dimers the Ti value corresponds to the temperature Ts of the beginning of the reverse transformation a®d:

(52)

In case of a-phase dimmers, the Ti value corresponds to the temperature Tf of the reverse transformation a®d finish:

(53)

Knowing the critical temperature T* of martensitic transition finish and the temperature Ts of the reverse transformation beginning at heating one can determine the temperature hysteresis Ts – T*.

A temperature factor in Formula (51) for an average size of the reverse transformation region is written as:

It is taken into consideration that at high temperatures the vibration energy in a crystal lattice per an atom is equal to (see [51] )

In the plane (F, T) there exists a special point, in which exponent index is equal to zero, and F(T) = 1, i.e., elastic vibrations of atoms are balanced by elasticity softening. With higher temperatures elasticity softening becomes predominant and inverse transformation “decelerates”. From here we obtain the following expression for the value TM corresponding to the decrease of average size growth rate with temperature for the given number of atoms M in cluster seed with vibrational spectrum of the solid:

(54)

It can be seen from relationship (54) that there exists a discrete spectrum of characteristic temperatures TM each of which corresponds to a particular number M. This discrete spectrum determines a stepwise character of particle growth during the heating due to the capture of clusters with different numbers of atoms.

Discrepancy in calculated and experimental data [64] relating to characteristic temperatures of the martensitic transition beginning and finish and of the reverse transformation is 2% - 24% [60] . The calculated, using Formula (54), discrete set of temperatures TM = 318, 335, 350, 364 K at which rate of growth of a particle average size with atom numbers in the seed M from 6 to 9 decreases, corresponds to the temperature interval [65] 328 - 360 K. Within this interval one can observe “steps” on dilatometric curves for relative extension of samples and the corresponding “bursts” on curves for velocities of sample extension in experiments [65] .

The established correspondence gives foundation for determining the least grain size of the d-phase initial matrix at which martensitic transformation into a-phase is still possible. As a necessary condition for martensitic transformation one can accept that an average size of a martensite particle should be no less than a characteristic longitudinal size of a twin llong: . The proper equation for the minimum grain size of d-phase matrix lower which the martensitic transformation does not occur, since a typical twin size is not reached at cooling down to reasonably low temperatures about 100 K obtained from (48), looks like:

Hence at llong = 1 µm [9] we obtain lmin = 2.2 µm. This grain size is in compliance with the calculated temperatures of the beginning T0 and finish T* of martensitic transition equal to 140 and 110 K, respectively. Values lmin = 2.2 µm, T0 = 140 K and T* = 110 K can be considered as limiting for this mechanism of martensitic particle growth. From Formula (48) it follows that. This means that the grain size less than lmin =

2.2 µm can be reached as a result of thermal cycling, i.e., with sequential cycles of “cooling down to a-phase transition®heating back up to d-phase®the following cooling down to a-phase transition” and so on. Since with d-phase recovery a typical matrix grain size decreases to the value 0.45 l, with the following cooling, the martensitic particle size will be equal to 0.2025l, hence at heating, a typical grain size of d-Pu matrix will be equal to 0.2025l and etc. in decreasing geometric progression. With the given initial size l, it is possible to obtain in just several cycles such size of d-phase matrix grain that will be lower than lmin. The fact that the size of a-phase martensitic particle is always less than that of d-phase initial matrix grain means that the martensite total volume at the following cycles will be always less than the martensite total volume at the previous phase. And at each further cycle the temperature of the martensitic transition beginning will be lower than that at the previous cycle. The results are in qualitative agreement with the experimental data [66] relating to d-Pu stabilization at thermal cycling.

Under long annealing of samples the homogenizing and increase of their grain sizes occur after each cycle that brings to martensite transformations in their former volume. Sample homogenizing can take place as a result of coalescence of particles undergone the a®d inverse transformation. With a point contact of such particles the parameter N in Equation (25) is equal to number of atoms in a channel formed by linear chain of elementary crystal cells arranged within the length 2: where Nc is a number of atoms in a cell, and ac is a typical cell size. Then from (25) we obtain the following equation for the growth of average size of the homogenized region:

(55)

Estimates in accordance with this formula show that with the temperature T = 600 K at the time t = 30 min the size of the region of homogenization reaches a value of about 2 cm. Exactly this explanation can be applied to the data [67] on plutonium sample annealing.

In [15] the nature of d-Pu®a-Pu transition with the sample cooling down to the lower temperatures than those used in [64] is determined. In this case an a-phase particle with the minimum size = 2 nm (the particle having a short-range crystal ordering) can take out of the surface of d-phase surrounding matrix groups of atoms (clusters) with the solid body vibration spectrum. In the course of the interaction of the cluster with the particle the vibration spectrum of the cluster is transformed into the vibration spectrum of the a-phase.

With temperatures the Equation (48) for the limiting average size of an a-phase particle is written as:

Here, it is taken into consideration that in Debye approximation [51] at minor temperatures the “specific” internal energy of a solid per an atom is equal to:

The value basically can reach a typical size of d-phase initial matrix (for example, the size of a homogenized grain), i.e., full transition from d-phase into a-phase occurs. The temperature Tft of full d®a transformation is determined from the condition. At very low temperatures, when in previous equations one can neglects the terms proportional to T3, we obtain the following approximated equation for the temperature of full transformation:

In accordance with Formula (47), the time when a martensitic particle achieves the size l, can be written as:

Hence we obtain that with the grain size d-Pu l = 30 µm the temperature of full d ® a transformation is Tft = 13.5 K, and the time of this transformation is about 9 × 10–2 s. If as a result of preliminary thermocycling with moderate cooling the d-Pu grain size is decreased to the critical value lmin = 2.2 µm, then the temperature of d®a full transformation is Tft = 17 K, and the corresponding time is about 1.3 × 10–3 s. With the given temperatures critical sizes of a-phase particles estimated according to Formula (40) significantly exceed the sizes of initial d-phase grains. For example, at Tft = 13.5 K we have = 220 µm which is significantly higher than the initial grain size l = 30 µm. This means that in the a-phase particle the long-range ordering is not reached. Such particle does not have macroscopic properties and represents a mesostructure with a short-range crystal ordering. At the same time sizes of such object are large as compared to those of a martensitic particle formed at cooling down to reasonably low temperatures about 100 K, when.

An important issue is the one relating to growing of macroscopic size plutonium crystals for investigation of their mechanical properties [67] . In the proposed models one can use as seeds nanoparticles with a critical size a*, at which a particle has quasi-long-range crystal ordering. With small seed flux, when each seed has time to occupy its energy advantageous place on a cluster surface before it interacts with the following nanoparticle, the size of a phonon excitation region can be approximately equal to a*. Then a number of excited atoms is equal to, where a0 is a characteristic size of an atomic cluster (crystal structure seed), M is a number of atoms in cluster. From (25) we obtain that the law of growth looks like:

(56)

Equation (56) determines growth rate of Pu crystals in the mode similar to a free-molecular flux of seeds. In this case nanoparticles play a role of seeds. From considerations of translation symmetry we will accept, for example, that the minimum size of d-Pu initial nanoparticle is equal to a* = 2ac = 0.9274 nm. Then with T = 700 K, t = 10 days we obtain that. The estimate corresponds in order of magnitude to the experimental data given in [67] on Pu large crystal growth by the “strain-annealing” method. Formula (56) can be useful when choosing conditions for growing large crystals of any material with strong interatomic bonds by the method of nanopowder sintering.

5. Discussion of Results

The above obtained approximate analytical expressions for characteristics of object formation include fundamental physical constants, physical properties of substances, and phenomenological parameters of objects. We will consider an issue about relation (32) between time-space scales and mass scales in subatomic physics. Formally relation (32) is exactly fulfilled for Planck values of mass, length and time which proves validity of this relation. Using this relation it was shown in Section 3 that in quark-gluon substance least time unit corresponding to the minimum value of the length unit that can be considered as the upper estimate of a current quark size with, is. A time unit for current quarks in the state of asymptotic freedom with a spatial scale is. Based on relation (32) one can obtain the following estimate for the unit time in processes of interaction with participation of constituent quarks in the state of confinement:, Þ. Time of transition from asymptotic freedom to confinement is evaluated in Section 3 as that corresponds to the time scale of strong interaction, and the time of stable hadron formation is evaluated as. Thus the relation (32) in principle shows that processes in microcosm are characterized by a spectrum of time units of “measurement”. It should be noted that close to the upper boundary this spectrum is overlapped with the “lower” region of the spectrum of typical nuclear times 10–23 - 10–22 s in case of direct reactions. In this way Formula (32), demonstrating that up to distances of about 10–18 - 10–16 m and times of about 10–29 - 10–26 s the “ordinary” space-time relations are valid, is in compliance with generally accepted concepts about space-time scales in microcosm [68] .

One can try to estimate the value of fundamental mass, if one accepts that the least space unit (fundamental length) is the value afund ~ 10–18 m, and that this value is in compliance with the time scale. Then from relation (32) we obtain the following expression for fundamental mass:

Hence we obtain that [25] . This value fits into the expected (indicated in book [29] ) experimental range in which neutral Higgs H0-boson mass should fall in:. In addition this value is close to the “critical” mass 180 - 200 GeV above which H0-boson can disintegrate into couples of Wand Z-bosons [69] .The reason for this compliance is in the fact that the proposed model of “quantum” compact cluster formation is scalar, it does not include spin parameters or object charges, and hence it can be applied to the scalar neutral boson H0 with zero spin.

It should be interesting to perform, using relation (32), a qualitative analysis of a wide spectrum of processes with the participation of fundamental and elementary particles. For example, one can try to estimate the state of the heaviest of available quarks: t-quark (tq) having conditional mass 176 GeV [19] . Using Formula (32) the minimum time unit for t-quark, corresponding to the minimum length unit estimated in [40] as, is evaluated as. Then it follows from Formula (21) that for formation of a bound state with a characteristic size of about 10–15 m, the required time would be. The resultant estimate corresponds to the lifetime of unstable resonances; therefore t-quark cannot form stable hadrons. The result does not contradict to the known concepts according to which t-quark is the only quark that comes into being and dies free [19] .

Within the nuclear scales the proposed method allows one to estimate the time of the process of approaching equilibrium based on characteristic mass numbers and the corresponding nuclide sizes. Such task was solved in [70] to fit the study into kinetics of deep inelastic interaction between a beam of heavy copper ions and a target of gold with the collision energy 365 MeV [39] . It has been demonstrated that the times of formation of interaction products, having average mass numbers = 60 and 100 typical for the process under consideration, are in the order of magnitude 10–20 s. The resultant values are significantly lower than the lifetime 10–16 - 10–14 s of an intermediate compound nucleus. This allows us to speak about relatively fast (explosive) process progress. It should be noted that the established in [39] product distribution behavior in the reaction corresponds qualitatively to asymptotic properties of a distribution density function established in Sections 2 and 3. Thus for large mass numbers law (4) is approximately fulfilled, and against the background of a smooth curve local maximums take place at A = 56 and 63. Calculations based on Formula (28) show that the most probable mass numbers of cluster nuclides in the interregion between cluster radioactivity and spontaneous fission are equal to 51, 58, 65, 67.

In case of mesostructures the developed qualitative models include the constants, physical characteristics of the material phenomenological parameters a0, lcov, and operate with the minimum necessary number of variable in the experiment parameters (observed values): the environment temperature T, the grain initial size l, the time t of object growth process and their average size reached by the end of the process. It should be noted that the constant is omitted from final expressions for an average size of volume-packed mesoobjects. This factor is connected with a large number of elementary acts of object interactions, i.e., with quasi-classical nature of growth processes of nanostructured diamond particles and light actinides. And yet, the expressions for a typical diameter and an average length of a single-wall carbon nanotube include constant explicitly. This makes it possible to consider carbon nanotube as a purely quantum object.

The proposed in Ref. [15] [53] mechanism of macroscopic diamond particle formation from nanodiamonds covered all the available at that time data related to both artificial diamonds produced under static and dynamic synthesis, and natural diamonds that could have been formed as a result of cataclysm, e.g., impact of meteorite. In addition we will consider an issue about mesoscopic “limits” based on Formula (26) for the maximum size of a crystal cluster. As it is reasonable to take the residence time of seed at one place on the surface of a particle-cluster, i.e., characteristic oscillations period. One can see that the above limits significantly depend on the seed size. By the example of diamond we will consider the following types of seeds: carbon skeleton of cyclohexane with the characteristic size 0.25 nm, as well as nanodiamonds with the “critical” size at which quasi-long range crystal ordering originates (see Section 4). For diamond the value is 2.12 × 10–14 s. If we consider synthesis from molecules of carbon skeleton of cyclohexane, then the maximum size turns out to be amax = 100 nm. This value can be considered as the “theoretical” limit of detonation nanodiamonds. If a seed size is equal to 2 nm, then the maximum size at is evaluated as amax ≈ 3 cm.

Like a theory of fuzzy martensitic transitions developing according to diffusionless mechanism [10] , the model [15] [60] of the d-Pu®a-Pu transformation, as diffusion formation of crystal particles, shows that with the decreased grain size l of the initial phase, the temperature T0 of the martensitic transition beginning decreases, and with l < lmin martensitic transition is “blocked”. Yet, as l decreases in diffusion model the temperature T* of martensitic transition finish increases and “fuzzing of temperature transformation” decreases, which is determined in [25] as the difference T0 – T*. It is connected with phonon effect (in particular, “zero” crystal lattice oscillations). The established in section 4 agreement between the calculated and experimental data for plutonium provides grounds for using the model in describing behavior of other light actinides with strong approach of atoms [9] (Pa, U, Np). In addition the model [15] [60] predicts the possible reverse a-Pu®d-Pu phase transformation under pulsed heating of the sample preliminary cooled down to temperatures of about 100 K and the absence of the reverse pulsed transformation in the sample preliminary cooled down to temperatures of about 15 K.

With all variety of above processes one can single out the following common features of growth dynamics of compact objects having strong bonds and significantly different space scales:

Ÿ At the end of the aggregation process the function of distribution density of large clusters with sizes (a0 is seed size) is in inverse proportion to a3 (to cluster masses), and an average size of clusters is much lower than the maximum size amax:;

Ÿ There are many local maximums of function connected with the most probable cluster sizes, and there is a set of zeros of this function reflecting cluster instability at certain sizes;

Ÿ The obtained approximate laws of average size growth of large clusters in time with small and large seed fluxes describe aggregation processes both in microcosm and in mesoscopic scales.

Every power exponent Z with t obtained in this work are positive which is indicative of suppression of an inverse damage wave and formation of resultant direct wave of object growth even in case of cluster interaction between themselves. This result is in agreement with the concept of “persistence” in cluster-cluster aggregation [8] . By analogy with the theory of dispersive waves [71] one can try to introduce the concept of “phase” and “group” velocities in dimensionless space. Phase velocity of growth process is infinitely large:. This means extremely fast spreading of any initial cluster distribution and “instant” formation of density wave in space. Group velocity of density wave in the majority of growth mechanisms decreases with time:

This reflects broadening of wave packet and slowing down of its propagation. At the final stage of the process of irreversible cluster aggregation in closed system such wave packet of a non-Gaussian shape has a decaying sudden change in a wave “front”, associated with the maximum possible size:

These considerations bring to a conclusion that the developed asymptotic method for investigation into kinetics of formation of compact objects having strong internal bonds has a certain degree of generality to use it in solving problems in physics of high density energy and high-intensive processes.

In particular, in [25] the consideration is given to a problem of nuclide “stability islands” [29] [43] , and the fact of such islands existence is explained in terms of formation of clusters with the most probable sizes from nuclear matter seeds. As an example the element roentgenium with Z = 111 has the mass number of A = 273. In the proposed model, this value of A obtained from Formula (28) corresponds to tritons as the seeds. At Z = 114, in compliance with 2β-stability, the number of neutrons is equal to 184 and the mass number is equal to A = 298 [43] . In the model [25] , the mass numbers closed to it are 299 and 302. These values obtained from the Formula (28) correspond to tritons and deuterons as the seeds, respectively. Work [72] , which is dedicated to the production and the decay of element 114, mentions its isotopes with A = 288 and 289. In the model [25] , the values A = 290 and 293 are the closest mass numbers corresponded to the tritons and deuterons as the seeds, respectively. Using various calculation options of the β-stability band the islands of stability are predicted also at Z = 164 and with neutron numbers of 272 or 318 [43] . Respective mass numbers A = 436 and 482 are close to the approximate value of the mass number of the final nuclide Aend » 470 calculated in [14] . In the model [25] , mass number A =436 corresponds to deuterons as the seeds. The closest to the mass number of 482 are the following values calculated with the help of Formula (28): A = 483 for the tritons as the seeds, and A = 485 for deuterons as the seeds. Thus, one can presume that the developed asymptotic model of clusters formation in the nuclear matter complements the method for predicting nuclear mass with the help of radial basis function [73] , which makes it possible to find mass numbers in the range from 20 to 260.

As for high-intensive nuclear processes, when a target is affected by fast, including relativistic, particles, then to describe them one should use Equation (25). From this equation as well as from Formula (27) one can derive the following expression for the time of formation of cluster nuclides with an average mass number:

Here, mnucl is nucleon mass, which can be accepted as an average value between proton and neutron masses, Eexc is nucleus excitation energy per a nucleon. If we accept the average kinetic nucleon energy in degenerate ideal Fermi-gas as Eexc, then from the previous equation we will have:

Here, c0 is average thermal nucleon velocity (see Section 3). In such “acoustic” (ac) approximation the time of light nuclide formation with ≈ 10 is equal to 2 × 10–22 s, for heavy nuclides with ≈ 100 this time is 5 × 10–22 s. If one takes the maximum possible energy mnuclc2 as Eexc, the time of nuclide cluster formation is expressed as:

In such “ultrarelativistic” (ur) approximation the times of nuclide formation are 2 × 10–23 s and 4 × 10–23 s. The obtained estimates do not contradict to the generally known notions about typical time of direct nuclear reactions [29] [35] [36] .

Another example of using the developed method is explanation and description of diamond synthesis in conditions of cavitation experimentally studied in [74] [75] . In [76] a solid stage of the process is considered, at this stage diamond particle growth in extreme condition occurs. Asymptotic models of nanostructured mesoobjects with strong interatomic bonds make the basis of this consideration [15] . These models include description of crystal growth models at high temperatures T, when. For obtaining laws of diamond particle growth in conditions of cavitation Equation (25) was used. It has been shown that nanodiamonds are produced at the first stage of synthesis. Hence production of diamond particles in a solid phase occurs as a result of nanodiamond coalescence and consolidation of the produced objects into larger particles both at the final stage of cavitation bubble collapse, and in the process of collapse “products” expansion as well. 

The next example is the problem of determining the width of detonation zone in solid high explosives (HE) [77] . In this case the subject to be investigated is not formation of an object with strong internal bonds, but failure of strong covalent bonds between carbon atoms, resulting in relaxation processes of new bonds formation and energy release. The detonation zone is considered as quantum mechanical system, being in an excited state. For example, it can be excitation of quasi-particles, phonons, caused by shock wave. As a result of system phonon excitations, there occur atom oscillations with relatively high amplitude with which bonds between atoms break and new bonds are formed that are typical for decomposition products of the original substance: various chemical compounds (molecules, radicals), carbon clusters, and others.

Relaxation processes in a detonation are considered using the concept of distribution density wave j(x, t) of objects with new interatomic bonds, formed when former bonds broke as a result of “excitation” of oscillations. The wave propagates with time t within a quantum mechanical system along the axis of HE charge (x-axis) from the shock wave front towards the surface, where chemical reactions of energy release finish. The motion of the wave within the detonation zone goes up to the time t = t, wheret is duration time of detonation-related chemical reactions. The estimates obtained for the average (in terms of its charge width section) detonation zone fit into the range of generally known data for typical solid HE, and indicate that it is possible to produce explosives with nanometer energy release regions [77] .

One can try to apply the proposed asymptotic method to estimate grain sizes formed as a result of recrystallization process in aluminum alloys at an impact of a group of solids with a rate of 1200 - 1300 m/s [78] . Numerical calculations have shown that the material temperature of crater walls formed, when solid impactors penetrate aluminum alloy targets, is T = 500 K, and a full time of impactor penetration up to its stopping is t = 40 µs. One can accept this calculated time as an estimate of life duration of intensive phonon excitations in the system “impactor-target material”. As in a previous example, the process of grain formation and growth can be represented as a process of coalescence and growth of nanocrystals formed under failure of target original material structure as a result of solid shock action and stress wave interference generated by neighboring impactors.

Since the calculated temperature is comparable with Debye aluminum parameter θD = 443 K [79] , then with small flux of nanocrystal seeds the time dependence of an average size of crystal particles (grains) is described by Formula (45). From here, we get that at t = 40 µs, the grain average size is 0.4 µm. This value can be accepted as a grain size estimate at primary recrystallization. This estimate corresponds to grain sizes at a crater bottom specified in [78] using scanning electronic microscope.

In general, the proposed phenomenological approach seems to be useful when considering the following problems from the list given in review [68] : mesoscopics, clusters, nanotubes, super-heavy elements, mass spectrum, quarks and gluons, fundamental length, particle interaction at high energies, neutron stars, and supernovas.

It should be noted that the estimated in [25] fundamental mass 196 GeV does not differ strongly from the mass 125 - 126 GeV [80] of the new particle with Higgs boson properties discovered on the Big Hadron collider in CERN. It could not be ruled out that the above new particle is a product of disintegration of another heavier particle. It is also important to mark approximate accordance of obtained evaluation for to the value of the upper limit of the Higgs mass identified in [81] as 170 GeV. These factors argue in favor of validity of relation (32) in terms of fundamental processes in microcosm.

Hence the developed models can be useful both for theoretical studies and experiments as well as for primary analysis of their first data.

6. Conclusions

The proposed phenomenological approach allows one to consider formation and growth of subatomic and mesoscopic objects under single formalism.

Based on linearization and qualitative analysis of solutions of Fokker-Planck kinetic equation the basic properties of a size-distribution density function of volume-packed bodies were established and approximate relations between average object sizes and their growth times were obtained.

An approximate differential relation for determining object growth in high-intensive processes with high energy was obtained.

The proposed asymptotic model provides adequate estimates of mass characteristics of typical intranuclear processes.

The calculated data obtained for formation of nanostructured mesoobjects having strong interatomic bonds agree well with generally known experimental data.

References

  1. Morokhov, I.D., Petinov, V.P., Turusov, L.P. and Petrunin, V.F. (1981) Structure and Properties of Fine Metallic Particles. Soviet Physics Uspekhi, 24, 295-334. http://dx.doi.org/10.1070/PU1981v024n04ABEH004800
  2. Kadomtsev, B.B. (1994) Dynamics and Information. Soviet Physics Uspekhi, 37, 425-499. http://dx.doi.org/10.1070/PU1994v037n05ABEH000109
  3. Von Oppen, G. (1994) Objects and Environment. Physics-Uspekhi, 39, 617-620.http://dx.doi.org/10.1070/PU1996v039n06ABEH000152
  4. Smirnov, B.M. (2011) Processes Involving Clusters and Small Particles in Buffer Gas. Physics-Uspekhi, 54, 691-721.
  5. Pool, C.P. and Owens, F.J. (2003) Introduction to Nanotechnology. John Wiley & Sons Inc., London.
  6. Imry, Y. (2002) Introduction to Mesoscopic Physics. Oxford University Press, Oxford.
  7. Bording, J.K. and Taft, J. (2000) Molecular-Dynamics Simulation of Growth of Nanocrystals in Amorphous Matrix. Physical Review B, 62, 8098-8103. http://dx.doi.org/10.1103/PhysRevB.62.8098
  8. Hellen, E.K. and Alava, M.J. (2002) Persistence in Cluster-Cluster Aggregation. Physical Review E, 66, 026120.
  9. Hecker, S.S. (2000) Plutonium and Its Alloys: From Atoms to Microstructures. In: Cooper, N.G., Buican, I.G. and Schecker, J.A., Eds., Challenges in Plutonium Science, Los Alamos Science, Los Alamos, 290-334.
  10. Malygin, G.A. (2008) Nanoscope Size Effects on Martensitic Transformation in Shape Memory Alloys. Physics of the Solid State, 50, 1538-1543. http://dx.doi.org/10.1134/S1063783408080258
  11. Lin, E.E. (1993) Kinetics of the Formation of Compact Quantum Clusters in Stochastic Conservative Systems. Technical Physics Letters, 19, 165-166.
  12. Suzdalev, I.P. (2006) Nanotechnology: Physics-Chemistry of Nanoparticles, Nanostructures and Nanomaterials. Synergetics: From the Past to Future, Komkniga, Moscow. (in Russian)
  13. Vanossi, A., Manini, N., Urbakh, M., Zapetti, S. and Tosatti, E. (2013) Modeling Friction: From Nanoscale to Mesoscale. Review of Modern Physics, 85, 529-562. http://dx.doi.org/10.1103/RevModPhys.85.529
  14. Lin, E.E. (2012) Calculating the Mass Numbers of Cluster Nuclides on the Basis of the Asymptotic Model. Bulletin of the Russian Academy of Sciences: Physics, 76, 881-883. http://dx.doi.org/10.3103/S1062873812080199
  15. Lin, E.E. (2010) Mesokinetics of Growth of Nanostructured Objects with Strong Interatomic Bonds. In: Taylor, J., Ed., Advances in Chemistry Research, Nova Science Publishers, New York, 5, 171-190.
  16. Fermi, E. (1960) Notes on Quantum Mechanics. The University of Chicago Press, Chicago.
  17. Man’ko, V.I. (1990) Coherent State. In: Prokhorov, A.M., Ed., Physical Encyclopedia, Soviet Encyclopedia, Moscow, 392-394. (in Russian)
  18. Wagner, P. and Zhong, Y.M. (1995) Cluster Formation in Disordered Systems and Nuclear Fragmentation. Nuclear Physics A, 592, 385-412.
  19. Okun, L.B. (1998) Current Status of Elementary Particle Physics. Physics-Uspekhi, 41, 553-557.http://dx.doi.org/10.1070/PU1998v041n06ABEH000403
  20. Dremin, I.M. (2010) The Quark-Gluon Medium. Physics-Uspekhi, 53, 123-1149.
  21. Lozovik, Y.E. and Popov, A.M. (1997) Formation and Growth of Carbon Nanostructures: Fullerens, Nanoparticles, Nanotubes and Cones. Physics Uspekhi, 40, 717-737. http://dx.doi.org/10.1070/PU1997v040n07ABEH000253
  22. Mukhin, K.N. and Patarakin, O.O. (2000) Exotic Processes in Nuclear Physics. Physics-Uspekhi, 43, 799-839. http://dx.doi.org/10.1070/PU2000v043n08ABEH000802
  23. Lin, E.E. (2010) Asymptotical Model of the Formation of Nuclear Matter Clusters. In: Il’kaev, R.I., Ed., Proceedings of International Conference 12th Khariton Thematic Scientific Readings, Problems of High Density Energy Physics, 19-23 April 2010, Sarov, 262-267. (in Russian)
  24. Panov, I.V. and Thielemann, F.K. (2004) Fission and the r-Process: Competition between Neutron-Induced and Beta-Delayed Fission. Astronomic Letters, 30, 647-655. http://dx.doi.org/10.1134/1.1795953
  25. Lin, E.E. (2011) Qualitative Kinetic Models of Formation of Compact Objects with Strong Internal Bonds. Monograph, Russian Federal Nuclear Center-VNIIEF, Sarov. (in Russian)
  26. Pugachev, V.C. and Sinitsyn, I.N. (1990) Stochastic Differential Systems: Analysis and Filtration. Nauka, Moscow. (in Russian)
  27. Reissland, J.A. (1973) The Physics of Phonons. John Wiley and Sons LTD., London; New York; Sydney; Toronto.
  28. Kulakov, V.M. (1991) Nuclear Properties of Nuclides. In: Grigor’ev, S. and Meilikhov, E.Z., Eds., Physical Quantities Handbook, Energoatomizdat, Moscow, 993-1054. (in Russian)
  29. Ishanov, B.S., Kapitonov, I.M. and Yudin, N.P. (2007) Particles and Atomic Nuclei. LKI, Moscow. (in Russian)
  30. Kadmenskii, S.G., Kurgalin, S.D. and Tchuvil’sky, Y.M. (2012) Cluster States in Atomic Nuclei and Cluser-Decay Processes. Physics of Particles and Nuclei, 38, 699-742.
  31. Zamyatnin, Y.S., Kadmensky, S.S., Kurgalin, S.D. and Tchuvil’sky, Y.M. (1994) Where the New Examples of Nuclei Cluster Decay Can Be Seek? Physics of Atomic Nuclei, 57, 1905-1918.
  32. Vandenbosch, R. and Huisenga, J.R. (1973) Nuclear Fission. Academic Press, New York; London.
  33. Obukhov, A.I. and Grigor’ev, I.S. (1991) Fissionof Nuclei. In: Grigor’ev, I.S. and Meilikhov, E.Z., Eds., Physical Quantities Handbook, Energoatomizdat, Moscow, 1087-1098. (in Russian)
  34. Panov, I.V., Korneev, I.Y. and Thielemann, F.K. (2009) Superheavy Elements and r-Process. Physics of Atomic Nuclei, 72, 1026-1033. http://dx.doi.org/10.1134/S1063778809060155
  35. Olkhovsky, V.S. (1984) On Investigation of Nuclear Reactions and Decays with the Help of Analysis of Theirs Durations. Soviet Journal of Particles and Nuclei, 15, 130-148.
  36. Olkhovsky, V.S., Recami, E. and Maydanuyk, S.P. (2012) Time as Quantum Observable, Canonical Conjugated to Energy. In: Pahlavani, M.R., Ed., Measurements in Quantum Mechanics, In Tech, Shaghai, 18-56.
  37. Segre, E. (1964) Nuclei and Particles, an Introduction to Nuclear and Subnuclear Physics. W. A. Benjamin Inc., New York.
  38. Segre, E. (1977) Nuclei and Particles. 2nd Edition, W. A. Benjamin Inc., New York.
  39. Valentin, L. (1982) Physique Subatomique: Noyaux Et Particles. Nouvelle Edition Entirement Refondue, Hermann, Paris.
  40. Jacob, M. and Landshoff, P. (1982) The Inner Structure of Proton. Scientific American, 242, 66-75. http://dx.doi.org/10.1038/scientificamerican0380-66
  41. Kane, G. (1987) Modern Elementary Particle Physics. University of Michigan, Addison-Wesley Publishing Company Inc., Michigan.
  42. Neufeld, R.B. and Mueller, B. (2009) Sound Produced by a Fast Parton in a Guark-Gluon Plasma Is a “Crescendo”. Physical Review Letters, 103, 043201.
  43. Selinov, I.P. (1990) Atomic Nuclei: Structure and Systematic. Nauka, Moscow. (in Russian)
  44. Nadezin, D.K. (1992) Neutron Stars. In: Prokhorov, A.M., Ed., Physical Encyclopedia, Big Russian Encyclopedia, Moscow, 280-283. (in Russian)
  45. Utrobin, V.P. (1994) Supernovae. In: Prokhorov, A.M., Ed., Physical Encyclopedia, Big Russian Encyclopedia, Moscow, Vol. 4, 433-435. (in Russian)
  46. Ponomarev, L.I. (1992) Uncertainty Relations. In: Prokhorov, A.M., Ed., Physical Encyclopedia, Big Russian Encyclopedia, Moscow, 321-322. (in Russian)
  47. Vereshchagin, A.L., Sakovich, G.V., Brylyakov, P.M., Zolotukhina, I.I., Petrova, L.A. and Novoselov, N.N. (1990) The Structure of Carbon Diamond-Like Phase of Detonated Synthesis. Soviet Physics Doklady, 35, 851-852.
  48. Lin, E.E. (2003) On the Growth Kinetics of a Single-Wall Carbon Nanotube. Doklady Physics, 48, 180-181. http://dx.doi.org/10.1134/1.1574375
  49. Badiali, J.P. (1999) Fractal Behavior in Quantum Statistical Physics. Physical Review E, 60, 2533-2539. http://dx.doi.org/10.1103/PhysRevE.60.2533
  50. Vollmer, M. (1983) Kinetics of New Phase Formation. Plenum, New York. (Kinetic der Phasenbildung. T. Steinkopf, Dresden, 1939)
  51. Syue-sen’, T. (1965) Physical Mechanics. Mir, Moscow. (Translated from Chinese)
  52. Yakubov, T.S. (1990) On Specific Heat of Solids Revealing Fractal Character. Doklady Academii Nauk SSSR, 310, 145-149.
  53. Lin, E.E. (2000) On the Cluster Mechanism of Diamond Synthesis from Different Solid Carbon Forms. Physics of the Solid State, 42, 1946-1951. http://dx.doi.org/10.1134/1.1318892
  54. Lin, E.E. (1994) Aggregation of Crystalline Clusters in the Shock Wave Front Propagating through Condensed Species. Soviet Journal of Chemical Physics, 12, 404-409.
  55. Lin, E.E. (1997) Shock-Induced Growth of Crystals in a Porous Medium of Diamond Nanoparticles. Chemical Physics Reports, 16, 2241-2244.
  56. Lin, E.E. (2000) Shock-Induced Coalescence of Nanodiamonds. Bulletin of the Russian Academy of Sciences: Physics, 64, 1215-1216.
  57. Lin, E.E. (2005) Pulsed Loading of Objects during Intense Expansion of Products of Solid Explosives (Review). Combustion, Explosion, and Shock Waves, 41, 241-263. http://dx.doi.org/10.1007/s10573-005-0030-2
  58. Lin, E.E. (2009) Determining the Growth Rate of Nanostructured Particles of Light Actinide at High and Moderate Temperatures. Technical Physics Letters, 35, 418-420. http://dx.doi.org/10.1134/S1063785009050101
  59. Van Bueren, H.G. (1960) Imperfections in Crystals. North-Holland Publishing Company, Amsterdam.
  60. Lin, E.E. (2010) Influence of Size Factors on the Character of Size Transformations in Light Actinides. Physics of the Solid State, 52, 153-157. http://dx.doi.org/10.1134/S1063783410010270
  61. Janke, E., Emde, F. and Losch, F. (1960) Tafeln Hoherer Funktionen. (Neubearbeitet von F. Losch), B. G. Verlagsgesellschaft, Stuttgart.
  62. Louson, E.S., Martinez, B., Roberts, J.A., Richardson, J.U. and Benne, B.I. (2000) Oscillations of Atoms and Plutonium Melting. In: Cooper, N.G., Buican, I.G. and Schecker, J.A., Eds., Challenges in Plutonium Science, Los Alamos Science, Los Alamos, 190-199.
  63. Dwight, H.B. (1961) Tables of Integrals and other Mathematical Data. The Macmillan Company, New York.
  64. Hecker, S.S., Harbur, D.R. and Zocco, T.G. (2004) Phase Stability and Phase Transformations in Pu-Ga Alloys. Progress in Materials Science, 49, 429-485. http://dx.doi.org/10.1016/S0079-6425(03)00032-X
  65. Kitching, S., Planterose, P.G. and Gill, D.C. (2003) Stabilized Plutonium. In: Jarvinen, G.D., Ed., Plutonium Futures— The Science, American Institute of Physics, New York, 79-81.
  66. Mitchell, J.N., Hecker, S.S., Freibert, F.J., Schwarts, D.S. and Bange, M.E. (2008) Inconventional Delta-Phase Stabilization in Pu-Ga Alloys. In: Fundamental Plutonium Properties, Proceedings of 8th International Workshop, 8-12 September 2008, Snezhinsk, 5-6.
  67. Leshli, D.S., Blau, M.S. and Moment, R.L. (2000) Manufacturing of Single Crystals of Plutonium. In: Cooper, N.G., Buican, I.G. and Schecker, J.A., Eds., Challenges in Plutonium Science, Los Alamos Science, Los Alamos, 233-245.
  68. Ginzburg, V.L. (1999) What Problems of Physics and Astrophysics Seem Now to Be Especially Important and Interesting (Thirty Years Later, Already on the Verge of XXI Century)? Physics-Uspekhi, 42, 353-373.http://dx.doi.org/10.3367/UFNr.0169.199904d.0419
  69. Okun, L.B. (1985) Particle Physics: The Guest for the Substance of Substance. Hardwood Academic Publishers, Chur.
  70. Lin, E.E. (1993) Kinetics of Deep-Inelastic Heavy Ions Interaction. Technical Physics Letters, 19, 669-670.
  71. Whitham, G.B. (1974) Linear and Nonlinear Waves. A Wiley-Interscience Publication, Hoboken.
  72. Düllmann, C.E., Schadell, M. and Yakushev, A. (2010) Production and Decay of Element 114: High Cross Sections and New Nucleus 277Hs. Physical Review Letters, 104, 252701-252705. http://dx.doi.org/10.1103/PhysRevLett.104.252701
  73. Wang, N. and Liu, M. (2011) Nuclear Mass Predictions with a Radial Basis Function Approach. Physical Review C, 84, 051303.
  74. Galimov, E.M., Kudin, A.M., Skorobogatskii, V.N., Plotnichenko, V.G., Bondarev, O.L., Zarubin, B.G., Strazdovskii, V.V., Aronin, A.S., Fisenko, A.V. and Bubov, I.V. (2004) Experimental Confirmation of the Diamond Synthesis during Cavitation. Doklady Physics, 395, 187-191.
  75. Khachatryan, A., Aloyan, S.S., May, P.W., Sargsyan, R., Khachatryan, V.A. and Bagdasaryan, V.S. (2008) Grafite-to-Diamond Transformation Induced by Ultrasound Cavitation. Diamond and Related Materials, 17, 931-936. http://dx.doi.org/10.1016/j.diamond.2008.01.112
  76. Lin, E.E. (2011) Cluster Mechanism of Diamond Synthesis under Severe Cavitation Conditions. Technical Physics Letters, 37, 593-595. http://dx.doi.org/10.1134/S1063785011070121
  77. Lin, E.E. (2011) Detonation Zone Width Determination Based on Uncertainty Principle. Technical Physics Letters, 37, 449-450. http://dx.doi.org/10.1134/S1063785011050257
  78. Auchev, A.A., Ilyushkina, N.Y., Lin, E.E., Popov, N.N. and Tanakov, Z.V. (2013) Recrystallization in Aluminum Alloyes under Impact of Solids. In: Proceedings of International Conference 15th Khariton Topical Scientific Readings “Extreme States of Substance. Detonation. Shock Waves”, 18-22 March 2013, Sarov, 300-301. (in Russian)
  79. Hlopkin, M.N. (1991) Heat Capacity. In: Grigor’ev, S. and Meilikhov, E.Z., Eds., Physical Quantities Handbook, Energoatomizdat, Moscow, 205. (in Russian)
  80. Rubakov, V.A. (2012) Large Hadron Collider’s Discovery of a New Particle with Higgs Boson Properties. Physics-Uspekhi, 55, 949-957.
  81. Elsayed, A., Khalil, S. and Moretti, S. (2012) Higgs Mass Corrections in the SUSY B-L Model with Inverse Seesaw. Physics Letters B, 715, 208-213. http://dx.doi.org/10.1016/j.physletb.2012.07.066