Signatures of Topological Superconductors

<p>Topological superconductors are particularly interesting in light of the active ongoing experimental efforts for realizing exotic physics such as Majorana zero modes. These systems have excitations with non-Abelian exchange statistics, which provides a path towards topological quantum information processing. Intrinsic topological superconductors are quite rare in nature. However, one can engineer topological superconductivity by inducing effective p-wave pairing in materials which can be grown in the laboratory. One possibility is to induce the proximity effect in topological insulators; another is to use hybrid structures of superconductors and semiconductors.</p> <p>The proposal of interfacing s-wave superconductors with quantum spin Hall systems provides a promising route to engineered topological superconductivity. Given the exciting recent progress on the fabrication side, identifying experiments that definitively expose the topological superconducting phase (and clearly distinguish it from a trivial state) raises an increasingly important problem. With this goal in mind, we proposed a detection scheme to get an unambiguous signature of topological superconductivity, even in the presence of ordinarily detrimental effects such as thermal fluctuations and quasiparticle poisoning. We considered a Josephson junction built on top of a quantum spin Hall material. This system allows the proximity effect to turn edge states in effective topological superconductors. Such a setup is promising because experimentalists have demonstrated that supercurrents indeed flow through quantum spin Hall edges. To demonstrate the topological nature of the superconducting quantum spin Hall edges, theorists have proposed examining the periodicity of Josephson currents respect to the phase across a Josephson junction. The periodicity of tunneling currents of ground states in a topological superconductor Josephson junction is double that of a conventional Josephson junction. In practice, this modification of periodicity is extremely difficult to observe because noise sources, such as quasiparticle poisoning, wash out the signature of topological superconductors. For this reason, We propose a new, relatively simple DC measurement that can compellingly reveal topological superconductivity in such quantum spin Hall/superconductor heterostructures. More specifically, We develop a general framework for capturing the junction's current-voltage characteristics as a function of applied magnetic flux. Our analysis reveals sharp signatures of topological superconductivity in the field-dependent critical current. These signatures include the presence of multiple critical currents and a non-vanishing critical current for all magnetic field strengths as a reliable identification scheme for topological superconductivity.</p> <p>This system becomes more interesting as interactions between electrons are involved. By modeling edge states as a Luttinger liquid, we find conductance provides universal signatures to distinguish between normal and topological superconductors. More specifically, we use renormalization group methods to extract universal transport characteristics of superconductor/quantum spin Hall heterostructures where the native edge states serve as a lead. Interestingly, arbitrarily weak interactions induce qualitative changes in the behavior relative to the free-fermion limit, leading to a sharp dichotomy in conductance for the trivial (narrow superconductor) and topological (wide superconductor) cases. Furthermore, we find that strong interactions can in principle induce parafermion excitations at a superconductor/quantum spin Hall junction.</p> <p>As we identify the existence of topological superconductor, we can take a step further. One can use topological superconductor for realizing Majorana modes by breaking time reversal symmetry. An advantage of 2D topological insulator is that networks required for braiding Majoranas along the edge channels can be obtained by adjoining 2D topological insulator to form corner junctions. Physically cutting quantum wells for this purpose, however, presents technical challenges. For this reason, I propose a more accessible means of forming networks that rely on dynamically manipulating the location of edge states inside of a single 2D topological insulator sheet. In particular, I show that edge states can effectively be dragged into the system's interior by gating a region near the edge into a metallic regime and then removing the resulting gapless carriers via proximity-induced superconductivity. This method allows one to construct rather general quasi-1D networks along which Majorana modes can be exchanged by electrostatic means.</p> <p>Apart from 2D topological insulators, Majorana fermions can also be generated in other more accessible materials such as semiconductors. Following up on a suggestion by experimentalist Charlie Marcus, I proposed a novel geometry to create Majorana fermions by placing a 2D electron gas in proximity to an interdigitated superconductor-ferromagnet structure. This architecture evades several manufacturing challenges by allowing single-side fabrication and widening the class of 2D electron gas that may be used, such as the surface states of bulk semiconductors. Furthermore, it naturally allows one to trap and manipulate Majorana fermions through the application of currents. Thus, this structure may lead to the development of a circuit that enables fully electrical manipulation of topologically-protected quantum memory. To reveal these exotic Majorana zero modes, I also proposed an interference scheme to detect Majorana fermions that is broadly applicable to any 2D topological superconductor platform.</p>

