Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “TUNNEL EFFECT”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

629 records · Page 35

Stochastic thermodynamic cycles of a mesoscopic thermoelectric engine

Here, we analyze a steady-state thermoelectric engine, whose working substance consists of two capacitively coupled quantum dots. One dot is tunnel-coupled to a hot reservoir serving as a heat source, the other one to two electrically biased reservoirs at a colder temperature, such that work is extracted under the form of a steady-state current against the bias. In single realizations of the dynamics of this steady-state engine autonomous, four-stroke cycles can be identified. The cycles are purely stochastic, in contrast to mechanical autonomous engines which exhibit self-oscillations. In particular, these cycles fluctuate in direction and duration and occur in competition with other spurious cycles. Using a stochastic thermodynamic approach, we quantify the cycle fluctuations and relate them to the entropy produced during individual cycles. We identify the cycle mainly responsible for the engine performance and quantify its statistics with tools from graph theory. We show that such stochastic cycles are made possible because the work extraction mechanism is itself stochastic instead of the periodic time dependence in the working-substance Hamiltonian which can be found in conventional mechanical engines. Our investigation brings new perspectives about the connection between cyclic and steady-state engines.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Spontaneous Localization at a Potential Saddle Point from Edge State Reconstruction in a Quantum Hall Point Contact

Quantum point contacts (QPCs) are an essential component in mesoscopic devices. Here, in this work, we study the transmission of quantum Hall edge modes through a gate-defined QPC in monolayer graphene. We observe resonant tunneling peaks and a nonlinear conductance pattern characteristic of Coulomb-blockaded localized states. The in-plane electric polarizability reveals the states are localized at a classically unstable electrostatic saddle point. We explain this unexpected finding within a self-consistent Thomas-Fermi model, finding that localization of a zero-dimensional state at the saddle point is favored whenever the applied confinement potential is sufficiently soft compared to the Coulomb energy. Our results provide a direct demonstration of Coulomb-driven reconstruction at the boundary of a quantum Hall system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Effects of Different Types of Entrances on Natural Ventilation in a Subway Station

In this study, the natural ventilation of a horizontal entrance of a typical subway station is investigated based on numerical simulations and experiments. In addition, a renormalization group k-e model (RNG k-e model) is applied to compute both the internal and external airflow patterns. Computational fluid dynamics (CFD) simulations are validated based on the experimental results. Furthermore, a multiple variable regression model is employed to study how the different parameters affect the internal airflow rates of the subway station statistically, according to the experimental results. For this typical model, the parameter importance can be ranked as follows: (1) outdoor wind speed; (2) flow resistance of the subway station; (3) height of the wind catcher; and (4) length of the wind catcher. To understand the detailed pressure distribution of the horizontal entrance of the subway station with and without the wind catcher, 3D numerical simulations are conducted for different scenarios. We attempt to alter the size of the wind catcher (including the length and height) to study the characteristics of the pressure on the horizontal entrance of the subway station under outdoor wind-driven conditions. The pressure distributions for the entrance for different scenarios are compared and analyzed. Finally, the interactions between the internal and external flows are investigated by changing the resistance of the subway station. When the internal flow resistance is changed, the pressure coefficient (Cp) of the entrance is different as well. Hence, the Cp is not only affected by the outdoor environment, but is also influenced by the internal airflows.

CFD↗

On the scaling of wind turbine rotors

