Stationary Fokker-Planck equation applied to fission dynamics

By use of both analytical and numerical techniques, we study the relaxation of time-dependent solutions of the Fokker-Planck equation for an inverted oscillator to Kramers' stationary solution. This is done by integrating over all time the time-dependent solutions for given initial conditions at the saddle point to obtain stationary solutions, whose densities and higher velocity moments are compared as functions of the coordinate with the corresponding quantities calculated from Kramers' stationary solution. For large values of the coordinate an a symptotic expansion of the density is obtained, but for general values of the coordinate and for higher velocity moments the time integration must be done numerically. With increasing dissipation the relaxation to Kramers' stationary solution occurs at successively smaller values of the coordinate. By use of Kramers' stationary solution, we derive analytical expressions as functions of nuclear temperature and dissipation strength for several quantities of interest in fission dynamics, including the mean time from the saddle point to scission, the mean fission-fragment kinetic energy at the scission point and the contribution to the variance in the fission-fragment kinetic energy resulting from fluctuations in the fission degree of freedom. We apply these results to some examples that have been studied experimentally, including the mean saddle-to-scission time for the heavy-ion-induced fission of the compound nucleus ¹⁶⁸Yb and the mean fission-fragment kinetic energy at scission and the contribution to its variance for the α-particle-induced fission of the compound nucleus ²¹³At.     

曲率能对热力学驱动力的影响

为了研究曲率能对核裂变热力学驱动力（TDF）的影响,首先利用包含曲率能的截断版小液滴模型计算了200Pb和224Th的位垒和熵垒,对比液滴模型的计算结果表明:曲率能未改变224Th的位垒鞍点,却将200Pb的位垒鞍点向后推移。能级密度参数的形变关系越强则两系统的熵垒鞍点越靠近基态。为了进一步探究曲率能如何通过位势和熵势影响TDF,以断前中子多重性（PNM）为探针,通过两种方案进行了模拟,结果表明:曲率能降低了两系统的位势驱动力,而增强了其熵势驱动力。结合PNM的计算表明,前一种效应要比后一种效应明显,因此曲率能总体减弱了200Pb和224Th的TDF,进而延缓了两系统的核裂变进程。In order to study the effect of curvature energy on the thermodynamic driving force (TDF) of nuclear fission, the potential and entropy barrier of 200Pb and 224Th systems are calculated by using the truncated droplet model including curvature energy, respectively. Compared with the liquid drop model, the results show that curvature energy does not affect the saddle point of 224Th, but pushes the saddle point of 200Pb backwards the ground state. The stronger the deformation dependence of the level density parameter is, the closer the saddle point of entropy barrier for these systems is to the ground state. In order to further investigate how curvature energy affects TDF through nuclear potential and entropy, respectively, the prescission neutron multiplicity (PNM) is selected as the probe, some simulations based on two schemes are carried out. The results show that curvature energy reduces the potential driving force of 200Pb and 224Th, and enhances the entropy potential driving force. Combined with the calculations and analyses of PNM, the former effect is more obvious than the latter, so curvature energy weakens TDF of two systems on whole, thus delaying the nuclear fission process of two systems.     

多维Kramers公式的研究

本文研究了鞍点附近多维的位能曲面及鞍点所在位置;并用Werner Wheeler及无旋液体等两种方法计算了多维的质量系数与粘滞系数,然后用多维Kramers公式计算了裂变速率.发现裂变速率随着维数的增加而适当增大.不同的形变参量以及不同的计算质量和粘滞系数方法对计算核裂变速率影响不大.从结果看,采用三维计算裂变几率已足够准确.     

Fokker-Planck方程的数值解

采用扩散模型研究核裂变,需要求解Fokker-Planck方程。本文提出一个数值计算方法-平均隐式差分方法。对具有粘滞性的核体系的有关裂变动力学量,如几率分布、裂变率、断点处的平均动能以及鞍点到断点的平均扩散时间等一系列物理量做了计算,并与适合大粘滞性的Kramers的解析解做了比较。通过与解析解的比较及对归一常数的检验,证明计算结果精确可靠。     