Wave uplift pressures on horizontal platforms

<p>The major objective of the study has been to investigate in detail the rapidly-varying peak uplift pressure and the slowly-varying positive and negative uplift pressures that are known to be exerted by waves against the underside of a horizontal pier or platform located above the still water level, but not higher than the crests of the incident waves.</p> <p>In a "two-dimensional" laboratory study conducted in a 100-ft long by 15-in.-wide by 2-ft-deep wave tank with a horizontal smooth bottom, individually generated solitary waves struck a rigid, fixed, horizontal platform extending the width of the tank. Pressure transducers were mounted flush with the smooth soffit, or underside, of the platform. The location of the transducers could be varied.</p> <p>The problem of a d equate dynamic and spatial response of the transducers was investigated in detail. It was found that unless the radius of the sensitive area of a pressure transducer is smaller than about one-third of the characteristic width of the pressure distribution, the peak pressure and the rise-time will not be recorded accurately. A procedure was devised to correct peak pressures and rise-times for this transducer defect.</p> <p>The hydrodynamics of the flow beneath the platform are described qualitatively by a si1nple analysis, which relates peak pressure and positive slowly-varying pressure to the celerity of the wave front propagating beneath the platform, and relates negative slowly-varying pressure to the process by which fluid recedes from the platform after the wave has passed. As the wave front propagates beneath the platform, its celerity increases to a maximum, then decreases. The peak pressure similarly increases with distance from the seaward edge of the platform, then decreases.</p> <p>Measured peak pressure head, always found to be less than five times the incident wave height above still water level, is an order of magnitude less than reported shock pressures due to waves breaking against vertical walls; the product of peak pressure and rise-time, considered as peak impulse, is of the order of 20% of reported shock impulse due to waves breaking against vertical walls. The maximum measured slowly-varying uplift pressure head is approximately equal to the incident wave height less the soffit clearance above still water level. The normalized magnitude and duration of negative pressure appears to depend principally on the ratio of soffit clearance to still water depth and on the ratio of platform length to still water depth.</p>

Electromagnetic wave propagation and scattering in a randomly-inhomogeneous dielectric sphere

<p>Electromagnetic wave propagation and scattering in a sphere composed of an inhomogeneous medium having random variations in its permittivity are studied by utilizing the Born approximation in solving the vector wave equation. The variations in the permittivity are taken to be isotropic and homogeneous, and are spatially characterized by a Gaussian correlation function. Temporal variations in the medium are not considered.</p> <p>Two particular problems are considered: i) finding the far-zone electric field when an electric or magnetic dipole is situated at the center of the sphere, and ii) finding the electric field at the sphere's center when a linearly polarized plane wave is incident upon it. Expressions are obtained for the mean-square magnitudes of the scattered field components; it is found that the mean of the product of any two transverse components vanishes. The cases where the wavelength is much shorter than correlation distance of the medium and where it is much longer than it are both considered.</p>

An Optically recombined laser interferometer for gravitational wave detection

<p>Kilometer scale interferometers for the detection of gravitational waves are currently under construction by the LIGO (Laser Interferometer Gravitational-wave Observatory) and VIRGO projects. These interferometers will consist of two Fabry-Perot cavities illuminated by a laser beam which is split in half by a beam splitter. A recycling mirror between the laser and the beam splitter will reflect the light returning from the beam splitter towards the laser back into the interferometer. The positions of the optical components in these interferometers must be controlled to a small fraction of a wavelength of the laser light. Schemes to extract signals necessary to control these optical components have been developed and demonstrated on the tabletop. In the large scale gravitational wave detectors the optical components must be suspended from vibration isolation platforms to achieve the necessary isolation from seismic motion. These suspended components present a new class of problems in controlling the interferometer, but also provide more exacting test of interferometer signal and noise models.</p> <p>This thesis discusses the first operation of a suspended-mass Fabry-Perot-Michelson interferometer, in which signals carried by the optically recombined beams are used to detect and control all important mirror displacements. This interferometer uses an optical configuration and signal extraction scheme that is planned for the full scale LIGO interferometers with the simplification of the removal of the recycling mirror. A theoretical analysis of the performance that is expected from such an interferometer is presented and the experimental results are shown to be in generally good agreement.</p>

