Search NASA⌕ Search

SEARCH · Search NASA

Results for “radiative properties”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Representing Fine‐Scale Topographic Effects on Surface Radiation Balance in Hyper‐Resolution Land Surface Models

Land surface models are increasingly used to simulate land surface processes at hyper-spatial resolutions (e.g., ∼1 km). As model resolution increases, grid-scale topographic effects on surface radiation fluxes and their interactions between adjacent grids become more pronounced. However, current land surface models routinely neglect the fine-scale topographic effects on surface radiation balance. This study developed physically-based and computationally-efficient parameterizations (fineTOP) that explicitly resolve fine-scale topographic effects on downward shortwave and longwave radiation as well as land surface radiative properties. The newly developed parameterizations were implemented and tested in the Energy Exascale Earth System Model (E3SM) Land Model (ELM). Multi-decadal km-resolution ELM simulations over the California Sierra Nevada show that fine-scale topography significantly impacts the surface energy balance and snow processes across seasons. Slope determines the magnitude of topographic effects, while aspect controls their sign. For slopes larger than 30°, topography-induced change in annual surface temperature can be as large as 3.3 K. Regionally, the mean value and standard deviation of topography-induced changes in annual surface temperature are −0.22 ± 0.38 K and +0.25 ± 0.37 K over north-facing and south-facing slopes, respectively. Topography-induced changes in surface radiative properties account for 3.5% ± 13.8% of total topographic effects on annual net radiation. With fineTOP, ELM captures the aspect-dependence of snow cover fraction, snow water equivalent, and land surface temperature found in MODIS satellite observations and a snow reanalysis data set, while the default ELM fails to capture this phenomenon. The enhanced capability to represent fine-scale topographic effects on surface radiation balance can be used to advance understanding of the role of fine-scale topography in land surface processes and land-atmosphere interactions over mountainous regions.

