How to Build a Laser Death Ray

Material Response to Laser Radiation

Rough empirical values used for laser machining and welding are described below (Lasers and their Applications, A. Sona, ed, chapter "Machining with lasers" by Dieter Roess)

For a sufficiently low power flux the thin surface layer will be heated to the fluid state but will stay beyond the evaporation temperature, while the solid-fluid interface will slowly progress into the bulk material by heat conduction. In iron the typical progress rate is about 10-2 cm ms-1, and the power flux for this situation may be around 106 W cm-2. [...]

At a somewhat higher power flux, between 106 to 108 Wcm2 [sic], the thin absorbing surface layer is heated up to its evaporation temperature before the solid-fluid interface has progressed appreciably into the material by heat conduction. Thus, with continuing laser power, a less than μ-thick layer of material will continually evaporate, with the material-air interface progressing into the material. Typically a gas jet develops [...]. The gas jet ejects also part of the molten material, thus that the progress rate of the hole is faster than with evaporation of all the material. [...]

At a still higher flux of 109 to 1010 Wcm-2, after initial evaporation of the surface layer, the gas jet is thermally ionized and absorbs most of the incident radiation, which is such blocked away from the material. The surface layer explodes with an ultrasonic jet, its temperature may rise beyond 10 105 ° K [sic], its pressure beyond 103 at.

Below is a calculator to estimate the effects of a heat ray incident on a surface of known material properties (or select from a set of pre-defined materials).

Beam Power, P
Beam Energy, E
choose one


Beam Duration (optional)
Beam Diameter at Target, S
Number of Pulses:
Time between Pulses:
Initial (Ambient) Temperature, T0
Ambient Pressure, p0
Geometry Thickness, X
Surface properties Absorptivity, a
Input custom material properties
Density, ρ
Heat of Fusion, Hf
Heat of Vaporization, Hv
Specific Heat, C
Melting Temperature, Tm
Boiling Temperature, Tb
Thermal Conductivity, κ
Viscosity Coefficient, μ @ Tμ
Surface Tension, λ @ Tλ
Molar mass
Speed of sound
Ultimate strength
Shear modulus
Adiabatic index
(adiabatic index is usually 1.666... for stuff that evaporates into single atoms, or 1.2 for stuff that evaporates into complex molecules)
Choose pre-defined material

The physics used in the calculator is described in the following sections. Keep in mind that many of these calculations are very approximate, more accurate treatments of heat flow into a surface generally require complex numerical simulations.

Light contains energy. When that light is absorbed by a material, its energy goes toward heating the material. However, the faster the material gains heat, the faster it sheds heat. If the rate of heat loss gets high enough to equal the rate of heat gain, the death ray can no longer raise the temperature of the surface. Otherwise, the surface continues to gain heat and its temperature rises. The temperature to which a beam of light can bring a material depends on the intensity of the light and the rates and mechanisms by which the material responds to the heat.

There are generally considered to be three mechanisms by which heat moves around: heat conduction, heat convection, and radiation. Heat conduction is what you get when heat moves through a material, transfered from atom to atom by the atoms bumping into each other and exchanging energy. Convection is the movement of heat by moving the material it is stored in. Radiation transfers heat by the emission or absorption of electromagnetic waves. When a surface is exposed to a death ray, it is being heated by radiation.

To estimate the temperature after illumination by a death ray, first note some of the basic parameters of the beam. The section on diffraction allows to to find the diameter of the spot, S, into which you can focus your beam. If you know the time average power, P, of the beam, the total intensity (power divided by area) of the beam is

I = P / (π S2 / 4).

Find the absorptivity, a, of the illuminated material. This is a number between 0 and 1 that determines the fraction of light that is absorbed. The intensity absorbed by the surface is

Ia = I × a.
Experience with laser machining tells us that the absorptivity of materials under intense laser irradiation is typically high. From this reference we find
  • "At room temperature, most metals are highly reflective of infrared energy, the initial absorptivity can be as low as 0.5% to 10%. But the focused laser beam quickly melts the metal surface and the molten metal can have an absorption of laser energy as high as 60~80%."
  • "Non-metallic materials are good absorbers of infrared energy. They also have lower thermal conductivity and relatively low boiling temperatures. Thus the laser energy can almost totally transmitted into the material at the spot and instantly vaporize the target material."

