text:
# 1 Introduction

## 1  Introduction

Data on nuclear reactions at the energies of hundred MeV and GeV are required for multiple purposes such as long-lived radioactive waste transmutation, material analysis, nuclear medicine as well as research of the cosmic ray effects on spaceships and astronauts. Experiments on measuring the data are costly to be carried out and there is a limited number of facilities to make them. Therefore, reliable computer models for a simulation of the reactions are created to provide the necessary data. Most of them are using the ideas of the cascade-evaporation model (CEM) (see  \- ).

The cascading of nucleons and $`\pi `$-mesons was considered only by the first variants of the model. A good description of hadron-nucleus interactions was reached at the studies. The best result was obtained by S. Mashnik . Though, the application of the model for description of nucleus-nucleus collisions has shown that the model gives a satisfactory yield of the nucleon but overestimates the meson production. Taking into account meson and baryon resonances production is one of the possible ways to solve the problem.

A number of authors were trying to do this, and they usually obtained a decreasing yield of mesons and baryons. It is natural because the effective decreasing of the multiplicity of the produced particles leads to a less powerful cascading. So, a problem of a self consistent description of the meson and baryon yields in hadron-nucleus and nucleus-nucleus interactions appeared.

Let us note that the nuclear destruction mechanism and the procedure of the excitation energy calculation were not changed in the mentioned approaches. Maybe, they led to the unsatisfactory results. In this paper we consider a synthesis of the FRITIOF model  which takes into account the resonances production and the reggeon theory inspired model of the nuclear destruction .

The FRITIOF code , that is a program of Monte Carlo simulation of the inelastic hadron-hadron, hadron-nucleus and nucleus-nucleus interactions, is very popular in high energy experimental physics. It is explained by its access, its physical ideas simplicity and a defined beauty of the code by itself. It is easy to use.

The FRITIOF model  assumes that an excitation of hadrons into continuum mass spectra takes place in the inelastic hadron-hadron collisions. In the case of the hadron-nucleus or the nucleus-nucleus interactions the excited hadrons can suffer an additional collisions with the nuclear nucleons and go into more excited states, or de-excite. The excited hadrons are considered as quark strings, and the corresponding quark model  is used for description of their decay. The probabilities of the multiple collisions are calculated within the Glauber model . The inelastic collisions are usually considered only . In order to reproduce the baryon yield, we introduce elastic re-scattering too.

It is assumed that the program can not be used at the relatively low energies as the hypothesis about the creation and decay of the quark strings is not valid. Attempts of compelling the program to operate at the energies below 5-10 GeV/nucleon for AA-interactions usually failed. Though the analysis of the code operation shows that the program cycles due to its quite simple-hearted interpretation of the Fermi-motion of the nucleons. A change of the Fermi-motion simulation algorithm  allowed one to decrease a formal limit of the model application region.

Now in order to use the model at the intermediate energies, one needs to correct a scheme of the fragmentation of the quark string with low masses because in this case the excited hadrons have low masses, too. Below the corresponding changes of the model will be presented which allow one to reach a defined success. To solve the problem, we have used the experimental data on the neutron-proton interactions at momentum 1.25 – 5.1 GeV/c  given amiable by the Yu. A. Troyan’s group of the neutron-proton interaction study.

One of the disadvantages of the FRITIOF code is an omission of the slow particle cascading into the nuclei. Under the ”cascading” one usually understands the standard intranuclear cascade scenario (see, e.g., ), which neglects the quantum mechanical effects. The quantum mechanical description of the particles cascading in the nuclei can be achieved in the framework of the reggeon theory. According to the reggeon approach , a consideration of the cascade interactions necessitates a calculation of the yields of the so-called enhanced diagrams to the elastic scattering amplitude. Using the procedure for the calculation proposed in  and Abramovski-Gribov-Kancheli cutting rules , one can obtain the positive defined cross-sections of the inelastic processes. At first glance, this returns one to the classical cascade picture of the interactions. However, there is an essential difference. The cascade model assumes that the cascade is developed in a three-dimensional space of a target nucleus. According to the reggeon approach , the ”cascade” of the reggeon exchanges occurs in a two-dimensional space of projected radius-vectors of nucleons on the plane perpendicular to the momentum of projectile particle (on a plane of impact parameter), with a ”cascade power” independent of the multiplicity of the produced particles and defined by the reggeon vertex constants and the size of nucleus. We shell give a corresponding algorithm in Sec. 3. The main calculation results are presented in Sec. 4. We are starting with a short description of the main assumptions of the FRITIOF model.

## 2 Theses of the FRITIOF model

The FRITIOF model assumes that the hadron-hadron interactions are of two-particle character

$$a+ba^{}+b^{},$$
(1)
where $`a^{}`$ and $`b^{}`$ are $`a`$ and $`b`$ hadrons in excited states. The kinematics of the reaction is determined as follows: in the center of mass of colliding hadrons the energy-momentum conservation low has the form:

$`E_a+E_b`$ $`=`$ $`E_a^{}+E_b^{}=\sqrt{s_{ab}},`$
$`p_{az}+p_{bz}`$ $`=`$ $`p_{a^{}z}+p_{b^{}z}=0,`$ (2)
$`0`$ $`=`$ $`\stackrel{}{p}_a^{}+\stackrel{}{p}_b^{},`$
where $`E_a`$ and $`E_b(E_a^{},E_b^{})`$\- energies of the initial (final) hadron $`a`$ and $`b(a^{},b^{})`$, $`p_{az}`$ and $`p_{bz}`$ \- longitudinal momentum components (projected momenta on the interaction axis).

Adding and subtracting the first two equations from (2), we get

$`P_a^++P_b^+`$ $`=`$ $`P_a^{}^++P_b^{}^+`$
$`P_a^{}+P_b^{}`$ $`=`$ $`P_a^{}^{}+P_b^{}^{}`$ (3)
$`0`$ $`=`$ $`\stackrel{}{p}_a^{}+\stackrel{}{p}_b^{},`$
where $`P^+=E+p_z,P^{}=Ep_z.`$

At high energies: $`P_a^{}^{}m_a^{}^2/2p_{a^{}z},P_b^{}^+m_b^{}^2/2p_{b^{}z}.`$ Thus, the $`P_a^{}^{}`$ and $`P_b^{}^+`$ distributions used in the code have the form

$`dW`$ $``$ $`dP_a^{}^{}/P_a^{}^{}dm_a^{}^2/m_a^{}^2,`$
$`dW`$ $``$ $`dP_b^{}^+/P_b^{}^+dm_b^{}^2/m_b^{}^2.`$ (4)
The limits of $`P_a^{}^{}`$ and $`P_b^{}^+`$ are defined as

$$[P_a^{},P_b^{}],[P_b^+,P_a^+].$$
(5)
The distributions (4) are typical for the so-called high-mass diffraction dissociation processes.

Knowing $`P_a^{}^{}`$, $`P_b^{}^+`$, $`\stackrel{}{p}_a^{}`$, $`\stackrel{}{p}_b^{}`$ and determining $`P_a^{}^+`$, $`P_b^{}^{}`$ from Eqs. (3 ), one can find the masses of the excited hadrons $`a^{}`$ and $`b^{}`$. $`\stackrel{}{p}_a^{}`$ and $`\stackrel{}{p}_b^{}`$ are sampled according to the low

$$dWexp(\stackrel{}{p}_a^{}^2/<p_{}^2>)d^2p_a^{}.$$
(6)
In case of hadron-nucleons interactions the kinematics governed by Eqs. (3), (4), (6) is applied to the first collision of the projectile nucleon with one of the target nucleons ($`a+N_1a^{}+N_1^{}`$). For the second collision ($`a^{}+N_2a^{\prime \prime }+N_2^{}`$), analogous relations are used, but (5) is replaced for

$$[P_a^{}^{},P_{N_2}^{}],[P_{N_2}^+,P_a^{}^+].$$
(7)
As a result, the consequent collisions involve a systematic increasing of the mass of hadron $`a`$ if transfers of the transverse momentum are small.

A similar approach is also applied to simulate the nucleus-nucleus interactions. Here the reactions $`a^{}+b^{}a^{\prime \prime }+b^{\prime \prime }`$ are considered. The above distributions on $`P_a^{}^{}`$ and $`P_b^{}^+`$ are replaced by those on $`P_{a^{\prime \prime }}^{}`$ and $`P_{b^{\prime \prime }}^+`$ and the limits of $`P_{a^{\prime \prime }}^{}`$ and $`P_{b^{\prime \prime }}^+`$ are redefined as

$$[P_a^{}^{},P_b^{}^{}],[P_b^{}^+,P_a^{}^+].$$
(8)
At relatively lower energies of the order of 5–10 GeV/nucleon, the FRITIOF model without taking into account hadron de-excitation overestimates the multiplicity of the produced particles in hA- and AA-interactions. To this end we have changed the condition (8) by the following giving an allowance for considering the excitation process with the increasing mass and the de-excitation process with decreasing mass:

$$[P_a^{},\sqrt{s_{a^{}b^{}}}m_b],[P_b^+,\sqrt{s_{a^{}b^{}}}m_a].$$
(9)
$`P_a^{}`$ and $`P_b^+`$ are calculated at $`s_{ab}=s_{a^{}b^{}},m_a^{}=m_a,m_b^{}=m_b`$. The minimal values of $`P_a^{}`$ and $`P_b^+`$, being the lower limits in (9), are obviously reached in the reaction $`a^{}+b^{}a+b`$. The maximal values are achieved at $`p_{a^{\prime \prime }z}=p_{b^{\prime \prime }z}=0`$ when hadrons come to rest in the c. m. frame of $`NN`$-collision.

The reactions (1), or

$`a+ba^{}+b,`$ (10)
$`a+ba+b^{},`$ (11)
are the so-called diffraction dissociation reactions. The reactions (10), (11) are one-vertex diffractions, the reaction (1) is a double vertex diffraction. It is obvious that a minimal mass of the excited nucleon in NN-interactions can be equal to $`m_a^{}=m_N+N_\pi =1080`$ MeV. In the FRITIOF model the minimal mass is equal to 1.2 GeV. Due to this, the one vertex diffraction with excitation only one hadron is possible at $`\sqrt{s_{NN}}2.4`$ GeV, or at $`P_{lab}1.91`$ GeV/c in NN-interactions. There can be a two vertex diffraction at higher energies. A relation between the cross-sections of the processes is determined by the hadrons $`a^{}`$ and $`b^{}`$ mass distributions.

The excited hadrons $`a^{}`$ and $`b^{}`$ are considered as quark strings, and the corresponding quark model is used for the simulation of their decay . It is assumed that the quark model can be used at sufficient large string masses what can be created at high energies. Thus the FRITIOF model was used mainly at high energies. One can expect that the FRITIOF model predictions fall into a contradiction with experimental data on the processes where the states with low masses can appear. In order to study the situation, we turn to the data on $`\pi ^{}`$-meson and proton distributions in $`np`$-interactions at $`P_n=`$ 1.25–5.1 GeV/c .

Fig. 1 shows the experimental and calculated $`\pi ^{}`$-meson rapidity distributions, $`y=\frac{1}{2}\mathrm{ln}(E+P_z)/(EP_z),`$ where $`E`$ and $`P_z`$ are laboratory energy and longitudinal momentum of $`\pi ^{}`$-meson, respectively. As seen, the distributions calculated according to the original code (dashed lines) are of two-bump structure more pronounced at low energies. It seems like it is a circumstance of the assumed diffractive character of the interactions: the bump at large rapidities is caused by a projectile hadron diffraction, the bump at low rapidities is connected with a target nucleon diffraction. Though, at the neutron momentum of 1.25 GeV/c a diffraction system with mass of 1.2 GeV (the minimal mass of excited nucleon assumed by the FRITIOF model) must be in a rest in the center of mass system, and there must not be a subdivision of the fragmentation regions. Thus, we conclude that the two-bump structure of the calculated distributions is not a circumstance of the diffraction dissociation of the hadrons. It only reflects an anysotropy of the low mass string decay. A direct simulation of the decay of the strings with low masses has shown that the bi-module structure is typical for the decay of the strings with masses lower than 1.7 GeV. Two-particle channel is dominating in the decay of such strings. At higher masses the multi-particle channel gets more probable, and the calculated distributions get more regular.