Analysis of fission excitation functions and the determination of shell effects at the saddle point

A method is proposed to deduce the shell correction energy corresponding to the fission transition state shape of nuclei in the mass region around 200, from an analysis of the first chance fission values of the ratio of fission to neutron widths, (Γ _f /Γ _n )₁. The method is applied to the typical case of the fissioning nucleus²¹²Po, formed by alpha bombardment of²⁰⁸Pb. For the calculation of the neutron width, the level densities of the daughter nucleus after neutron emission were obtained from a numerical calculation starting from shell model single particle energy level scheme. It is shown that with the use of standard Fermi gas expression for the level densities of the fission transition state nucleus in the calculation of the fission width, an apparent energy dependence of the fission barrier height is required to fit the experimental data. This energy dependence, which arises from the excitation energy dependence of shell effects on level densities, can be used to deduce the shell correction energy at the fission transition state point. It is found that in the case of²¹²Po, the energy of the actual transition state point is higher than the energy of the liquid drop model (LDM) saddle point by (3 ± 1) MeV, implying significant positive shell correction energy at the fission transition state. Further, the liquid drop model value of level density parametera is found to be a few per cent smaller for the saddle point shape as compared to its spherical shape.     

Dynamical effects on the way to fusion of very heavy nuclei

The central collision of ⁴⁰Ar and ²⁰⁸Pb is studied considering the ellipsoidal deformations and isovector dipole mode of motion in the approaching phase. The collective energy dissipation is suggested to originate from the Fermi surface deformation which is treated as a kinematically independent mode of motion within the canonical Lagrange-Rayleigh dynamics. The possible extensions of the approach are discussed. The potential energy surface, calculated using the generalized (folded) surface potential, is studied. The saddle point in the potential energy surface lying at the border of strongly deformed compact configurations is located. The potential energy at this point is about 10MeV smaller than that of the ions touching each other in the spherical shape. The examination of trajectories followed by the system in its evolution shows that the inertia forces strongly hinder the motion of ions along the potential energy valley. The collective energy dissipated during the approach is found to be smaller than the difference in the potential energies at saddle point and at the touching configuration of unpolarized ions.     

Fission Barrier for ^240Pu in the Quadrupole Constrained Relativistic Mean Field Approach

The fission barrier for ^240pu is investigated beyond the second saddle point in the potentiM energy surface by the constrained relativistic mean field method with the newly proposed parameter set PK1. The microscopic correction for the centre-of-mass motion is essentiM to provide the correct potential energy surface. The shell effects that stabilize the nuclei against the fission is also investigated by the Strutinsky method. The shapes for the ground state, fission isomer and saddle-points, etc, are studied in detail.     

d-Mott phases in one and two dimensions

We use exact diagonalization to determine the spectrum of reduced Hamiltonians based on renormalization group flows to strong coupling. For the half-filled two-leg Hubbard ladder we reproduce the known insulating d-Mott ground state with spin and charge gaps. For the saddle point regions of the two-dimensional Hubbard model near half filling we find a crossover to a similar strong coupling state, which truncates the Fermi surface near the saddle points. At lower scales d-wave superconductivity appears on the remaining Fermi surface.     

Two-Dimensional Saddle Point Equation of Ginzburg-Landau Hamiltonian for the Diluted Ising Model

The saddle point equation of Ginzburg-Landau Hamiltonian for the diluted Ising model is developed. The ground state is solved numerically in two dimensions. The result is partly explained by the coarse-grained approximation.     

Coloured noise controlled dynamics of nonlinear polaritons in semiconductor microcavity