This paper formulates laws for scaling wind turbine rotors. Although the analysis is general, the article primarily focuses on the subscaling problem, i.e., on the design of a smaller-sized model that mimics a full-scale machine. The present study considers both the steady-state and transient response cases, including the effects of aerodynamic, elastic, inertial, and gravitational forces. The analysis reveals the changes to physical characteristics induced by a generic change of scale, indicates which characteristics can be matched faithfully by a subscaled model, and states the conditions that must be fulfilled for desired matchings to hold. Based on the scaling laws formulated here, the article continues by considering the problem of designing scaled rotors that match desired indicators of a full-scale reference. To better illustrate the challenges implicit in scaling and the necessary tradeoffs and approximations, two different approaches are contrasted. The first consists in a straightforward geometric zooming. An analysis of the consequences of zooming reveals that, although apparently simple, this method is often not applicable in practice, because of physical and manufacturing limitations. This motivates the formulation of scaling as a constrained optimal aerodynamic and structural matching problem of wide applicability. Practical illustrations are given considering the scaling of a large reference 10 MW wind turbine of about 180 m in diameter down to three different sizes of 54, 27, and 2.8 m. Results indicate that, with the proper choices, even models characterized by very significant scaling factors can accurately match several key performance indicators. Additionally, when an exact match is not possible, relevant trends can at least be captured.

17 WIND ENERGY↗

Nanoscale Imaging of Magnetotransport around a Circular 𝑝−𝑛 Junction in Graphene

Magnetoresistance studies of 2D systems are often shaped by the motion of electrons that occupy spatially confined wave functions, such as topological edge modes and disorder-induced bound states. Directly probing how such states form and behave in situ , under applied currents, provides a clear way of connecting microscopic physics to the macroscopic transport response. Here, in this Letter, scanning tunneling potentiometry is used to probe the local, current-induced electrochemical potential of carriers in graphene near circular 𝑝−𝑛 junctions in an out-of-plane magnetic field ranging from 0 to 1.4 T. These measurements provide detailed information about the motion of carriers at the nanometer scale, revealing how it evolves with increasing field. The electrochemical potential displays distinct patterns, such as dipoles, spirals, and concentric disks in weak, moderate, and high fields, respectively. The size and orientation of these patterns can be used to understand how local carrier dynamics change in different transport regimes as well as to directly extract physical parameters such as the electron mean free path and cyclotron diameter.

36 MATERIALS SCIENCE↗

Transmon Qubit Using Sn as a Junction Superconductor

Superconductor qubits typically use aluminum–aluminum oxide tunnel junctions to provide the nonlinear inductance. Junctions with semiconductor barriers make it possible to vary the superconductor material and explore beyond aluminum. We use InAs semiconductor nanowires coated with thin superconducting shells of β-Sn to realize transmon qubits. By tuning the Josephson energy with a gate voltage, we adjust the qubit frequency over a range of 3 GHz. The longest energy relaxation time, T 1 = 27 μs, is obtained at the lowest qubit frequencies, while the longest echo dephasing time, T 2 = 1.8 μs, is achieved at higher frequencies. Finally, we assess the possible factors limiting coherence times in these devices and discuss steps to enhance performance through improvements in materials fabrication and circuit design.

B-Sn↗

Time-Resolved Spectroscopy of Insulator-Metal Transitions: Exploring Low-Energy Dynamics in Strongly Correlated Systems. Final Report

