HOME

TheInfoList



OR:

Rutherford backscattering spectrometry (RBS) is an
analytical technique Analytical technique is a method used to determine a chemical or physical property of a chemical substance, chemical element, or mixture. There is a wide variety of techniques used for analysis, from simple weighing to advanced techniques using high ...
used in materials science. Sometimes referred to as high-energy ion scattering (HEIS) spectrometry, RBS is used to determine the structure and composition of materials by measuring the backscattering of a beam of high energy ions (typically
proton A proton is a stable subatomic particle, symbol , H+, or 1H+ with a positive electric charge of +1 ''e'' elementary charge. Its mass is slightly less than that of a neutron and 1,836 times the mass of an electron (the proton–electron mass ...
s or
alpha particles Alpha particles, also called alpha rays or alpha radiation, consist of two protons and two neutrons bound together into a particle identical to a helium-4 nucleus. They are generally produced in the process of alpha decay, but may also be prod ...
) impinging on a sample.


Geiger–Marsden experiment

Rutherford backscattering spectrometry is named after Lord Rutherford, a
physicist A physicist is a scientist who specializes in the field of physics, which encompasses the interactions of matter and energy at all length and time scales in the physical universe. Physicists generally are interested in the root or ultimate cau ...
sometimes referred to as the father of
nuclear physics Nuclear physics is the field of physics that studies atomic nuclei and their constituents and interactions, in addition to the study of other forms of nuclear matter. Nuclear physics should not be confused with atomic physics, which studies t ...
. Rutherford supervised a series of experiments carried out by
Hans Geiger Johannes Wilhelm "Hans" Geiger (; ; 30 September 1882 – 24 September 1945) was a German physicist. He is best known as the co-inventor of the detector component of the Geiger counter and for the Geiger–Marsden experiment which discover ...
and
Ernest Marsden Sir Ernest Marsden (19 February 1889 – 15 December 1970) was an English-New Zealand physicist. He is recognised internationally for his contributions to science while working under Ernest Rutherford, which led to the discovery of new theories ...
between 1909 and 1914 studying the scattering of
alpha particles Alpha particles, also called alpha rays or alpha radiation, consist of two protons and two neutrons bound together into a particle identical to a helium-4 nucleus. They are generally produced in the process of alpha decay, but may also be prod ...
through metal foils. While attempting to eliminate "stray particles" they believed to be caused by an imperfection in their alpha source, Rutherford suggested that Marsden attempt to measure backscattering from a gold foil sample. According to the then-dominant
plum-pudding model The plum pudding model is one of several historical scientific models of the atom. First proposed by J. J. Thomson in 1904 soon after the discovery of the electron, but before the discovery of the atomic nucleus, the model tried to explain two p ...
of the atom, in which small negative electrons were spread through a diffuse positive region, backscattering of the high-energy positive alpha particles should have been nonexistent. At most small deflections should occur as the alpha particles passed almost unhindered through the foil. Instead, when Marsden positioned the detector on the same side of the foil as the alpha particle source, he immediately detected a noticeable backscattered signal. According to Rutherford, "It was quite the most incredible event that has ever happened to me in my life. It was almost as incredible as if you fired a 15-inch shell at a piece of tissue paper and it came back and hit you." Rutherford interpreted the result of the Geiger–Marsden experiment as an indication of a
Coulomb collision A Coulomb collision is a binary elastic collision between two charged particles interacting through their own electric field. As with any inverse-square law, the resulting trajectories of the colliding particles is a hyperbolic Keplerian orbit. This ...
with a single massive positive particle. This led him to the conclusion that the atom's positive charge could not be diffuse but instead must be concentrated in a single massive core: the
atomic nucleus The atomic nucleus is the small, dense region consisting of protons and neutrons at the center of an atom, discovered in 1911 by Ernest Rutherford based on the 1909 Geiger–Marsden gold foil experiment. After the discovery of the neutron ...
. Calculations indicated that the charge necessary to accomplish this deflection was approximately 100 times the charge of the electron, close to the atomic number of gold. This led to the development of the Rutherford model of the atom in which a positive nucleus made up of ''N''e positive particles, or
protons A proton is a stable subatomic particle, symbol , H+, or 1H+ with a positive electric charge of +1 ''e'' elementary charge. Its mass is slightly less than that of a neutron and 1,836 times the mass of an electron (the proton–electron m ...
, was surrounded by ''N'' orbiting electrons of charge -e to balance the nuclear charge. This model was eventually superseded by the
Bohr atom In atomic physics, the Bohr model or Rutherford–Bohr model, presented by Niels Bohr and Ernest Rutherford in 1913, is a system consisting of a small, dense nucleus surrounded by orbiting electrons—similar to the structure of the Solar Syste ...
, incorporating some early results from
quantum mechanics Quantum mechanics is a fundamental theory in physics that provides a description of the physical properties of nature at the scale of atoms and subatomic particles. It is the foundation of all quantum physics including quantum chemistry, ...
. If the energy of the incident particle is increased sufficiently, the
Coulomb barrier The Coulomb barrier, named after Coulomb's law, which is in turn named after physicist Charles-Augustin de Coulomb, is the energy barrier due to electrostatic interaction that two nuclei need to overcome so they can get close enough to undergo a ...
is exceeded and the
wavefunctions A wave function in quantum physics is a mathematical description of the quantum state of an isolated quantum system. The wave function is a complex-valued probability amplitude, and the probabilities for the possible results of measurements ...
of the incident and struck particles overlap. This may result in nuclear reactions in certain cases, but frequently the interaction remains
elastic Elastic is a word often used to describe or identify certain types of elastomer, elastic used in garments or stretchable fabrics. Elastic may also refer to: Alternative name * Rubber band, ring-shaped band of rubber used to hold objects togethe ...
, although the scattering cross-sections may fluctuate wildly as a function of energy and no longer be calculable analytically. This case is known as "Elastic (non-Rutherford) Backscattering Spectrometry" (EBS). There has recently been great progress in determining EBS scattering cross-sections, by solving Schrödinger's equation for each interaction. However, for the EBS analysis of matrices containing light elements, the utilization of experimentally measured scattering cross-section data is also considered to be a very credible option.


Basic principles

We describe Rutherford backscattering as an
elastic Elastic is a word often used to describe or identify certain types of elastomer, elastic used in garments or stretchable fabrics. Elastic may also refer to: Alternative name * Rubber band, ring-shaped band of rubber used to hold objects togethe ...
, hard-sphere collision between a high kinetic energy particle from the incident beam (the ''projectile'') and a stationary particle located in the sample (the ''target''). ''Elastic'' in this context means that no energy is transferred between the incident particle and the stationary particle during the collision, and the state of the stationary particle is not changed. (Except that for a small amount of momentum, which is ignored.) Nuclear interactions are generally not elastic, since a collision may result in a nuclear reaction, with the release of considerable quantities of energy.
Nuclear reaction analysis Nuclear may refer to: Physics Relating to the nucleus of the atom: *Nuclear engineering * Nuclear physics * Nuclear power * Nuclear reactor * Nuclear weapon *Nuclear medicine * Radiation therapy * Nuclear warfare Mathematics * Nuclear space * ...
(NRA) is useful for detecting light elements. However, this is not Rutherford scattering. Considering the
kinematics Kinematics is a subfield of physics, developed in classical mechanics, that describes the motion of points, bodies (objects), and systems of bodies (groups of objects) without considering the forces that cause them to move. Kinematics, as a fiel ...
of the collision (that is, the conservation of momentum and kinetic energy), the energy E1 of the scattered projectile is reduced from the initial energy E0: :E_1 = k \cdot E_0, where k is known as the ''kinematical factor'', and :k = \left(\frac\right)^2,Oura ''et al''. (2003) p. 110 where particle 1 is the projectile, particle 2 is the target nucleus, and \theta_1 is the scattering angle of the projectile in the laboratory
frame of reference In physics and astronomy, a frame of reference (or reference frame) is an abstract coordinate system whose origin, orientation, and scale are specified by a set of reference points― geometric points whose position is identified both math ...
(that is, relative to the observer). The plus sign is taken when the mass of the projectile is less than that of the target, otherwise the minus sign is taken. While this equation correctly determines the energy of the scattered projectile for any particular scattering angle (relative to the observer), it does not describe the probability of observing such an event. For that we need the ''differential cross-section'' of the backscattering event: :\frac = \left(\frac\right)^2 \frac, where Z_1 and Z_2 are the atomic numbers of the incident and target nuclei. This equation is written in the
centre of mass In physics, the center of mass of a distribution of mass in space (sometimes referred to as the balance point) is the unique point where the weighted relative position of the distributed mass sums to zero. This is the point to which a force may ...
frame of reference In physics and astronomy, a frame of reference (or reference frame) is an abstract coordinate system whose origin, orientation, and scale are specified by a set of reference points― geometric points whose position is identified both math ...
and is therefore not a function of the mass of either the projectile or the target nucleus. The scattering angle in the laboratory frame of reference \theta_1 is ''not'' the same as the scattering angle in the centre of mass frame of reference \theta (although for RBS experiments they are usually very similar). However, heavy ion projectiles can easily
recoil Recoil (often called knockback, kickback or simply kick) is the rearward thrust generated when a gun is being discharged. In technical terms, the recoil is a result of conservation of momentum, as according to Newton's third law the force r ...
lighter ions which, if the geometry is right, can be ejected from the target and detected. This is the basis of the
Elastic Recoil Detection Elastic recoil detection analysis (ERDA), also referred to as forward recoil scattering (or, contextually, spectrometry), is an ion beam analysis technique in materials science to obtain elemental concentration depth profiles in thin films. This ...
(ERD, with synonyms ERDA, FRS, HFS) technique. RBS often uses a He beam which readily recoils H, so simultaneous RBS/ERD is frequently done to probe the hydrogen isotope content of samples (although H ERD with a He beam above 1 MeV is not Rutherford: see http://www-nds.iaea.org/sigmacalc). For ERD the scattering angle in the lab frame of reference is quite different from that in the centre of mass frame of reference. Heavy ions cannot ''back''scatter from light ones: it is kinematically prohibited. The kinematical factor must remain real, and this limits the permitted scattering angle in the laboratory frame of reference. In ERD it is often convenient to place the recoil detector at recoil angles large enough to prohibit signal from the scattered beam. The scattered ion intensity is always very large compared to the recoil intensity (the Rutherford scattering cross-section formula goes to infinity as the scattering angle goes to zero), and for ERD the scattered beam usually has to be excluded from the measurement somehow. The singularity in the Rutherford scattering cross-section formula is unphysical of course. If the scattering cross-section is zero it implies that the projectile never comes close to the target, but in this case it also never penetrates the electron cloud surrounding the nucleus either. The pure Coulomb formula for the scattering cross-section shown above must be corrected for this
screening effect In physics, screening is the damping of electric fields caused by the presence of mobile charge carriers. It is an important part of the behavior of charge-carrying fluids, such as ionized gases (classical plasmas), electrolytes, and charge carr ...
, which becomes more important as the energy of the projectile decreases (or, equivalently, its mass increases). While large-angle scattering only occurs for ions which scatter off target nuclei, inelastic small-angle scattering can also occur off the sample electrons. This results in a gradual decrease in the kinetic energy of incident ions as they penetrate into the sample, so that backscattering off interior nuclei occurs with a lower "effective" incident energy. Similarly backscattered ions lose energy to electrons as they exit the sample. The amount by which the ion energy is lowered after passing through a given distance is referred to as the
stopping power Stopping power is the ability of a weapon – typically a ranged weapon such as a firearm – to cause a target (human or animal) to be incapacitated or immobilized. Stopping power contrasts with lethality in that it pertains only to a weapo ...
of the material and is dependent on the electron distribution. This energy loss varies continuously with respect to distance traversed, so that stopping power is expressed as : S(E) = -. For high energy ions stopping power is usually proportional to \frac; however, precise calculation of stopping power is difficult to carry out with any accuracy. Stopping power (properly, ''stopping force'') has units of energy per unit length. It is generally given in thin film units, that is eV /(atom/cm2) since it is measured experimentally on thin films whose thickness is always measured absolutely as mass per unit area, avoiding the problem of determining the density of the material which may vary as a function of thickness. Stopping power is now known for all materials at around 2%, see http://www.srim.org.


Instrumentation

An RBS instrument generally includes three essential components: * An ion source, usually
alpha particle Alpha particles, also called alpha rays or alpha radiation, consist of two protons and two neutrons bound together into a particle identical to a helium-4 nucleus. They are generally produced in the process of alpha decay, but may also be prod ...
s (He2+ ions) or, less commonly,
proton A proton is a stable subatomic particle, symbol , H+, or 1H+ with a positive electric charge of +1 ''e'' elementary charge. Its mass is slightly less than that of a neutron and 1,836 times the mass of an electron (the proton–electron mass ...
s. * A
linear particle accelerator A linear particle accelerator (often shortened to linac) is a type of particle accelerator that accelerates charged subatomic particles or ions to a high speed by subjecting them to a series of oscillating electric potentials along a linear b ...
capable of accelerating incident ions to high energies, usually in the range 1-3 MeV. * A detector capable of measuring the energies of
backscatter In physics, backscatter (or backscattering) is the reflection of waves, particles, or signals back to the direction from which they came. It is usually a diffuse reflection due to scattering, as opposed to specular reflection as from a mirror, a ...
ed ions over some range of angles. Two common source/acceleration arrangements are used in commercial RBS systems, working in either one or two stages. One-stage systems consist of a He+ source connected to an acceleration tube with a high positive potential applied to the ion source, and the ground at the end of the acceleration tube. This arrangement is simple and convenient, but it can be difficult to achieve energies of much more than 1 MeV due to the difficulty of applying very high voltages to the system. Two-stage systems, or "tandem accelerators", start with a source of He ions and position the positive terminal at the center of the acceleration tube. A stripper element included in the positive terminal removes electrons from ions which pass through, converting He ions to He++ ions. The ions thus start out being attracted to the terminal, pass through and become positive, and are repelled until they exit the tube at ground. This arrangement, though more complex, has the advantage of achieving higher accelerations with lower applied voltages: a typical tandem accelerator with an applied voltage of 750 kV can achieve ion energies of over 2 MeV. Detectors to measure backscattered energy are usually
silicon Silicon is a chemical element with the symbol Si and atomic number 14. It is a hard, brittle crystalline solid with a blue-grey metallic luster, and is a tetravalent metalloid and semiconductor. It is a member of group 14 in the periodic ...
surface barrier detectors, a very thin layer (100 nm) of P-type silicon on an N-type substrate forming a p-n junction. Ions which reach the detector lose some of their energy to
inelastic scattering In chemistry, nuclear physics, and particle physics, inelastic scattering is a fundamental scattering process in which the kinetic energy of an incident particle is not conserved (in contrast to elastic scattering). In an inelastic scattering proces ...
from the electrons, and some of these electrons gain enough energy to overcome the
band gap In solid-state physics, a band gap, also called an energy gap, is an energy range in a solid where no electronic states can exist. In graphs of the electronic band structure of solids, the band gap generally refers to the energy difference ( ...
between the semiconductor valence and
conduction band In solid-state physics, the valence band and conduction band are the bands closest to the Fermi level, and thus determine the electrical conductivity of the solid. In nonmetals, the valence band is the highest range of electron energies in ...
s. This means that each ion incident on the detector will produce some number of
electron-hole pairs In the solid-state physics of semiconductors, carrier generation and carrier recombination are processes by which mobile charge carriers ( electrons and electron holes) are created and eliminated. Carrier generation and recombination processes a ...
which is dependent on the energy of the ion. These pairs can be detected by applying a voltage across the detector and measuring the current, providing an effective measurement of the ion energy. The relationship between ion energy and the number of electron-hole pairs produced will be dependent on the detector materials, the type of ion and the efficiency of the current measurement; energy resolution is dependent on thermal fluctuations. After one ion is incident on the detector, there will be some ''
dead time For detection systems that record discrete events, such as particle and nuclear detectors, the dead time is the time after each event during which the system is not able to record another event. An everyday life example of this is what happens when ...
'' before the electron-hole pairs recombine in which a second incident ion cannot be distinguished from the first. Angular dependence of detection can be achieved by using a movable detector, or more practically by separating the surface barrier detector into many independent cells which can be measured independently, covering some range of angles around direct (180 degrees) back-scattering. Angular dependence of the incident beam is controlled by using a tiltable sample stage.


Composition and depth measurement

The energy loss of a backscattered ion is dependent on two processes: the energy lost in scattering events with sample nuclei, and the energy lost to small-angle scattering from the sample electrons. The first process is dependent on the scattering cross-section of the nucleus and thus on its mass and atomic number. For a given measurement angle, nuclei of two different elements will therefore scatter incident ions to different degrees and with different energies, producing separate peaks on an N(E) plot of measurement count versus energy. These peaks are characteristic of the elements contained in the material, providing a means of analyzing the composition of a sample by matching scattered energies to known scattering cross-sections. Relative concentrations can be determined by measuring the heights of the peaks. The second energy loss process, the stopping power of the sample electrons, does not result in large discrete losses such as those produced by nuclear collisions. Instead it creates a gradual energy loss dependent on the electron density and the distance traversed in the sample. This energy loss will lower the measured energy of ions which backscatter from nuclei inside the sample in a continuous manner dependent on the depth of the nuclei. The result is that instead of the sharp backscattered peaks one would expect on an N(E) plot, with the width determined by energy and angular resolution, the peaks observed trail off gradually towards lower energy as the ions pass through the depth occupied by that element. Elements which only appear at some depth inside the sample will also have their peak positions shifted by some amount which represents the distance an ion had to traverse to reach those nuclei. In practice, then, a compositional depth profile can be determined from an RBS N(E) measurement. The elements contained by a sample can be determined from the positions of peaks in the energy spectrum. Depth can be determined from the width and shifted position of these peaks, and relative concentration from the peak heights. This is especially useful for the analysis of a multilayer sample, for example, or for a sample with a composition which varies more continuously with depth. This kind of measurement can only be used to determine elemental composition; the chemical structure of the sample cannot be determined from the N(E) profile. However, it is possible to learn something about this through RBS by examining the crystal structure. This kind of spatial information can be investigated by taking advantage of blocking and channeling.


Structural measurements: blocking and channeling

To fully understand the interaction of an incident beam of nuclei with a crystalline structure, it is necessary to comprehend two more key concepts: ''blocking'' and '' channeling''. When a beam of ions with parallel trajectories is incident on a target atom, scattering off that atom will prevent collisions in a cone-shaped region "behind" the target relative to the beam. This occurs because the repulsive potential of the target atom bends close ion trajectories away from their original path, and is referred to as ''blocking''. The radius of this blocked region, at a distance L from the original atom, is given by :R = 2\sqrt When an ion is scattered from deep inside a sample, it can then re-scatter off a second atom, creating a second blocked cone in the direction of the scattered trajectory. This can be detected by carefully varying the detection angle relative to the incident angle. ''Channeling'' is observed when the incident beam is aligned with a major symmetry axis of the crystal. Incident nuclei which avoid collisions with surface atoms are excluded from collisions with all atoms deeper in the sample, due to blocking by the first layer of atoms. When the interatomic distance is large compared to the radius of the blocked cone, the incident ions can penetrate many times the interatomic distance without being backscattered. This can result in a drastic reduction of the observed backscattered signal when the incident beam is oriented along one of the symmetry directions, allowing determination of a sample's regular crystal structure. Channeling works best for very small blocking radii, i.e. for high-energy, low-atomic-number incident ions such as He+. The tolerance for the deviation of the ion beam angle of incidence relative to the symmetry direction depends on the blocking radius, making the allowable deviation angle proportional to :\sqrt While the intensity of an RBS peak is observed to decrease across most of its width when the beam is channeled, a narrow peak at the high-energy end of larger peak will often be observed, representing surface scattering from the first layer of atoms. The presence of this peak opens the possibility of surface sensitivity for RBS measurements.


Profiling of displaced atoms

In addition, channeling of ions can also be used to analyze a crystalline sample for lattice damage. If atoms within the target are displaced from their crystalline lattice site, this will result in a higher backscattering yield in relation to a perfect crystal. By comparing the spectrum from a sample being analyzed to that from a perfect crystal, and that obtained at a random (non-channeling) orientation (representative of a spectrum from an amorphous sample), it is possible to determine the extent of crystalline damage in terms of a fraction of displaced atoms. Multiplying this fraction by the density of the material when amorphous then also gives an estimate for the concentration of displaced atoms. The energy at which the increased backscattering occurs can also be used to determine the depth at which the displaced atoms are and a defect depth profile can be built up as a result.


Surface sensitivity

While RBS is generally used to measure the bulk composition and structure of a sample, it is possible to obtain some information about the structure and composition of the sample surface. When the signal is channeled to remove the bulk signal, careful manipulation of the incident and detection angles can be used to determine the relative positions of the first few layers of atoms, taking advantage of blocking effects. The surface structure of a sample can be changed from the ideal in a number of ways. The first layer of atoms can change its distance from subsequent layers ( relaxation); it can assume a different two-dimensional structure than the bulk (
reconstruction Reconstruction may refer to: Politics, history, and sociology * Reconstruction (law), the transfer of a company's (or several companies') business to a new company *''Perestroika'' (Russian for "reconstruction"), a late 20th century Soviet Unio ...
); or another material can be
adsorbed Adsorption is the adhesion of atoms, ions or molecules from a gas, liquid or dissolved solid to a surface. This process creates a film of the ''adsorbate'' on the surface of the ''adsorbent''. This process differs from absorption, in which a ...
onto the surface. Each of these cases can be detected by RBS. For example, surface reconstruction can be detected by aligning the beam in such a way that channeling should occur, so that only a surface peak of known intensity should be detected. A higher-than-usual intensity or a wider peak will indicate that the first layers of atoms are failing to block the layers beneath, i.e. that the surface has been reconstructed. Relaxations can be detected by a similar procedure with the sample tilted so the ion beam is incident at an angle selected so that first-layer atoms should block backscattering at a diagonal; that is, from atoms which are below and displaced from the blocking atom. A higher-than-expected backscattered yield will indicate that the first layer has been displaced relative to the second layer, or relaxed. Adsorbate materials will be detected by their different composition, changing the position of the surface peak relative to the expected position. RBS has also been used to measure processes which affect the surface differently from the bulk by analyzing changes in the channeled surface peak. A well-known example of this is the RBS analysis of the premelting of lead surfaces by Frenken, Maree and van der Veen. In an RBS measurement of the Pb (110) surface, a well-defined surface peak which is stable at low temperatures was found to become wider and more intense as temperature increase past two-thirds of the bulk melting temperature. The peak reached the bulk height and width as temperature reached the melting temperature. This increase in the disorder of the surface, making deeper atoms visible to the incident beam, was interpreted as pre-melting of the surface, and computer simulations of the RBS process produced similar results when compared with theoretical pre-melting predictions. RBS has also been combined with nuclear microscopy, in which a focused ion beam is scanned across a surface in a manner similar to a scanning electron microscope. The energetic analysis of backscattered signals in this kind of application provides compositional information about the surface, while the microprobe itself can be used to examine features such as periodic surface structures.Hobbs ''et al''. (1988)


See also

*
Collision cascade In condensed-matter physics, a collision cascade (also known as a displacement cascade or a displacement spike) is a set of nearby adjacent energetic (much higher than ordinary thermal energies) collisions of atoms induced by an energetic par ...
*
Elastic recoil detection Elastic recoil detection analysis (ERDA), also referred to as forward recoil scattering (or, contextually, spectrometry), is an ion beam analysis technique in materials science to obtain elemental concentration depth profiles in thin films. This ...
* Geiger–Marsden experiment * Ion beam analysis * Nuclear microscopy *
Nuclear reaction analysis Nuclear may refer to: Physics Relating to the nucleus of the atom: *Nuclear engineering * Nuclear physics * Nuclear power * Nuclear reactor * Nuclear weapon *Nuclear medicine * Radiation therapy * Nuclear warfare Mathematics * Nuclear space * ...
* Particle induced X-ray emission *
Rutherford scattering In particle physics, Rutherford scattering is the elastic scattering of charged particles by the Coulomb interaction. It is a physical phenomenon explained by Ernest Rutherford in 1911 that led to the development of the planetary Rutherford model ...
*
Secondary ion mass spectrometry Secondary-ion mass spectrometry (SIMS) is a technique used to analyze the composition of solid surfaces and thin films by sputtering the surface of the specimen with a focused primary ion beam and collecting and analyzing ejected secondary ions. ...
*
Stopping power (particle radiation) In nuclear and materials physics, stopping power is the retarding force acting on charged particles, typically alpha and beta particles, due to interaction with matter, resulting in loss of particle kinetic energy. Its application is important in ...
*
Surface science Surface science is the study of physical and chemical phenomena that occur at the interface of two phases, including solid– liquid interfaces, solid– gas interfaces, solid–vacuum interfaces, and liquid– gas interfaces. It includes th ...


References

* * * * * * * *{{cite journal , author=Frenken, J.W.M., author2=Maree, P.M.J., author3=van der Veen, J.F. , date=1986 , title=Observation of surface-initiated melting , journal= Phys. Rev. B , volume=34 , issue= 11, pages=7506–7516 , doi=10.1103/PhysRevB.34.7506 , pmid=9939429, bibcode = 1986PhRvB..34.7506F , hdl=1887/71635, s2cid=27028634 , hdl-access=free Materials science Scientific techniques Spectroscopy Ernest Rutherford