The absorptivity becomes more complicated once the laser intensity is high enough to flash part of the irradiated surface into plasma. The plasma will absorb some or all of the laser light, which then conducts or re-radiates the laser energy to the surface. This is good when the material is highly reflective, since now the laser light is not diectly interacting with the surface. Unfortunately, some of the laser energy can also be re-radiated elsewhere. Further, the plasma can heat the air in front of it to a plasma, and then the air plasma can absorb the beam which will insulate the surface of the material from the beam. You can get interesting phenomena, such as laser-supported detonation or radiation waves rushing as plasma waves through the air away from the target in the direction of the laser source, but this is bad for getting energy to the target. Normally, for laser intensities high enough to cause propagating plasma waves away from the target, it is best to use short pulses where the laser light turns off before the plasma wave can start and then repeat once the plasma has cleared.

You will also want to know some properties of the material being illuminated that deal with heat movement, specifically:
  • ρ, the density. This is the mass of material divided by its volume.
  • κ, the thermal conductivity. This is the power of heat flow (heat energy crossing a surface per unit time) per unit area divided by the temperature gradient (the temperature gradient is the change in temperature ΔT for a given change in distance ΔX). If you know the temperture gradient, the thermal conductivity tells you how much power is flowing P = κ A ΔT/ΔX. Conversely, if you know how much heat power is flowing through a material, you know how fast the temperature changes. Note that this gives all the fundamentals of heat conduction - heat power flows from high temperature to low temperature at a rate proportional to the temperature change per unit distance and proportioanl to the area of the surface across which the heat can flow.
  • Cp, the specific heat. This is the heat energy E a material gains divded by its mass M and divided by its change in temperature ΔT. If you know how much the temperature has changes in a material, you can find out the change in heat energy in it by E = M Cp ΔT. Similarly, if you know how much energy a material has gained, you can determine the change in temperature ΔT = E / (M Cp). Note - if the material reaches its melting point or boiling point, you can't use this equation! Heat will go into melting or vaporizing the material rather than changing its temperature (or, if cooling, the condensing or solidifying material will release heat). Use the specific heat to find how much energy goes into the material up to the melting point, use the heat of fusion (below) to find out how much energy is required to melt the material, and then use the specific heat again to describe the material's rise in temperature once all of it has melted.
  • Hv, the heat of vaporization. This is the heat energy E needed to vaporize a material divided by its mass M, Hv = E/M.
  • Hf, the heat of fusion. This is the heat energy E needed to melt a material divided by its mass M, Hv = E/M.
  • Tm, the melting temperature.
  • χ = κ/(ρ Cp), the thermal diffusivity.
  • λ, the surface tension, at a known temperature Tλ.
  • μ, the coefficient of viscosity, at a known temperature Tμ.
Also note the ambient pressure p0 and ambient temperature T0. If the target of the death ray is of finite thickness, note the thickness X of the slab being heated.


One of the first significant means heat is removed from the illuminated spot is via conduction. The heat from the beam immediately begins to diffuse into the surrounding colder material. As described in the definition of κ heat flow is driven by temperature gradients - the faster the temperature changes, the more heat power flows across a given area. As before, the equation is

P = κ A ΔT/ΔX.
After being illuminated for a time t, heat will have penetrated a distance D into the material approximately given by
D = √ 4 χ t .
The front of the heated area will advance at a speed of
vD = √ χ/t 
= 2 χ / D.
Note that the heated area initially advances rapidly, but eventually gets large enough to insulate the illuminated spot from the surrounding cold material. Also note that these equations break down once the material starts melting (or evaporating), since when this happens heat energy is being used to melt or evaporate the material rather than raising the temperature.

Conduction from a stationary source in 3D

One limiting case deserves mention. For a point source of heat on the surface of a half-infinite slab, the temperature distribution at a distance r and time t is