The coloured noise induced escape rate from the lower energy stable state of a driven nonlinear microcavity oscillator has been investigated by means of quasi-classical kinetic equations. We show that for coloured, i.e. narrow-band, relatively intense noise, the escape time is controlled by the interplay of two mechanisms: the noise induced drift and adiabatic regular shift of the oscillator state towards unstable saddle point. The cross-over between these mechanisms takes place in a particular range of the driving field intensity values, depending on the ratio between the oscillator damping and the coloured noise spectrum width. The dependence of the transition rate on the noise correlation time is analyzed for wide range of correlation time values. The article is published in the original.     

Dynamics of the universe with disformal coupling between the dark sectors

We use a dynamical analysis to study the evolution of the universe at late time for the model in which the interaction between dark energy and dark matter is inspired by a disformal transformation. We extend the analysis in the existing literature by assuming that the disformal coefficient depends both on the scalar field and its kinetic terms. We find that the dependence of the disformal coefficient on the kinetic term of scalar field leads to two classes of the scaling fixed points that can describe the acceleration of the universe at late time. The first class exists only for the case where the disformal coefficient depends on the kinetic terms. The fixed points in this class are saddle points unless the slope of the conformal coefficient is sufficiently large. The second class can be viewed as the generalization of the fixed points studied in the literature. According to the stability analysis of these fixed points, we find that the stable fixed point can take two different physically relevant values for the same value of the parameters of the model. These different values of the fixed points can be reached for different initial conditions for the equation of state parameter of dark energy. We also discuss the situations in which this feature disappears.     

Study of Fission Barrier Heights of Uranium Isotopes by the Macroscopic-Microscopic Method

Potential energy surfaces of uranium nuclei in the range of mass numbers 229 through 244 are investigated in the framework of the macroscopic-microscopic model and the heights of static fission barriers are obtained in terms of a double-humped structure. The macroscopic part of the nuclear energy is calculated according to Lublin-Strasbourg-drop (LSD) model. Shell and pairing corrections as the microscopic part are calculated with a folded-Yukawa single-particle potential. The calculation is carried out in a five-dimensional parameter space of the generalized Lawrence shapes. In order to extract saddle points on the potential energy surface, a new algorithm which can effectively find an optimal fission path leading from the ground state to the scission point is developed. The comparison of our results with available experimental data and others' theoretical results confirms the reliability of our calculations.     

The influence of the anion vibrational temperature on the fs dynamics in a NeNePo experiment

A negative-neutral-positive (NeNePo) femtosecond charge-reversal experiment studying the temperature-dependent relaxation dynamics of linear Ag₃ is described. A wavepacket is prepared on the potential-energy surface of the electronic ground state of the neutral trimer by photodetachment from an anion ensemble held at different, well-defined temperatures between 20 K and 350 K. The wavepacket dynamics is probed by resonant two-photon ionization. The cooled octopole ion trap developed to prepare the cold anions is described. The relaxation dynamics of the initially linear Ag₃ from a saddle point of the potential-energy surface into the triangular equilibrium configuration shows a significant dependence on the anion temperature, which is rationalized in terms of a simple model. We demonstrate that a low anion vibrational temperature results in the generation of “narrower” wavepackets on the neutral potential-energy surface. This allows us to probe coherent effects more sensitively. Received: 10 January 2000 / Published online: 24 July 2000     

Leading large N modification of QCD2 on a cylinder by dynamical fermions

We consider two-dimensional QCD on a cylinder, where space is a circle. We find the ground state of the system in case of massless quarks in a 1/Nexpansion. We find that coupling to fermions nontrivially modifies the large N saddle point of the gauge theory due to the phenomenon of “decompactification” of eigenvalues of the gauge field. We calculate the vacuum energy and the vacuum expectation value of the Wilson loop operator both of which show a nontrivial dependence on the number of quarks flavours at the leading order in 1/N.     

Non-trivial saddle points and band structure of bound states of the two-dimensional O(N) vector model