Transition metal oxides (TMOs) have attracted much attention in recent years. The strong coupling between charge, spin, and orbital degrees of freedom creates novel functional properties. The research at the interplay between ferromagnetism, ferroelectricity and conductivity provides a promising way to the development of future spintronic devices. The interface characterization of multifunctional oxide heterostructures is crucial to understand deeply the magnetic response, spin injection and spin transport. In year one, we characterized optically the magnetic property with the method of magnetization-induced second-harmonic generation (MSHG), which is a sensitive technique to selectively probe the interface where the spatial inversion symmetry is broken. The interface of n-type BaTiO 3 (BTO)/La 0.67 Sr 0.33 MnO 3 (LSMO) forms a Schottky barrier and displays a good diode effect. By applying an external electric field across the sample, we observe the magnetoelectric (ME) coupling effect to modulate the interface magnetization. The voltage dependent magnetic contrast sharply vanishes at positive voltage (applied to LSMO contact) and reveals a ferromagnetic (FM)-to- antiferromagnetic (AFM) transition occurring at the interface, while the magneto-optic Kerr effect (MOKE) indicates that the bulk property is not affected. The strain mediated mechanism and ferroelectric (FE) polarization induced charge medicated mechanism are both discussed, but they do not play an important role in the observed magnetization transition. A new mechanism is proposed - the injected minority spins through strong Hund's interaction with the local magnetic moments at positive voltage weaken the double-exchange interaction of nearby Mn e g electrons, thereby reducing the ferromagnetic (FM) ordering. The dominant superexchange interaction of localized t 2g electrons finally leads to the AFM configuration. By reducing the electron carrier concentration of BTO, a positive voltage shift of the transition is observed. The oxygen richer sample with lower carrier concentration requires a higher voltage to increase the hole doping level of LSMO to induce the interfacial AFM phase. By inserting a thin (1nm) La 0.5 Ca 0.5 MnO 3 layer between BTO and LSMO, the interface magnetic response from LSMO is enhanced. In year two, a non-multiferroic heterostructure SrTiO 3 (STO)/ La 0.5 Ca 0.5 MnO 3 (LCMO)/La 0.67 Sr 0.33 MnO 3 (LSMO) was studied. In this system, the interface FM-AFM transition is observed at negative voltage. The STO layer may act as a hole-donating layer, which induces the hole transfer to LSMO making it overdoped, while the LCMO interlayer at the 0.5 hole doping level displays complicated CE type AFM phase. Thus the AFM phase is dominant at the interface and electron accumulation is required to induce the FM phase by applying a negative voltage. For oxygen rich sample, the magnetic MSHG loop reappears at high positive voltage (around +2.5 V). At high doping level (x > 0.5), LSMO exhibits A-type AFM order, consisting of AFM superexchange coupling between adjacent layers and intralayer FM double-exchange coupling. The less tensile strain of STO and LSMO due to a smaller lattice mismatch as compared to BTO would result in a weaker interlayer coupling . Thus, interfacial A-type AFM order can occur in STO/LCMO/LSMO heterostructure at high hole doping concentration. The study of ME coupling of multifunctional oxide heterostructures is definitely important to the development and performance of future spintronic devices. The results presented here regarding the interfacial magnetic transition controlled by an external electric field offer some new perspective to study the spin transport of magnetic tunnel junctions.

36 MATERIALS SCIENCE↗

Integrated Circuit Design for Radiation-Hardened Charge-Sensitive Amplifier Survived up to 2 Mrad

According to the continuous development of metal-oxide semiconductor (MOS) fabrication technology, transistors have naturally become more radiation-tolerant through steadily decreasing gate-oxide thickness, increasing the tunneling probability between gate-oxide and channel. Unfortunately, despite this radiation-hardened property of developed transistors, the field of nuclear power plants (NPPs) requires even higher radiation hardness levels. Particularly, total ionizing dose (TID) of approximately 1 Mrad could be required for readout circuitry under severe accident conditions with 100 Mrad around a reactor in-core required. In harsh radiating environments such as NPPs, sensors such as micro-pocket-fission detectors (MPFD) would be a promising technology to be operated for detecting neutrons in reactor cores. For those sensors, readout circuits should be fundamentally placed close to sensing devices for minimizing signal interferences and white noise. Therefore, radiation hardening ability is necessary for the circuits under high radiation environments. This paper presents various integrated circuit designs for a radiation hardened charge-sensitive amplifier (CSA) by using SiGe 130 nm and Si 180 nm fabrication processes with different channel widths and transistor types of complementary metal-oxide-semiconductor (CMOS) and bipolar CMOS (BiCMOS). These circuits were tested under γ–ray environment with Cobalt-60 of high level activity: 490 kCi. The experiment results indicate amplitude degradation of 2.85%–34.3%, fall time increase of 201–1730 ns, as well as a signal-to-noise ratio (SNR) of 0.07–11.6 dB decrease with irradiation dose increase. These results can provide design guidelines for radiation hardening operational amplifiers in terms of transistor sizes and structures.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Surface facet dependence of competing alloying mechanisms