T = [P/(4 π κ r)][1-erf(r/sqrt(2 χ t))] + T0.
Here, erf(X) is known as the error function of X. The main detail we need to know is that for large X, erf(X)→1, and for small X, erf(X)→0. This means that for small distances or large times, the temperature approaches a constant
T = P/(2 π κ r) + T0.
Conversely, if we know the temperature at some distance r, we know the power being conducted away
P = 2 π κ r (T - T0).
These relations only hold if the spot size S is much smaller than the thickness of the slab X, and while the heat diffusion distance D remains less than X. When D is not much larger than S, these relations are only approximate. In particular, for thin slabs, power cannot be conducted away through the bottom of the slab so the temperature will keep rising indefinitely (neglecting radiation). Still, for thick slabs, you can estimate the maximum temperature due to conduction
T = Pa/(2 π κ S/2) + T0
and if you know the temperature, the conducted power is very approximately
Pc = 2 π κ (T - T0) S/2.

Conduction from a moving source in 1D

Another limiting case will be useful later on. If a moving plane of heat is passing through the material at a uniform velocity of v, adding a constant intensity of heat Ic to the material as it passes, the temperature at a distance x in front of the heat plane is

T(x) = T0 + (Ic/(Cp ρ v)) exp(-v x / χ)
This will come in useful when the laser death ray is boring into the surface, drilling a hole at constant velocity while adding enough heat to raise the temperature of the surface high enough to remove the surface.

Conduction from a moving source in 3D

Related to both of the above, a point heat source or spherically symmetric heat source of constant power P passing through a bulk volume of material at velocity v will create a temperature field around it of

T(r,x) = T0 + (P/(4 π κ r)) × exp(-(r+x)v/(2 χ))
where r is the distance from the heat source and x is the distance in front of the heat source along the direction of motion. Again, if the laser is drilling into a material this lets us estimate how much heat is leaking out in all directions.


One of the methods by which hot surfaces dispose of their heat is by radiating it as electromagnetic waves. The total radiated intensity (power per unit area) is given by

Ir = e σ T4,
where T is the temperature in Kelvins, e is the emissivity (a number between 0 and 1 that determines how well the surface can radiate), and σ is a constant of nature
σ = 5.67×10-8 W m-2 K-4.
Note that σ is a constant - the letters after it just tell you want units to use in the equation. If you use temperature in kelvins, for example, it gives you intensity in watts per square meter.

For this rough estimate, you can set the absorptivity equal to the emssivity (this is strictly always true for absorption and emission at the same wavelength, but both a and e vary with wavelength). If we assume that radiation is the only mechanism through which a surface sheds its heat, then the surface will gain heat until the intensity of radiation equals the intensity of beam energy being absorbed. In other words, Ia = Ir. This allows us to solve for the maximum temperature the surface will reach

T = (I/σ)1/4.

This estimate ignores other mechanisms of shedding heat, such as heat conduction from the surface into the bulk material, or transfer of heat to the surrounding air or water which is then carried away by a rising plume of this warmed fluid. When a surface gets very hot, radiation is the dominant method of disposing of heat. However, for cooler temperatures (like "just" the heat of a stove burner on high), the transfer of heat by conduction to cooler matter in contact with the material can be faster than radiation.


Once a material reaches its melting temperature, any further input of heat goes toward melting the material rather than raising its temperature (although keep in mind that if not all of the material is at the melting temperature, heat will be conducted in to the melting zone from hotter areas and away from the melting zone toward cooler areas). It takes an energy

Emelt = Hf ρ Vmelt
to melt a volume Vmelt of the material. Thus the energy needed to melt a volume V of material initially at temperature T0 will be
E = Emelt + (Tm - T0) Cp ρ V.

If the melt is not removed from the hole, the heat from the laser will need to conduct through the melt to reach the melt front. This will raise the temperature of the surface above the melting temperature. When as much heat is being conducted away from the melt front as is being conducted in to it, the melt front will no longer advance. A rough idea of this condition can be found using the relations for conduction into a thick slab, giving a maximum melt distance

rmelt = Pa/(2 π κ (Tm-T0)).


