• Volume/Page
  • Keyword
  • DOI
  • Citation
  • Advanced
   
 
 
 

Flickr Twitter iResearch App Facebook

Year Range: 
Search Issue | RSS Feeds RSS
Previous Issue Next Issue

1 Jul 1962

Volume 33, Issue 7, pp. 2145-2400

Page 1 of 3 Pages Next Page | Jump to Page

Isotope Effect in an Electric Field and Jump Frequencies for Diffusion in Ionic Crystals

John R. Manning

J. Appl. Phys. 33, 2145 (1962); http://dx.doi.org/10.1063/1.1728919 (7 pages) | Cited 4 times

Online Publication Date: 11 June 2004

Full Text: | Download PDF

Show Abstract
The simultaneous diffusion of two isotopes of a given ionic impurity in an electric field is considered, and an improved method of analyzing the experimental data to obtain ΔD∕D is discussed. (Here D is the diffusion coefficient of one of the isotopes, and ΔD the difference between those of the two isotopes.) As noted by Chemla, an electric field increases the accuracy with which ΔD∕D can be determined. Equations are found giving this increase in accuracy as a function of field and diffusion time. If the ratio ww of the jump frequencies of the two isotopes is known, relative values of the vacancy or interstitialcy jump frequencies near an impurity can be calculated from ΔD∕D and μ∕D, where μ is the drift mobility. As an example, data in the literature for Na22 and Na24 diffusing in KCl are considered. For self‐diffusion, ww may be calculated from ΔD∕D and known values of the correlation factor. When data in the literature are used to estimate ww for Na22 and Na24 diffusing in NaCl, a value smaller than that predicted by classical rate theory is found.

Interaction of Sliding Metal Surfaces

M. Cocks

J. Appl. Phys. 33, 2152 (1962); http://dx.doi.org/10.1063/1.1728920 (10 pages) | Cited 52 times

Online Publication Date: 11 June 2004

Full Text: | Download PDF

Show Abstract
The friction of sliding metal surfaces is often attributed to the formation and shearing of welded junctions between the surfaces. This paper describes a more detailed investigation of the phenomena which occur in the junctions. Experiments were performed with a hemispherically ended copper rider sliding against the surface of a large cylindrical copper drum. Observation of the sliding contact through the microscope shows that most of the time during sliding, the drum and rider are held apart by a wedge of displaced metal from the drum, which becomes trapped between them. Whenever a wedge escapes from between the drum and rider, another wedge forms rapidly to replace it. The wedges are also observed with several other combinations of metals, and with pairs of flat copper surfaces. Sections cut through the wedges and adjacent copper contact members were examined. It is found that the shearing of the metal by which the sliding proceeds occurs in a direction slightly inclined to the drum surface. This causes the formation and growth of the wedges. The shearing is primarily a continuous plastic deformation process somewhat similar to the shearing process by which the chip is generated in the cutting of ductile metals with machine tools. It is possible that the adhesion between the rider and wedge may often be weaker than the junctions usually assumed to be responsible for friction. Wear debris is observed to be produced in a number of different ways.

Alternating‐Current Generation with Moving Conducting Fluids

H. A. Haus

J. Appl. Phys. 33, 2161 (1962); http://dx.doi.org/10.1063/1.1728921 (12 pages) | Cited 5 times

Online Publication Date: 11 June 2004

Full Text: | Download PDF

Show Abstract
A one‐dimensional, small‐amplitude analysis is presented of alternating‐current generation by a stationary electromagnetic circuit coupling to waves in a moving, perfectly conducting fluid. In a uniform applied dc magnetic field the waves in the fluid are either compressional waves if the time‐average velocity v0 is transverse to the dc magnetic field, or Alfvén waves if it is along the dc magnetic field. It is shown that, for velocities v0 greater than the velocity of propagation c of the compressional or Alfvén waves, self‐excited oscillations may be produced in a resistively loaded coil coupling to the fluid. In a traveling‐wave circuit and in a resistively loaded distributed circuit coupling to the fluid, growing waves are found which result from a transformation of the fluid energy into electromagnetic energy. There is a close analogy between these power‐generation mechanisms and those employed in microwave electron‐beam tubes.
The results of two‐ and three‐dimensional analyses are summarized in order to show that the phenomena predicted by the simplified one‐dimensional analysis occur in less idealized geometries. Possible applications are briefly discussed with reference to ac power generation by using combustion gases or a flow of a high‐temperature plasma.

