Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “peak finding”

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.

At least 217 records · Page 12

Multiple Peaks and a Long Precursor in the Type IIn Supernova 2021qqp: An Energetic Explosion in a Complex Circumstellar Environment

Abstract We present optical photometry and spectroscopy of the Type IIn supernova (SN) 2021qqp. Its unusual light curve is marked by a long precursor for ≈300 days, a rapid increase in brightness for ≈60 days, and then a sharp increase of ≈1.6 mag in only a few days to a first peak ofM r ≈ −19.5 mag. The light curve then declines rapidly until it rebrightens to a second distinct peak ofM r ≈ −17.3 mag centered at ≈335 days after the first peak. The spectra are dominated by Balmer lines with a complex morphology, including a narrow component with a width of ≈1300 km s −1 (first peak) and ≈2500 km s −1 (second peak) that we associate with the circumstellar medium (CSM) and a P Cygni component with an absorption velocity of ≈8500 km s −1 (first peak) and ≈5600 km s −1 (second peak) that we associate with the SN–CSM interaction shell. Using the luminosity and velocity evolution, we construct a flexible analytical model, finding two significant mass-loss episodes with peak mass loss rates of ≈10 and ≈5M ⊙ yr −1 about 0.8 and 2 yr before explosion, respectively, with a total CSM mass of ≈2–4M ⊙ . We show that the most recent mass-loss episode could explain the precursor for the year preceding the explosion. The SN ejecta mass is constrained to be ≈5–30M ⊙ for an explosion energy of ≈(3–10) × 10 51 erg. We discuss eruptive massive stars (luminous blue variable, pulsational pair instability) and an extreme stellar merger with a compact object as possible progenitor channels.

Astronomy & Astrophysics↗

Quasi-Periodic Variability in NGC 5408 X-1

We report the discovery with XMM-Newton of quasiperiodic variability in the 0.2 - 10 keV X-ray flux from the ultraluminous X-ray source NGC 5408 X-1. The average power spectrum of all EPIC-pn data reveals a strong 20 mHz QPO with an average amplitude (rms) of 9%, and a coherence, Q identical with nu(sub 0)/sigma approximately equal to 6. In a 33 ksec time interval when the 20 mHz QPO is strongest we also find evidence for a 2nd QPO peak at 15 mHz, the first indication for a close pair of QPOs in a ULX source. Interestingly, the frequency ratio of this QPO pair is inconsistent with 3:2 at the 3 sigma level, but is consistent with a 4:3 ratio. A powerlaw noise component with slope near 1.5 is also present below 0.1 Hz with evidence for a break to a flatter slope at about 3 mHz. The source shows substantial broadband variability, with a total amplitude (rms) of about 30% in the 0.1 - 100 mHz frequency band, and there is strong energy dependence to the variability. The power spectrum of hard X-ray photons (greater than 2 keV) shows a "classic" flat-topped continuum breaking to a power law with index 1.5 - 2. Both the break and 20 mHz QPO are detected in the hard band, and the 20 mHz QPO is essentially at the break. The QPO is both strong and narrow in this band, having an amplitude (rms) of 15%, and Q approx. equal to 25. The energy spectrum is well fit by three components, a "cool" disk with kT = 0.15 keV, a steep power law with index 2.56, and a thermal plasma at kT = 0.87 keV. The disk, power law, and thermal plasma components contribute 35, 60, and 5% of the 0.3 - 10 keV flux, respectively. Both the timing and spectral properties of NGC 5408 X-1 are strikingly reminiscent of Galactic black hole systems at high inferred accretion rates, but with its characteristic frequencies (QPO and break frequencies) scaled down by a factor of 10 - 100. We discuss the implications of these findings in the context of models for ULXs, and their implications for the object's mass.

Strohmayer, Tod E.↗

The Activity Profiles and Peak Flux of Radar Meteor Showers