Metal alloys are ubiquitous in many branches of heterogeneous catalysis, and it is now fairly well established that the local atomic structure of an alloy can have a profound influence on its chemical reactivity. While these effects can be difficult to probe in nanoparticle catalysts, model studies using well defined single crystal surfaces alloyed with dopants enable these structure–function correlations to be drawn. The first step in this approach involves understanding the alloying mechanism and the type of ensembles formed. In this study, we examined the atomic structure of RhCu single-atom alloys formed on Cu(111), Cu(100), and Cu(110) surfaces. Our results show a striking difference between Rh atoms alloying in Cu(111) vs the more open Cu(100) and Cu(110) surface facets. Unlike Cu(111) on which Rh atoms preferentially place-exchange with Cu atoms in the local regions above step edges leaving the majority of the Cu surface free of Rh, highly dispersed, homogeneous alloys are formed on the Cu(100) and (110) surfaces. These dramatically different alloying mechanisms are understood by quantifying the energetic barriers for atomic hopping, exchange, swapping, and vacancy filling events for Rh atoms on different Cu surfaces through theoretical calculations. Density functional theory results indicate that the observed differences in the alloying mechanism can be attributed to a faster hopping rate, relatively high atomic exchange barriers, and stronger binding of Rh atoms in the vicinity of step edges on Cu(111) compared to Cu(110) and Cu(100). Furthermore, these model systems will serve as useful platforms for examining structure sensitive chemistry on single-atom alloys.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Surface dependence of electronic growth of Cu(111) on MoS 2

Scanning tunneling microscopy shows that copper deposited at room temperature onto a freshly exfoliated MoS 2 surface forms Cu(111) clusters with periodic preferred heights of 5, 8, and 11 atomic layers. These height intervals correlate with Fermi nesting regions along the necks of the bulk Cu Fermi surface, indicating a connection between physical and electronic structures. Density functional theory calculations of freestanding Cu(111) films support this as well, predicting a lower density of states at the Fermi level for these preferred heights. This is consistent with other noble metals deposited on MoS 2 that exhibit electronic growth, in which the metal films self-assemble as nanostructures minimizing quantum electronic energies. Here, we have discovered that it is critical for the metal deposition to begin on a clean MoS 2 surface. If copper is deposited onto an already Cu coated surface, even if the original film displays electronic growth, the resulting Cu film lacks quantization. Instead, the preferred heights of the Cu clusters simply increase linearly with the amount of Cu deposited upon the surface. We believe this is due to different bonding conditions during the initial stages of growth. Newly deposited copper would bond strongly to the already present copper clusters, rather than the weak bonding, which exists to the van der Waals terminated surface of MoS 2 . The stronger bonding with previously deposited clusters hinders additional Cu atoms from reaching their lowest quantum energy state. The interface characteristics of the van der Waals surface enable surface engineering of self-assembled structures to achieve different applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Photothermal alternative to device fabrication using atomic precision advanced manufacturing techniques

The attachment of dopant precursor molecules to depassivated areas of hydrogen-terminated silicon templated with a scanning tunneling microscope (STM) has been used to create electronic devices with subnanometer precision, typically for quantum physics experiments. This process, which we call atomic precision advanced manufacturing (APAM), dopes silicon beyond the solid-solubility limit and produces electrical and optical characteristics that may also be useful for microelectronic and plasmonic applications. However, scanned probe lithography lacks the throughput required to develop more sophisticated applications. In this work, we demonstrate and characterize an APAM device workflow where scanned probe lithography of the atomic layer resist has been replaced by photolithography. An ultraviolet laser is shown to locally and controllably heat silicon above the temperature required for hydrogen depassivation on a nanosecond timescale, a process resistant to under- and overexposure. STM images indicate a narrow range of energy density where the surface is both depassivated and undamaged. Modeling that accounts for photothermal heating and the subsequent hydrogen desorption kinetics suggests that the silicon surface temperatures reached in our patterning process exceed those required for hydrogen removal in temperature-programmed desorption experiments. A phosphorus-doped van der Pauw structure made by sequentially photodepassivating a predefined area and then exposing it to phosphine is found to have a similar mobility and higher carrier density compared with devices patterned by STM. Lastly, it is also demonstrated that photodepassivation and precursor exposure steps may be performed concomitantly, a potential route to enabling APAM outside of ultrahigh vacuum.