On the Conductivity of Cadmium Sulfide Following Electron Bombardment

R. G. Schulze and B. A. Kulp

J. Appl. Phys. 33, 2173 (1962); http://dx.doi.org/10.1063/1.1728922 (3 pages) | Cited 5 times

Online Publication Date: 11 June 2004

Full Text: | Download PDF

Show Abstract
Fine rod and whisker‐type crystals of undoped CdS whose dark conductivity was of the order of 10−10 (Ω‐cm)−1 were bombarded with 100‐keV electrons. Following bombardment dark conductivity of the order of 10−3 (Ω‐cm)−1 was observed. Heat treatment to about 200°C reduced the conductivity to near its original value. Because of the high efficiency with which the electron beam raised the conductivity of the crystals, an explanation based on the redistribution of electrons and holes over existing levels rather than the production or removal of defects is required.

Dislocations and Plastic Flow in NaCl Single Crystals. I

S. Mendelson

J. Appl. Phys. 33, 2175 (1962); http://dx.doi.org/10.1063/1.1728923 (7 pages) | Cited 34 times

Online Publication Date: 11 June 2004

Full Text: | Download PDF

Show Abstract
Dislocation multiplication and motion in NaCl was studied using etch pit and birefringence techniques. The observations show that dislocations in glide bands of NaCl are in the form of large concentric loops with sources staggered in the center of the bands. They propagate and multiply by cross‐slip and are braked through intersection jogs. The slip distance parallel to b is estimated to be 2–3 mm and that normal to b 0.5–1 mm, at a stress of 30 g∕mm2 in crystals having 5×104 dislocations∕cm2. ``Deformation bands'' form normal to the operating system during the latter part of stage I by ``glide polygonization'' of edge dislocations on closely spaced slip planes. These limit dislocation motion to shorter distances in stage II. Stage III appears when significant flow occurs on a secondary system oblique to the primary operating one.

Role of Surfaces in Plastic Flow of NaCl Single Crystals. II

S. Mendelson

J. Appl. Phys. 33, 2182 (1962); http://dx.doi.org/10.1063/1.1728924 (5 pages) | Cited 8 times

Online Publication Date: 11 June 2004

Full Text: | Download PDF

Show Abstract
The effects of surface sources, ozonized surface films, and surface dissolution, on the nature of plastic flow in NaCl single crystals were studied using etch pit and birefringence techniques. Slip steps increase in depth and become rounded as the slip band broadens. This is also often accompanied by a bending of the steps as much as 3° off 〈100〉. These effects are explained in terms of both internal and surface cross‐slip mechanisms. Ozonization of NaCl for short periods restricts the action of surface sources and for longer times provides barriers to the egress of dislocations. The ``Suzuki effect'' (formation of a work‐hardened surface layer) is found to depend upon the existence of a multiplicity of surface sources such as is produced by polishing on wet silk. Cross‐slip of screw dislocations at the surface is found to be greater and flow more homogeneous when the surface is dissolved during deformation. This is proposed as a likely mechanism to explain the ``Joffe effect'' when atmospheric variables are eliminated.

Properties of Bi2Te3‐Sb2Te3 Alloys

Mac J. Smith, R. J. Knight, and C. W. Spencer

J. Appl. Phys. 33, 2186 (1962); http://dx.doi.org/10.1063/1.1728925 (5 pages) | Cited 17 times

Online Publication Date: 11 June 2004

Full Text: | Download PDF

Show Abstract
The phase diagram for the Bi2Te3☒Sb2Te3 pseudo‐binary system is of the solid‐solution type, where the distribution coefficient k is equal to unity at 33.3 and 66.7 mole % Sb2Te3. The c lat tice parameter remains essentially constant across the diagram at 30.49 Å for both slowly crystallized and quenched alloys. For quenched alloys the a lattice parameter decreases almost linearly, from a value of 4.487 Å, for pure Bi2Te3, to a value of 4.275 Å, for pure Sb2Te3; however, a significant contraction from linear variation is found in slowly crystallized materials. EG diminishes in an essentially linear fashion from 0.16 eV, for pure Bi2Te3, to 0.12 eV at 24.2 mole % Sb2Te3 in both slowly crystallized and quenched materials. EG remains approximately constant from 24.2 to 66.7 mole % Sb2Te3 for slowly crystallized materials but continues to drop for quenched materials.