We discuss O(N) invariant scalar field theories in 0 + 1 space-time dimensions (quantum mechanics) and in 1 + 1 space-time dimensions (field theory). Combining ordinary “Large N” saddle point techniques and simple properties of the diagonal resolvent of one-dimensional Schrödinger operators we find non-trivial (non-constant) solutions to the saddle point equations of these models in addition to the saddle point describing the ground state of the theory. In the “Large N” limit these saddle points are exact for the quantum mechanical case, but only approximate in the two-dimensional theory. In the latter case they are the leading contributions to the time evolution kernel at short times, or equivalently, the leading contribution to the high temperature expansion of partition function stemming from space dependent static configurations in case of the Euclidean theory. We interpret these novel saddle points as collective O(N) singlet excitations of the field theory, each embracing a host of finer quantum states arranged in O(N) multiplets, in an analogous manner to the band structure of molecular spectra.     

The Rubik’s cube problem revisited: a statistical thermodynamic approach

Inspired by the protein folding problem, we propose a Rubik’s cube model and study its thermodynamic and kinetic behavior. We find that the energy landscape contains a tiny funnel-like region, as the dynamics towards the native state is mostly diffusive. In particular, from Monte Carlo simulations we observe exponential kinetics in the first-passage-time distribution towards the native state at all temperatures considered, while the complexity of the energy landscape is exhibited through a stretched-exponential relaxation of the energy autocorrelation function. The rollover feature in the mean first passage time, as observed in many protein-folding dynamics studies, is captured again in our model and discussed under the statistical energy landscape approach.     

LEVEL DENSITY AND FINITE-TEMPERATURE SPECIFIC HEAT OF NUCLEUS 104Pd UNDER MICROSCOPIC IBM

By using the microscopic sdgIBM-F_max approach, procedure of canonical ensemble average and saddle point approximation, thermodynamics of nucleus is established under microscopic IBM. Calculations of spectrum, level density and finite-temperature specific heat for the nucleus ¹⁰⁴Pd are carried out. The calculated values are coincident with the experimental data reported recently. The results predict that the shape phase transition in ground state band appears at about T≈0.230MeV and the phase transition of thermal excitation mode takes place at T≈0.630MeV for nucleus ¹⁰⁴Pd.     

POX(X=1,2)的光谱常数与从头算势能曲线

采用从头算的多种方法对PO、PO₂自由基的基态结构进行优化计算，结果表明：使用密度泛函（DFT）方法计算的结果最接近实验值.对PO双原子分子优选6-311G（3df）基组进行计算、扫描并拟合.对PO₂三原子分子优选出6-311+G（3df）方法计算结构参数、谐振频率、离解能及力常数，借助多体项展式理论导出PO₂自由基的势能函数并绘制等值势能图.发现：PO₂自由基的对称伸缩振动和旋转势能图中，在O+PO→OPO反应通道上都有鞍点出现，O原子需要越过0.55 eV的能量才能生成稳定的PO₂自由基.要形成PO₂自由基只能通过两等价的通道越过势垒才能形成.     

Funnels in energy landscapes

Local minima and the saddle points separating them in the energy landscape are known to dominate the dynamics of biopolymer folding. Here we introduce a notion of a “folding funnel” that is concisely defined in terms of energy minima and saddle points, while at the same time conforming to a notion of a “folding funnel” as it is discussed in the protein folding literature.     

Decay of 1+ states as a new probe of the structure of 0+ shape isomers

The nuclides 98Mo and 100Mo have been studied in photon-scattering experiments by using bremsstrahlung produced from electron beams with kinetic energies from 3.2 to 3.8 MeV. Six electromagnetic dipole transitions in 98Mo and 19 in 100Mo were observed for the first time in the energy range from 2 to 4 MeV. A specific feature in the two nuclides is the de-excitation of one state with spin J = 1 to the 0+ ground state as well as to the first excited 0+ state, which cannot be explained in standard models. We present a model that allows us to deduce the mixing coefficients for the two 0+ shape-isomeric states from the experimental ratio of the transition strengths from the J = 1 state to the 0+ ground state and to the 0+ excited state.