A material can evaporate at any temperature (although if it is at a temperature less than the melting temperature, it is called sublimation rather than evaporation). If the temperature is low, the rate of evaporation will be very low, but the vaporization rate increases rapidly with temperature, until, at high intensities, most of the energy might be going in to evaporation. The energy needed to evaporate a volume of material Vvap is

Eevap = Hv ρ Vvap

From this, we can see that the energy needed to completely evaporate a volume V of material at temperature T will be

E = Eevap + Emelt + (T - T0) Cp ρ V.

The difficulty is in finding the final temperature T of the evaporating material.

The amount of material that evaporates depends on the vapor pressure of the material at the temperature to which the laser heats the top layer. The vapor pressure is, as its name implies, the pressure of the vapor emitted at that temperature. If the vapor pressure inside the molten material ever exceeds the pressure on the material at that point, it can push the material aside as a vapor bubble spontaneously grows. When this happens with a pot of water on the stove, we call it boiling. This tells us that the vapor pressure is one atmosphere (1 × 105 Pa) at the boiling temperature Tb. However, at different pressures, the boiling temperature will be different - this is why it is not sufficient to say that the laser has to heat the material to its boiling temperature to vaporize it. Instead, we can use the Clausius-Clapeyron relation to find the vapor pressure p at any temperature T if you know the vapor pressure pref at any other temperature Tref.

p = pref exp[Hv Mmol (1/Tref - 1/T) / R]
p = 1 × 105 Pa exp[Hv Mmol (1/Tb - 1/T) / R]
or equivalently
T = 1/[1/Tref - (R / (Hv Mmol)) ln(p/pref)]
T = 1/[1/Tb - (R / (Hv Mmol)) ln(p / 1 × 105 Pa)]
where R is the gas constant
R = 8.314472 J/(Mol K)
and Mmol is the mass of one mole of the material (be careful - this is usually given in grams, and if you work in SI units, you will need to divide by 1000 to get kg).

If the vapor pressure is less than the ambient pressure p0, the vapor cannot push the surrounding air out of the way. At most, it slowly diffuses away, remaining close to equilibrium and re-absorbing almost as much of the material as evaporates. As soon as the vapor pressure exceeds p0, the evaporate will be blown away by its excess pressure at speed vv.

When the vapor jet develops, it will take energy to get the jet going. This energy comes from the internal energy of the hot gas, resulting in a temperature drop as the jet expands and accelerates. For vapors escaping from a laser-heated surface, this acceleration, temperature drop, and expansion happen rapidly across a discontinuous boundary called the Knudsen layer. The temperature and density drop across this layer is

Tv / T = (√(1+q2) - q)2
ρv / ρs = √(T / Tv) ((m2+0.5) e erfc(m) - m/√π) + 0.5 (T / Tv) (1 - √π m e erfc(m))
q = √π (m/2) ((γ-1)/(γ+1))
m = ηv √(γ/2)
ρs = p Mmol / (R T).
ρs is the density of the vapor at the surface due to the vapor pressure, before the density drop across the Knudsen layer

γ is the adiabatic index of the gas. Many evaporated vapors are made up of single atoms, in which case γ=5/3. If the evaporate is composed of molecules, γ will be smaller but larger than 1 - for molecules made of two atoms γ is approximately 1.4 and for atoms of three or more atoms γ is approximately 1.3 to 1.2, although gamma tends to decrease as the temperature increases.

The function erfc(x) is known as the complimetary error function of x

erfc(x) = 1 - erf(x).
There are a number of approximations of the error function you can use if you do not have this capability on your calculator
erfc(x) ≈ (0.3480242 t - 0.0958798 t2 + 0.7478556 t3) exp(-x2)
t = 1/(1 + 0.47047 x)

ηv is the Mach number of the vapor (how fast the jet is traveling with respect to the speed of sound in the jet). This needs to be found self-consistently using the relation