Contribution of Divacancies to the Magnetic Aftereffect in Nickel

Miles V. Klein and Helmut Kronmüller

J. Appl. Phys. 33, 2191 (1962); http://dx.doi.org/10.1063/1.1728926 (7 pages) | Cited 7 times

Online Publication Date: 11 June 2004

Full Text: | Download PDF

Show Abstract
This paper presents a theoretical discussion of the results to be expected from magnetic aftereffect measurements on Ni containing a sufficient number of divacancies. Calculations are made of the stabilization energy for 109°−〈001〉 and 180°−〈110〉 Bloch walls, the walls expected to be most mobile at low magnetic fields. The results show features characteristic of the 〈110〉 symmetry of the divacancies that are not shared by centers with 〈100〉 symmetry: (1) The existence of two relaxation modes having different relaxation times and different interaction energies with the magnetization vector, and (2) the nonvanishing of the stabilization force on 109° walls at large Bloch wall displacements.

Some Properties of Dirty Contacts on Semiconductors and Resistivity Measurements by a Two‐Terminal Method

George G. Harman and Theodore Higier

J. Appl. Phys. 33, 2198 (1962); http://dx.doi.org/10.1063/1.1728927 (9 pages) | Cited 4 times

Online Publication Date: 11 June 2004

Full Text: | Download PDF

Show Abstract
The surface and bulk properties of semiconductors have been studied by a two‐terminal method using dirty contacts. These contacts are defined as ones that are easily applied and removed and that are separated from the bulk by surface states, oxides, adhered gasses, and chemical films. The method essentially involves measuring the resistance‐voltage (R‐V) characteristics from the millivolt range up to about 100 V, and from these data the sample resistivity can be obtained directly. The effect of work damaging or etching the semiconductor surface can be readily evaluated. By observation of the shape of the R‐V curves it is possible to separate the bulk from the surface effects, calculate the surface barrier height and thickness from tunneling equations, and determine whether the barrier is a surface film or due to metal‐semiconductor contact potential difference. An application of the theory of electric field tunneling of Holm to the data of silicon carbide gives values of about 2.4 eV for the barrier height and 13 Å for the barrier width of the film on the surface. Sufficient information is included on electrode techniques and pitfalls so that the experimentalist can readily make the measurements. Efforts were concentrated on silicon carbide and silicon, but the techniques are applicable to all types of semiconductors.

Some Electrical Properties of the Porous Graphite Contact on p‐Type Silicon

George G. Harman, Theodore Higier, and Owen L. Meyer

J. Appl. Phys. 33, 2206 (1962); http://dx.doi.org/10.1063/1.1728928 (3 pages) | Cited 2 times

Online Publication Date: 11 June 2004

Full Text: | Download PDF

Show Abstract
Some unusual properties of the porous graphite contact on p‐type silicon are described. Current through the sample reaches a saturation level that is inversely proportional to the amount of adhered water vapor. Other gases such as ammonia, H2S, and HCl modify the shape or amplitude of the saturation current. There is a hysteresis loop in the low voltage region which is similar in appearance to that of a ferroelectric. The general method of measurement can be applied to studying the semiconductor surface as well as the contact phenomena. The possible applications include such devices as current regulators, humidity detectors, and surface‐barrier radiation detectors.
An electronic band model, which includes a trap‐dominated inversion layer, is presented to explain the phenomena. This model also integrates various conflicting theories of metal‐semiconductor contacts.

Fluctuations in Bistable Tunnel Diode Circuits

R. Landauer

J. Appl. Phys. 33, 2209 (1962); http://dx.doi.org/10.1063/1.1728929 (8 pages) | Cited 122 times

Online Publication Date: 11 June 2004

Full Text: | Download PDF

Show Abstract
In a series of preceding papers connection has been made between the physics of the activated jump over a static potential barrier and related computing devices. This paper treats a more complicated system, namely, the tunnel diode, which stores information in one of two possible dissipative states. The activated jump between dissipative states is analyzed. An idealized physical model is used to find which of the two states is the really preferred one, and to also evaluate the rate at which fluctuations bring about an approach to this preferred distribution. To simplify the analysis, the case where the tunnel diode is in series with a vacuum diode, rather than a resistor, is emphasized. For typical germanium Esaki diodes, activated jumps are improbable if the junction cross section exceeds 10−11 cm2.