Taking into account the character of the experimental distributions at $`P_n`$ = 1.25, 1.73 GeV/c, it seems reasonable to simulate an isotropic decay of the strings with low masses in the case of the two-particle decay channel. A boundary value of the string mass of 1.7 GeV below of which we have used the proposed procedure, was chosen requiring a good description of the data. As seen in the Fig. 1 (solid curves), this allows us to describe the experimental data quite well.

A more complete situation appears in the description of the $`\pi ^{}`$-meson distributions in transverse momentum, $`P_T`$. The original model predicted the average transverse momentum larger than the experimental ones (see the dashed curves in Fig. 2), though we had changed the character of the low mass string decay. Since at $`P_n=1.25`$, we simulate an isotropic decay of the excited hadrons, the maximum of the distribution in $`P_T`$ is determined by mass of the decayed system. Thus, in order to decrease the average transverse momentum, the minimal mass of the excited nucleon state was decreased to the value of 1.1 GeV. This, as seen from Fig. 2, gave a better result at $`P_n<2`$ GeV/c. At larger energies one needs to take into account the transverse momentum transferred by the colliding nucleons.

An analysis of proton spectra gives some additional information, and allows one to determinate the model parameters more exactly. For the analysis we have used the data on the following reactions

$`np`$ $``$ $`pp\pi ^{};`$ (12)
$`np`$ $``$ $`pp\pi ^{}\pi ^0;`$ (13)
$`np`$ $``$ $`np\pi ^+\pi ^{}.`$ (14)
In the reaction (12) according to the model the diffraction dissociation processes of the projectile particle is a dominant one. In the reaction (13) the two-vertex diffraction processes do the same. At last, in the reaction (14) we have equal yields of the one-vertex and two-vertex diffraction. Thus, the study of the reactions allows one to check the different components of the model.

Fig. 3 gives the experimental and calculated rapidity distributions of the protons. The experimental data are presented by the histograms, the original model calculations - by the dashed curves, and the last calculations - by the solid curves. As seen, there is a two-bump structure of the experimental distributions on the reactions (12). The bump at small rapidities is caused by the saved target protons. The bump at large rapidities is connected with the protons created in the projectile particle diffraction. The distribution of the reaction (13) has no structure. At last, in the reactions (14) there is a dominant production of the protons in the target fragmentation region.

The original model calculations are in agreement with the data on the reaction (14). At the same time, the calculated distribution for the reaction (12) has a bump in the region $`y1.7`$ what does not observe at the experiment. For its elimination we introduce a charge exchange between the nucleons in 50 % of the two-vertex diffraction. This allowed us to improve the proton spectra description in part.

Enumeration of the changes made in the FRITIOF code

1. The minimal mass of the excited nucleon decreases from 1.2 GeV to 1.1 GeV;
2. In case of the two-particle decay channel of a string with a mass lower than 1.7 GeV the isotropic decay is simulated in the center of the mass system;
3. The charge exchange between the colliding nucleons is allowed in 50 % of the two-vertex diffraction;
4. The value of the average square of the transverse momentum which transferred between the colliding nucleons increases from 0.08 (GeV/c)2 to 0.15 (GeV/c)2.

## 3 Simulation of nuclear destruction at the fast stage of interaction.

### 3.1 Determination of the number of knocked-out nucleons

In the last few years there have been some successful attempts to describe the hadron-hadron elastic scattering at low and intermediate energies (below $`12GeV`$) within the quark-gluon approach (see Refs.  \- ). In Ref.  - the amplitudes of $`\pi \pi `$, $`K\pi `$ and $`NN`$ scattering were found and an agreement of the theoretical calculations with corresponding experimental data was reached at the assumption that in the elastic hadron scattering one-gluon exchange with the following quark interchange between hadrons takes place (see Fig. 4a). At high energies two-gluon exchange appropriation (Fig.4b) works quite well (see Ref. ,  and ). What kind of exchanges can dominate in hadron-nucleus and nucleus-nucleus interactions?

The simplest possible diagrams of the processes with three nucleons are given in Fig. 5. Calculation of their amplitudes according to Refs. - is a serious mathematical problem. It can be simplified if one takes into account an analogy between the quark-gluon diagrams and the reggeon diagrams: the quark diagram of Fig. 4a corresponds to a one-nonvacuum-reggeon exchange diagram; the diagram of Fig. 4b describes the pomeron exchange in the $`t`$ channel; the diagram of Fig. 5a is in a correspondence with the enhanced reggeon diagram of the pomeron splitting into two non-vacuum reggeons. The three pomeron diagram (Fig. 5d) represents a more complicated process. It is rather hard to find a correspondence between the reggeon diagrams and the diagrams of Fig. 5b, 5c.

The reggeon parameters and the functional forms of the amplitudes of 3-reggeon processes are well known. The constants of the reggeon interaction vertexes are poor determined. The 3-pomeron vertex constant $`G_{PPP}`$ is well established ($`G_{PPP}=1.35^2(GeV)^2`$, Ref. ). There are only old data  and the estimations of Ref.  on the values of other constants - $`G_{PRR}`$ and $`G_{RRR}`$, which are large. Nevertheless, we believe that the properties of the reggeon amplitudes must be taken into account in consideration of the nuclear destruction.

It is obvious that the processes like that in Fig. 5d can not dominate in the elastic hadron-nucleus scattering because they are accompanied by production of a high mass diffraction beam of the particles in the intermediate state. Thus, their yields are dumped by a nuclear form-factor. According to the same reason, the yields of the processes like ones in Figs. 5a, 5b can be small, too. If it is not so, one will expect large corrections to Glauber’s cross-sections. The practice shows that the corrections to the hadron-nucleus cross-sections must be lower than 5 – 7 $`\%`$.

The yield of the diagram of Fig. 5c gives a correction to Glauber’s one-scattering amplitude. There must be analogous corrections to the other terms of Glauber’s series. The sum of the corrections must lead to small effects in the elastic small angle scattering because the corrections are large at small impact parameters. So, they can manifest themselves at large scattering angles. We assume that they have a big influence on the inelastic process characteristics, too.

According to the reggeon theory, a description of the inelastic reactions can be reached in a consideration of the different cuts of the reggeon diagrams. Here the Abramovski - Gribov - Kancheli cutting rules  are often used. The corrections to them were discussed in Ref.  in the application to the problem of a particle cascading on the nucleus. As was shown in Ref. , summation of the yields of enhanced diagrams allows one to describe increasing of the one-particle spectra in the target fragmentation region. At the same time, the authors of Ref.  did not take into account the shadowing effects considered in Ref. .

Here we have to note that the yields of the diagrams like that shown in Fig. 5c have no shadowing corrections. The yield of the enhanced diagram of Fig. 5a has a form

$$Y_aexp[(\stackrel{}{b}_1\stackrel{}{b}_2)^2/3r_a^2(\stackrel{}{b}_1\stackrel{}{b}_3)^2/3r_a^2(\stackrel{}{b}_2\stackrel{}{b}_3)^2/3r_a^2]$$

where $`\stackrel{}{b}_1,\stackrel{}{b}_2`$ and $`\stackrel{}{b}_3`$ are the impact coordinates of the nucleons. At the same time, the yield of the diagram of Fig. 5c according to Refs.  \-  is given by

$$Y_cexp[(\stackrel{}{b}_1\stackrel{}{b}_2)^2/r_c^2]exp[(\stackrel{}{b}_2\stackrel{}{b}_3)^2/r_c^2].$$

In the limit of $`r_a^2,r_c^2R_A^2`$, where $`R_A`$ is a nucleus radius, the yields coincide. Thus, we can save the results of Ref.  considering them as a summation of the yields of the quark-gluon diagrams.

Let us note that neither $`Y_a`$, nor $`Y_c`$ depend on the longitudinal coordinates or on the multiplicity of produced particles. It is the main difference between the ”reggeon cascading” and the ”usual” cascading.

As well known, the intranuclear cascade model (-) assumes that in a hadron-nucleus collision the secondary particles are produced due to an inelastic interaction of the projectile particle with a target nucleon. The produced particles can interact with other target nucleons. A distribution on distance $`l`$ between the first interaction and the second one has a form

$$W(l)dl\frac{n}{<l>}exp(\frac{n}{<l>}l),$$

where $`<l>=1/\sigma \rho _A`$, $`\sigma `$ is a hadron-nucleon cross-section, $`n`$ is the multiplicity of the produced particles and $`\rho _A0.15fm^3`$ is the nuclear density. At the same time the amplitudes or the cross-sections of the processes shown in Fig.5 have no dependence on $`l`$ or $`n`$. Thus, we expect that in the quark-gluon or reggeon approach the ”cascade” will be more restricted than in the cascade model. The difference between approaches can lead to the different predictions for the light nuclei destruction (an effect of the limited volume) and for the characteristics of the heavy nuclei interactions (an influence of a large multiplicity of the produced particles).

To show this, we use a simple method to estimate the nuclear destruction in the framework of the quark-gluon approach.

1. As it was said above, the ”reggeon cascade” is developed in the space of the impact parameter. Thus, for its description it is needed to determinate a probability to involve a nucleon into the ”cascade”. It is obvious that the probability depends on a difference of the impact coordinates of the new and the previously involved nucleons. Looking at the yield of the diagram of Fig. 5c, we choose the functional form of the probability as

$$P(\stackrel{}{b}_i\stackrel{}{b}_j)=C_{nd}exp((\stackrel{}{b}_i\stackrel{}{b}_j)^2/r_nd^2).$$
(15)
Here $`\stackrel{}{b}_i`$ and $`\stackrel{}{b}_j`$ are projections of the radiuses of $`i^{th}`$ and $`j^{th}`$ nucleons on the impact parameter plane.
2. The ”cascade” is initiated by the primary involved, wounded nucleons. If the constant $`C_{nd}`$ is small, we can use the Glauber theory for their determination.
3. We assume that all the involved and wounded nucleons are ejected from the nucleus.

The ”cascade” looks as follows: a projectile particle interacts with some of the intranuclear nucleons. They are called ”wounded” nucleons. The wounded nucleons initiate the ”cascade”. A wounded nucleon can involve a spectator nucleon into the ”cascade” with the probability (15). The latter can involve a second nucleon. The second nucleon can involve a third one, and so on.

A Monte Carlo algorithm for estimation of the nuclear destruction in the nucleus-nucleus interactions which corresponds to the model formulation, includes the following steps:

1. The calculation of the impact parameter distribution in the framework of the Glauber theory ;
2. The sampling of the impact parameter and the nucleon coordinates;
3. The determination of the wounded nucleons (see Ref. );
4. The determination of the spectator nucleons involved in the ”cascade” by the wounded nucleons. If the number of the involved nucleons is equal to zero - exit;
5. If the number of the involved nucleons is not equal to zero, a possibility is considered to involve the other spectators nucleons by the involved ones. If the number of the new involved nucleons is equal to zero - exit. In other case - it is needed to repeat the step 5 taking into account only the new involved nucleons.

The first step is performed only once at the given mass numbers of the projectile and target nuclei. The steps 2 – 5 are repeated until the needed statistics is reached. The steps 4, 5 are applied to the nucleons of projectile and target nuclei.

It is suggested that all the newly involved participants and the ”wounded” nucleons are knocked-out from the nucleus.

### 3.2 Fermi-motion of nucleons