Meteor showers are transient phenomena; a limited duration is what distinguishes a shower from a sporadic source. They do not, however, have clear start and end times; instead, shower activity gradually (or rapidly) increases over time, reaches a peak, and then decays. In this report, we examine the activity profiles of 38 meteor showers using 20 years of single-station shower fluxes from the Canadian Meteor Orbit Radar (CMOR). We find that a number of showers, such as the Southern delta Aquariids (SDAs), exhibit a broad maximum in activity rather than a sharp peak. We also find that approximately one-third of showers exhibit asymmetric activity profiles in which the rate in which activity rises prior to the peak differs from the rate at which it decays after the peak. We provide a functional form for meteor shower activity profiles that incorporates these features and use these profile fits to estimate each shower’s peak flux. This report also includes an uncertainty analysis for CMOR fluxes. This includes the first (to our knowledge) characterization of systematic uncertainty in radar meteor fluxes.

A V Moorhead↗

What can the GRACE seasonal cycle tell us about lake-aquifer interactions?

Lake-aquifer interactions have been the subject of investigation and debate for decades. Traditional investigation methods include measurement of water flux across the groundwater-surface water interface, application of heat and environmental tracer methods, conducting numerical simulations of the water flow, and mass balance-based approaches. We first review, evaluate and describe the different approaches that have been applied for examining lake-aquifer interactions and propose an additional complementary approach. While the leakage from lakes and its slow movement through porous media (laminar flow) is well characterized by existing methods, modeling the rapid and turbulent flow through preferred pathways (faults/karst) remains a challenging task. To better understand the nature, and the full scale, of the lake aquifer interactions, and using Lake Nasser and the underlying Dakhla subbasin of the Nubian Sandstone Aquifer System (NSAS) in the Western Desert of Egypt (area: 0.66 x 10 6 km 2 ) as a test site, we examine, model, and correlate temporally and spatially, the variations in Gravity Recovery and Climate Experiment terrestrial water storage (GRACE TWS ) solutions, with precipitation, lake levels, area, and water volume. We review current GRACE applications in hydrology and present our novel approach. Findings include: firstly, large seasonal fluctuations (peak: Nov./Dec.; trough: July/Aug.) in Lake Nasser's surface water levels (average annual fluctuations: from 3.7 m to 7.5 m), area (average area: 3622 km 2 but up to 4530 km 2 ), and volume (annual average: 13.4 km 3 , and up to 34 km 3 ) are observed throughout years 2006 to 2015. These fluctuations are accompanied by an increase in GRACE TWS (average: 50 ± 13 mm/yr, up to 77 ± 18 mm/yr) over Lake Nasser and by a progression of a front of increasing GRACE TWS values (> 50 ± 13 mm) away from the lake reaching distances of up to 700 km some 3 to 5 months following peak lake level periods. The areas witnessing the seasonal increase in GRACE TWS display a progressive increase in phase and decrease in amplitude with distance from the lake. Secondly, the negligible precipitation over the Dakhla subbasin cannot account for the observed seasonal GRACE TWS patterns and neither can the leakage signal from Lake Nasser. Thirdly, overall similarities in the distribution of GRACE TWS seasonal spatial patterns are observed. During high lake level periods (e.g., Dec. 2007; Feb. 2015) additional, and more distant, areas from Lake Nasser saw increases in GRACE TWS and vice versa during low lake level periods (e.g., July 2006; Aug. 2010). These observations are consistent with Lake Nasser being the main source of modern recharge for the Dakhla subbasin and suggest a new conceptual model for the subbasin: a slow groundwater flow through a porous matrix and a rapid groundwater flow along a network of faults, fractures, and karst topography across the matrix. We suggest that the proposed conceptual model for the interaction between Lake Nasser and the Dakhla subbasin could be applicable to aquifers of similar geologic, climatic, and hydrologic settings worldwide and that approaches similar to those advocated here could be used to investigate the validity of this suggestion.