p‐n Photovoltaic Effect in Cadmium Sulfide

H. G. Grimmeiss and R. Memming

J. Appl. Phys. 33, 2217 (1962); http://dx.doi.org/10.1063/1.1728930 (6 pages) | Cited 25 times

Online Publication Date: 11 June 2004

Full Text: | Download PDF

Show Abstract
So‐called front‐wall cells have been produced from CdS by inwards diffusion and subsequent selective etching away of the metal (Cu, Ag, Ni) producing the barrier layer. In particular, measurements were carried out of the dark current as a function of the voltage applied at different temperatures, of the shortcircuit current and the open‐circuit voltage as a function of tempering time, illumination and temperature; and the spectral sensitivity of the short‐circuit current was compared with the photoconductivity. It is shown that the photovoltaic effect can be explained not by a photoionization of electrons from the free metal but only by a p‐n photoeffect. On the assumption of an impurity conduction in the p region the properties of the p‐n photovoltaic effect are discussed.

Oxygen Precipitation Effects on the Deformation of Dislocation‐Free Silicon

J. R. Patel and A. R. Chaudhuri

J. Appl. Phys. 33, 2223 (1962); http://dx.doi.org/10.1063/1.1728931 (2 pages) | Cited 12 times

Online Publication Date: 11 June 2004

Full Text: | Download PDF

Show Abstract
The yield phenomenon observed in dislocation‐free crystals of silicon is critically dependent on prior heat treatment at 1000°C. Since oxygen is always present in crystals grown from quartz crucibles, heat treatment at 1000°C is believed to result in silicon‐oxygen clusters. It is shown that these clusters are effective sites for the generation of dislocations. Since the velocities of dislocations in heat‐treated and unheat‐treated crystals is not altered, the lowered maximum stress is largely due to the increased numbers of dislocations generated at a relatively early stage of the deformation process.

Electrification of the Dust of NaCl and KCl Whiskers

Stanisław Jakubiszyn, Jan Nikliborc, Henryk Wolniewicz, and Anna Szaynok

J. Appl. Phys. 33, 2224 (1962); http://dx.doi.org/10.1063/1.1728932 (3 pages)

Online Publication Date: 11 June 2004

Full Text: | Download PDF

Show Abstract
The dust of NaCl and KCl ordinary monocrystals shows greater charge per particle surface area than the dust of whiskers of these materials. The main reason is a greater cleavage surface of monocrystals than the cleavage surface of whiskers. The electrification of the dust is also influenced by the presence of contaminations in the crystal.

Effect of Crystal Perfection and Polarity on Absorption Edges Seen in Bragg Diffraction

H. Cole and N. R. Stemple

J. Appl. Phys. 33, 2227 (1962); http://dx.doi.org/10.1063/1.1728933 (7 pages) | Cited 39 times

Online Publication Date: 11 June 2004

Full Text: | Download PDF

Show Abstract
When a single crystal is used in an x‐ray spectrometer to display the continuous spectrum of an x‐ray tube by Bragg diffraction, absorption edges characteristic of the atomic species in the crystal may be observed. The change in diffracted power across an edge is fundamentally due to a change in the atomic scattering, especially the imaginary part of the scattering, of the species whose edge it is. This produces a sharp anomaly in both the structure factor and the absorption factor of the reflection being used. The explicit effect of the absorption is a function of the perfection of the crystal face; the change in the structure factor depends uniquely on (hkl) for noncentrosymmetric crystals. Thus, the absorption edges may be used as a tool to study perfection of crystal faces and polarity of polar axes. The K edges of Ga and As seen in reflection from (111) and (111) faces of perfect and mosaic crystals of GaAs are discussed both theoretically and experimentally in some detail. The linear absorption coefficient of GaAs and Ge have been measured in the course of the work.

Two‐Wave Shock Structures in the Ferroelectric Ceramics Barium Titanate and Lead Zirconate Titanate

C. E. Reynolds and G. E. Seay

J. Appl. Phys. 33, 2234 (1962); http://dx.doi.org/10.1063/1.1728934 (8 pages) | Cited 34 times

Online Publication Date: 11 June 2004

Full Text: | Download PDF