Wave induced oscillations in harbors of arbitrary shape

<p>Theoretical and experimental studies were conducted to investigate the wave induced oscillations in an arbitrary shaped harbor with constant depth which is connected to the open-sea.</p> <p>A theory termed the ���arbitrary shaped harbor��� theory is developed. The solution of the Helmholtz equation, ���²f + k²f = 0, is formulated as an integral equation; an approximate method is employed to solve the integral equation by converting it to a matrix equation. The final solution is obtained by equating, at the harbor entrance, the wave amplitude and its normal derivative obtained from the solutions for the regions outside and inside the harbor.</p> <p>Two special theories called the circular harbor theory and the rectangular harbor theory are also developed. The coordinates inside a circular and a rectangular harbor are separable; therefore, the solution for the region inside these harbors is obtained by the method of separation of variables. For the solution in the open-sea region, the same method is used as that employed for the arbitrary shaped harbor theory. The final solution is also obtained by a matching procedure similar to that used for the arbitrary shaped harbor theory. These two special theories provide a useful analytical check on the arbitrary shaped harbor theory.</p> <p>Experiments were conducted to verify the theories in a wave basin 15 ft wide by 31 ft long with an effective system of wave energy dissipators mounted along the boundary to simulate the open-sea condition.</p> <p>Four harbors were investigated theoretically and experimentally: circular harbors with a 10�� opening and a 60�� opening, a rectangular harbor, and a model of the East and West Basins of Long Beach Harbor located in Long Beach, California.</p> <p>Theoretical solutions for these four harbors using the arbitrary shaped harbor theory were obtained. In addition, the theoretical solutions for the circular harbors and the rectangular harbor using the two special theories were also obtained. In each case, the theories have proven to agree well with the experimental data.</p> <p>It is found that: (1) the resonant frequencies for a specific harbor are predicted correctly by the theory, although the amplification factors at resonance are somewhat larger than those found experimentally,(2) for the circular harbors, as the width of the harbor entrance increases, the amplification at resonance decreases, but the wave number bandwidth at resonance increases, (3) each peak in the curve of entrance velocity vs incident wave period corresponds to a distinct mode of resonant oscillation inside the harbor, thus the velocity at the harbor entrance appears to be a good indicator for resonance in harbors of complicated shape, (4) the results show that the present theory can be applied with confidence to prototype harbors with relatively uniform depth and reflective interior boundaries.</p>

The resolution of the thermodynamic paradox and the theory of guided wave propagation in anisotropic media

<p>The resolution of the so-called thermodynamic paradox is presented in this paper. It is shown, in direct contradiction to the results of several previously published papers, that the cutoff modes (evanescent modes having complex propagation constants) can carry power in a waveguide containing ferrite. The errors in all previous ���proofs��� which purport to show that the cutoff modes cannot carry power are uncovered. The boundary value problem underlying the paradox is studied in detail; it is shown that, although the solution is somewhat complicated, there is nothing paradoxical about it.</p> <p>The general problem of electromagnetic wave propagation through rectangular guides filled inhomogeneously in cross-section with transversely magnetized ferrite is also studied. Application of the standard waveguide techniques reduces the TM part to the well-known self-adjoint Sturm Liouville eigenvalue equation. The TE part, however, leads in general to a non-self-adjoint eigenvalue equation. This equation and the associated expansion problem are studied in detail. Expansion coefficients and actual fields are determined for a particular problem. </p>

Steady surface wave pattern in a shear flow