To take into account the energy-momentum conservation law by simulating compound system, let us consider a reaction of the compound system $`(1,2)`$ with hadron $`h`$: $`(1,2)+h1+2+h`$. Neglecting the transverse momenta, a final state of the reaction will be fully characterized by a value of merely one independent kinematical variable. As the variable, let us take

$$x_1^+=(E_1+p_1)/(E_1+E_2+p_1+p_2).$$

It is useful to introduce the analogous quantity $`x_2^+`$ which satisfies obviously the relation $`x_1^++x_2^+=1`$. At given $`x_1^+`$ and $`x_2^+`$, the rest of the kinematical variables can be determined by the energy-momentum conservation law.

In case of dissociation of two compound systems $`A`$ and $`B`$ containing $`A`$ and $`B`$ constituents respectively, let us introduce for the $`i`$-th constituent of system $`A`$

$$x_i^+=(E_{Ai}+p_{iz})/W_A^+and\stackrel{}{p}_i,$$

and for the $`j`$-th constituent of system $`B`$

$$y_j^{}=(E_{Bj}q_{jz})/W_B^{}and\stackrel{}{q}_i,$$

where, $`E_{A_i}(E_{B_i})`$ and $`\stackrel{}{p}_i(\stackrel{}{q}_i)`$ are energy and momentum of $`i`$-th constituent from $`A(B)`$,

$$W_A^+=\underset{i=1}{\overset{A}{}}(E_{Ai}+p_{iz}),W_B^{}=\underset{i=1}{\overset{B}{}}(E_{Bi}q_{iz}).$$

One can find $`W_A^+`$ and $`W_B^{}`$ using the energy-momentum conservation low at given $`\{x_i^+,\stackrel{}{p}_i\}`$, $`\{y_i^{},\stackrel{}{q}_i\}`$, and determine all kinematical variables for all constituents .

According to the experimental evidence , the average transverse momentum of spectator fragments obeys the parabolic law:

$$<P_{}^2>=\frac{A(AF)}{A}<p_{}^2>,\sqrt{<p_{}^2>}=0.07GeV/c.$$

To reproduce this result, the values of $`\stackrel{}{p}_i`$ for knocked-out nucleons are simulated according to the distribution

$$dWexp(\stackrel{}{p}_i^2/<p_{}^2>)d^2p_i,\sqrt{<p_{}^2>}=0.07.$$
(16)
The sum of the transverse momenta (with sign ”minus”) was ascribed for the residual nucleus.

The choice of $`x_i^+`$ is carried out by

$$dWexp[(x_i^+1/A)^2/(d_x/A)^2]dx_i^+,0x_i^+<1,d_x=0.07.$$
(17)
$`x^+`$ of the residual nucleus is determined as $`1x_i^+`$ .

It was assumed that the knocked-out nucleons changed their characteristics again. The new values of $`x_i^+`$ and $`p_i`$ were simulated using the distributions (16) and (17) at $`<p_{}^2>=0.385`$ $`(GeV/c)^2`$ and $`d_x=0.2`$. The results from  were used for determination of the fitting parameters.

### 3.3 Excitation energy of residual nucleus.

The change of the nuclear destruction mechanism requires a change of the nuclear residual excitation energy calculation procedure. In a self-consistent diagram approach one needs to consider the more complete diagrams and describes bound states of the nucleons. The last problem is not solved in the presented quark-gluon approach which allows one only to calculate an repulsive part of the NN-potential. It is expected that a taking into account the terms with high order vertex constant allows one to calculate the attractive part of the potential too. Until this is not made, we are to use a phenomenological approach. Here we follow the Ref. .

In paper  proton-nucleus interactions at intermediate energies were analyzed. The first stage of the interactions was considered within the Glauber approach. It was assumed that the projectile hadron undergoes successive interactions (elastic and inelastic) with target nucleon. In each of the collisions a part of the energy of the projectile hadron, $`E`$, is transferred to the target nucleon. The distribution on the transfer energy was chosen in the form:

$$F_1(E)=\frac{1}{E}e^{E/E}.$$
(18)
The authors of Ref.  supposed that the excitation energy of residual nucleus was the sum of the recoil nucleon energies. As a result, they described experimental data on neutron multiplicity dependence upon the excitation energy of the residual nuclei.

To apply this approximation to $`AA`$-interactions at high energies, one needs to evaluate the number of re-scatterings of each of the knocked-out nucleons. Taking into account that most of $`AA`$-interactions are of peripheral nature, we assume that the slow nucleons knocked-out from the peripheral parts of nucleus can not penetrate deeply inside the nucleus because of the large $`NN`$ cross-section. Their re-scatterings thus occur in the nearest environment. We consider the environment as the spectator nucleons being inside a sphere with radius $`r_0=2`$ fm surrounding a wounded or involved nucleon in its initial state. It is assumed that each spectator nucleon may acquire recoil energy distributed according to Eq. (18). Remember, that the nucleon coordinates were chosen randomly and independent according to the Saxon-Woods distribution. If a spectator nucleon was a neighbor of two ”wounded” ones, or involved nucleons, it received the sum of two energies chosen according to Eq. (18), and so on. The sum of energies transferred to all spectator nucleons was considered as the excitation energy of the residual nucleus.

Unlike the CEM, this method obviously will lead to the zero-excitation energy when all nucleons are be ejected.

It also seems evident that the boundary between the spectator part of the nucleus and the part affected by the first stage of the interaction which governs the excitation energy of the residual nucleus, depends on the impact parameter. For example, with heavy projectile nuclei the excitation energy rises with decreasing impact parameter from $`R_A+R_B`$ to 0, approaching a maximum and then falls. According to the CEM, it must gradually increase.

We have used $`E=8`$ MeV and the standard evaporation model  (see, also Ref. ) for simulation of the residual nucleus de-excitation.

## 4 Description of nucleus-nucleus interactions

Fig. 6 shows the proton and $`\pi ^{}`$-meson rapidity distributions in different nucleus-nucleus interactions. They were calculated by the modified FRITIOF code and by the code of the cascade-evaporation model333 )The model takes into account the trailing effect, Pauli principle, the dependence of the Fermi momentum on the local nuclear density, the pre-equilibrium emission and the evaporation of the nuclei. . Events with at least one inelastic NN-collision were selected at the simulations. As seen, for $`dd`$-interactions the models give close results. Since in this case we can neglect the cascade interactions, the coincidence of the results tells us about correctness of the NN-interactions description in the CEM (the CEM code used by us does not allow a direct simulation of NN-interactions). We can mark only a little enhanced baryon production in the central rapidity region in the CEM.

In $`\alpha \alpha `$-interactions where the influence of the cascade interactions is sufficiently large, we observe a difference of the predicted $`\pi ^{}`$-meson spectra. As expected, taking into account the resonances in the FRITIOF model leads to a decreased meson yield. The enhanced meson production in CEM becomes more pronounced in $`CC`$-interactions.

To describe the proton production in the nuclear fragmentation regions in the framework of the FRITIOF model, we take into account both inelastic interactions of the nucleons considered above and elastic re-scatterings. As seen, we reproduce the baryon yields in the fragmentation regions at $`y0,2.2`$ for the $`dd`$\- and $`\alpha \alpha `$-interactions. The CEM calculations for $`CC`$-collisions were used for determination of the nuclear destruction model parameters, $`C_{nd}=1,r_{nd}^2=1.4`$ (fm2). Then we introduced the experimental criteria for the proton registration.

The calculations for heavy nuclei look most interesting. In Fig. 7 the proton distributions are presented on total and transverse momentum in $`n+Ta`$ interactions at $`P_n=4.2`$ GeV/c. As seen, the model predictions are close to each other. The same closeness we observe for $`C+Ta`$ interactions at $`P=4.2`$ GeV/c/nucleon presented in Fig. 8. Though, there is a wide difference between the predicted meson spectra. The FRITIOF model calculations are close to the experimental data. At the calculation we use the following values of the parameters $`C_{nd}=0.2,r_{nd}^2=1.1`$ (fm2) to reproduce the $`Ta`$ nuclei destruction.

Summing up, we can conclude that we have reached a satisfactory description of the meson and nucleon production in the nucleus-nucleus interactions at the energy of 3.3 GeV/nucleon in the framework of the sufficiently simple FRITIOF model. The model can be applied for practical calculation of nucleus-nucleus interaction characteristics.

The authors of the paper are thankful to Prof. V.S. Barashenkov, Zh.Zh. Musulmanbekov and B.F. Kostenko for useful discussions. One of the authors (V.V.U.) thanks RFBR (grand No. 00-01-00307) for its financial support.


format: text

text:
# Wedge filling, cone filling and the strong-fluctuation regime

## 1 Introduction

Following the seminal studies of Cahn  and Ebner and Saam , two fundamental questions concerning the nature of the wetting transition have drawn an enormous amount of attention : Firstly, is the transition first-order or continuous? In particular, how does the order depend on the range and balance of the competing intermolecular forces? Secondly, for continuous (critical) wetting transitions, what role do the interfacial fluctuations of the unbinding fluid interface play in determining the critical exponents characterising the diverging lengthscales?

The answers to both these questions reveal a host of reasons why the laboratory observation of fluctuation effects at critical wetting is beset with difficulties. In brief, the well developed fluctuation theory of wetting concludes  that three dimensional critical wetting requires a fine tuning of the effective Hamaker constants appearing in the binding potential. This sensitivity means that critical wetting is a rather rare phenomena of which only very few examples are known for fluid-fluid interfaces . Secondly, for systems with long-ranged forces, the critical exponents at critical wetting are mean-field-like (and depend on the range of the forces) and the interfacial roughness induced by fluctuations is very small. Even for short-ranged forces, where renormalization group (RG) theory based on effective interfacial models famously predicts dramatic fluctuation induced effects , no appreciable deviations from simple mean-field like behaviour have been observed in either Ising model simulation  or recent experimental studies . At the very least, this indicates that the asymptotic critical regime for critical wetting is much smaller than initially expected . The conclusion inevitably forced upon us appears to be that fluctuation effects at wetting transitions are of negligible practical importance.

The purpose of the present article is to compare and contrast these features of critical wetting transitions with the properties of continuous (critical) filling transitions occuring in both wedge and cone shaped geometries. The thermodynamics of the filling transition has been discussed independently by a number of authors  but only recently have the analogues of the above questions been addressed . Here we present full details of a novel interfacial theory of wedge filling and also extend this approach to study fluctuation effects occuring at the filling transition of a fluid adsorbed in a cone shaped geometry. These turn out to be completely different to those predicted to occur in the wedge. In particular we wish to draw attention to two striking features of conic filling which have not been reported before. First, that the conditions for critical conic filling, even more than conditions for critical wedge filling, do not rely on any fine tuning of the binding potential structure and should be observable in cone geometries made from substrates that exhibit even strong first-order wetting. Secondly we point out what appears to be an unexpected connection between the geometry-induced aspects of critical conic filling and the theory of the universal strong-fluctuation (SFL) regime of critical wetting . Specifically, from analysis of effective interfacial Hamiltonian models, we show that for the physically relevant dimensions $`d=2`$ and $`d=3`$, and also for almost all choices of the intermolecular interactions, the influence of the conic geometry on the asymptotic critical behaviour and the detailed scaling structure of the interfacial height probability distribution function (PDF) as one approaches a critical cone filling transition is to mimic the way in which interfacial fluctuations determine the analogous exponents and the PDF for critical wetting transitions for planar wall-fluid (or fluid-fluid) interfaces belonging to the SFL scaling regime. This remarkable connection between conic filling and wetting and has already been noted in $`d=2`$ from the explicit results of exact transfer matrix and replica trick studies of two-dimensional filling transitions (where, of course, the cone is analogous to a wedge) in both pure and impure bulk systems . Extending this theory our analysis shows that the fluctuation properties of three-dimensional conic filling are quite distinct from those occuring at three-dimensional wedge filling  and are instead related to some of the most dramatic predictions of the theory of wetting for planar wall-fluid interfaces. Thus in $`d=3`$ the divergence of the equilibrium filling layer thickness and roughness at cone filling are universal and the same as the renormalization group predictions of Brezin, Halperin and Leibler  for critical wetting with short-ranged forces in the large wetting parameter ($`\omega >2`$) regime.