Show Abstract
Two ferroelectric ceramics, Pb (Zr0.52, Ti0.48)O3 with 1 wt % Nb2O5 and pure BaTiO3, were studied under shock wave loading. Plane shock waves having pressures ranging from 2 to 175 kbar were produced by high‐explosive driving systems. Measurements were made with a streak camera and shock‐detecting ferroelectric wafers. Both ceramic materials showed evidence of a two‐wave structure over a wide range of shock pressures. Results indicate that these two‐wave structures are characteristic of the elastic properties of the materials and are not due to phase transitions from ferroelectric to paraelectric. Pressure‐volume Hugoniot curves are presented for both materials. Electrical charge released to an external load by Pb (Zr0.52, Ti0.48)O3 with 1 wt % Nb2O5 was measured as a function of the pressure of a shock which was moving parallel to the axis of polarization. The primary mechanism in this reduction of polarization by shocks was found to be, not a phase transition or domain switching, but merely a reduction of the dipole moment by axial compression.

Critical Currents of Superconducting Nb 25 at.% Zr in High Magnetic Fields

Paul R. Aron and Harley C. Hitchcock

J. Appl. Phys. 33, 2242 (1962); http://dx.doi.org/10.1063/1.1728935 (3 pages) | Cited 7 times

Online Publication Date: 11 June 2004

Full Text: | Download PDF

Show Abstract
Measurements of critical current as a function of applied magnetic field were made on wire samples of Nb 25% at. Zr alloy at 4.22°K and 2.24°K. Curves are plotted for the cases for which the applied magnetic field is parallel to or perpendicular to the current direction. At 4.22°K, Hc(I → 0) for H perpendicular to I is found to be 70 kG, while for H parallel to I, it is found to be 82 kG. Reduction of the temperature to 2.24°K results in an increase of 13 kG in both cases. For H perpendicular to I, current densities of 4.4×104 A∕cm2 are observed at 60 kG and 4.22°K, and at 73 kG and 2.24°K. The effective demagnetization factor of the current carrying elements is 0.6 π. All measurements were made in a dc magnetic field with quasidirect currents by observation of the first resistive voltage to appear across the sample.

Maximum Field of Magnetized Disks and Cylinders

W. L. Mammel and S. P. Morgan

J. Appl. Phys. 33, 2244 (1962); http://dx.doi.org/10.1063/1.1728936 (3 pages)

Online Publication Date: 11 June 2004

Full Text: | Download PDF

Show Abstract
Tables and plots are given of certain definite integrals which facilitate the calculation of an upper bound on the field attainable with a given configuration of magnetic material, when the volume occupied by magnetic material consists of a combination of solid and hollow circular cylinders. Such estimates are useful in the design of powerful electromagnets.

Axial Quasi‐Brillouin Streams

M. H. Miller

J. Appl. Phys. 33, 2247 (1962); http://dx.doi.org/10.1063/1.1728937 (2 pages) | Cited 1 time

Online Publication Date: 11 June 2004

Full Text: | Download PDF

Show Abstract
A nonlaminar model for axially uniform electron streams is postulated and conditions for which the model is consistent with physical requirements are derived. The character of the current density profiles permitted by the model are discussed.

Influence of Crystal Growth on Domain Behavior in Yttrium Iron Garnet

R. A. Lefever, A. B. Chase, and K. A. Wickersheim

J. Appl. Phys. 33, 2249 (1962); http://dx.doi.org/10.1063/1.1728938 (2 pages) | Cited 2 times

Online Publication Date: 11 June 2004

Full Text: | Download PDF

Show Abstract
The behavior of magnetic domains in certain sections of single crystal yttrium iron garnet is related to growth history. Discontinuities in domain patterns coincide with lines of intersection of vicinal faces on the growth surfaces of the crystals, and regions delineated by the discontinuities show differences in magnetic behavior. The relationship of these effects to oriented strains resulting from growth processes is considered.

A Note on Logarithmic Rate Laws

P. T. Landsberg

J. Appl. Phys. 33, 2251 (1962); http://dx.doi.org/10.1063/1.1728939 (3 pages) | Cited 8 times

Online Publication Date: 11 June 2004

Full Text: | Download PDF