<p>The subject under investigation concerns the steady surface wave patterns created by small concentrated disturbances acting on a non-uniform flow of a heavy fluid. The initial value problem of a point disturbance in a primary flow having an arbitrary velocity distribution (U(y), 0, 0) in a direction parallel to the undisturbed free surface is formulated. A geometric optics method and the classical integral transformation method are employed as two different methods of solution for this problem. Whenever necessary, the special case of linear shear (i.e. U(y) = 1+��y)) is chosen for the purpose of facilitating the final integration of the solution.</p> <p>The asymptotic form of the solution obtained by the method of integral transforms agrees with the leading terms of the solution obtained by geometric optics when the latter is expanded in powers of small �� r.</p> <p>The overall effect of the shear is to confine the wave field on the downstream side of the disturbance to a region which is smaller than the wave region in the case of uniform flows. If U(y) vanishes, and changes sign at a critical plane y = y_cr (e.g. ��y_cr = -1 for the case of linear shear), then the boundary of this asymmetric wave field approaches this critical vertical plane. On this boundary the wave crests are all perpendicular to the x-axis, indicating that waves are reflected at this boundary. </p> <p>Inside the wave field, as in the case of a point disturbance in a uniform primary flow, there exist two wave systems. The loci of constant phases (such as the crests or troughs) of these wave systems are not symmetric with respect to the x-axis. The geometric optics method and the integral transform method yield the same result of these loci for the special case of U(y) = U_o(1 + ��y) and for large Kr (��r ���� 1 ���� Kr).</p> <p>An expression for the variation of the amplitude of the waves in the wave field is obtained by the integral transform method. This is in the form of an expansion in small ��r. The zeroth order is identical to the expression for the uniform stream case and is thus not applicable near the boundary of the wave region because it becomes infinite in that neighborhood. Throughout this investigation the viscous terms in the equations of motion are neglected, a reasonable assumption which can be justified when the wavelengths of the resulting waves are sufficiently large. </p>

High frequency wave propagation in the earth: theory and observation

<p>The wave-theoretical analysis of acoustic and elastic waves refracted by a spherical boundary across which both velocity and density increase abruptly and thence either increase or decrease continuously with depth is formulated in terms of the general problem of waves generated at a steady point source and scattered by a radially heterogeneous spherical body. A displacement potential representation is used for the elastic problem that results in high frequency decoupling of P-SV motion in a spherically symmetric, radially heterogeneous medium. Through the application of an earth-flattening transformation on the radial solution and the Watson transform on the sum over eigenfunctions, the solution to the spherical problem for high frequencies is expressed as a Weyl integral for the corresponding half-space problem in which the effect of boundary curvature maps into an effective positive velocity gradient. The results of both analytical and numerical evaluation of this integral can be summarized as follows for body waves in the crust and upper mantle:</p> <p>1) In the special case of a critical velocity gradient (a gradient equal and opposite to the effective curvature gradient), the critically refracted wave reduces to the classical head wave for flat, homogeneous layers.</p> <p>2) For gradients more negative than critical, the amplitude of the critically refracted wave decays more rapidly with distance than the classical head wave.</p> <p>3) For positive, null, and gradients less negative than critical, the amplitude of the critically refracted wave decays less rapidly with distance than the classical head wave, and at sufficiently large distances, the refracted wave can be adequately described in terms of ray-theoretical diving waves. At intermediate distances from the critical point, the spectral amplitude of the refracted wave is scalloped due to multiple diving wave interference.</p> <p>These theoretical results applied to published amplitude data for P-waves refracted by the major crustal and upper mantle horizons (the Pg, P*, and Pn travel-time branches) suggest that the 'granitic' upper crust, the 'basaltic' lower crust, and the mantle lid all have negative or near-critical velocity gradients in the tectonically active western United States. On the other hand, the corresponding horizons in the stable eastern United States appear to have null or slightly positive velocity gradients. The distribution of negative and positive velocity gradients correlates closely with high heat flow in tectonic regions and normal heat flow in stable regions. The velocity gradients inferred from the amplitude data are generally consistent with those inferred from ultrasonic measurements of the effects of temperature and pressure on crustal and mantle rocks and probable geothermal gradients. A notable exception is the strong positive velocity gradient in the mantle lid beneath the eastern United States (2 x 10^-3 sec^-1), which appears to require a compositional gradient to counter the effect of even a small geothermal gradient.</p> <p>New seismic-refraction data were recorded along a 800 km profile extending due south from the Canadian border across the Columbia Plateau into eastern Oregon. The source for the seismic waves was a series of 20 high-energy chemical explosions detonated by the Canadian government in Greenbush Lake, British Columbia. The first arrivals recorded along this profile are on the Pn travel-time branch. In northern Washington and central Oregon their travel time is described by T = ��/8.0 + 7.7 sec, but in the Columbia Plateau the Pn arrivals are as much as 0.9 sec early with respect to this line. An interpretation of these Pn arrivals together with later crustal arrivals suggest that the crust under the Columbia Plateau is thinner by about 10 km and has a higher average P-wave velocity than the 35-km-thick, 62-km/sec crust under the granitic-metamorphic terrain of northern Washington. A tentative interpretation of later arrivals recorded beyond 500 km from the shots suggests that a thin 8.4-km/sec horizon may be present in the upper mantle beneath the Columbia Plateau and that this horizon may form the lid to a pronounced low-velocity zone extending to a depth of about 140 km.</p>