54 ENVIRONMENTAL SCIENCES↗

Investigating the Electronic Structure of Prospective Water-Splitting Oxide BaCe 0.25 Mn 0.75 O 3-δ before and after Thermal Reduction

BaCe 0.25 Mn 0.75 O 3-δ (BCM), a non-stoichiometric oxide with a layered perovskite-like crystal structure, has recently emerged as a prospective contender for application in renewable energy harvesting by solar thermochemical hydrogen generation. Using solar-thermal energy and a reducing environment, oxygen vacancies can be created in high-temperature BCM, and the reduced crystal so obtained can, in turn, produce H 2 by stripping oxygen from H 2 O. Therefore, a first step toward understanding the working mechanism and optimizing the performance of BCM is a thorough and comparative analysis of the electronic structure of the pristine and the reduced material. In this paper, we probe the electronic structure of BCM using the combined effort of first-principles calculations and experimental O K-edge X-ray absorption spectroscopy (XAS). The computed projected density of states (PDOS) and orbital plots are used to propose a simplified model for orbital mixing between the oxygen and metal atoms. With the help of state-of-the-art simulations, we are able to find the origins of the XAS peaks and categorize them on the basis of contribution from Ce and Mn. For the reduced crystal, the calculations show that the change in electron density resulting from the reduction is strongly localized around the oxygen vacancy. Experimental measurements reveal a marked lowering of the first O K-edge peak in the reduced crystal. Using theoretical analysis, this is shown to result from lifting of spin degeneracy in the absorption peaks as well as from a diminished O 2p contribution to the frontier unoccupied orbitals, in accordance with the tight binding scheme. The simulated results serve as a reference for the extent of spectral change as a function of the percentage of oxygen vacancies in the reduced crystal. Here, our study paves the way for the investigation of the working mechanism of BCM and for computational and experimental efforts aimed at design and discovery of efficient water-splitting oxides.

08 HYDROGEN↗

Transient Mesoscopic Immiscibility, Viscosity Anomaly, and High Internal Pressure at the Semiconductor–Metal Transition in Liquid Ga 2 Te 3

Gallium tellurides appear to be promising phase-change materials (PCMs) of the next generation for brain-inspired computing and reconfigurable optical metasurfaces. They are different from the benchmark PCMs because of sp 3 gallium hybridization in both cubic Ga 2 Te 3 and amorphous pulsed laser deposition (PLD) films. Liquid Ga 2 Te 3 also shows a viscosity η(T) anomaly just above melting when η(T) first increases and only then starts decreasing. We used high-energy X-ray diffraction to observe a transient mesoscopic immiscibility that suggested dense metallic liquid droplets in a semiconducting melt. The η(T) shape was consistent with this finding. A vanishing first sharp diffraction peak that also shifts to a higher Q indicates a high internal pressure in the metallic melt, which produces a remarkable asymmetry of the Ga–Te nearest neighbor distances and is reminiscent of high-pressure rhombohedral Ga 2 Te 3 . Finally, the observed phenomena provide a realistic scenario for a fast, multilevel SET–RESET response, which also unravels similar trends in the purported density-driven liquid polyamorphism of water, phosphorus, sulfur, and other materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Theory of Photoluminescence Spectral Line Shapes of Semiconductor Nanocrystals

Single-molecule photoluminescence (PL) spectroscopy of semiconductor nanocrystals (NCs) reveals the nature of exciton–phonon interactions in NCs. Understanding the homogeneous spectral line shapes and their temperature dependence remains an open problem. Here, we develop an atomistic model to describe the PL spectrum of NCs, accounting for excitonic effects, phonon dispersion relations, and exciton–phonon couplings. We validate our model using single-NC measurements on CdSe/CdS NCs from T = 4 to 290 K, and we find that the slightly asymmetric main peak at low temperatures is comprised of a narrow zero-phonon line (ZPL) and acoustic phonon sidebands. Furthermore, we identify the specific phonon modes that give rise to the optical phonon sidebands. At temperatures above 200 K, the spectral line width shows a stronger dependence upon the temperature, which we demonstrate to be correlated with higher order exciton–phonon couplings. We also identify the line width dependence upon reorganization energy, NC core sizes, and shell thicknesses.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