Show Abstract
An earlier theory for the rate law ψ(t) = β ln(1+t/τ) for strain, or q(t) = b−1 ln(1+t/τ) for chemisorption, is extended to cover a law for creep in materials which has been of interest recently, and takes the form ψ(t) = (β/ν)[(1+t/τ)ν−1]. The theory is formal in that it utilises only general concepts such as probability of success and attempt frequency.

Thermal Generation of Vacancies and Substitutional Sites in Crystalline Polymers

R. K. Eby

J. Appl. Phys. 33, 2253 (1962); http://dx.doi.org/10.1063/1.1728940 (4 pages) | Cited 21 times

Online Publication Date: 11 June 2004

Full Text: | Download PDF

Show Abstract
It is proposed that differences between the thermal expansion of the macroscopic and crystal lattice volumes of crystalline polymers might be attributed to thermal generation of new substitutional sites. Consideration of the molecular arrangements suggests that vacancies can be generated within the lamellae at molecular folds and ends and at any point along some helical molecules. Part of the displaced molecule can be transported by point defects to the crystal surface where new sites are formed. Analysis of the available data for polyethylene shows the concentration of new sites (Δn∕n) to be 2.2×10−2 at 100°C with an energy of creation of 0.2 eV and an entropy of creation of 5.3×10−4 eV∕°K per site. It is pointed out that the phenomenon could account for the so‐called glass transitions and is, therefore, closely related to the ``amorphous'' relaxations observed in crystalline polymers. The need for additional experiments is discussed.

A Simple Approximation for Classical Scattering at Large Angles

G. Leibfreid and O. S. Oen

J. Appl. Phys. 33, 2257 (1962); http://dx.doi.org/10.1063/1.1728941 (6 pages) | Cited 19 times

Online Publication Date: 11 June 2004

Full Text: | Download PDF

Show Abstract
The classical scattering for repulsive potentials is treated approximately by a matching method. The matching potential is a cut‐off potential which matches the value and the slope of the actual potential at the distance of closest approach. The method is thought to be successfully applicable to interatomic potentials of the screened Coulomb type. Comparison is made with machine results for an exponentially screened Coulomb potential. The differential cross section and the stopping power are compared for various energies. The agreement in the differential cross section is fair except for very low energy transfers. The stopping powers agree over a wide energy range with an error of only a few percent. Further comparison is made with another approximation given by Lehman and Shapiro.

Statistical Stability of the Preisach Diagram for Particles of γ‐Fe2O3

G. Bate

J. Appl. Phys. 33, 2263 (1962); http://dx.doi.org/10.1063/1.1728942 (7 pages) | Cited 33 times

Online Publication Date: 11 June 2004

Full Text: | Download PDF

Show Abstract
An experiment is described that enables the statistical stability of the Preisach diagram for interacting ferromagnetic particles to be tested. The samples used were of γ—Fe2O3 dispersed in a nonferromagnetic matrix to a nominal packing density of 20% by volume. Both oriented and nonoriented samples were examined. In every case the diagram was found to be stable over most of the area bounded by the major remanence loop. The region of departure from stability occurred at the higher levels of remanent intensity. A method of calculating the diagrams using only dc measurements is given. They are then used to predict accurately the ascending remanence curve after ac demagnetization for each sample. Finally, it is shown how the experimental ascending and descending remanence curves, which could not be explained by Wohlfarth's relations [J. Appl. Phys. 29, 595 (1958)] based on noninteracting particles, can be completely described in terms of the Preisach diagram.

Kinetic Interpretation of Non‐Newtonian Flow

D. A. Denny and Robert S. Brodkey

J. Appl. Phys. 33, 2269 (1962); http://dx.doi.org/10.1063/1.1728943 (6 pages) | Cited 16 times

Online Publication Date: 11 June 2004

Full Text: | Download PDF

Show Abstract
A new theoretical approach to the problem of describing the nonlinear behavior of non‐Newtonian materials is presented. The theory is based upon the assumption that the nonlinear characteristics can be associated with some structural change of the material, and that reaction kinetics can be applied to this breakdown. The theory allows treatment of time‐dependent (thixotropic) materials, and in a limiting form (equilibrium) can be applied to time‐independent materials. The theory has been applied to a variety of materials, including polymer latex solutions, polyethylene polymer melts, and a heavy thixotropic mineral oil. The results are gratifying.
Page 1 of 3 Pages Next Page | Jump to Page
Close
Google Calendar
ADVERTISEMENT

close