Spin wave resonance in ferromagnetic films

<p>The objective of this investigation has been a theoretical and experimental understanding of ferromagnetic resonance phenomena in ferromagnetic thin films, and a consequent understanding of several important physical properties of these films. Significant results have been obtained by ferromagnetic resonance, hysteresis, torque magnetometer, He ion backscattering, and X-ray fluorescence measurements for nickel-iron alloy films.</p> <p>Taking into account all relevant magnetic fields, including the applied, demagnetizing, effective anisotropy and exchange fields, the spin wave resonance condition applicable to the thin film geometry is presented. On the basis of the simple exchange interaction model it is concluded that the normal resonance modes of an ideal film are expected to be unpinned. The possibility of nonideality near the surface of a real film was considered by means of surface anisotropy field, inhomogeneity in demagnetizing field and inhomogeneity of magnetization models. Numerical results obtained for reasonable parameters in all cases show that they negligibly perturb the resonance fields and the higher order mode shapes from those of the unpinned modes of ideal films for thicknesses greater than 1000 ��. On the other hand for films thinner than 1000 �� the resonance field deviations can be significant even though the modes are very nearly unpinned. A previously unnoticed but important feature of all three models is that the interpretation of the first resonance mode as the uniform mode of an ideal film allows an accurate measurement of the average effective demagnetizing field over the film volume. Furthermore, it is demonstrated that it is possible to choose parameters which give indistinguishable predictions for all three models, making it difficult to uniquely ascertain the source of spin pinning in real films from resonance measurements alone.</p> <p>Spin wave resonance measurements of 81% Ni-19% Fe coevaporated films 30 to 9000 �� thick, at frequencies from 1 to 8 GHz, at room temperature, and with the static magnetic field parallel and perpendicular to the film plane have been performed. A self-consistent analysis of the results for films thicker than 1000 ��, in which multiple excitations can be observed, shows for the first time that a unique value of exchange constant A can only be obtained by the use of unpinned mode assignments. This evidence and the resonance behavior of films thinner than 1000 �� strongly imply that the magnetization at the surfaces of permalloy films is very weakly pinned. However, resonance measurements alone cannot determine whether this pinning is due to a surface anisotropy, an inhomogeneous demagnetizing field or an inhomogeneous magnetization. The above analysis yields a value of 4��M=10,100 Oe and A = (1.03 �� .05) x 10^-6 erg/cm for this alloy. The ability to obtain a unique value of A suggests that spin wave resonance can be used to accurately characterize the exchange interaction in a ferromagnet.</p> <p>In an effort to resolve the ambiguity of the source of pinning of the magnetization, a correlation of the ratio of magnetic moment and X-ray film thickness with the value of effective demagnetizing field 4��NM as determined from resonance, for films 45 to 300 �� has been performed. The remarkable agreement of both quantities and a comparison with the predictions of five distinct models, strongly imply that the thickness dependence of both quantities is related to a thickness dependent average saturation magnetization, which is far below 10,100 Oe for very thin films. However, a series of complementary experiments shows that this large decrease of average saturation magnetization cannot be simply explained by either oxidation or interdiffusion processes. It can only be satisfactorily explained by an intrinsic decrease of the average saturation magnetization for very thin films, an effect which cannot be justified by any simple physical considerations.</p> <p>Recognizing that this decrease of average saturation magnetization could be due to an oxidation process, a correlation of resonance measurements, He ion backscattering, X-ray fluorescence and torque magnetometer measurements, for films 40 to 3500 �� thick has been performed. On basis of these measurements it is unambiguously established that the oxide layer on the surface of purposefully oxidized 81% Ni-19% Fe evaporated films is predominantly Fe-oxide, and that in the oxidation process Fe atoms are removed from the bulk of the film to depths of thousands of angstroms. Extrapolation of results for pure Fe films indicates that the oxide is most likely ��-Fe₂O₃. These conclusions are in agreement with results from old metallurgical studies of high temperature oxidation of bulk Fe and Ni-Fe alloys. However, X-ray fluorescence results for films oxidized at room temperature, show that although the preferential oxidation of Fe also takes place in these films, the extent of this process is by far too small to explain the large variation of their average saturation magnetization with film thickness.</p>