0 = (p / p0) (Tv / T) (ρv / ρ) - 1 - γ ηv (cv/c0)(r + √(1 + r2))
r = ((γ + 1)/4) ηv (cv/c0).
cv = √(γ R Tv / Mmol)
c0 = √(γ0 R T0 / Mmol 0)
This can be done by, for example, plotting the above function and finding where it crosses zero (more advanced root finding techniques, such as bisection or Newton's method, will allow determining the solution faster). The Mach number has a maximum value of 1, if you find a value larger than 1 set ηv = 1. The quantity cv has the physical meaning of the speed of sound in the vapor, and c0 is the speed of sound in the ambient air. γ0 is the adiabatic index of the air (1.4 for an atmosphere of diatomic molecules like that of Earth) and Mmol 0 is the molar mass of the air (0.0288 kg/mol for Earth's air).

In addition, by the ideal gas law we can determine the pressure of the vapor jet

pv = ρv R Tv / Mmol.
We also know the speed of the jet
vv = cv ηv.

To estimate the behavior of the vaporizing target material in the presence of a jet, we consider a sample of the target material moving into the beam at a speed vi such that the position of the melt-vapor interface remains at a fixed position x=0. If we consider the material moving into the interface to evaporate, it must be balanced by the material moving away as the fuly developed jet. In a time Δt, a length of target material vi Δt moves into the interface and a portion of the jet of length ve Δt leaves. The mass of material moving into the interface to be vaporized is thus ρ A vi Δt and the mass of material leaving as a jet is ρv A vv Δt. Since these masses must be equal, we find

vi ρ = vv ρv
vi = vv ρv/ρ.

At this point, we can describe all relevant properties of the melt and the jet in terms of the temperature T of the interface between the melt and the vapor, but we do not yet know the interface temperature. To determine this, we now balance the energy flow between the interface and the jet, similar to the momentum and mass flow. The energy leaving the interface as a jet is the heat energy of the jet plus the kinetic energy of the jet

Ev = ρv A vv Δt (γ R Tv / Mmol + vv2/2)
where we used ρv vv2/2 as the kinetic energy per unit volume and ρv γ R Tv / Mmol as the heat energy per unit volume. Since the evaporation of a given mass M of material adds an energy Hv M, and since the resulting vapor, if at rest and at the same temperature, has a heat energy of γ (M/Mmol) R T, we determine that the heat energy in the liquid at melting temperature must be γ (M/Mmol) R T - Hv M. Therefore, the energy carried into the evaporating interface is
El = ρ A vi Δt (γ R T / Mmol - Hv + vi2/2)
= ρ A vi Δt (vi2/2 - Hv) + ρv A vv Δt γ R T / Mmol.
There is also energy being conducted away from the interface into the melt.
Ec = κ A Δt ΔT/ΔX.
At steady state, this must equal the energy needed to raise the solid at ambient temperature to a liquid at the vaporization temperature
Ec = ρ A vi Δt (Hf + (T-T0) Cp).
Finally, there is energy being added by the laser and energy being radiated away as heat from the high temperature interface
Ea = Ia A Δt
Er = e σ T4 A Δt.
When we require the energy in to equal the energy out, we find
El + Ea = Er + Ec + Ev
which evaluates to
Ia - e σ T4 - κ ΔT/ΔX + ρ vi (vi2/2 - Hv) + ρv vv (γ R (T - Tv) / Mmol - vv2/2) = 0.
For steady state drilling, this becomes
Ia - e σ T4 - ρ vi (Hf + (T-T0) Cp) + ρ vi (vi2/2 - Hv) + ρv vv (γ R (T - Tv) / Mmol - vv2/2) = 0.
All of the variables are either determined by the material properties, the incident laser radiation, or the temperature of the evaporation front. Thus, we can use this to solve for the temperature at which evaporation occurs, T. This equation has no analytic solution for T, but it can be solved numerically using a computational root finding algorithm or with pencil and paper by graphing the function and finding where it crosses 0. Once T is determined, we will have found along the way the speed at which the hole is being drilled vi, along with other useful properties such as the pressure p and jet velocity vv.

Melt ejection