Our article is arranged as follows. In the next section we briefly recall the well known fluctuation theory predictions for critical wetting detailing the critical behaviour predicted to occur in the SFL regime in $`d=2`$ and $`d=3`$. In §3 we present the main methods and results of the fluctuation theory of three dimensional wedge filling emphasising both the conditions for critical filling and the predictions for fluctuation dominated behaviour. A brief summary of these results has been given before . We finish this section by turning to the special case of filling in $`d=2`$. The role of a breather-like soft mode which increases the width of the interface as it translates up and down the wedge is shown to be directly responsible for the scaling behaviour of the PDF. In §4 we discuss three dimensional cone filling, at mean-field level and beyond, predicting universal geometry-dominated critical behaviour characterised by an appreciable interfacial roughness and critical exponents identical to those predicted for the planar wetting SFL regime. We argue that this equivalence of the PDF’s at cone filling and SFL wetting is special to $`d=2`$ and $`d=3`$ only. We finish our article with a summary of our main results.

## 2 Fluctuation regimes at critical wetting.

Consider the preferential adsorption of a thick liquid film at the interface between a planar wall and a bulk vapour phase at chemical potential $`\mu `$ and a sub-critical temperature $`T`$. We suppose that at bulk two phase coexistence, corresponding to chemical potential $`\mu =\mu _{\text{sat}}`$, the wall-vapour interface is completely wet by the liquid phase for temperatures $`T>T_\pi `$. Here we are interested in the values of the critical exponents describing the continuous divergence of the mean wetting film thickness $`l_\pi t^{\beta _s}`$, roughness $`\xi _{}t^\nu _{}`$ and transverse correlation length $`\xi _{}t^\nu _{}`$ as $`t^{}=(T_\pi T)/T_\pi 0`$. The subscript $`\pi `$ is used to emphasise that these quantities are pertinent to planar wetting transitions. We denote the reduced temperature $`t^{}`$ to avoid confusion with the filling reduced temperature (q.v.). The standard theory of interfacial fluctuation effects at the wetting transition is based on analysis of the effective interfacial Hamiltonian

$$H[l]=𝑑𝐱\left\{\frac{\mathrm{\Sigma }}{2}\left(l\right)^2+W(l)\right\}$$
(2.1)
where $`l`$ is a collective co-ordinate describing the local height of the interface above the wall, $`\mathrm{\Sigma }`$ is the surface stiffness (tension) of the unbinding liquid-vapour interface and $`W(l)`$ is the usual binding potential which accounts for the direct influence of the molecular forces. Note that throughout this article we shall suppose that the factor $`\beta (1/k_BT)`$ is absorbed into all the effective Hamiltonians, stiffnesses and potentials etc. Mean-field (MF) studies indicate that for critical wetting (with long-ranged forces) $`W(l)`$ is necessarily of the form

$$W(l)=\frac{a}{l^p}+\frac{b}{l^q}+\mathrm{};l>0$$
(2.2)
where $`a(T_\pi ^{MF}T)`$ and $`b>0`$ are effective Hamaker constants. Here the exponents $`p`$ and $`q`$ are specific to the particular ranges of the intermolecular forces. In $`d=3`$ and for van der Waals forces the appropriate values are $`p=2`$ and $`q=3`$. For fixed dimension $`d<3`$ the scaling theory of Lipowsky and Fisher  (which is supported by transfer-matrix and approximate non-linear renormalization group (RG) studies) predicts that the critical behaviour falls into three generic fluctuation classes labelled mean-field (MF), weak-fluctuation (WFL) and strong-fluctuation (SFL) scaling regimes. These regimes arise due to the interplay of the direct binding potential with the entropic repulsion originating from the interfacial wandering. We briefly recall the scaling argument here since the same ideas generalise to filling transitions in wedge geometries. First recall that for interfaces, the roughness $`\xi _{}`$ and the transverse correlation length $`\xi _{}`$ are related via the wandering exponent $`\xi _{}\xi _{}^{\zeta _\pi }`$. We write $`\zeta _\pi `$ keeping the sub-script to avoid confusion with different wandering exponents that emerge for filling in wedges and cones. For thermal forces and $`d<3`$ it is well known that

$$\zeta _\pi =\frac{3d}{2}$$
(2.3)
with $`\xi _{}`$ finite for $`d>3`$ and $`\xi _{}\sqrt{\mathrm{ln}}\xi _{}`$ for the marginal case. Thus we are led rather naturally  to deduce that the bending energy term in (2.1) and also the entropic effect of interface collisions with the wall gives rise to a fluctuation contribution to the effective potential of the form $`l^\tau `$ with $`\tau =2(d1)/(3d)`$ for $`d<3`$. The competition between these three terms gives rise to the three fluctuation regimes mentioned above. In the MF and WFL regimes the leading term in 2.2 is still dominant and the wetting transition temperature is unchanged by fluctuations. The critical exponents are non-universal and are given by

$$\beta _s^{MF}=\frac{1}{qp},\nu _{}^{MF}=\frac{q+2}{2(qp)},\nu _{}^{MF}=\zeta _\pi \nu _{}^{MF}$$
(2.4)
$$\beta _s^{WF}=\frac{1}{\tau p},\nu _{}^{WF}=\frac{\tau +2}{2(\tau p)},\nu _{}^{WF}=\zeta _\pi \nu _{}^{WF}$$
(2.5)
In the SFL, corresponding to $`p<\tau `$, and representing the critical behaviour in systems with short-ranged forces, the transition temperature is lowered by fluctuations and the critical exponents are universal. In $`d=2`$ the critical exponents and the interfacial PDF, $`P_\pi (l)`$, are very well known (see for example  and also references in ) and given by

$$\beta _s^{SF}=1,\nu _{}^{SF}=1,\nu _{}^{SF}=2\text{ and }P_\pi (l)=\frac{1}{l_\pi }e^{\frac{l}{l_\pi }}$$
(2.6)
Notice that in both the SFL and WFL the interfacial fluctuations are very strong and $`l_\pi \xi _{}`$ for $`d<3`$. The unusual bifurcation mechanism controlling the fixed point behaviour as $`d3`$ means that critical wetting in three-dimensions systems is rather different in the SFL regime, though for systems with long-ranged forces MF theory remains valid. Fluctuation effects are only predicted to occur for systems with short-ranged forces and are traditionally and most simply described by a bare binding potential of the form

$$W(l)=ae^{\kappa l}+be^{2\kappa l}+\mathrm{};l>0$$
(2.7)
together with the hard-wall restriction $`W(l)=\mathrm{}`$ for $`l<0`$. Here $`\kappa `$ denotes the inverse bulk correlation length of the adsorbed bulk phase whilst the Hamaker constants $`a`$ and $`b`$ have the same interpretation as for long-ranged forces discussed above. The mean-field predictions for short-ranged critical wetting are therefore $`\beta _s=0`$ (ln) and $`\nu _{}=1`$. The fluctuation theory of critical wetting with short-ranged forces was famously investigated by Brezin, Halperin and Leiblier  using a linearised RG scheme and shows that the critical behaviour is strongly non-universal and falls into three regimes (labelled I - III) depending on the value of the dimensionless wetting parameter $`\omega `$ where

$$\omega =\frac{k_BT_\pi \kappa ^2}{4\pi \mathrm{\Sigma }}$$
(2.8)
The three regimes correspond to the ranges $`\omega <1/2`$, $`1/2<\omega <2`$ and $`\omega >2`$ respectively. Usually attention focuses on the value of the transverse correlation length critical exponent which shows the strongest non-universal critical behaviour:

$$\text{(I)}\nu _{}=\frac{1}{1\omega },\text{(II)}\nu _{}=\frac{1}{(\sqrt{2}\sqrt{\omega })^2},\text{(III)}\nu _{}=\mathrm{}$$
(2.9)
where the last case corresponds to an essential singularity. In all three regimes the roughness $`\xi _{}\sqrt{l_\pi }`$ so that interfacial wandering has a much more pronounced influence on the equilibrium distribution of matter in the wetting film compared to systems with long-ranged forces. In regimes (I) and (II) the wetting temperature is unchanged from its mean-field value and the film thickness still grows logarithmically but with a different amplitude to mean-field theory. The most dramatic effect of interfacial wandering is manifest in the regime (III) which bears all the hallmarks of the SFL regime.Here the wetting temperature is depressed by the fluctuations, and the critical exponents are given by

$$\text{(III)}\beta _s=1,\nu _{}=\frac{1}{2}$$
(2.10)
with an equilibrium PDF which is a simple Gaussian.

## 3 Fluctuation Regimes at Critical Wedge Filling.

### 3.1 Three-dimensional wedge filling: phenomenology.

Consider a three-dimensional wedge formed by the junction of two walls at angles $`\pm \alpha `$ to the horizontal with the height of the wall above the plane is described by a height function $`z_w(x,y)\alpha |x|`$. Thus the wedge has a ‘V’ shaped cross-section in the $`x`$-direction with the wedge bottom oriented along the $`y`$ axis. Following our earlier discussion of wetting at planar walls we suppose that the wedge is in contact with a bulk vapour phase at saturation chemical potential $`\mu =\mu _{\text{sat}}`$. General thermodynamic considerations  indicate that the wedge is completely filled by liquid for temperatures greater then the filling temperature $`T_F`$, which is specified by the elegant equality

$$\theta _\pi (T_F)=\alpha $$
(3.1)
where $`\theta _\pi (T)`$ denotes the temperature dependent contact angle of the liquid drop at the planar wall-vapour interface. Thus the influence of the wedge geometry is to lower the temperature at which the liquid-vapour interface unbinds from the wedge bottom compared to the planar wall. The filling transition occurring as $`t=(T_FT)/T_F0`$ may be first-order or critical corresponding to the discontinuous or continuous divergence of the equilibrium interfacial height $`l_0l(x=0,y)`$ measured from the bottom of the wedge. Critical filling is also characterised  by the divergence of distinct correlation lengths describing fluctuations in the interfacial height along ($`\xi _y`$) and across ($`\xi _x`$) the wedge and also by the roughness $`\xi _{}`$ of the unbinding interface. Note that the cross-wedge correlation length essentially measures the width of the flat, filled region of the wedge. The flatness of the filled region simply reflects the absence of any macroscopic curvature of the meniscus at two phase coexistence. Recalling our earlier notation we define the wedge wetting critical exponents via

$$l_0t^{\beta _s^W},\xi _{}t^{\nu _{}^W},\xi _xt^{\nu _x},\xi _yt^{\nu _y}$$
(3.2)
where we have adopted a suitable super-script, $`W`$ (for wedge), for some of the exponents to avoid confusion with those pertinent to critical wetting and also for cone filling which we will consider later.

Wedge filling can also be studied using a simple modification of the effective interfacial Hamiltonian and for rather open wedges corresponding to small $`\alpha `$ (where $`\mathrm{tan}\alpha \alpha `$) the appropriate model is

$$H[l]=𝑑x\left\{\frac{\mathrm{\Sigma }}{2}\left(l\right)^2+W(l\alpha |x|)\right\}$$
(3.3)
where $`l(x,y)`$ denotes the height of the interface measured above the plane. As shown by Rejmer et al this model can be justified from analysis of a more general drum-head model and with a binding potential interaction depending on the local normal distance of the interface to the wall rather than the simple vertical height used in (3.3). The simplifying approximations implicit in the above model do not effect any of the universal physics associated with the filling phase transition.