An efficient frequency-domain algorithm for wave scattering problems with application to jet noise

This paper presents a pseudo-time-step method to calculate a (vector) Green function for the adjoint linearised Euler equations as a scattering problem in the frequency domain, for use as a jet-noise propagation prediction tool. A method of selecting the acoustics-related solution in a truncated spatial domain while suppressing any possible shear-layer-type instability is presented. Numerical tests for 3-D axisymmetrical parallel mean flows against semi-analytical reference solutions indicate that the new iterative algorithm is capable of producing accurate solutions with modest computational requirements.

Tsunami wave runup on the clustered islands using boussinesq-type model

This paper describes the implementation of the Boussinesq-type model and extends its application to the tsunami wave runup on the clustered islands (multiple adjacent conical islands), in turn, an extensively validated two-dimensional Boussinesq-type model is employed to examine the interaction between a propagating solitary wave and multiple idealised conical islands, with particular emphasis on a combination effect of two adjustable parameters for spacing interval/diameter ratio between the adjacent conical islands, S/D, and the rotating angle of the structural configuration,�� on maximum soliton runup heights. An extensive parameter study concerning the combination effect of altering�� and S/D on the maximum soliton runup with the multi-conical islands is subsequently carried out and the distributions of the maximum runup heights on each conical island are obtained and compared for the twin-island cases. The worst case study is performed for each case in respect of the enhancement in the maximum wave runup heights by the multi-conical islands. It is found that the nonlinear wave diffraction, reflection and refraction play a significant role in varying the maximum soliton runup heights on multiconical islands. The comparatively large maximum soliton runups are generally predicted for the merged and bottom mounted clusteredislands. Furthermore, the joints of the clustered-merged islands are demonstrated to suffer the most of the tsunami wave attack. The conical islands that position in the shadow regions behind the surrounding islands are found to withstand relatively less extreme wave impact. Although, these numerical investigations are considerable simplifications of the multi conical islands, they give a critical insight into certain important hydrodynamic characteristics of the interaction between an extreme wave event and a group of clustered conical islands, and thus providing a useful engineering guidance for extreme wave mitigation and coastal development. Copyright �� 2012 by the International Society of Offshore and Polar Engineers (ISOPE).

Measurement of the magnetic field inside the holes of a drilled bulk high-Tc superconductor

We use macroscopic holes drilled in a bulk YBCO superconductor to probe its magnetic properties in the volume of the sample. The sample is subjected to an AC magnetic flux with a density ranging from 30mT to 130mT and the flux in the superconductor is probed by miniature coils inserted in the holes. In a given hole, three different penetration regimes can be observed: (i) the shielded regime, where no magnetic flux threads the hole; (ii) the gradual penetration regime, where the waveform of the magnetic field has a clipped sine shape whose fundamental component scales with the applied field; and (iii) the flux concentration regime, where the waveform of the magnetic field is nearly a sine wave, with an amplitude exceeding that of the applied field by up to a factor of two. The distribution of the penetration regimes in the holes is compared with that of the magnetic flux density at the top and bottom surfaces of the sample, and is interpreted with the help of optical polarized light micrographs of these surfaces. We show that the measurement of the magnetic field inside the holes can be used as a local characterization of the bulk magnetic properties of the sample.

Investigation of DyBa2Cu3O7-d superconducting domains grown by the infiltration technique starting with small size Dy-211 particles

An infiltration and growth process is here used as an alternative to the classical top-seeded melt-textured growth process for the production of Dy-123 single-domains with finely dispersed small size Dy-211 particles. The starting materials are the 211-particles and a barium and copper rich liquid phase precursor. The infiltration and growth process allows for controlling both the spatial and size distribution of the 211-particles in the final superconducting 123-single-domain. The main parameters (set-ups, maximum processing temperature with respect to the peritectic temperature, nature of reactant, porosity of the 211-preform) of the infiltration and growth process are discussed. Moreover, different processes of chimie douce are shown in order to produce Dy-211 particles with controlled shape and size, particles that can be used as precursors for the infiltration and growth process. �� 2005 IOP Publishing Ltd.