The high pressure of the hot evaporated vapor jetting from the molten surface can act like a piston to squirt the melt layer off to the side, where it is then broken up into droplets and blown out of the hole by the high speed jet. This removes the molten material without evaporating it. Since melting takes much less energy than vaporization, the more melt ejection that occurs, the faster the laser can drill at a given energy.

Melt ejection can be the dominant effect in laser machining when the spot sizes are small and powers are in the multi-kilowatt range. For larger spot sizes or higher powers, melt ejection becomes less important compared to vaporization. For lower power lasers, the vapor pressure is insufficient to overcome the surface tension and melt ejection does not work.

To determine the rate of melt ejection, it is first necessary to estimate the thickness of the liquid layer. This can be determined from the equation for heat conduction in front of a moving heat source. The intensity being conducted into the surface has to be enough to raise its temperature from the initial temperature T0 to the temperature of the surface T when advancing at a speed of v

Ic = ρ v (Hf + (T-T0) Cp).
The thickness of the melt layer is then approximately
Dm = (χ/v) ln[Ic/((Tm-T0) Cp ρ v)]
= (χ/v) ln[(T-T0 + Hf/Cp)/(Tm-T0)]

Melt ejection is driven by the pressure ramming down onto the molten layer. It is resisted by the viscosity of the melt, which retards shear motion of the liquid - in this case, of the melt layer squirting out to the sides due to the pressure pushing down on it. Further, the central pressure of the beam is resisted by the ambient pressure and the surface tension from the liquid at the edge of the spot. The surface tension gives an effective pressure

pst = 4 λ / S
where λ is the surface tension.

It has been found that the rate at which the melt squirts out the bottom of the hole and up onto the sides can be reasonably approximated by assuming a rate of mass loss of

dM/dt = π Dm3 (p-p0-pst) / (3 μ)
where μ is the coefficient of viscosity. This adds a speed of vm to the interface speed
v = vm + vi
where vi is the interface speed from vaporization alone and
vm = 4 Dm3 (p-p0-pst) / (3 ρ μ S2).
This requires the interface speed and melt layer to be determined self-consistently
[4 (p-p0-pst) / (3 ρ μ S2)] Dm3 + vi Dm - χ ln[(T-T0 + Hf/Cp)/(Tm-T0)] = 0.
This gives a quartic equation, which can be solved analytically.

In a time Δt, this mass removal carries away a heat energy of

Ee = (π S2/4) vm Δt ρ (T-Tm) Cp
and thus an intensity (energy per time per area) of
Ie = vm ρ (T - Tm) Cp.
This can be directly added to the energy balance equation
Ia - e σ T2 - ρ vi (Hf + (T-T0) Cp) + ρ vi (vi2/2 - Hv) + ρv vv (γ R (T - Tv) / Mmol - vv2/2) - vm ρ (T - Tm) Cp = 0.

One complication is that the viscosity coefficient μ and surface tension λ depend on temperature. Since you have to solve for the rate of drilling numerically anyway, this isn't too much more of a problem. The viscosity for many materials is approximately given by

μ = μ0 exp[Hv Mmol (1/T - 1/Tμ) / (3 R)]
and the surface tension by
λ = λ0 + k (Mmol/ρ)2/3 [Tλ - T]
where k = 2.1×10-7 J / (K mol2/3).

These results are for steady state beams. In practice, it is found that melt ejection can be significantly although not overwhelmingly enhanced by illuminating the target surface with a relatively long but low intensity pulse to make a sizeable melt pool, and then deliver a short but very high intensity driver pulse to blast the melt pool away. This is not yet accounted for in the calculator.


The beam is initially incident on the target. A thin layer of target material is flashed to plasma, and the plasma absorbs the rest of the pulse. At this point, the plasma is at extremely high pressure, much greater pressure than the strength of the target material.
The expanding plasma launches a blast wave. A shell of hot vapor blasts out into the air, while a shell of super-hot, highly compressed material propagates into the target. Pressures are so high the target material flows like a fluid.
As the blast front propagates into the material, it slows down.
Eventually, the blast wave has done so much work pushing through the target material and has spread out so much that the pressure is no longer far in excess of the material strength. The flowing blast wave, no longer able to bull its way through the target material, is instead re-directed along the crater walls and splashes out into the air.
Once the fluidized material that made up the blast wave has squirted out of the crater along the crater sides and into the air, what is left is a permanent crater, surrounded by a layer of compressed and deformed material, and beyond that the mostly undamaged target material.