### 3.2 Mean-field Theory

At MF level the equilibrium interfacial profile follows from the solution to the Euler-Lagrange equation found from the minimization of (3.3)

$$\mathrm{\Sigma }\ddot{l}=W^{}(l\alpha |x|)$$
(3.4)
which is solved subject to the boundary conditions that $`\dot{l}(0)dl/dx|_0=0`$ and $`l(x)l_\pi \pm \alpha x`$ as $`x\pm \mathrm{}`$. This equation can be integrated once and leads to simple condition for the mid-point height :

$$\frac{\mathrm{\Sigma }\alpha ^2}{2}=W(l_0)W(l_\pi )$$
(3.5)
Note that in the limit $`l_0\mathrm{}`$, the RHS $`\mathrm{\Sigma }\theta _\pi ^2/2`$ which follows from the Youngs equation in the present small angle limit. The above equation has an elegant graphical interpretation  from which one can determine the order of the filling phase transition occurring as $`TT_F`$ at fixed $`\alpha `$ or equivalently as $`\alpha \theta _\pi `$ at fixed temperature. Thus the MF condition for first-order/critical wedge filling is that, at the filling temperature $`T_F`$, the binding potential does/does not possess a local maximum separating the minimum at $`l=l_\pi `$ from the extremum at infinity. Thus wedges made from walls that exhibit critical wetting, show critical filling whilst wedges made from walls that exhibit first-order wetting show both first order and critical filling depending on whether the filling temperature is above or below the spinodal temperature $`T_s`$ at which the first-order wetting binding potential develops a local maximum. Note that since the filling temperature can be lowered almost arbitarily by increasing the wedge angle $`\alpha `$ it follows that the wedge will typically show first-order filling if $`\alpha `$ is small and critical filling for larger values of $`\alpha `$. The tricritical value of at which the order of the filling transition changes character is given by $`\alpha ^{}=\theta _\pi (T_s)`$ the precise value of which depends on how weak the first-order wetting transition is. Of course this argument presupposes that a spinodal temperature for the binding potential exists. This is certainly true for systems with short-ranged forces where it follows from the Landau theory of wetting. For example in the global phase diagrams of Nakanishi and Fisher  the locus of spinodal temperatures (for different surface fields) corresponds to the analytic extension of the critical wetting curve. The prediction that wedges made from walls exhibiting first-order wetting show critical filling can certainly be tested in Ising model simulation studies. For experimental systems it may be that the mean-field condition for critical filling cannot be met because the spinodal temperature $`T_s`$ is too low. However as we shall see, fluctuation effects lead to a second mechanism for critical filling beyond the MF condition described above which may overcome this difficulty. Also, and perhaps of greater practical importance, we shall show that the mechanism for critical filling in a cone is rather different to that for a wedge and is much more easily fulfilled.

To continue, from (3.2) it immediately follows that at a critical filling transition the MF value of the order-parameter critical exponent is simply determined by the leading order decay of $`W(l)`$. The exponent $`q`$ plays no role and we find $`\beta _s^W=1/p`$ which is totally different to the corresponding MF critical exponent $`\beta _s1/(qp)`$ for critical wetting. The MF values of the other critical exponents follow from analysis of the height-height correlation function (with $`\stackrel{~}{y}y^{}y`$)

$$H(x,x^{};\stackrel{~}{y})\delta l(x,y)\delta l(x^{},y^{})$$
(3.6)
and its Fourier transform

$$S(x,x^{};Q)=𝑑\stackrel{~}{y}e^{iQ\stackrel{~}{y}}H(x,x^{};\stackrel{~}{y})$$
(3.7)
which exploits the translational invariance along the wedge. At MF level, $`S(x,x^{};Q)`$ follows from the solution to the Ornstein-Zernike equation

$$𝑑x^{\prime \prime }C(x,x^{\prime \prime };Q)H(x^{\prime \prime },x^{};Q)=\delta (xx^{})$$
(3.8)
where the direct-correlation function is defined by

$$C(x,x^{};Q)\frac{\delta ^2H[l]}{\delta l(x)\delta l(x^{})}$$
(3.9)
evaluated at the equilibrium height. Thus we arrive at the differential equation

$$(\mathrm{\Sigma }_x^2+\mathrm{\Sigma }Q^2+W^{\prime \prime }(l(x)\alpha |x|))S(x,x^{};Q)=\delta (xx^{})$$
(3.10)
To proceed it is best to consider the moment expansion

$$S(x,x^{};Q)\underset{n=0}{\overset{\mathrm{}}{}}Q^{2n}S_{2n}(x,x^{})$$
(3.11)
and first concentrate on the zeroth moment $`S_0(x,x^{})`$ which describes the position dependence of correlations in the fluctuations of the interfacial height across the wedge. This is found by using standard techniques to be

$`S_0(x,x^{})=\left(|\dot{l}(x)|\alpha \right)\left(|\dot{l}(x^{})|\alpha \right)\left\{{\displaystyle \frac{1}{2\alpha W^{}(l_0)}}+{\displaystyle \frac{H(xx^{})}{\mathrm{\Sigma }}}{\displaystyle _0^{\mathrm{min}(|x|,|x^{}|)}}{\displaystyle \frac{dx}{(\dot{l}(x)\alpha )^2}}\right\}`$
(3.12)
where $`H(x)`$ denotes the Heaviside step function $`(H(x)=1;x>0,H(x)=0;x<0)`$. Thus the structure of correlations across the wedge reflects the properties of the equilibrium height profile $`l(x)`$ found from the solution to the Euler-Lagrange equation (3.4), which is particularly simple. In essence, the interface is flat such that $`l(x)l_0`$, for $`xl_0/\alpha `$, while for $`xl_0/\alpha `$ the height decays exponentially quickly to its planar value $`l_\pi `$. Importantly, the lengthscale which controls the exponential decay is the wetting correlation length $`\xi _{}`$, which remains microscopic at the filling transition temperature. One consequence of this is that we can anticipate that the lateral width of the filled region of the wedge is trivally identified as $`\xi _x2l_0/\alpha `$ so that $`\nu _x=\beta _s^W`$. In this way we conclude that $`S_0(x,x^{})`$ is essentially position independent provided $`|x|,|x^{}|<\xi _x/2`$ implying that the whole filled portion of the wedge fluctuates coherently. Moreover that there is negligible correlation between height fluctuations inside the wedge with those outside and simlarly between positions on either side of the filled region. From (3.12) one can iteratively construct the position dependence of the higher moments using the integral equation

$$S_{2n+2}(x,x^{})=\mathrm{\Sigma }S_0(x,x^{\prime \prime })S_{2n}(x^{\prime \prime },x^{})𝑑x^{\prime \prime }$$
(3.13)
which follows from the Ornstein-Zernike equation above. From this it follows that in the asymptotic critical (scaling) regime $`S_{2n}(x,x^{})`$ is position independent for $`|x|,|x^{}|<\xi _x/2`$ with

$$S_{2n+2}(x,x^{})=\frac{\alpha }{2W(l_0)}\left(\frac{\mathrm{\Sigma }l_0}{W^{}(l_0)}\right)^n$$
(3.14)
and is negligible otherwise. This behaviour is consistent with the simple Lorentzian structure factor

$$S(x,x^{};Q)=\frac{S_0(x,x^{})}{1+\xi _y^2\xi _{}^2};|x|,|x^{}|<\xi _x/2$$
(3.15)
where $`S_0(0,0)\alpha /2W^{}(l_0)`$ and identifies the MF correlation length along the wedge $`\xi _y`$ as

$$\xi _y=\left(\frac{\mathrm{\Sigma }l_0}{W^{}(l_0)}\right)^{\frac{1}{2}}$$
(3.16)
Substituting for $`l_0`$ we find the desired MF critical exponent $`\nu _y`$ as $`\nu _y=1/p+1/2`$. The important observation here is that $`\nu _y\nu _x`$ so that the fluctuations are highly anisotropic and dominated by the modes along the wedge. Indeed as one approaches the bulk phase boundary $`\xi _y`$ becomes arbitrarily larger than $`\xi _x`$ implying that the fluctuations become pseudo-one dimensional. The last lengthscale to calculate within the MF analysis is the midpoint roughness defined by $`\xi _{}^2[l(0,y)l_0]^2`$ which can be identified from the relation

$$\xi _{}^2=H(0,0;0)$$
(3.17)
where the RHS can be obtained from the Fourier inversion of $`S(0,0;Q)`$. This leads to the roughness relation for the $`d=3`$ wedge

$$\xi _{}\left(\frac{\xi _y}{\mathrm{\Sigma }l_0}\right)^{\frac{1}{2}}$$
(3.18)
which is one of the central results of the MF analysis. Substituting for the divergence of $`\xi _y`$ and $`l_0`$ we find that all dependence on the binding potential exponent $`p`$ cancels leading to a remarkable universal divergence

$$\xi _{}t^{\frac{1}{4}}$$
(3.19)
valid for all ranges of intermolecular potential. In summary, the MF values of the lenghscale critical exponents at critical filling are

$$\beta _s^W=\frac{1}{p},\nu _y=\frac{1}{2}+\frac{1}{p},\nu _{}^W=\frac{1}{4}$$
(3.20)
Next we turn to the fluctuation theory of filling beyond MF approximation.

### 3.3 Fluctuation Theory: breather modes and their excitations

The MF analysis presented above is only valid if the fluctuations are small in some sense. Since the divergence of $`\xi _{}`$, even within MF, is algebraic we can simply use the ‘contact condition’  that $`\xi _{}`$ cannot be greater than $`l_0`$ to determine the Ginzburg criterion for the breakdown of MF theory. Thus we can assert that the above critical exponents are valid provided

$$\frac{\xi _{}}{l_0}t^{(\frac{1}{4}\frac{1}{p})}1$$
(3.21)
which in the limit $`t0`$ implies that $`p`$ must be less then $`4`$. For $`p>4`$, fluctuation effects dominate and we can anticipate (analogous to the WFL and SFL regime behaviour of critical wetting) that the roughness $`\xi _{}`$ is comparable with the interfacial height $`l_0`$. Unfortunately it is very difficult to develop a full fluctuation theory of filling transitions based on the effective interfacial Hamiltonian (3.3) since, as is clear from the above remarks, the upper critical dimensions for wetting and filling are quite different. The situation is rather similar to that encountered in the theory of wetting where one would ideally like to develop a fluctuation theory based on analysis of, say, a semi-infinite Ising or Landau-Ginzburg Wilson (LGW) model. Such a approach would afford a description of both bulk (and surface) critical behaviour and also of wetting. Since this has not proved possible one instead focuses on effective models constructed to capture the essential physics of the phase transition in question. Of course, all Hamiltonians are effective Hamiltonians but some are more effective than others. The LGW model based on a local magnetization order-parameter can (in principle) describe bulk and surface criticality in addition to wetting and filling. In turn the effective interfacial Hamiltonian (3.3) based on a collective co-ordinate $`l(x,y)`$ is a theory of wetting and filling. The task here is to develop a simpler model than (3.3) which just affords a description of filling.

The theory of wedge filling which we propose is constructed from the full interfacial model under the assumption that the only fluctuations that determine the asymptotic critical behaviour arise from the pseudo-one dimensional long-wavelength undulations in the local height of the filled region. We refer to fluctuations in the local height as breather modes since these increase/decrease the local width of the filled region (and the volume of liquid beneath it). Thus we seek to derive a model based on a collective co-ordinate $`l_0(y)`$ which only describes these modes (see Fig 3.1) and which is valid for sufficiently small wave-vectors $`Q1/\xi _x`$. Such a model can be easily constructed by generalising the method introduced by Fisher and Jin  for the derivation of wetting interfacial models from the LGW Hamiltonian. Assuming that for a given constrained configuration of the mid-point height $`l_0(y)`$ that all other fluctuations are small (relative to the soft mode), it follows that the desired filling Hamiltonian can be identified by

