Christopher J.
Crissman
*abc,
Mianzhen
Mo
c,
Zhijiang
Chen
c,
Jie
Yang
cd,
Diego A.
Huyke
ce,
Siegfried H.
Glenzer
c,
Kathryn
Ledbetter
fg,
J. Pedro
F. Nunes
g,
May Ling
Ng
c,
Hengzi
Wang
c,
Xiaozhe
Shen
c,
Xijie
Wang
c and
Daniel P.
DePonte
*c
aUnited States Military Academy, West Point, NY 10996, USA. E-mail: chris@crissman.us
bDepartment of Applied Physics, Stanford University, Stanford, California 94305, USA
cSLAC National Accelerator Laboratory, Menlo Park, California 94025, USA. E-mail: deponte@slac.stanford.edu
dDepartment of Entrepreneurship, Innovation, and Strategy, Tsinghua University, Beijing, China
eDepartment of Mechanical Engineering, Stanford University, Stanford, California 94305, USA
fStanford PULSE Institute, SLAC National Accelerator Laboratory, Menlo Park, California 94025, USA
gDepartment of Physics and Astronomy, University of Nebraska-Lincoln, Lincoln, Nebraska 68588, USA
First published on 16th February 2022
We report on the design and testing of glass nozzles used to produce liquid sheets. The sheet nozzles use a single converging channel chemically etched into glass wafers by standard lithographic methods. Operation in ambient air and vacuum was demonstrated. The measured sheet thickness ranges over one order of magnitude with the smallest thickness of 250 nm and the largest of 2.5 μm. Sheet thickness was shown to be independent of liquid flow rate, and dependent on the nozzle outlet area. Sheet surface roughness was dependent on nozzle surface finish and was on the order of 10 nm for polished nozzles. Electron transmission data is presented for various sheet thicknesses near the MeV mean free path and the charge pair distribution function for D2O is determined from electron scattering data.
Desired sheet qualities are application-dependent, though good temporal and spatial stability are generally required. Soft X-ray transmission and MeV electron scattering require thickness in the range of 0.1 μm to 1 μm while thicker sheets in the range of 10 μm are preferable for optimal signal generation when considering the higher transmission of hard X-ray probes. Damage to the nozzle can occur when using high power laser sources and so sheet length must also be considered to maintain a safe working distance between the nozzle and the interaction point. For static structural biology applications, large biomolecules or viruses will be easier to image using fast X-ray or electron probes if the carrier fluid is submicron thickness thereby contributing less to the background scatter. For time resolved or dynamics studies in which the target will first be optically pumped, such as in the study of photosynthetic molecules, sheet flatness is necessary to eliminate any lensing effect which might result in an uneven illumination of the sample. In addition to flatness, surface roughness is of particular importance to high energy density (HED) science. In HED science, targets have typically been solid state with very carefully prepared surfaces.6 These targets are heated and compressed by optical absorption of high-intensity laser pulses. Surface roughness plays a role in absorption and surfaces that are optically smooth or better are desired.
The need for liquid sample delivery and, perhaps more importantly, flat jet technology is increasing as pulsed optical and X-ray laser facilities move to higher repetition rates. In the case of free electron lasers (FEL), higher repetition rate implies an increased sample consumption for liquid phase studies unless the sample can be reused or recirculated. By spreading out the incident X-ray flux over a wider area the intensity can be lowered, reducing the sample damage. In the case of high-powered laser facilities, the replacement of solid targets several times per minute (or even per second), may not be feasible and liquid targets may provide a suitable alternative to enable the production of energetic particle sources of ions, electrons, neutrons, and X-rays.7
Converging nozzles operate in a manner similar to the more commonly used colliding jet nozzles.1,8–14 Both rely on a transfer of momentum transverse to the direction of the jet travel to generate a pressure gradient that spreads out the liquid into a sheet. In Fig. 1C, the outlet of a converging nozzle is shown with a resulting sheet flowing from the outlet. The momentum generated by the liquid flow within the converging channel results in a transfer of momentum flux to a plane perpendicular to the flow of the channel. Ultimately, surface tension causes the jet to collapse and again, this inward momentum flux generates another sheet orthogonal to the first. This can continue on as a “chain” of sheets until breakup. Converging nozzles have been previously made individually with either multipart assembly or 3D printing.7–9 In this study, the fabrication and characterization of the glass nozzles made by wet isotropic etch are reported together with their application for transmission electron scattering experiments.
Liquid sheet nozzles can be fabricated by several methods with various respective merits. One of the most accessible to fabricate is a colliding jet in which cylindrical jets emanating from two aligned tubes collide to form a sheet.11 In this case, alignment can be challenging and pulsation of the liquid flow can cause fluctuations in sheet position or thickness. Lithographically formed channels have been used to aid in alignment, but this comes with another set of alignment challenges related to the detachment of the liquid from the chip surface. Koralek et al., mitigated this issue by using a lithographically formed nozzle that employs three channels – a central liquid channel and two gas orifices on opposing sides of the liquid channel to flatten out the liquid flow.2 This device can produce very thin sheets, to the limit of their reported measuring capability of 16 nm. Sheet thickness and the width to length ratios can be tuned by changing the gas and liquid flow rates, however the thickness reliance on both means that pulsation in flow rate causes pulsation of the sheet thickness. Gas load on vacuum systems and scattering of X-ray or optical probes from the gas is another disadvantage. The difficulty in fabricating single aperture converging nozzles lies somewhere between colliding and gas accelerated nozzles. Using a single orifice, such as with the converging nozzles examined here, eliminates the need for alignment and decreases clogging by using a larger nozzle aperture for similar flow rate colliding jets.
Converging nozzles can be produced by 3D printing or other soft methods in which case the exit orifice thickness, which must be non-zero for structural stability affects the nozzle performance. A straight section at the nozzle exit will decrease the ratio of spanwise to streamwise momentum flux and will result in a narrower and thicker sheet. In the limit that the exit length becomes arbitrarily large, the sheet becomes a cylinder. A 3D printed nozzle reported by Galinis et al. fabricated using a two-photon polymerization method produced very smooth and flat sheets of micron thickness in vacuum and ambient air.8 Theses sheets and those produced by the laser cut polyimide channels reported by Ha et al. produced similar results with thickness and length/width ratios both being somewhat smaller for the polyimide channels due to having closer to the ideal zero-length, “knife edge” exit.10 Razor blades were used to form the exit aperture in a device reported by Runge and Rosenberg for dye lasers. Dye laser nozzles were, for decades, produced by compression of a steel tube. In their simplest form, a hammer was the only tool needed for fabrication of these relatively large nozzles. A compressed-tube, dye laser nozzle with a razor blade exit reported by Watanabe et al. is perhaps closest to the geometry presented here but is limited in minimum size.9 Although lithography involves a more complex manufacturing process than hammering, highly repeatable structures in a wide variety of sizes are more consistently achieved. This method is not as well suited as 3D printing to producing single prototypes but can scale up to large quantities with ease, bringing the per-part costs down considerably.10
Exit channel dimensions and convergence angle were chosen by practical considerations. The smallest nozzle dimension chosen was 20 μm depth to limit clogging. Channel width, prior to converging to outlet, was then chosen to minimize sheet thickness in the channel width to channel depth ratio of 2.25 as given by Ha et al.10 A convergence angle of 60 degrees was chosen as a trade-off between sheet thickness and ease of manufacturing. Larger angles produce thinner sheets at comparable distance from the nozzle but are more sensitive to the location of the chip edge with uncertainties in both etch depth and cleave location contributing to a larger variation in nozzle width. Once the smallest practical nozzle design was determined, a set of larger nozzles was also designed to provide a range of sheet sizes from 100's of nm to 10's of microns. As a finishing step, some of the nozzle exit faces were polished with 3 μm lapping film to decrease roughness around the fluid exit. As an unintended consequence, this step increased the nozzle width and added some variation to the range or nozzle widths tested. Sheet thickness data was collected prior to polishing, and the effects of polishing nozzle outlets is further discussed with regards to sheet surface characteristics.
Although sheet thickness does not depend on flow rate or fluid properties, as shown in Results and discussion – sheet thickness, the sheet length and width are directly impacted by flow rate changes that would result from varying the pressure. For the nozzles tested, flow rates ranged from 3 to 5 ml min−1. Flow rates could be increased beyond the point at which the first sheet would fray and a chain would not form. All sheet measurements were recorded below this limit.
For all measurements, the nozzle was oriented vertically with the liquid flowing downwards from the nozzle. Nozzles were held in an O-ring sealed interface supplied by vendor Neptune Fluid Flow Systems. The fluid was collected within a reservoir that also fed through a filter and back into the HPLC pump. For vacuum operation, a temperature and pressure controlled collector/reservoir was purchased from vendor Innovative Research Solutions.3 In the case of ambient operation, a beaker sufficed as the reservoir.
Given the sheet's similarity to a thin film, a thin film interferometer, Filmetrics model F40, was used to measure the sheet thickness, h. The illumination from the interferometer was focused into the object plane of an optical microscope. Once a sheet was generated, the sheet was maneuvered into the focal plane of the interferometer through use of micron precision X, Y, Z stages as depicted in Fig. 2. By moving the sheet upward, thickness measurements were taken in intervals that were selected for each individual sheet, as each sheet length varied with respect to its nozzle outlet dimensions. Sheet thickness measurements were made in ambient air only.
The sheet length, lj, and sheet width, wj, were measured via positioning in the microscope field of view. The nozzle exit and sheet apex positions were recorded from the stage positions. These positions were recorded for all ambient experiments and a subset was repeated in vacuum for comparison.
The data can be reasonably fit by the relation
(1) |
Fig. 3 Nozzle outlets formed the geometry shown as a result of isotropic etch and bonding of two glass wafers. Nozzle outlet width, w0, and depth, d, were measured at their maximum distances. |
Fig. 4 Sheet thickness as measured along the sheet centerline is shown for 11 different nozzles. The width of the nozzle exit aperture varied from 60 μm to 140 μm. Nozzle exit aperture depths ranged from 16–18.4 μm. Solid lines represent the thicknesses calculated by eqn (1). For each respective nozzle outlet width, w0. |
The work of Ha et al. used laser-cut polyimide channels which resulted in rectangular channel geometry. They showed that sheet thickness is directly related to the ratio of streamwise to spanwise momentum flux which is itself a function only of the nozzle converging angle.10
For comparison, flow simulations for nozzles of width 100 μm and 45 μm width and 20 μm depth were carried out for both rectangular channel geometry and the rounded channels due to isotropic etch. Flux momentum ratios are in agreement to within 2% for the two nozzle geometries. Nozzle shape did show one difference related to thickness – while the rectangular outlet geometry is dependent on the aspect ratio, that is not found to be true for geometry formed by isotropic etching.
There is no flow rate term in the above expression for sheet thickness. This can be demonstrated by white light interference. Repeating colored bands indicate a change in thickness along the sheet. These fringes were spatially stable with decreasing flow rate, as shown in Fig. 6. This is a valuable feature of sheets created by converging nozzles, as any changes to the pressure or flow rate of the fluid entering the nozzle will not affect the thickness of the targeted area.
The difference in thickness between any two points on a sheet can be estimated from the color at each location and therefore provides a measure of sheet roughness. The red band near the center of the sheet shown in Fig. 7 is 1.75 mm from the nozzle exit and corresponding to a thickness of 640 nm. Further downstream 100 nm in the orange band, the thickness is 600 nm thick. This implies orange-colored fluctuations on the red band at the sheet midpoint are under 40 nm. Between these ridges and valleys, the sheet surface roughness is too small to measure by this method. A separation of the red/blue/green channels for color data collected along the centerline of the sheet is shown in Fig. 7 and shows smoothly oscillating curves except at the points of surface disturbances. The largest of the red peaks shown was limited by detector saturation.
Fig. 7 The RGB channels from the centerline of the image in Fig. 9B with the detector saturation level set to 1 on the vertical axis. The curves are smoothly periodic except at distances corresponding to small ripples on the sheet surface. |
To explore the possibility that rapid fluctuations were being averaged out over the camera exposure time, several shorter exposure snapshots were compared. These images as shown in Fig. 5, were taken with an Andor ICCD camera with a 10 ns gated window and a HeNe laser centered at 632.8 nm. No additional short-timescale disturbances within the sheet were found. The use of monochromatic light also allowed for the direct calculation of thickness change, with each fringe period indicating a change in thickness of 257 nm. This is in line with the thickness data collected using thin film reflectometry and gives a higher resolution measurement of thickness change along the sheet.
Fig. 10 shows the transmission of the direct beam as a function of the sheet jet thickness which was obtained by scanning along the vertical direction of the same sheet jet. Note that the transmission data for thickness below 800 nm are consistent with the previously reported data for a thinner liquid sheet produced by the gas-accelerated jet.3 The experimental data were fit to a simple exponential function showing an 1/e attenuation length of approximately 0.94 μm. This implies that the elastic mean-free-path, (MFP) of MeV electrons in liquid D2O is approximately 0.94 μm.
Fig. 10 The transmission of the direct electron beam or main beam as a function of the heavy water sheet thickness. The red solid line is the exponential fit to the data. |
The effect of multiple scattering events for such a sheet jet can be evaluated with the ratio of the high-Q scattering signal over the low Q scattering signal. Here Q is the momentum transfer and is defined as Q = 4πsinθ/λ, with θ and λ represent the scattering angle and the de Broglie wavelength of the electrons. Fig. 11a shows the radial averaged intensity signal of the diffraction patterns taken for four different sheet jet thicknesses. The rapid drop in scattering intensity at low Q is due to a through hole in the detector center to pass through the direct beam. The scattering signal that is below 4 Å−1 decreases noticeably with increasing the sheet jet thickness, which is indicative of the multiple scattering effect. Fig. 11b shows the ratio of the scattering signal between Q regions of 6.75–9 Å−1 and 0.75–3 Å−1 as a function of sheet jet thickness. Note the ratio was normalized to the value at the thinnest position. As indicated, the data shows clearly a two-trend behavior with the inflection point located approximately at the MFP value. These results facilitate the design of future pump-probed experiments in choosing the appropriate jet thickness to avoid multiple elastic scattering effect. Assuming the total elastic scattered electron intensity Is decreases as
Is = I0(1 − e−x/l) | (2) |
The data were of sufficiently high scattering angle that a pair distribution function could be determined. Fig. 12 shows the charge pair distribution function as described by Yang.15 In the 1st Born approximation, the angle averaged electron scattering intensity I(Q) for a liquid is given by
(3) |
(4) |
(5) |
Fig. 12 The D2O charge pair distribution function is shown for data collected at four different sheet thickness and compared to simulation for H2O at 250 K. |
As a first demonstration of the nozzles' utility, they we used at SLAC's MeV-UED facility. MeV electrons were shot through a liquid sheet consisting of heavy water and electron transmission and scattering data were recorded. The collection of electron scattering data in water without noise from a container and for a target thickness that is otherwise difficult to consistently achieve allowed for the determination of the D2O pair distribution function. The thickness of the liquid sheets used was from 34% below to 36% above the MFP of electrons in D2O. Above the MFP, multiple scattering made data analysis more difficult but reasonably good CPDF fits were found for smaller values. The optimal thickness for sheets was not determined; vacuum limitations prevented the thinnest sheets from being probed. However, the sheet range presented in this study covers a range including from 1 MFP to well below the peak in single elastic scattering fraction. Extrapolating from the existing data we would expect 4% of the elastic scattering to be from multiple scattering events at 250 nm, the thinnest sheets reported in this study.
Wet etch lithography seems well suited to nozzle fabrication. While some deviation in nozzle width was found, these deviations can be kept to within a range that there is reasonable coverage of the intended thickness range for each nozzle. Most importantly, the sheet generated is consistent, per nozzle; once characterization is complete, each nozzle can be used to repeatedly re-create the same sheet as would be expected from a simple model for sheet thickness given as a function of nozzle distance and nozzle outlet cross sectional area. All nozzles fabricated by this method were able to produce liquid sheets within the expected range of sizes. Nozzle surface roughness was highly variable for unpolished sheets. For applications requiring 10's of nm surface roughness or better, the polishing and characterization steps carried out in this study may likely be eliminated by improved fabrication. Perhaps lithographically defined chip boundaries, laser dicing, use of crystalline silicon or quartz rather than amorphous glass, or some combination will suffice. Most important is that these nozzles can be produced in large batches and based on standard microfluidic chip layouts thereby reducing costs and simplifying interfacing.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/d1lc00757b |
This journal is © The Royal Society of Chemistry 2022 |