How Cooperatively Folding Are Homopolymer Molecular Knots?

Detailed thermodynamic analysis of complex systems with multiple stable configurational states allows for insight into the cooperativity of each individual transition. In this work, we derive a heat capacity decomposition comprising contributions from each individual configurational state, which together sum to a baseline heat capacity, and contributions from each state-to-state transition. We apply this analysis framework to a series of replica exchange molecular dynamics simulations of linear and 1-1 coarse-grained homo-oligomer models, which fold into stable, configurationally well-defined molecular knots, in order to better understand the parameters leading to stable and cooperative folding of these knots. We find that a stiff harmonic backbone bending angle potential is key to achieving knots with specific 3D structures. Tuning the backbone equilibrium angle in small increments yields a variety of knot topologies. Populations of different knotted states as a function of temperature can also be manipulated by tuning the backbone torsion stiffness or by adding side-chain beads. We find that sharp total heat capacity peaks for the homo-oligomer knots are largely due to a coil-to-globule transition rather than a cooperative knotting step. Furthermore, in some cases, the cooperativity of globule-to-knot and coil-to-globule transitions is comparable, suggesting that highly cooperative folding to knotted structures can be achieved by refining the model parameters or adding sequence specificity.

36 MATERIALS SCIENCE↗

Thermal decomposition of neptunyl ammonium nitrate: mechanistic insights and structural characterization of the Np 2 O 5 intermediate phase

Neptunium (Np) possesses a rich and unique chemistry that often diverges from other actinide elements yet remains relatively underexplored compared with the other light actinides. A resurgence of interest in Np has been spurred by the application of 237 Np for plutonium-238 ( 238 Pu) production for use in radioisotope thermoelectric generators (RTGs), necessitating evaluation of Np chemical reactions and materials. The work presented here studied the thermal decomposition of neptunyl ammonium nitrate (NH 4 Np VI O 2 (NO 3 ) 3 ) for synthesis of neptunium dioxide (NpO 2 ), which is the target material used for production of 238 Pu. Additionally, structural characterization of the intermediate solid Np pentoxide (Np 2 O 5 ) was performed. Advanced solid-state characterization techniques, including simultaneous thermal analysis (STA), powder X-ray diffraction (pXRD), Raman spectroscopy, and density functional theory (DFT) modeling have been combined to study the reaction pathways. Analysis revealed that NH 4 Np VI O 2 (NO 3 ) 3 thermally decomposes to a proposed neptunyl nitrate intermediate, followed by Np 2 O 5 and finally NpO 2 , all within the temperature range of 150 °C–600 °C. Further characterization of the pentoxide intermediate provided the first Raman spectra of pure-phase Np 2 O 5 and associated DFT modeling confirmed Raman peak assignments for this phase. These findings provide mechanistic information to advance production of the critical radioisotope 238 Pu and advance the state of knowledge on Np materials chemistry using modern characterization techniques.