$$\mathrm{exp}(H_W[l_0])=𝒯\text{r}\left\{\mathrm{exp}(H[l])\right\}$$
(3.22)
where the $`𝒯\text{r}`$ denotes a partial trace with respect to the generalised crossing-criterion (GCC) constraint that $`l(x=0,y)=l_0(y)`$. Since this trace is over non-critical fluctuations one can employ a saddle-point idenifiction

$$H_W[l_0(y)]=H[l^{}(x,y;l_0())]$$
(3.23)
where $`l^{}(x,y;l_0())`$ denotes the profile that minimises (3.3) subject to the GCC. As in the Fisher-Jin theory of wetting, the desired non-planar constrained configurations can be constructed perturbatively in terms of constrained profiles that are translationally invariant along the wedge (analogous to the planar constrained magnetization profiles) which we simply write as $`l^{}(x;l_0)`$ to avoid a proliferation of indices. These translationally invariant profiles satisfy the Euler-Lagrange (3.4) subject to the GCC $`l^{}(x=0;l_0)=l_0`$. The resulting wedge filling model has the form

$$H_W[l_0]=𝑑y\left\{\frac{\sigma (l_0)}{2}\left(\frac{dl_0}{dy}\right)^2+V_W(l_0)\right\}$$
(3.24)
where $`V_W(l_0)`$ denotes the wedge binding potential and $`\sigma (l_0)`$ is a position dependent wedge line stiffness. Formally these quantities are determined by the relations

$$V_W(l_0)=_{\mathrm{}}^{\mathrm{}}𝑑x\left\{\frac{\mathrm{\Sigma }}{2}\left(\frac{dl^{}(x;l_0)}{dx}\right)^2+W(l^{}(x;l_0)\alpha |x|)\right\}$$
(3.25)
(up to unimportant additive constants defined such that $`V_W(l_\pi )=0`$) and

$$\sigma (l_0)=\mathrm{\Sigma }_{\mathrm{}}^{\mathrm{}}𝑑x\left(\frac{l^{}(x;l_0)}{l_0}\right)^2$$
(3.26)
The planar constrained profiles $`l^{}(x;l_0)`$ have the same simple near-flat structure as the equilibrium MF profile discussed earlier. This means that the evaluation of (3.25, 3.26) is particularly simple so that for first-order and critical wedge filling we can approximate

$$H_W[l_0]=𝑑y\left\{\frac{\mathrm{\Sigma }l_0}{\alpha }\left(\frac{dl_0}{dy}\right)^2+V_W(l_0)\right\}$$
(3.27)
The qualitative and quantitative form of the filling potential depends on the order of the MF phase transition. For critical filling it has the expansion

$$V_W(l_0)=\frac{\mathrm{\Sigma }(\theta _\pi ^2\alpha ^2)l_0}{\alpha }+a_Fl_0^{1p}+\mathrm{}$$
(3.28)
and of course fluctuations are restricted to the region $`l_0>0`$ analogous to the hard-wall condition for wetting. Here the effective Hamaker constant $`a_Fa/(p1)\alpha `$ and note that the coefficient of the linear term is simply proportional to $`t`$, the relevant reduced temperature. The critical filling potential shows a single minimum (located precisely at the MF value of $`l_0`$) which continuously moves out to infinity as $`t0`$. The one-dimensional filling Hamiltonian can be studied using a number of different techniques including transfer-matrix theory  and approximate RG methods . However there is no need to go to these lengths to extract the fluctuation regimes and the values of the critical exponents since they follow from a simple extension of the Lipowsky-Fisher scaling theory considered earlier. To see this we suppose that the influence of fluctuations at critical filling follows from analysis of an effective potential

$$V_{\text{eff}}(l_0)V_W(l_0)+V_{\text{fl}}(l_0)$$
(3.29)
which accounts for the entropic repulsion of the unbinding interface from the wedge bottom. Again the fluctuation term can be estimated from both the form of the bending energy contribution to (3.27) and also the number of collisions the interface has with the wedge bottom. This we estimate to of order $`l_0^{32/\zeta _F}`$ where $`\zeta _F`$ denotes the thermal wandering exponent for wedge filling (in $`d=3`$) relating the roughness and pertinent correlation length $`\xi _{}\xi _y^{\zeta _F}`$. The value $`\zeta _F=1/3`$ follows from considering the invariance of the free part of the filling Hamiltonian under the scale transformation $`yy/b`$ and $`l_0b^{\zeta _F}l_0`$. Thus we arrive at an effective filling potential of the form

$$V_{\text{eff}}(l_0)V_W(l_0)+Cl_0^3$$
(3.30)
which has a single minimum at the equilibrium mid-point height. In this way we find two different fluctuation regimes for filling depending on whether the direct intermolecular potential or entropic contribution determines the next-to-leading order correction to the linear term in (3.28). The leading linear term is always relevant (in the RG sense) and reflects the influence of the wedge geometry on the statistical/thermodynamic properties of filling. This is rather similar to the theory of fluctuation effects at complete wetting transitions . In this way we predict there exists a filling MF regime (FMF) for $`p<4`$, with exponents given by (3.20) and a filling fluctuation dominated regime (FFL) for $`p>4`$ with

$$\beta _s^W=\frac{1}{4},\nu _y=\frac{3}{4},\nu _{}^W=\frac{1}{4};FFL$$
(3.31)
As a final remark we turn to the effective one-dimensional Hamiltonian description of first-order filling. This is based on the model (3.27) but with a different potential $`V_W(l_0)`$ the form of which is shown in Fig. 3.3. In fact the structure of $`V_W(l_0)`$ for this case can be qualitatively deduced by noticing that the extrema of $`V_W(l_0)`$ must correspond to the possible solutions of the MF equation (3.5). For $`\theta _\pi >\alpha `$, $`V_W(l_0)`$ has only a single minimum which remains finite at $`T=T_F`$. However for $`\theta _\pi <\alpha `$ a potential barrier appears separating the local minimum and the infinitely deep potential well at infinity. Note that exactly at the first-order filling transition temperature the minimum of $`V_W(l_0)`$ has a finite well-depth relative to $`V_W(\mathrm{})`$. This is quite unlike the behaviour of $`W(l)`$ at first-order wetting transitions. These remarks indicate even if the mean-field condition for critical filling is not fulfilled, the influence of the one-dimensional breather mode fluctuations may still drive the transition critical by causing the interface to tunnel away from the local minimum of the potential. Whether this happens or not depends on the range and depth of the minimum of $`V_W(l_0)`$ but will certainly be most important for systems with short-ranged forces (such as Ising models). Note this tunnelling mechanism does not alter the filling transition phase boundary (3.1) which simply reflects the relevance of the linear term in $`V_W(l_0)`$. One may also note that this tunnelling must definitely occur out of coexistence for temperatures $`T>T_F`$ where it is responsible for the washing out of the pre-filling line.

We mention only briefly here that it is trivial to investigate critical wedge filling from off-bulk co-existence using the effective model by allowing for a term of order $`l^2`$ which couples to the undersaturation in the filling potential. This clearly represents the contribution from the volume of undersaturated liquid in the filled region of the wedge. From this one can readily establish that the filling thickness, roughness and correlation lengths all exhibit scaling with a gap exponent $`\delta =5/4`$ in the FFL regime.Thus the divergence of the mean interface height along the critical wedge isotherm as one approaches saturation is characterised by a $`1/5`$ power law.

### 3.4 Wedge Filling in two dimensions

It is a simple matter to extend the wedge filling fluctuation theory to arbitrary dimensions $`d`$ and in particular $`d=2`$ where one can compare the predictions with the exact transfer matrix solution of the full interfacial model (3.3). Consider a wedge in bulk $`d`$ dimensions which is of ‘V’ shape cross-section in the $`x`$ direction but which is translationally invariant in the remaining $`d2`$ (surface) dimensions. Denoting the position of the interface above the wedge bottom by $`l(\stackrel{}{x}_{})`$ where $`\stackrel{}{x}_{}`$ denotes the $`(d2)`$ dimensional vector parallel to the wedge orientation. The filling Hamiltonian trivially generalises to

$$H_W[l_0]=𝑑\stackrel{}{x}_{}\left\{\frac{\mathrm{\Sigma }l_0}{\alpha }\left(l_0\right)^2+V_F(l_0)\right\}$$
(3.32)
which clearly describes the $`(d2)`$ dimensional thermal excitations of the breather interfacial mode. The scaling theory also generalises in a natural way. First we calculate the $`d`$ dimensional wedge wandering exponent $`\zeta _W`$ by considering the invariance of the free part of (3.32) under scale transformations. This determines $`\zeta _W`$ as

$$\zeta _W=\frac{4d}{3}$$
(3.33)
relating the roughness and correlation length measured along the wedge (which we here denote $`\xi _y`$) via $`\xi _{}\xi _y^{\zeta _W}`$. Similarly the fluctuation contribution to the effective potential arising from the bending energy and entropy terms scales as

$$V_F(l_0)l_0^{\frac{3(2d)}{4d}}l_0^{\tau _W}$$
(3.34)
where $`\tau _W`$ is the wedge fluctuation exponent analogous to $`\tau `$ in the Lipowsky-Fisher theory of wetting.

This heuristic scaling theory predicts that for general dimension $`d<4`$, the criticality at filling falls into a mean-field (FMF) and fluctuation dominated (FFL) regime depending on whether $`p`$ is less or greater than the marginal value $`p^{}=2(d1)/(4d)`$. The critical exponents are given by

$$\beta _s^W=\frac{1}{p},\nu _y=\frac{1}{p}+\frac{1}{2},\nu _{}^W=\frac{4d}{4}+\frac{(3d)}{2p};FMF$$
(3.35)
$$\beta _s^W=\nu _{}^W=\frac{(4d)}{2(d1)},\nu _y=\frac{3}{2(d1)};FFL$$
(3.36)
and notice that it is only in $`d=3`$ that the roughness exponent is universal independent of the range of the forces. These predictions for critical singularities at wedge filling are, in general very different to the predictions for critical wetting. However for the special case of $`d=2`$ there appears to be a connection between filling and critical wetting discussed below. In $`d=2`$ the breather fluctuation theory predicts that the marginal intermolecular range is $`p=1`$ and that the values of the FFL regime critical exponents are

$$\beta _s^W=\nu _{}^W=1;d=2FFL$$
(3.37)
with $`\nu _y`$ undefined. Note that these are pertinent to all physical examples of intermolecular forces. The first thing to observe here is that the numerical values of the filling critical exponents $`\beta _s^W`$ and $`\nu _{}^W`$ in the $`d=2`$ FFL are identical to the analogous critical wetting exponents in the ($`d=2`$) SFL regime. As we shall now see this coincidence extends to the full form of the PDF.

As well as predicting the critical exponents at filling the fluctuation theory also predicts the position dependence of the scaling form of the mid-point PDF, $`P_W(l_0)`$. To see this note that the case of wedge filling in $`d=2`$ is special because the breather mode is only subject to zero dimensional fluctuations since there is no direction along the wedge. As a result one can go beyond the above predictions for the values of the critical exponents and assert that in the asymptotic critical region, the mid-point PDF is given by

$$P_W(l_0)=Ne^{V_F(l_0)}$$
(3.38)
where here (and below) $`N`$ denotes a suitable normalization factor. Now from (3.28) we have $`V_F(l_0)\mathrm{\Sigma }(\theta _\pi ^2\alpha ^2)l_0/\alpha 2\mathrm{\Sigma }(\theta _\pi \alpha )l_0`$ so that in the critical regime the breather fluctuation theory leads to the scaling form of the PDF