nanoscale devices↗

Towards Low-Temperature CVD Synthesis and Characterization of Mono- or Few-Layer Molybdenum Disulfide

Molybdenum disulfide (MoS 2 ) transistors are a promising alternative for the semiconductor industry due to their large on/off current ratio (>10 10 ), immunity to short-channel effects, and unique switching characteristics. MoS 2 has drawn considerable interest due to its intriguing electrical, optical, sensing, and catalytic properties. Monolayer MoS 2 is a semiconducting material with a direct band gap of ~1.9 eV, which can be tuned. Commercially, the aim of synthesizing a novel material is to grow high-quality samples over a large area and at a low cost. Although chemical vapor deposition (CVD) growth techniques are associated with a low-cost pathway and large-area material growth, a drawback concerns meeting the high crystalline quality required for nanoelectronic and optoelectronic applications. This research presents a lower-temperature CVD for the repeatable synthesis of large-size mono- or few-layer MoS 2 using the direct vapor phase sulfurization of MoO 3 . The samples grown on Si/SiO 2 substrates demonstrate a uniform single-crystalline quality in Raman spectroscopy, photoluminescence (PL), scanning electron microscopy (SEM), atomic force microscopy (AFM), X-ray photoelectron spectroscopy (XPS), and scanning transmission electron microscopy. These characterization techniques were targeted to confirm the uniform thickness, stoichiometry, and lattice spacing of the MoS 2 layers. The MoS 2 crystals were deposited over the entire surface of the sample substrate. With a detailed discussion of the CVD setup and an explanation of the process parameters that influence nucleation and growth, this work opens a new platform for the repeatable synthesis of highly crystalline mono- or few-layer MoS 2 suitable for optoelectronic application.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effect of Controlled Artificial Disorder on the Magnetic Properties of EuFe 2 (As 1– x P x ) 2 Ferromagnetic Superconductor

Static (DC) and dynamic (AC, at 14 MHz and 8 GHz) magnetic susceptibilities of single crystals of a ferromagnetic superconductor, EuFe 2 ( As 1 – x P x ) 2 ( x = 0.23), were measured in pristine state and after different doses of 2.5 MeV electron or 3.5 MeV proton irradiation. The superconducting transition temperature, T c ( H ) , shows an extraordinarily large decrease. It starts at T c ( H = 0 ) ≈ 24 K in the pristine sample for both AC and DC measurements, but moves to almost half of that value after moderate irradiation dose. Remarkably, after the irradiation not only T c moves significantly below the FM transition, its values differ drastically for measurements at different frequencies, ≈16 K in AC measurements and ≈12 K in a DC regime. We attribute such a large difference in T c to the appearance of the spontaneous internal magnetic field below the FM transition, so that the superconductivity develops directly into the mixed spontaneous vortex-antivortex state where the onset of diamagnetism is known to be frequency-dependent. We also examined the response to the applied DC magnetic fields and studied the annealing of irradiated samples, which almost completely restores the superconducting transition. Overall, our results suggest that in EuFe 2 ( As 1 – x P x ) 2 superconductivity is affected by local-moment ferromagnetism mostly via the spontaneous internal magnetic fields induced by the FM subsystem. Another mechanism is revealed upon irradiation where magnetic defects created in ordered Eu 2 + lattice act as efficient pairbreakers leading to a significant T c reduction upon irradiation compared to other 122 compounds. On the other hand, the exchange interactions seem to be weakly screened by the superconducting phase leading to a modest increase of T m (less than 1 K) after the irradiation drives T c to below T m . Our results suggest that FM and SC phases coexist microscopically in the same volume.