If the vapor pressure exceeds the material strength of the target, the pressure of the vapor can drive an expanding cavity in the target material. We can approximate the effects of this pressure driven penetration by the beam by modifying the theory of long-rod tank penetrators. This starts with Bernoulli's equation for the dynamic pressure of a moving fluid

p = ½ ρ v2.
Consider a purely fluid target moving the target into the beam at a speed v chosen so that the point of deepest penetration remains the same. Then we have to balance the dynamic pressure of the fluid with the pressure of the vapor in order to keep the interface stationary. This corresponds to the laser boring a hole into the fluid at a speed of v.

Of course, the targets we are interested in are not fluids, but rather solids. It is found that a reasonable approximation is to add a constant to the dynamic pressure

p = ½ ρ v2 + Kcrater.
This approximation performs reasonably well for estimating the penetration of shape-sharge explosive jets, anti-tank long-rod kinetic penetrators, hypervelocity impacts with debris in orbit by satellites, and bullets from high powered smokeless powder centerfire rifles. In all cases, the penetration is largely due to the dynamic pressure of the projectile balancing the dynamic pressure of the target across the interface, with dynamic pressures in excess of the target's strength pushing aside the target material. Here, we replace the dynamic pressure of the projectile with the vapor pressure generated by the beam. This means that the beam will bore into the target at a speed of
v = √(2 (p - Kcrater) / ρ).

The material constant Kcrater has a simple physical interpretation. It is the energy needed to gouge a cavity in the material divided by the volume of the cavity. Kcrater is approximately 3 to 4 times the compressive strength Kcomp, or can be found as

Kcrater = (2/3) Kcomp × (1+ln(2 G/Kcomp))
where G is the shear modulus. From this definition, we know the volume of the crater excavated by the beam
V = E / Kcrater.
For a given cross-sectional area of the hole A and depth x of the hole, the volume is V = A x. In a time t we can then write V/t = A x/t = A v, where v is the speed of the deformation-driven hole into the target material. Since E = P t, we therefore have P/Kcrater = A v. Since the hole will be more or less cylindrical, the radius r of the hole will be approximately r = √(A/π). Thus
r = √(P/(π Kcrater v)).

More generally for pulses of finite duration, consider the depth of the hole made while the laser is on

h = v t.
This describes a cylinder directly illuminated by the laser of depth h and diameter S (where, again, S is the spot size). If the pressure is sufficient to expand the cavity by a distance rc, the volume of the final cavity is
E / Kcrater = π [ (4/3) rc3 + (π S / 2 + h) rc2 + (S2/2 + S h) rc + h S2/4 ]
This requires solving a cubic equation, which has analytical solutions. The total depth D is then
D = h + rc
and the width of the hole is
r = S + 2 rc.
This latter version has the advantage of smoothly transitioning from the steady state case for long pulses to the circular cratering limit of very short pulses.

Using deformation to drive a hole into your target has an advantage in that it doesn't matter if most of the beam gets absorbed by the plasma and little of that energy goes into heating the target. The high pressure of the plasma itself can act to deform the target just as well as evaporated vapors. However, you do have to avoid laser supported plasma waves rushing away from the target where the front of the wave is too far away from the target for its pressure to support the target. Consequently a beam composed of a series of pulses, where each pulse is too short to initiate a plasma wave, would be preferred. Each pulse is incident on the back of the previous crater, making a deeper and deeper hole.

If subsequent pulses are incident on the bottom of the hole before the crater from the previous pulse has had time to fully expand, the pressure can pile up, leading to a shallower but wider hole. This is accounted for in the calculatr by using the time average power over the entire pulse train to compute the depth of the hole if this results in a wider hole than adding together the individual pulses. Otherwise, the depth is the sum of the depths of each individual pulse.