$$P_W(l_0)=Ne^{2\mathrm{\Sigma }(\theta _\pi \alpha )l_0}$$
(3.39)
or equivalently

$$P_W(l_0)=\frac{1}{l_0}e^{\frac{l}{l_0}}$$
(3.40)
which is of course consistent with the values of the critical exponents quoted in (3.36). This prediction is remarkable for two reasons: (I) it is in exact agreement with the transfer matrix result for the PDF based on the ($`d=2`$) interfacial model (3.3). For quite arbitrary choices of binding potential the mid-point interface height distribution function is given by

$$P_W(l_0)=Ne^{2\mathrm{\Sigma }\alpha l_0}|\mathrm{\Psi }_0^{(\pi )}(l_0)|^2$$
(3.41)
where $`\mathrm{\Psi }_0^{(\pi )}(l)`$ is the ground state eigenfunction appearing in the transfer matrix spectrum for $`d=2`$ wetting at a planar wall. This satisfies the well known Schrödinger equation

$$\frac{1}{2\mathrm{\Sigma }}\frac{^2}{l^2}\mathrm{\Psi }_0^{(\pi )}(l)+W(l)\mathrm{\Psi }_0^{(\pi )}(l)=E_0\mathrm{\Psi }_0^{(\pi )}(l)$$
(3.42)
where $`E_0`$ is the corresponding groundstate eigenvalue which is related to the contact angle by $`E_0=\mathrm{\Sigma }\theta _\pi ^2/2`$. From the properties of $`\mathrm{\Psi }_0^{(\pi )}(l)`$ it follows that the filling transition is indeed characterised by the scaling regimes and corresponding critical singularities predicted by the general fluctuation theory. Moreover for $`p>1`$, and far from the wall the wave function decays as

$$\mathrm{\Psi }_0^{(\pi )}(l)e^{\mathrm{\Sigma }\theta _\pi l}$$
(3.43)
so that (3.39) and (3.41) are identical. The equivalence of the scaling form of the PDF with the breather mode fluctuation theory and the full effective interfacial model completely supports our conjecture that it is the overall breather mode which translates the interface up and down the sides of the wedge, rather than the thermal wandering of the interface within the filled region that determines the criticality at filling.

(II) The derived expression for the filling PDF in the FFL regime has an identical structure to the analogous interfacial height PDF, $`P_\pi (l)`$ for critical wetting in the SFL. Recall that this quantity is determined as $`P_\pi (l)|\mathrm{\Psi }_0^{(\pi )}(l)|^2`$ with the transfer matrix approach so that we in the SFL we have the well known result, quoted earlier

$$P_\pi (l)=\frac{1}{l_\pi }e^{\frac{l}{l_\pi }}$$
(3.44)
which is identical to the result for filling. In fact in $`d=2`$ the connection between filling and critical wetting for purely thermal disorder extends to the case marginal case $`p=1`$ for filling where the structure of the FMF/FFL PDF has the same form as the PDF for critical wetting at a WFL/MF borderline for which the specific heat exponent vanishes . The PDF for this case is not simply a pure exponential function.

The equivalence of the film thickness (and roughness) critical exponents and PDF’s for two dimensional, fluctuation-dominated filling and wetting has also been established for the case of filling/wetting with random bond disorder described by the interfacial model

$$H[l]=\left\{\frac{\mathrm{\Sigma }}{2}(\frac{l}{x})^2+V_r(x,l(x))+W(l\alpha |x|)\right\}$$
(3.45)
with disorder averages satisfying

$$\overline{V}(x,l(x))=0$$
(3.46)
and

$$\overline{V_r(x,l(x))V_r(x^{},l^{}(x))}=\mathrm{\Delta }\delta (xx^{})\delta (l(x)l^{}(x^{}))$$
(3.47)
where $`\mathrm{\Delta }`$ is a measure of the disorder. The study of filling transitions in this model allows us to test the possible connection between fluctuation dominated two dimensional filling and critical wetting in a system where the interfacial wandering exponent is different to the thermal case. It transpires that for systems with purely short-ranged forces the model (3.45) can be solved exactly using the transfer matrix method by extending Kardar’s Bethe Ansatz approach  used in his analysis of wetting. This again confirms the thermodynamic prediction for the location of the filling transition and identifies the universal filling critical exponents as

$$\beta _s^W=2,\nu _{}^W=2$$
(3.48)
which are identical to the analogous critical exponents $`\beta _s`$ and $`\nu _{}`$, determined by Kardar for critical wetting with random bonds with short-ranged forces corresponding to the SFL regime . The expression for the scaling form of the mid-point interfacial height PDF $`P_F(l_0)`$ for $`d=2`$ random bond filling is highly non-trivial but again analysis shows that it is identical to the PDF for random bond wetting .

## 4 Conic Filling.

### 4.1 Phenomenology.

Consider a non-planar wall-fluid interface in the shape of an infinite cone which makes an angle $`\alpha `$ to the horizontal. Thus, for open wedges (small $`\alpha `$) the height of the wall above the plane is described by a height function $`z_C(\stackrel{}{r})=\alpha \stackrel{}{r}`$ with $`\stackrel{}{r}`$ the displacement vector along the two-dimensional displacement vector along the horizontal plane with origin at the cone apex. As before the wall is supposed to be in contact with a vapour phase at two-phase coexistence with thermodynamics indicating that the cone is completely filled for $`T>T_F`$, where the conic filling temperature $`T_F`$ satisfies the same condition $`\theta _\pi (T_F)=\alpha `$ as filling in the $`d=3`$ wedge (3.1). We allow for the possibility of first-order and critical cone filling corresponding to the discontinuous or continuous divergence of the mid-point filling height $`l_0`$ as $`t(T_FT)/T_F0`$. Also of interest is the roughness $`\xi _{}`$ arising from the interfacial fluctuations of the unbinding liquid vapour interface. The radius of the near flat, filled region of the cone, which also characterises height fluctuations in the central region is trivially related to the height via $`\xi _rl_0/\alpha `$. As with wedge filling we anticipate that for radial positions $`r\xi _r`$ the interface will simply take the microscopic value $`l_\pi `$ it would have on a planar wall. Thus we define two lengthscale critical exponents for critical cone filling by

$$l_0t^{\beta _s^C},\xi _{}t^{\nu _{}^C}$$
(4.1)
in an obvious notation. Here we seek to understand both the conditions under which critical cone filling can occur and the values of the critical exponents at MF level and beyond. As we shall show these are quite distinct from the situation for wedge filling.

Following our earlier treatment of wedge filling we anticipate that interfacial fluctuations manifest themselves in two different ways at conic filling. Firstly there is the usual wandering of the interface within the filled region determined by the value of the standard wandering exponent $`\zeta _\pi `$ which is of course marginal in $`d=3`$ (concentrating on thermal fluctuations). Such fluctuations would occur even if the edges of the filled region were fixed to a given height, $`l_0`$ say, and would contribute terms of order $`\sqrt{l_0}`$ and $`\sqrt{\mathrm{ln}l_0}`$ to the interfacial roughness in $`d=2`$ and $`d=3`$ respectively. However much more important than this is the breather mode that translates the entire flat region up and down the cone. The fundamental distinction between three-dimensional wedge and cone filling arises because the breather mode is thermally excited by quasi-one dimensional and zero-dimensional excitations respectively. In this respect three dimensional conic filling is the higher dimensional version of two-dimensional filling considered in the last section. Note that the effective dimensionality of the breather mode does not destroy the phase transition because the fluctuating field diverges as the filling temperature is approached.

### 4.2 Mean-field theory.

Denoting the local interfacial height above the plane by $`l(\stackrel{}{r})`$, the natural starting point for the study of conic filling for open cones (small $`\alpha `$) is the standard interfacial model

$$H[l]=𝑑\stackrel{}{r}\left\{\frac{\mathrm{\Sigma }}{2}\left(l\right)^2+W(l\alpha r)\right\}$$
(4.2)
where, again, $`\stackrel{}{r}`$ denotes the radial vector from the apex. To begin we chose binding potentials $`W(l)`$ of the type (2.2) corresponding to walls that exhibit critical wetting and set the exponent parameter $`q=p1`$. Numerical minimization of the model Hamiltonian for different ranges of the intermolecular potential shows that the filling transition is critical (for the present choice of critical wetting walls) and located precisely at $`\theta _\pi =\alpha `$. However the numerics also indicates that the MF value of the critical exponent characterising the divergence of $`l_0`$ appears to fall into two distinct regimes. For $`p>1`$ the critical behaviour is universal and dominated by the geometry such that $`\beta _s^C=1`$. On the other hand if $`p<1`$ the critical exponent is consistent with the identification $`\beta _s^C=1/p`$, identical to that found within the MF theory of wedge filling. This intriguing MF critical behaviour follows from a very simple and highly accurate variational theory of conic filling. Noting the simplicity of the equilibrium interfacial profiles obtained in our numerical minimization we parameterise the possible height function by