Hao, Dalei [Pacific Northwest National Laboratory ↗
DOE OSTI2022

Platform for probing radiation transport properties of hydrogen at conditions found in the deep interiors of red dwarfs

We describe an experimental concept at the National Ignition Facility for specifically tailored spherical implosions to compress hydrogen to extreme densities (up to ∼800× solid density, electron number density ne∼4×1025 cm−3) at moderate temperatures (T∼200 eV), i.e., to conditions, which are relevant to the interiors of red dwarf stars. The dense plasma will be probed by laser-generated x-ray radiation of different photon energy to determine the plasma opacity due to collisional (free–free) absorption and Thomson scattering. The obtained results will benchmark radiation transport models, which in the case for free–free absorption show strong deviations at conditions relevant to red dwarfs. This very first experimental test of free–free opacity models at these extreme states will help to constrain where inside those celestial objects energy transport is dominated by radiation or convection. Moreover, our study will inform models for other important processes in dense plasmas, which are based on electron–ion collisions, e.g., stopping of swift ions or electron–ion temperature relaxation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Geometrically Tunable Beamed Light Emission from a Quantum‐Dot Ensemble Near a Gradient Metasurface

Abstract Optical metasurfaces have been widely investigated in recent years as a means to tailor the wavefronts of externally incident light for passive device applications. At the same time, their use in active optoelectronic devices such as light emitters is far less established. This work explores their ability to control the radiation properties of a nearby continuous ensemble of randomly oriented incoherent dipole sources via near‐field interactions. Specifically, a film of colloidal quantum dots is deposited on a plasmonic metasurface consisting of a 1D array of metallic nanoantennas on a metal film. The array is designed to introduce a linear phase profile upon reflection, and a bi‐periodic nanoparticle arrangement is introduced to ensure adequate sampling of the desired phase gradient. Highly directional radiation patterns are correspondingly obtained from the quantum dots at an enhanced emission rate. The underlying radiation mechanism involves the near‐field excitation of surface plasmon polaritons at the metal film, and their selective diffractive scattering by the metasurface into well‐collimated beams along predetermined geometrically tunable directions. These results underscore the distinctive ability of metasurfaces to control radiation properties directly at the source level, which is technologically significant for the continued miniaturization and large‐scale integration of optoelectronic devices.

Wang, Xiaowei↗
DOE OSTI2022

The High-Frequency Signature of Slow and Fast Laboratory Earthquakes

Tectonic faults fail through a spectrum of slip modes, ranging from slow aseismic creep to rapid slip during earthquakes. Understanding the seismic radiation emitted during these slip modes is key for advancing earthquake science and earthquake hazard assessment. In this work, we use laboratory friction experiments instrumented with ultrasonic sensors to document the seismic radiation properties of slow and fast laboratory earthquakes. Stick-slip experiments were conducted at a constant loading rate of 8 μm/s and the normal stress was systematically increased from 7 to 15 MPa. We produced a full spectrum of slip modes by modulating the loading stiffness in tandem with the fault zone normal stress. Acoustic emission data were recorded continuously at 5 MHz. We demonstrate that the full continuum of slip modes radiate measurable high-frequency energy between 100 and 500 kHz, including the slowest events that have peak fault slip rates <100 μm/s. The peak amplitude of the high-frequency time-domain signals scales systematically with fault slip velocity. Stable sliding experiments further support the connection between fault slip rate and high-frequency radiation. Experiments demonstrate that the origin of the high-frequency energy is fundamentally linked to changes in fault slip rate, shear strain, and breaking of contact junctions within the fault gouge. Our results suggest that having measurements close to the fault zone may be key for documenting seismic radiation properties and fully understanding the connection between different slip modes.

58 GEOSCIENCES↗
DOE OSTI2020

Reconciling Contrasting Relationships Between Relative Dispersion and Volume-Mean Radius of Cloud Droplet Size Distributions

Cloud droplet spectral relative dispersion is critical to parameterizations of cloud radiative properties, warm-rain initiation, and aerosol-cloud interactions in models; however, there is no consistent relationship between relative dispersion and volume-mean radius in literature, which hinders improving relative dispersion parameterization and calls for physical explanation. Here we show, by analyzing aircraft observations of cumulus clouds during Routine AAF [Atmospheric Radiation Measurement (ARM) Aerial Facility] Clouds with Low Optical Water Depths (CLOWD) Optical Radiative Observations, that the correlation between relative dispersion and volume-mean radius changes from positive to negative as volume-mean radius increases. With the new observation, we postulate that the sign of the correlation is determined by whether or not condensation (evaporation) occurs simultaneously with significant new activation (deactivation). The hypothesis is validated by simulations of both an adiabatic cloud parcel model and a parcel model accounting for entrainment-mixing. A new quantity, first bin strength, is introduced to quantify this new observation. Theoretical analysis of truncated gamma and modified gamma size distributions further supports the hypothesis and reconciles the contrasting relationships between relative dispersion and volume-mean radius, including the results in polluted fog observations. The results could shed new light on the so-called “twilight zone” between cloudy and cloud-free air, which in turn affects evaluation of aerosol-cloud interactions and retrieval of aerosol optical depth.

Lu, Chunsong↗
DOE OSTI2025

The diurnal susceptibility of subtropical clouds to aerosols

The diurnal susceptibility of clouds and their radiative properties to aerosols is examined during their Lagrangian transition from subtropical stratocumulus to shallow cumulus regimes. Using large-eddy simulations, we analyze the six-day evolution of an air mass along a 3800 km observed trajectory from the coast of Peru toward the equator. Pristine and polluted scenarios are simulated with forcing imposed from weather reanalysis. The polluted scenario exhibits stronger diurnal variations in cloud water, cloud fraction, and albedo, with enhanced nighttime entrainment and suppressed precipitation. The overall response of cloud properties and outgoing shortwave radiation to droplet number concentration follows a distinct diurnal pattern: strong positive cloud adjustments dominate at night and in the morning, while weak negative adjustments prevail in the afternoon. This cycle is driven by the competition between precipitation suppression, which enhances cloud water and coverage, and entrainment drying, which depletes them. In polluted conditions, enhanced entrainment leads to a deeper and more decoupled boundary layer that cannot be sustained by surface fluxes in the afternoon, resulting in negative cloud adjustments. The enhanced entrainment rate under polluted conditions is caused by the reduced sedimentation of cloud and precipitation water from the entrainment zone. While the Twomey effect dominates the diurnal average albedo response, the diurnal variation in the competing cloud adjustments lead to a near-neutral net susceptibility in the afternoon, highlighting the critical role of diurnally varying processes in aerosol-cloud interactions.

Geosciences↗
DOE OSTI2022

Thermal radiation analysis of a pressurized oxy-coal combustor with discrete ordinate model

Concerns over climate change have led to numerous efforts in developing low-carbon energy technologies. Pressurized oxy-combustion (POC) is a promising candidate to reduce carbon emissions in power generation. Pressure strongly impacts thermal radiation, leading to a substantial difference in heat transfer between POC and conventional atmospheric pressure combustion. A good understanding of thermal radiation in POC is needed to aid new combustor development. The present work is a step in this direction, initiating a systematic analysis of thermal radiation and heat transfer in a pilot-scale POC combustor, which has been developed at Washington University in St. Louis (WUSTL). In a POC process, pulverized coal is burned under elevated pressure and O2-CO2 environment. Few researchers have investigated the thermal radiation of pressurized oxy-combustion. This work adopted the discrete ordinate model to simulate the radiation of a 1-D cylinder with predetermined conditions (temperature, pressure, gas compositions) to predict the influence of radiative properties of flue gas on the radiation in the combustor.

Li, Lei↗
DOE OSTI2025

Contact Freezing of Water Droplets by Crystalline Organic Acids

The ability of water to freeze into ice crystals in mixed-phase clouds affects physical properties, including particle size, precipitation rates, and radiative properties. The presence of an insoluble particle at the surface of water droplets can promote ice nucleation at temperatures higher than that of pure water, even in the absence of a collision. However, contact freezing remains an underexplored mode of ice nucleation. Here, we present a study of atmospherically relevant organic acids and their role as effective ice-nucleating particles (INP) in contact mode using a Raman-microscope-equipped environmental chamber. We determined contact freezing temperatures induced by solid crystals of docosanol, adipic acid, cis-pinonic acid, fumaric acid, 4-hydroxybenzoic acid, palmitic acid, phthalic acid, sebacic acid, stearic acid, terephthalic acid, and vanillic acid. All solids except fumaric acid promoted contact freezing of water droplets at significantly higher temperatures than pure water in the chamber (−15.0 to −18.5 °C vs −21.3 °C). Physical and chemical properties were identified which correlate with greater effectiveness of INPs in the contact mode, including crystal lattice mismatch with ice, carbon number, and insolubility in water. In conclusion, we suggest that the presence of these organic solids in atmospheric aerosols may promote atmospheric ice nucleation at warm temperatures.

Clouds↗
DOE OSTI2019

Tunable High Spatio-Spectral Purity Undulator Radiation from a Transported Laser Plasma Accelerated Electron Beam

Undulator based synchrotron light sources and Free Electron Lasers (FELs) are valuable modern probes of matter with high temporal and spatial resolution. Laser Plasma Accelerators (LPAs), delivering GeV electron beams in few centimeters, are good candidates for future compact light sources. However the barriers set by the large energy spread, divergence and shot-to-shot fluctuations require a specific transport line, to shape the electron beam phase space for achieving ultrashort undulator synchrotron radiation suitable for users and even for achieving FEL amplification. Proof-of-principle LPA based undulator emission, with strong electron focusing or transport, does not yet exhibit the full specific radiation properties. We report on the generation of undulator radiation with an LPA beam based manipulation in a dedicated transport line with versatile properties. After evidencing the specific spatio-spectral signature, we tune the resonant wavelength within 200–300 nm by modification of the electron beam energy and the undulator field. We achieve a wavelength stability of 2.6%. We demonstrate that we can control the spatio-spectral purity and spectral brightness by reducing the energy range inside the chicane. We have also observed the second harmonic emission of the undulator.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Measurements and Analyses to Enable Science for the Eastern Pacific Cloud Aerosol Precipitation Experiment (EPCAPE)

Coastal cities provide the opportunity to characterize the substantial effects of manmade particles on marine cloud properties and processes. La Jolla lies to the north of San Diego, California, but it is often about a day directly downwind of the major pollution sources located in the ports of Los Angeles and Long Beach. The large dynamic range of aerosol particle concentrations combined with the multi-hour to multi-day persistence of stratocumulus cloud layers makes the site ideal for investigating the seasonal changes in cloud and aerosol properties as well as the quantitative relationships between cloud and aerosol properties. The Eastern Pacific Cloud Aerosol Precipitation Experiment (EPCAPE) characterized the extent, radiative properties, aerosol interactions, and precipitation characteristics of stratocumulus clouds in the Eastern Pacific across all four seasons at two coastal sites in La Jolla. This project was designed to enhance and expand the scientific uses of the ARM AMF1 measurements from EPCAPE. The goal was to ensure meaningful observations were collected that would enable science. The project included the following objectives: (1) Reviewing the ARM AMF1 measurements [ARM, 2021a; b] and circulating a summary of the measurements each week, (2) Comparing AMF1 measurements to those provided by collaborators (including filter measurements at the pier), (3) Collecting and analyzing filter samples from Scripps Pier by Fourier Transform Infrared spectroscopy (FTIR) and X-ray Fluorescence (XRF), (4) Assisting in operations of instruments provided by Guest PIs when possible, and (5) Providing an initial compilation of EPCAPE aerosol and cloud seasonal differences. The expected outcomes of these objectives were enhanced proposals and publications using EPCAPE measurements by helping to identify instrument issues, expanded data access and awareness by distributing weekly plots and related summaries, improved source-related attribution of aerosols with elemental tracers, additional observations provided by Guest PIs, and accelerated ACI studies enabled by the compiled seasonal summaries of aerosol and cloud properties. Three examples of the science enabled by this project are findings that (i) aerosol and cloud aqueous production contributes more than half of sulfate particle mass concentration, (ii) upwind sources make chemical composition very similar at nearby sites despite local differences in meteorology, and (iii) most of the large mass concentration of semi-volatile organic components is co-emitted and co-evaporated with nitrate. Together these findings illustrate how ARM extended field campaigns in coastal regions can be used to constrain ACI processes with direct observations. By using the unique ARM suite of cloud radiative products in addition to the measured aerosol properties at a coastal location, we were able to address the more specific question of which aerosol particles cause how much of the effects on clouds. Identifying this signature in coastal areas provides an opportunity to test the representation of aerosol sources by global models in a range of clean and urban-influenced conditions.

Russell, Lynn [Univ. of California, San Diego, CA ↗
DOE OSTI2021

Effects of turbulence-radiation coupling in pressurized oxy-combustion

Concerns over climate change have led to numerous efforts in developing low-carbon energy technologies. Pressurized oxy-combustion (POC) is a promising candidate to reduce carbon emission in power generation. Due to the strong impact of pressure on thermal radiation, heat transfer in POC differs significantly from the situation of conventional atmospheric pressure combustion. Thus thermal radiation in POC needs to be investigated to aid new combustor development. The present computational work is a step in this direction, initiating a systematic analysis of thermal radiation and heat transfer in a pilot-scale POC combustor, which has been developed at Washington University in St. Louis (WUSTL). In a POC process, pulverized coal is burned under elevated pressure and O2-CO2 environment. While most previous computational works on flame radiation focused on atmospheric pressure condition, this work considered a pressurized flame. Specifically, a 15-bar POC combustor of power 50 kW is modeled employing the Ansys FLUENT commercial platform, using both Reynolds-averaged Navier-Stokes (RANS) modeling and large-eddy simulation (LES). A recently published and validated global radiation model is used to predict the radiative property of the flue gas. The discrete ordinates (DO) radiation model is chosen to solve the radiative transport equation. Incident radiation on the walls of the combustor is identified and investigated. It is revealed that for this pilot-scale, pressurized combustor, thermal radiation exhibits similar patterns in the RANS and LES models. However, the magnitudes of radiation are different in both models, presumably because interactions between radiation and turbulence are embedded in the LES model but not in the RANS model. This difference thereby underlines the major goal of this research: investigating the impact of such turbulence-to-radiation coupling on thermal radiation in a pilot-scale POC setting. Most previous studies on turbulence-to-radiation interactions focused on bench-scale, atmospheric pressure flames. This work is the first effort to extend this study to a pilot-scale, pressurized flame.

Li, Lei↗
DOE OSTI2023

Relationships among diurnal variations of polar night cloud, precipitation, surface temperatures, and the fair-weather return current of the global electric circuit (GEC)

Over the past five years (June 2017-current), the vertical electric field (E z ) as well as numerous cloud, precipitation and radiation properties have been monitored at the Department of Energy-Atmospheric Radiation Measurement (DOE ARM) North Slope of Alaska (NSA) field site. Comparisons between the composite diurnal averaged fair-weather E z , and composite ceilometer derived cloud base height during the polar night, reveal a significant correlation between the parameters (r = 0.62), supporting previous studies that there is high correlation between local electric field and cloud properties. With the use of extensive instrumentation at the site, such as the Micro Pulse Lidar (MPL), Ka-band Zenith Radar (KAZR), ceilometer, SKYRAD, Precipitation Imaging Package (PIP), among others, this study provides a more comprehensive examination of the diurnal cycle of cloud and precipitation properties along with the localized fair-weather return current of the larger Global Electric Circuit (GEC) system. Comparisons between the composite diurnal averaged fair-weather E z , and cloud thickness, maximum column backscatter, and precipitation particle counts all show similar diurnal variability during the polar night, indicating that during the largest magnitude fair weather E z time periods, clouds bases tend to be higher, clouds are thicker, have a larger column backscatter, and display more precipitating particles. Furthermore, a slight diurnal variability in the polar night surface temperature was found to be highly correlated (r = 0.87) to the longwave downwelling irradiance measured by SKYRAD, indicating that the variations in the physical properties of local clouds could modulate the diurnal polar night surface temperature variability on the order of 0.5 °C/day. In conclusion, these findings emphasize the importance and global nature of the GEC system, with the global aggregate of thunderstorms and electrified clouds potentially influencing polar night cloud properties as well as diurnal wintertime polar surface temperature variation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Structural and optical properties of Sm{sup 3+}-doped borate glasses for luminescent applications

In this work, most interesting optical investigations are carried out with the glass composition (20 − x) CdCO{sub 3} + 15 Li{sub 2}CO{sub 3} + 15 Na{sub 2}CO{sub 3} + 50 H{sub 3}BO{sub 3} + xSm{sub 2}O{sub 3} (where x = 0.1, 0.5, 1.0, 3.0, 4.0 and 5.0 mol%) prepared with conservative melt-quenching technique. DSC verifies the annealing temperature and thermal stability of the prepared glass. X-ray diffraction technique is employed to inspect the amorphous structure of the glasses. FTIR technique helps to view the vibrational band positions of the host material. For optical studies, optical absorption and emission measurements are used. From the absorption spectra, Judd–Ofelt (JO) parameters (Ω{sub 2}, Ω{sub 4}, Ω{sub 6}), rms value of oscillator strength (δ{sub rms}), nephelauxetic parameter ( β-bar ), bonding parameter (δ) and optical band gaps are evaluated. Radiative properties such as radiative transition probabilities (A{sub T}), branching ratios (β{sub cal}, β{sub exp}) and stimulated emission cross section (σ{sub emi}) are calculated using JO theory. The higher stimulated emission cross section (14.39 × 10{sup −22} cm{sup 2}), optical gain (36.69 × 10{sup −25} cm{sup 2} s) indicate that 1.0 mol% of Sm{sup 3+}-doped borate glass would be a potential gain medium for solid state laser application at 603 nm. Color coordinates of the glasses located in the reddish-orange region could be useful for white light-emitting application.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

A Machine Learning Framework for Predicting Microphysical Properties of Ice Crystals From Cloud Particle Imagery

The microphysical properties of ice crystals are important because they significantly alter the radiative properties and spatiotemporal distributions of clouds, which in turn strongly affect Earth's climate. However, it is challenging to measure key properties of ice crystals, such as mass or morphological features. Here, we present a proof-of-concept framework for predicting three-dimensional (3D) microphysical properties of ice crystals from in situ two-dimensional (2D) imagery. First, we computationally generated synthetic ice crystals using 3D modeling software along with geometric parameters estimated from the 2021 Ice Cryo-Encapsulation Balloon (ICEBall) field campaign. Then, we used synthetic crystals to train machine learning (ML) models to predict effective density ($ρ_e$), effective surface area ($A_e$), and number of bullets ($N_b$) from synthetic rosette imagery. On unseen synthetic images, our ML models accurately predicted ice crystal properties. ResNet-18 performed best, achieving $R^2$ values of 0.99 and 0.98 for $ρ_e$ and $A_e$, respectively, and MAE of 0.10 for mathematical equation in single view tasks. Stereo view ResNet-18 further reduced RMSE by 40% for $ρ_e$ and $A_e$ and reduced MAE by 0.08 for $N_b$. This work provides a novel ML-driven framework for estimating ice microphysical properties from in situ imagery, which will allow for downstream constraints on microphysical parameterizations, such as the mass-size relationship.

Ko, J. [Columbia Univ., New York, NY (United State↗
DOE OSTI2023

Observations of Fog-Aerosol Interactions Over Central Greenland

Supercooled fogs can have an important radiative impact at the surface of the Greenland Ice Sheet, but they are difficult to detect and our understanding of the factors that control their lifetime and radiative properties is limited by a lack of observations. This study demonstrates that spectrally resolved measurements of downwelling longwave radiation can be used to generate retrievals of fog microphysical properties (phase and particle effective radius) when the fog visible optical depth is greater than ~0.25. For 12 cases of fog under otherwise clear skies between June and September 2019 at Summit Station in central Greenland, nine cases were mixed-phase. The mean ice particle (optically-equivalent sphere) effective radius was 24.0 ± 7.8 µm, and the mean liquid droplet effective radius was 14.0 ± 2.7 µm. These results, combined with measurements of aerosol particle number concentrations, provide evidence supporting the hypotheses that (a) low surface aerosol particle number concentrations can limit fog liquid water path, (b) fog can act to increase near-surface aerosol particle number concentrations through enhanced mixing, and (c) multiple fog events in quiescent periods gradually deplete near-surface aerosol particle number concentrations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Blowing Snow at McMurdo Station, Antarctica During the AWARE Field Campaign: Multi-Instrument Observations of Blowing Snow

In polar regions, blowing snow (BLSN) can play a substantial role in regional thermodynamics and radiation properties. The extent of the impact depends upon the depth of the BLSN layer, a property that is difficult to measure and not well understood. The West Antarctic Radiation Experiment (AWARE) Field Campaign saw the deployment of a large suite of observational instruments to McMurdo Station, Antarctica, allowing for in-depth investigation of BLSN. During the year of the campaign, BLSN occurred ~7.4% of the time. In this study, additional remote sensing observations are used to supplement an existing ceilometer-based BLSN characterization algorithm creating a multi-instrument (MI) depth estimation algorithm to yield layer depths with more certainty. Backscatter coefficient and linear depolarization ratios from the micropulse and high spectral resolution lidar (HSRL) are incorporated along with color ratio derived from HSRL backscatter and Ka-band radar reflectivity observations. This provided estimates of BLSN depth that were less impacted by noise/artifacts in any one set of observations. When applied to the data from AWARE, the MI algorithm yielded an average depth of 168.2 m with ~60% of plumes less than 200 m in depth. Events occurring with precipitation saw the most value in the additional instruments, with the suite of observations allowing for more separation of BLSN particles from fallstreaks.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Nuclear Structure and Decay Data for A=35 Isobars

Here, this work presents a comprehensive and critical evaluation of experimental nuclear spectroscopic data from reactions and decays for all 11 known nuclides with mass number 35 (Ne, Na, Mg, Al, Si, P, S, Cl, Ar, K, Ca). Recommended values are produced for level energies, spins and parities, half-lives, and radiation properties including energies, branching ratios, and multipolarities of γ rays, as well as characteristics of β radiation decays, based on a rigorous assessment of all available experimental data. Discrepancies among existing results are carefully addressed. This work supersedes earlier full evaluations of A=35 published by 2011Ch48, 1990En08 (also 1998En04 update) and 1978En02.

Sun, Lijie [Michigan State University, East Lansin↗