Since these relations are approximate, at edge cases we might end up with inconsistent results, such as a hole that is narrower than the beam. In this case we enforce consistency by expanding the hole to the beam diameter and reducing the rate at which the hole propagates into the target so that E = V Kcrater.


Thanks to Jez Weston and Anthony Jackson for stimulating discussions and useful information.

Further Reading

  1. M. von Allmen, "Laser drilling velocity in metals", J. Appl. Phys. 47, 5460 (1976)
  2. S. Basu and T. DebRoy, "Liquid metal expulsion during laser irradiation", J. Appl. Phys. 72, 3317 (1992)
  3. C. L. Chan and J. Mazumder, "One-dimensional steady-state model for damage by vaporization and liquid expulsion due to laser-material interaction", J. Appl. Phys. 62, 4579 (1987)
  4. B. N. Chichkov, C. Momma, S. Nolte, F. von Alvensleben, and A. Tünnermann, "Femptosecond, picosecond and nanosecond laser ablation of solids", Appl. Phys. A 63, 109-115 (1996)
  5. A. Luft, U. Franz, A. Emsermann, J. Kaspar, "A study of thermal and mechanical effects on materials induced by pulsed laser drilling", Appl. Phys. A 63, 93-101 (1996)
  6. R. K. Ganesh and A. Faghri, "A generalized thermal modeling for laser drilling process-I. Mathematical modeling and numerical methodology", Int. J. Heat mass Transfer 40, 3351-3360 (1997)
  7. A. Kar and J. Mazumder, "Two-dimensional model for material damage due to melting and vaporization during laser irradiation", J. Appl. Phys. 68, 3884 (1990)
  8. C. Körner, R. Mayerhofer, M. Hartmann, and H. W. Bergmann, "Physical and material aspects in using visible laser pulses of nanosecond duration for ablation", Appl. Phys. A 63, 123-131 (1996)
  9. Ladislav Grad and Janez Mozina, "Laser pulse shape influence on optically induced dynamic processes", Appl. Surf. Sci. 127-129, 999-1004 (1998)
  10. Leon J. Radziemski and David A. Cremers, Ed. Laser-Induced Plasmas and Applications, Marcel Dekker, New York, 1989
  11. Matter, Energy, and Radiation Hydrodynamics
  12. Sven-Olov Roos, "Laser drilling with different pulse shapes", J. Appl. Phys. 51, 5061 (1980)
  13. Pablo Solana, Phiroze Kapadia, John Dowden, William S. O. Rodden, Sean S. Kudesia, Duncan P. Hand and Julian D. C. Jones, "Time dependent ablation and liquid ejection processes during the laser drilling of metals", Optics Communications 191, 97-112 (2001)
  14. A. Sona, Ed. Lasers and their Applications, Gordon and Breach, New York, London, Paris, 1976
  15. B. Wu, Y. C. Shin, H. Pakhal, N. M. Laurendeau, and R. P. Lucht, "Modeling and experimental verification of plasmas induced by high-power nanosecond laser-aluminum interactions in air", Phys. Rev. E, 76, 026405 (2007)
  16. A. Tate, "A Theory for the Deceleration of Long Rods After Impact", J. Mech. Phys. Solids 15, 387-399 (1967)
  17. A. Tate, "Further Results in the Theory of Long Rod Penetration", J. Mech. Phys. Solids 17, 141-150 (1969)
  18. Y. Laurence Yao, Hongqiang Chen, and Wenwu Zhang, "Time scale effects in laser material removal: a review", Int. J. Adv. Manuf. Technol. (2004)
  19. Ya. B. Zel'dovich and Yu. P. Raizer, Physics of Shock Waves and High-Temperature Hydrodynamic Phenomena, Dover, Mineola N.Y., 2002
  20. W. Zhang, Y. L. Yao, and K. Chen, "Modelling and Analysis of UV Laser Micromachining of Copper", Int. J. Adv. Manuf. Technol., 18, 323-331 (2001)
  21. A. D. Zweig, "A thermo-mechanical model for laser ablation", J. Appl. Phys. 70, 1684 (1991)

Back to main death ray page