$$l^{}(\stackrel{}{r})=\{\begin{array}{cc}l_0,\hfill & r<\frac{l_0l_\pi }{\alpha }\hfill \\ l_\pi +\alpha r,\hfill & r>\frac{l_0l_\pi }{\alpha }\hfill \end{array}$$
(4.3)
and determine the mid-point height dependence of the conic filling potential defined by

$$V_c(l_0)=H[l^{}]H[l^{}]|_{l^{}(0)=l_\pi }$$
(4.4)
Note that this variational approximation becomes increasingly accurate for as one approaces a critical filling transition since the profile tends to a macroscopically flat meniscus. Minimization of the conic filling potential leads to the following integral equation for the mid-point height

$$\frac{\mathrm{\Sigma }(\alpha ^2\theta _\pi ^2)(l_0l_\pi )}{2}=_{l_\pi }^{l_0}W(l)𝑑l$$
(4.5)
This equation has an simple graphical interpretation which recovers the numerical results quoted above and leads to an elegant criterion for the order of the filling transition. The filling transition is continuous provided

$$_{l_\pi }^{\mathrm{}}W(l)𝑑l<0$$
(4.6)
corresponding to an integral measure of the overall attractive/repulsive nature of the binding potential. By allowing for pertubations to the simple variational ansatz it is easy to show that this is the exact MF condition for criticality within the effective interfacial model to order $`𝒪(\alpha ^4)`$. The condition for criticality of conic filling is much less severe than that for both critical wetting and critical wedge filling and implies that even substrates which exhibit rather strong first order-wetting will exhibit a continuous filling transition in a cone geometry. A comparison of the respective shapes of the binding potential need for criticality for three dimensional wetting, wedge filling and cone filling is shown in Fig. 4.1. Provided the condition for critical cone filling is met the MF critical exponent falls into two regimes:

$$\beta _s^C=\frac{1}{p}\text{ for }p<1\text{ and }\beta _s^C=1\text{ for }p>1$$
(4.7)
The universality of the critical behaviour in the geometry dominated critical regime is reminiscient of the way in which fluctuation effects lead to universal critical behaviour at wetting transitions. This connection is explored further below.

### 4.3 Fluctuation theory: breather modes and the PDF.

The effective zero-dimensionality of the breather mode at conic filling suggests that we should analyse fluctuation effects in the same manner as we did for filling of a wedge/cone in $`d=2`$. Recall that this simple approach gave the precise scaling form of the mid-point interfacial height PDF (and hence the critical exponents) at filling as that found from the exact transfer-matrix theory. Accordingly we propose that the critical behaviour at conic filling is described by a PDF for the mid-point interfacial height is of the form

$$P_C(l_0)=Ne^{V_C(l_0)}$$
(4.8)
since the role of the cone potential in $`d=3`$ is analogous to the wedge potential $`V_W(l_0)`$ in $`d=2`$. It is a simple matter to show that $`V_c(l_0)`$ has the form

$$V_C(l)=\frac{\mathrm{\Sigma }(\theta _\pi ^2\alpha ^2)l^2}{\alpha }\sigma _Cl+a_Cl_0^{2p}+\mathrm{}$$
(4.9)
consistent with the MF values of the critical exponents discussed above. The three terms in this expression represent the purely geometrical contributions arising form the area and circumference of the filled region and also the direct influence of the specific intermolecular forces. Both $`\sigma _C`$ and $`a_c`$ are positive effective Hamaker constants. For $`p>1`$ the final term in (4.8) is irrelevant, geometry dominates and the PDF can be written

$$P_C(l_0)=Ne^{\frac{1}{2\xi _{}^2}(l_0l_0)^2}$$
(4.10)
and identifies the universal values of the critical exponents as

$$\beta _s^C=1\text{ and }\nu _{}^C=\frac{1}{2}\text{ for }p>1$$
(4.11)
On the other hand if the second term in (4.9) is irrelevant and close to its most probable position the PDF is Gaussian with critical exponents given by

$$\beta _s^C=\frac{1}{p}\text{ and }\nu _{}^C=\frac{1}{2}\text{ for }p<1$$
(4.12)
Equations (4.6) and (4.8-4.12) are the main predictions of this section. Here we make two pertinent remarks:

(I) the MF values of the filling height critical exponent $`\beta _s^C`$ is unchanged by fluctuations. However the thermal fluctuations of the breather mode do lead to an appreciable interfacial roughness. The predictions (4.11) are appropriate for both experimental systems with long-ranged dispersion forces and simulation studies of Ising-like systems with short-ranged forces. For discrete lattice based Ising studies, it is probably more convenient to study adsorption in an inverted pyramid shaped geometry. The critical exponents are the same as the cone for this case.

(II) the scaling form of the PDF and the values of the critical exponents in the universal geometry dominated regime of cone filling are identical to the analogous critical singularities characteristic of the SFL regime of critical wetting short-ranged forces i.e. regime (III) corresponding to $`\omega >2`$. This is directly analogous to the observations made earlier for wedge/cone filling in $`d=2`$.

The above predictions point to a remarkable coincidence between the universality of conic filling in $`d=2`$ and $`d=3`$ and the SFL regime of critical wetting in the same dimension which is encapsulated in the equivalence of the film thickness and roughness critical exponents as well as the interfacial height PDF’s. Thus in the important geometry dominated regimes of cone filling and for the physically relevant dimensions $`d=2`$ and $`d=3`$ one can identify

$$\beta _s^C=\beta _s^{SFL},\nu _{}^C=\nu _{}^{SFL}\text{ and }P_C(l_0)=P_\pi ^{SFL}(l_0)$$
(4.13)
### 4.4 Universal conic filling and the SFL regime.

The above predictions beg the question, why does the influence of the cone geometry on the (universal) unbinding at filling mimic the fluctuations characteristic of critical wetting in the SFL? Usually, of course, if two different systems exhibit the same critical exponents one concludes, after identifying the appropriate order-parameter, that they belong to the same universality class. We do not believe that this is the correct explanation because in this case there does not appear to be an analogue of the wetting correlation length critical exponent $`\nu _{}`$ for filling. In other words the equivalence of the one-point PDF’s does not extend to the two point function. In view of this we forward a simple plausibility argument suggesting that the equivalence of the PDF’s is not simply coincidence but is special to the physical dimensions $`d=2`$ and $`d=3`$. We emphasise that the specific results of our model calculations for the values of the critical exponents stand quite apart from this discussion.

Let us focus on the structure of the PDF for critical wetting at a planar wall in $`d=2`$ and $`d=3`$. The distribution function $`P_\pi (l)`$ is found by evaluating the (normalised) Boltzmann sum over all interfacial configurations that pass through the specified height $`l`$ at an arbitarily chosen position $`\stackrel{}{r}=0`$ (say), equivalent to the identification $`P_\pi (l)=Ne^{F^{}[l]}`$ where $`F^{}=\mathrm{ln}Z_1^{}`$ and $`Z_1^{}`$ is the partial partition function summing over this class of configurations. In turn the constrained free-energy $`F^{}[l]`$ is a functional of the constrained height profile corresponding to the (partial) average over all profiles that satisfy the appropriate condition $`l(x=0)=l_0`$. For large $`ll_\pi `$ the constrained profile determining the PDF naturally forms the shape of a triangle (in $`d=2`$) and cone (in $`d=3`$) characterised by the equilibrium contact angle $`\theta _\pi `$. Thus the droplet shape of the constrained profile determining the asymptotics of $`P_\pi (l)`$ near critical wetting is essentially identical to the equilibrium height profile measured with respect to the local wall height near critical cone filling (since $`\theta _\pi \alpha `$). This observation hints at a possible connection between the PDF’s at wetting and cone filling. One may also reason that if there is any connection between the PDF’s of filling and wetting then it will only be valid when the direct effect of intermolecular interactions is unimportant since the non-universal influence of these forces manifests itself in very different ways. This is only true in the geometry dominated regime of cone filling and the SFL regime of critical wetting.

On the other hand it is easy to see that any connection between the two PDF’s is not valid for all dimensions. The generalisation of our simple fluctuation theory of cone filling presented to arbitary dimension $`d`$ is immediate. Using the same variational approach described above we arrive at a cone binding potential of the form (neglecting unimportant constants)

$$V_C(l)=tl^{d1}l^{d2}+\mathrm{}$$
(4.14)
representing the leading order surface and line contributions. The influence of intermolecular forces is present through a term $`𝒪(l^{d1p})`$ which is higher-order for $`p>1`$ independent of the dimensionality. Using this cone potential to evaluate the PDF arising from the (zero-deimensional) thermal fluctuations of the breather mode leads to the following prediction for the universal, geometry dominated cone filling critical singularities in arbitrary dimension $`2<d<4`$,

$$\beta _s^C=1\text{ and }\nu _{}^C=\frac{4d}{2}$$
(4.15)
which of course recovers our results for wedge/cone filling in $`d=2`$ and cone filling in $`d=3`$ as special cases. For $`d>4`$ the roughness remains finite. Only for the physically important dimensions $`d=2,3`$ to the numerical values of these critical exponents coincide with the predictions for SFL regime critical wetting. Nevertheless we find it remarkable that for these special dimensions not only the values of the critical exponents but also the detailed form of the PDF’s are identical for filling and critical wetting. Finally we mention that the exponents quoted above allow us to identify the value of the wandering exponent for critical cone filling defined by $`\xi _{}\xi _r^{\zeta _C}`$. Noting that interfacial height is trivially related to the radial correlation length (width) of the filled region we are led to the identification

$$\zeta _C=\frac{4d}{2}$$
(4.16)
which shows the influence of the zero-dimensional breather mode on the distribution of matter. Obviuously in two dimensions $`\zeta _C=\zeta _W=1`$ since the cone and wedge geometries are identical.

## 5 Conclusions.

In this article we have presented details of a general fluctuation theory of filling transitions at wedges and cones. The main results and conclusions of our study are summarised below:

(I) The dominant interfacial fluctuations at wedge and cone filling correspond to breather-like motion that translates the flat, filled region up and down the sides of the confining geometry.

(II) For three dimensional wedge filling the thermal excitations of the breather mode are described by a novel one-dimensional effective Hamiltonian which predicts a universal value of $`1/4`$ for the roughness critical exponent for all ranges of the intermolecular forces and also large scale interfacial roughness and universal critical singularities for $`p>4`$.

(III) For filling in $`d=2`$, the simple picture of an effective zero-dimensional fluctuating breather mode correctly captures the critical exponents and precise scaling behaviour of the mid-point height PDF as found from the full transfer-matrix solution to the full interfacial model (of wetting and filling).

(IV) For two-dimensional filling with both thermal and random bond disorder the scaling of the mid-point height PDF is identical to the corresponding results for SFL regime critical wetting in this dimensionality.

(V) Conic filling shows two different types of critical behaviour including a universal geometry dominated regime for $`p>1`$ pertinent to all physical types of intermolecular force. Applying the simple breather-like fluctuation picture to this phase transition yields a scaling form for the PDF and associated interfacial height and roughness critical exponents identical to the SFL regime prediction for critical wetting in this dimensionality. We have argued that this equivalence of the PDF’s is special to $`d=2`$ and $`d=3`$.

(VI) Even at MF level the conditions under which three dimensional wedge and cone filling are continuous are different to the fine tuning required for critical wetting and are sensitive to the particular geometry. For critical wedge wetting one requires that the binding potential shape at $`T=T_F`$ shows no local maximum. On the other hand cone-filling is critical provided the integrated weight of the potential remains negative (attractive).

These conclusions, regarding both the conditions for criticality and the manifestations of interfacial fluctuations at wedge and cone filling, are in sharp contrast to the situation for critical wetting. Concentrating on the values of the thickness and roughness critical exponents as pertinent to realistic three-dimensional systems with non-retarded van der Waals forces ($`p=2`$) we have

$$\beta _s=1\text{ and }\nu _{}=0(\mathrm{ln})\text{ for critical wetting}$$
(5.1)
$$\beta _s^W=\frac{1}{2}\text{ and }\nu _{}^W=\frac{1}{4}\text{ for critical wedge filling}$$
(5.2)
$$\beta _s^C=1\text{ and }\nu _{}^C=\frac{1}{2}\text{ for critical cone filling}$$
(5.3)
indicating that the geometry has a profound influence on the roughness of the interface. This reflects the special role of the breather mode at filling transitions. Consequently if critical wetting at a wall-fluid interface were ever seen experimentally, the observation of wedge and cone filling in this system would certainly provide a means of observing large scale interfacial fluctuations. It should certainly be possible to test the predictions of universal and large scale interfacial fluctuation related behaviour in Ising model systems. We have also emphasised that the geometry dramatically alters the conditions for criticality and certainly relaxes the requirement of fine tuning the values of the Hamaker constants necessary for critical wetting. The condition for inducing continuous filling transitions is easiest to fulfil for the cone geometry and is highly encouraging that the observation of critical filling in cones made from walls exhibiting first-order wetting is experimentally feasible. It may even be possible to further manipulate the conditions for criticality by considering filling in other geometries and by introducing chemical heterogeniety. Throughout this article we have focused exclusively on the filling transition occuring as $`TT_F`$ at bulk two-phase coexistence. Also of interest is the influence of the geometry on complete wetting adsorption isotherms corresponding to the continuous divergence of the filling thickness for temperatures $`T>T_F`$ as $`\mu \mu _{\text{sat}}`$ . These transitions have also received much recent attention and show great sensitivity to the shape of the confining geometry. For cones and wedges in all dimensions the divergence of the filling height is super-universal and given by

$$l_0\frac{1}{\delta \mu }\text{ and }\xi _{}\frac{1}{\delta \mu ^{\zeta _\pi }}$$
(5.4)
where $`\delta \mu `$ denotes the undersaturation. However unlike critical filling the breather mode does not change the nature of the roughness so that $`\xi _{}\delta \mu ^{\zeta _\pi }`$. It would be very interesting to explore the possible connection between conic filling and SFL critical wetting in two and three dimensions in more detail. Whilst for $`d=2`$ the equivalence of the critical exponents and PDF’s has been establised exactly using transfer-matrix methods it would be extremely instructive to understand this from a RG perspective. This would perhaps throw light on a more detailed theory of three dimensional cone filling based on the full interfacial model. Also our analysis of wedge and conic filling in $`d=3`$ was limited to pure thermal disorder unlike our analysis in $`d=2`$. Studies of three dimensional wedge and cone filling in impure systems would be very illuminating. It may be that for impure systems, filling and SFL critical wetting only share common features in $`d=2`$.

AJW and CR acknowledge support from the EPSRC and the E.C. under contract ERBFMBICT983229 respectively.

## References