36 MATERIALS SCIENCE↗

Mechanical sensors for ultraheavy dark matter searches via long-range forces

Dark matter candidates with masses around the Planck scale are theoretically well-motivated, and it has been suggested that it might be possible to search for dark matter solely via gravitational interactions in this mass range. In this work, we explore the pathway towards searching for dark matter candidates with masses around the Planck scale using mechanical sensors while considering realistic experimental constraints, and develop analysis techniques needed to conduct such searches. These dark matter particles are expected to leave tracks as their signature in mechanical sensor arrays, and we show that we can effectively search for such tracks using statistical approaches to track-finding. We analyze a range of possible experimental setups and compute sensitivity projections for searches for ultraheavy dark matter coupling to the Standard Model via long-range forces. We find that while a search for Planck-scale dark matter purely via gravitational couplings would be exceedingly difficult, requiring ∼80 dB of quantum noise reduction with a 100 3 array of devices, there is a wide range of currently unexplored dark matter candidates which can be searched for with already existing or near-term experimental platforms.

Qin, Juehang [Rice Univ., Houston, TX (United Stat↗

Novel and effective surface passivation for high efficiency n- and p-type Silicon solar cell

The project objective was to develop a novel Si surface passivation method using chalcogens, sulfur (S) and/or selenium (Se), as passivating elements, to withstand industry-standard high temperature contacting and metallization schemes for p-type Si based passivated emitter and rear contact (p-PERC) solar cells. The back surface passivation of PERC cells has been improved drastically with the invention and successful application of an Al 2 O 3 passivation layer. However, the front n + diffused junction surface is still poorly passivated by the standard amorphous silicon nitride (SiNx) anti-reflection coating (ARC) layer. This project sought to address the passivation challenges of both front n+ emitter and undiffused p-Si back surface. Improved p-PERC solar cell performance with open circuit voltage (V OC ) > 680 mV and efficiency of 22% were targeted to validate superior defect passivation properties as compared to standard SiO 2 / Al 2 O 3 passivation. During this project, we systematically investigated process-structure-properties-performance relationships of this novel advanced defect passivation approach. The S/Se passivation was carried out by reacting industrial Czochralski (Cz) Si wafers in H 2 S and H 2 Se gases in a chemical vapor deposition (CVD) reactor at temperatures up to 700°C. After an exhaustive optimization of the process parameters (temperature, time, and gas concentration), we established an optimized process and demonstrated extremely low surface recombination velocities (SRVs) of 1.5 cm/s and 8 cm/s on n-type and p-type Si, respectively, by S-passivation. In-depth surface and interface characterization were performed using soft x-ray and photoelectron spectroscopies (XPS, UPS, XES), combined with capacitance-voltage-frequency (C-V-f) measurements, to decipher the surface chemical/electronic structure and interface defect state densities. These measurements provided critical understanding of the defect passivation mechanism and elucidated the presence of surface S-Si bonds, a reduction of surface dipoles, and low interface state densities (D it ) < 10 11 cm -2 ev -1 . We also found that the Se-passivation is inferior to the S-passivation (by at least one order of magnitude in SRV). Application of the optimized S-passivation to the n+ diffused emitter surface led to a low surface recombination current density, J0 ≈ 40 fA/cm 2 (~ 1/4 of the industry-standard SiNx-passivation), and high implied V OC (686 mV) in p-PERC solar cell structures. The S-passivation process also was found to improve the bulk quality of the p-type Si, better than the SiO 2 or Al 2 O 3 passivation processes. After successful demonstration of efficient passivation of Si surface defects by S, we extensively studied the air, thermal, and illumination stability of the passivation structure. S-passivation itself degrades in air due to competing reactions with moisture and oxygen to form oxides, which can be eliminated by a SiNx capping layer (also acting as a anti-reflective coating). After SiNx process optimization, we demonstrated illumination and thermally stable S-passivation with SRV < 5 cm/s and J 0 < 80 fA/cm 2 . These enhancements in Si passivation, incorporated into p-PERC cells, achieved an efficiency of 19.93% with V OC = 649 mV, using manufacturing metallization and contacting schemes. The low cell performance (cell V OC is much less than the implied V OC = 686 mV, anticipated from surface passivation) was identified due to degradation of S-passivation during the metal firing step (out-diffusion of S from the Si interface to the SiNx surface). The S-passivation of Si surfaces shows significant promise with excellent passivation quality, essential for high performance (high V OC , high efficiency) solar cells. Integration of this innovative defect passivation into devices, however, demands further development of the capping layer, low temperature (<700°C) metallization process, and/or engineering of advanced device structures. Surface passivation-dominated advanced Si solar cells, such as tunnel oxide passivated contacts and Si heterojunctions, are increasingly of interest due to their high-performance potential and will have a growing photovoltaic market share in the near future.

14 SOLAR ENERGY↗

Electronic spectroscopy and ionization potentials for YbOH and YbOCH 3

The polyatomic molecules YbOH and YbOCH 3 have been recognized as being of potential value for spectroscopic experiments that explore charge-parity and time reversal symmetry violation effects. These measurements require very high precision, which, in turn, will necessitate that the molecules be manipulated at ultra-cold temperatures. Both YbOH and YbOCH 3 have electronic transitions that appear suitable for laser cooling ($\tilde{A}$ 2 Π 1/2 -$\tilde{X}$ 2 Σ + and $\tilde{A}$2E 1/2 -$\tilde{X}$ 2 A 1 , respectively) but the currently available spectroscopic data is not sufficient to determine the extent to which population leaks may compromise the optical cooling processes. A further complication is that the quantum states of interest for these measurements will need to be selectively populated. The $\tilde{A}$-$\tilde{X}$ band systems of both YbOH and YbOCH 3 show evidence of vibronic perturbations, such that there are unassigned vibronic features at energies that are just above the origin bands. Here, in the present study we have recorded spectra for the $\tilde{A}$ 2 Π 1/2 -$\tilde{X}$ 2 Σ + transition of jet-cooled YbOD to facilitate the vibronic assignments. In addition, spectra for the $\tilde{B}$ 2 Σ + - $\tilde{X}$ 2 Σ + transition of YbOH were recorded, establishing the origin band at 20473.8 cm -1 . Previously, the reaction of Yb with CH 3 OH has been used to generate gas-phase YbOCH 3 . As this reaction also yields YbOH, there have been complications in spectroscopic studies of YbOCH 3 due to overlap of the $\tilde{A} -\tilde{X}$ band systems. To identify specific regions of overlap, resonantly enhanced two-photon ionization spectra were recorded using mass-resolved detection of the YbOH + and YbOCH 3 + ions. These data confirmed the overlap of vibronic bands near 17640 and 17680 cm -1 . Two-photon ionization spectroscopy also provided accurate ionization energies of IE(YbOH)=45788(10) and IE(YbOCH 3 )=45283(10) cm -1 . The IE for YbOH is relevant to problems encountered in previous attempts to determine the bond dissociation energy of YbOH + .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dynamic Behavior of Natural Seep Vents: Analysis of Field and Laboratory Observations and Modeling (Final Scientific/Technical Report)

In this project, we have analyzed data collected by the U.S. Department of Energy (DOE), National Energy Technology Laboratory (NETL) in a high pressure water tunnel (HPWT) and data from two research cruises to natural seeps in the Gulf of Mexico to adapt and validate a numerical model to predict the dynamics of natural seeps in the deep oceans. The HPWT data include video observations of the shrinkage rate of individual methane and natural gas bubbles under simulated deep-water conditions. Field data were collected during two cruises by the Gulf Integrated Spill Research (GISR) Consortium led by Texas A&M University and funded by the Gulf of Mexico Research Initiative (GoMRI). These data included in situ observations from a remotely operated vehicle (ROV) of gas bubbles at two natural seep sites in the Gulf and acoustic observations of the natural seep bubble flares in the ocean water column. The acoustic data were from multibeam echosounders, one mounted in a forward-looking orientation on the ROV and another mounted down-looking in the haul of the ship. All of these laboratory and field data were focused on the dynamics of natural gas bubbles at temperatures and pressures favorable for clathrate hydrate formation between the gas and water. Our analyses of this data focused on understanding the mechanisms responsible for gas bubble dissolution within the hydrate stability zone (HSZ) of the oceans. Ice-like hydrate shells may form on the bubble-water interface under these conditions, and it was unknown how this might affect the mass transfer of gas into the ocean. We were able to extract bubble shrinkage rates from the HPWT datasets. Using this data we determined that mass transfer coefficients with and without a hydrate shell match empirical values for bubbles in contaminated systems (so-called dirty bubbles contaminated by naturally occurring surfactants). We also showed that free gas, and not gas hydrate, is the dominant dissolving phase when the hydrate sub-cooling is below 11 degree Celsius (temperature difference between hydrate the hydrate formation temperature and ambient temperator) or the pressure is reducing as bubbles rise through the ocean water column. Using this mass transfer model, our numerical model of bubble dissolution matched the over 200 HPWT experiments with an average error of 10% for predicting the bubble size at the end of an experiment. From field data in the literature, we also observed that gas bubbles dissolve faster when they are initially released, following mass transfer coefficients for so-called clean-bubbles (those not yet contaminated by surfactants). Shortly after release within the HSZ, a hydrate shell forms on the bubble-water interface, and the mass transfer reduces to rates matching those of dirty bubbles. We correlated this transition time from clean to dirty bubble behavior with the initial bubble surface area and the hydrate sub-cooling. With this model for hydrate formation time and using the mass transfer coefficients deduced from the HPWT data, we validated our numerical model for predicting the rise heights of natural seep flares in the oceans. Flare heights are commonly observed in haul-mounted acoustic multibeam data. The numerical model predicts bubbles to rise high in the ocean water column owing to the slower mass transfer rates for dirty bubbles that accompany the majority of their rise time. We found that the numerical model predictions matched the observed flare heights within 5% to 10% accuracy when we compared the rise heights of the largest bubbles released from the seafloor with the bubbles acoustically visible in the multibeam data. Bubbles become acoustically transparent as they shrink to sizes of order 1 mm in diameter for the multibeam frequencies used in the field. The forward-looking multibeam on the ROV also provided data on the lateral spreading of bubbles in natural seep flares. Our analysis of this data showed that spreading follows a diffusion process, with the effective diffusivity correlating with the wobbling length scale of these ellipsoidal bubbles. When we apply this diffusivity in a random displacement model of bubble spreading, our numerical simulations match closely the lateral spread observed by the M3 in the ocean water column. Finally, we compared the seep model predictions for the acoustic properties of these natural seep plumes with that observed by the acoustic instruments in the field. The M3 and EM 302 observations were converted to relative values of target strength using a calibration we obtained in the laboratory for the M3 and using an algorithm from the manufacturer for the EM 302. Comparing the numerical seep model to these data, we obtain good agreement over the whole height of rise of these bubble flares. This further validates the numerical model. Overall, our validated seep model captures the key dynamics of gas bubbles released from natural seeps in the oceans and helps to predict the fate of methane in the water column.

03 NATURAL GAS↗