Lawson, Kathryn M. [Oak Ridge National Laboratory ↗

Evidence for sub-Chandrasekhar Type Ia supernovae from the last major merger

ABSTRACT We investigate the contribution of sub-Chandrasekhar mass Type Ia supernovae to the chemical enrichment of the Gaia Sausage galaxy, the progenitor of a significant merger event in the early life of the Milky Way. Using a combination of data from Nissen & Schuster, the GALactic Archaeology with HERMES (GALAH) Data Release 3 [with 1D non-local thermal equilibrium (NLTE) abundance corrections], and the Apache Point Observatory Galactic Evolution Experiment (APOGEE) Data Release 16, we fit analytic chemical evolution models to a nine-dimensional chemical abundance space (Fe, Mg, Si, Ca, Cr, Mn, Ni, Cu, and Zn) in particular focusing on the iron-peak elements, Mn and Ni. We find that low [Mn/Fe] $\sim -0.15\, \mathrm{dex}$ and low [Ni/Fe] $\sim -0.3\, \mathrm{dex}$ Type Ia yields are required to explain the observed trends beyond the [α/Fe] knee of the Gaia Sausage (approximately at [Fe/H] $=-1.4\, \mathrm{dex}$). Comparison to theoretical yield calculations indicates a significant contribution from sub-Chandrasekhar mass Type Ia supernovae in this system (from ${\sim} 60\, \mathrm{per\, cent}$ to $100\, \mathrm{per\, cent}$ depending on the theoretical model with an additional ${\pm} 10\, \mathrm{per\, cent}$ systematic from NLTE corrections). We compare to results from other Local Group environments including dwarf spheroidal galaxies, the Magellanic Clouds, and the Milky Way’s bulge, finding the Type Ia [Mn/Fe] yield must be metallicity dependent. Our results suggest that sub-Chandrasekhar mass channels are a significant, perhaps even dominant, contribution to Type Ia supernovae in metal-poor systems, whilst more metal-rich systems could be explained by metallicity-dependent sub-Chandrasekhar mass yields, possibly with additional progenitor mass variation related to star formation history, or an increased contribution from Chandrasekhar mass channels at higher metallicity.

Sanders, Jason L.↗

Tunable Layer Circular Photogalvanic Effect in Twisted Bilayers

In this work, we develop a general theory of the layer circular photogalvanic effect (LCPGE) in quasi-two-dimensional chiral bilayers, which refers to the appearance of a polarization-dependent, out-of-plane static dipole moment induced by circularly polarized light. We elucidate the geometric origin of the LCPGE as two types of interlayer coordinate shift weighted by the quantum metric tensor and the Berry curvature, respectively. As a concrete example, we calculate the LCPGE in twisted bilayer graphene, and find that it exhibits a resonance peak whose frequency can be tuned from visible to infrared as the twisting angle varies. The LCPGE thus provides a promising route toward frequency-sensitive, circularly polarized light detection, particularly in the infrared range.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Optimizing Simulation Parameters for Weak Lensing Analyses Involving Non-Gaussian Observables

We performed a series of numerical experiments to quantify the sensitivity of the predictions for weak lensing statistics obtained in ray-tracing dark matter (DM)-only simulations, to two hyper-parameters that influence the accuracy as well as the computational cost of the predictions: the thickness of the lens planes used to build past light cones and the mass resolution of the underlying DM simulation. The statistics considered are the power spectrum (PS) and a series of non-Gaussian observables, including the one-point probability density function, lensing peaks, and Minkowski functionals. Counterintuitively, we find that using thin lens planes (< 60 h {sup −1} Mpc on a 240 h {sup −1} Mpc simulation box) suppresses the PS over a broad range of scales beyond what would be acceptable for a survey comparable to the Large Synoptic Survey Telescope (LSST). A mass resolution of 7.2 × 10{sup 11} h {sup −1} M {sub ⊙} per DM particle (or 256{sup 3} particles in a (240 h {sup −1} Mpc){sup 3} box) is sufficient to extract information using the PS and non-Gaussian statistics from weak lensing data at angular scales down to 1′ with LSST-like levels of shape noise.

79 ASTRONOMY AND ASTROPHYSICS↗

The Most Metal-poor Stars in the Inner Bulge

The bulge is the oldest component of the Milky Way. Since numerous simulations of Milky Way formation have predicted that the oldest stars at a given metallicity are found on tightly bound orbits, the Galaxy’s oldest stars are likely metal-poor stars in the inner bulge with small apocenters (i.e., R {sub apo} ≲ 4 kpc). In the past, stars with these properties have been impossible to find due to extreme reddening and extinction along the line of sight to the inner bulge. We have used the mid-infrared metal-poor star selection of Schlaufman & Casey (2014) on Spitzer/Galactic Legacy Infrared Mid-Plane Survey Extraordinaire data to overcome these problems and target candidate inner bulge metal-poor giants for moderate-resolution spectroscopy with Anglo-Australian Telescope/AAOmega. We used those data to select three confirmed metal-poor giants ([Fe/H] = −3.15, −2.56, −2.03) for follow-up high-resolution Magellan/Magellan Inamori Kyocera Echelle spectroscopy. A comprehensive orbit analysis using Gaia DR2 astrometry and our measured radial velocities confirms that these stars are tightly bound inner bulge stars. We determine the elemental abundances of each star and find high titanium and iron-peak abundances relative to iron in our most metal-poor star. We propose that the distinct abundance signature we detect is a product of nucleosynthesis in the Chandrasekhar-mass thermonuclear supernova of a CO white dwarf accreting from a helium star with a delay time of about 10 Myr. Even though chemical evolution is expected to occur quickly in the bulge, the intense star formation in the core of the nascent Milky Way was apparently able to produce at least one Chandrasekhar-mass thermonuclear supernova progenitor before chemical evolution advanced beyond [Fe/H]∼−3.

79 ASTRONOMY AND ASTROPHYSICS↗

Simulating Continuum-based Redshift Measurement in the Roman’s High Latitude Spectroscopic Survey

We investigate the capability of the Nancy Grace Roman Space Telescope’s (Roman) Wide-Field Instrument G150 slitless grism to detect red, quiescent galaxies based on the current reference survey. We simulate dispersed images for Roman reference High-Latitude Spectroscopic Survey (HLSS) and analyze two-dimensional spectroscopic data using the grism Redshift and Line Analysis (Grizli) software. This study focus on assessing Roman grism’s capability for continuum-level redshift measurement for a redshift range of 0.5 ≤ z ≤ 2.5. The redshift recovery is assessed by setting three requirements of: σ z = $\frac{|z–z_{true}|}{1+z}$ ≤ 0.01, signal-to-noise ratio≥ 5 and the presence of a single dominant peak in redshift likelihood function. We find that, for quiescent galaxies, the reference HLSS can reach a redshift recovery completeness of ≥50% for F158 magnitude brighter than 20.2 mag. We also explore how different survey parameters, such as exposure time and the number of exposures, influence the accuracy and completeness of redshift recovery, providing insights that could optimize future survey strategies and enhance the scientific yield of the Roman in cosmological research.

Astronomical simulations↗

Monitoring the Morphology of M87* in 2009-2017 with the Event Horizon Telescope

The Event Horizon Telescope (EHT) has recently delivered the first resolved images of M87*, the supermassive black hole in the center of the M87 galaxy. These images were produced using 230 GHz observations performed in 2017 April. Additional observations are required to investigate the persistence of the primary image feature—a ring with azimuthal brightness asymmetry—and to quantify the image variability on event horizon scales. To address this need, we analyze M87* data collected with prototype EHT arrays in 2009, 2011, 2012, and 2013. While these observations do not contain enough information to produce images, they are sufficient to constrain simple geometric models. We develop a modeling approach based on the framework utilized for the 2017 EHT data analysis and validate our procedures using synthetic data. Applying the same approach to the observational data sets, we find the M87* morphology in 2009–2017 to be consistent with a persistent asymmetric ring of ∼40 μas diameter. The position angle of the peak intensity varies in time. In particular, we find a significant difference between the position angle measured in 2013 and 2017. These variations are in broad agreement with predictions of a subset of general relativistic magnetohydrodynamic simulations. We show that quantifying the variability across multiple observational epochs has the potential to constrain the physical properties of the source, such as the accretion state or the black hole spin.

79 ASTRONOMY AND ASTROPHYSICS↗

Anisotropic Energy Transfer and Conversion in Magnetized Compressible Turbulence

We present a spatial filtering (or coarse-graining) analysis on 3D magnetized magnetohydrodynamic (MHD) turbulence simulations. The filtered compressible MHD formulae show transfer of kinetic and magnetic energies from large to small scales, as well as energy conversion between kinetic, magnetic, and thermal energies. The anisotropic filtering enables separate analyses of the energy flows perpendicular and parallel to the global mean magnetic field. Anisotropy in energy cascade is demonstrated by the larger perpendicular energy cascade rate and also the larger perpendicular wavenumbers associated with the peak energy transfer rate. We also find that the “inertial range” along the parallel (perpendicular) direction in the anisotropic energy cascade formulation is no longer strictly dissipation-free, because it includes the dissipation in the perpendicular (parallel) direction. A change in the driving force (kinetic only versus kinetic and magnetic) affects the energy conversion between kinetic and magnetic energies. While the compressibility of the driving force changes the partition of different channels of energy transfer and conversion, and also increases the total energy transfer rate, the global energy flow remains unaffected by compressibility qualitatively. Our analysis can be applied to multispacecraft observations of turbulence in the solar wind or a planetary magnetosphere.

79 ASTRONOMY AND ASTROPHYSICS↗

What Sets the Star Formation Rate of Molecular Clouds? The Density Distribution as a Fingerprint of Compression and Expansion Rates

We use a suite of 3D simulations of star-forming molecular clouds, with and without stellar feedback, magnetic fields, and driven turbulence, to study the compression and expansion rates of the gas as functions of density. We show that, around the mean density, supersonic turbulence promotes rough equilibrium between the amounts of compressing and expanding gas, consistent with continuous gas cycling between high- and low-density states. We find that the inclusion of protostellar jets produces rapidly expanding and compressing low-density gas. We find that the gas mass flux peaks at the transition between the lognormal and power-law forms of the density probability distribution function (PDF). This is consistent with the transition density tracking the post-shock density, which promotes an enhancement of mass at this density (i.e., shock compression and filament formation). At high densities, the gas dynamics are dominated by self-gravity: the compression rate in all of our runs matches the rate of the run with only gravity, suggesting that processes other than self-gravity have little effect at these densities. The net gas mass flux becomes constant at a density below the sink formation threshold, where it equals the star formation rate. The density at which the net gas mass flux equals the star formation rate is one order of magnitude lower than our sink threshold density, corresponds to the formation of the second power-law tail in the density PDF, and sets the overall star formation rates of these simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

Recovery of motion parameters from distortions in scanned images

Scanned images, such as those produced by the scanning-laser ophthalmoscope (SLO), show distortions when there is target motion. This is because pixels corresponding to different image regions are acquired sequentially, and so, in essence, are slices of different snapshots. While these distortions create problems for image registration algorithms, they are potentially useful for recovering target motion parameters at temporal frequencies above the frame rate. Stetter, Sendtner and Timberlake measured large distortions in SLO images to recover the time course of rapid horizontal saccadic eye movements. Here, this work is extended with the goal of automatically recovering small eye movements in two dimensions. Eye position during the frame interval is modeled using a low dimensional parametric description, which in turn is used to generate predicted distortions of a reference template. The input image is then registered to the distorted template using normalized cross correlation. The motion parameters are then varied, and the correlation recomputed, to find the motion which maximizes the peak value of the correlation. The location and value of the correlation maximum are determined with sub-pixel precision using biquadratic interpolation, yielding eye position resolution better than 1 arc minute. This method of motion parameter estimation is tested using actual SLO images as well as simulated images. Motion parameter estimation might also be applied to individual video lines in order to reduce pipeline delays for a near real-time system.

Mulligan, Jeffrey B.↗