Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dl”

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 307 records · Page 17

IoT Intrusion Detection Taxonomy, Reference Architecture, and Analyses

This paper surveys the deep learning (DL) approaches for intrusion-detection systems (IDSs) in Internet of Things (IoT) and the associated datasets toward identifying gaps, weaknesses, and a neutral reference architecture. A comparative study of IDSs is provided, with a review of anomaly-based IDSs on DL approaches, which include supervised, unsupervised, and hybrid methods. All techniques in these three categories have essentially been used in IoT environments. To date, only a few have been used in the anomaly-based IDS for IoT. For each of these anomaly-based IDSs, the implementation of the four categories of feature(s) extraction, classification, prediction, and regression were evaluated. We studied important performance metrics and benchmark detection rates, including the requisite efficiency of the various methods. Four machine learning algorithms were evaluated for classification purposes: Logistic Regression (LR), Support Vector Machine (SVM), Decision Tree (DT), and an Artificial Neural Network (ANN). Therefore, we compared each via the Receiver Operating Characteristic (ROC) curve. The study model exhibits promising outcomes for all classes of attacks. The scope of our analysis examines attacks targeting the IoT ecosystem using empirically based, simulation-generated datasets (namely the Bot-IoT and the IoTID20 datasets).

97 MATHEMATICS AND COMPUTING↗

A Deep Learning Modeling Framework to Capture Mixing Patterns in Reactive-Transport Systems

Prediction and control of chemical mixing are vital for many scientific areas such as subsurface reactive transport, climate modeling, combustion, epidemiology, and pharmacology. Due to the complex nature of mixing in heterogeneous and anisotropic media, the mathematical models related to this phenomenon are not analytically tractable. Numerical simulations often provide a viable route to predict chemical mixing accurately. However, contemporary modeling approaches for mixing cannot utilize available spatial-temporal data to improve the accuracy of the future prediction and can be compute-intensive, especially when the spatial domain is large and for long-term temporal predictions. To address this knowledge gap, in this work we will present in this paper a deep learning (DL) modeling framework applied to predict the progress of chemical mixing under fast bimolecular reactions. This framework uses convolutional neural networks (CNN) for capturing spatial patterns and long short-term memory (LSTM) networks for forecasting temporal variations in mixing. By careful design of the framework—placement of non-negative constraint on the weights of the CNN and the selection of activation function, the framework ensures non-negativity of the chemical species at all spatial points and for all times. Our DL-based framework is fast, accurate, and requires minimal data for training. The time needed to obtain a forecast using the model is a fraction (≈ O(-6)) of the time needed to obtain the result using a high-fidelity simulation. To achieve an error of 10% (measured using the infinity norm) for capturing local-scale mixing features such as interfacial mixing, only 24% to 32% of the sequence data for model training is required. To achieve the same level of accuracy for capturing global-scale mixing features, the sequence data required for model training is 64% to 70% of the total spatial-temporal data. Hence, the proposed approach—a fast and accurate way to forecast long-time spatial-temporal mixing patterns in heterogeneous and anisotropic media—will be a valuable tool for modeling reactive-transport in a wide range of applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Thermodynamics-guided machine learning model for predicting convective boundary layer height and its multi-site applicability

Accurate estimation of convective boundary layer height (CBLH) is vital for weather, climate, and air quality modeling. Machine learning (ML) shows promise in CBLH prediction, but input parameter selection often lacks physical grounding, limiting generalizability. This study introduces a novel ML framework for CBLH prediction, integrating thermodynamic constraints and the diurnal CBLH cycle as an implicit physical guide. Boundary layer growth is modeled as driven by surface heat fluxes and atmospheric heat absorption represented with the low tropospheric stability, using the diurnal cycle as input and output. TPOT and AutoKeras are employed to select optimal models, validated against Doppler lidar-derived CBLH data, achieving an R 2 of 0.84 across untrained years. Comparisons of eddy covariance (ECOR) and energy balance Bowen ratio (EBBR) flux measurements show the same prediction capability. Models trained on the ARM SGP C1 site with ECOR data and tested at E37 and E39 yield R 2 values of 0.79 and 0.81, respectively, demonstrating their adaptability. The ML model trained with all sites' data slightly enhances the performance compared with ML models trained over single-site data. The interquartile range for predicted CBLH is consistently narrower than that for DL-derived CBLH, reflecting lower variability in predicted CBLH compared to DL-derived CBLH, which is influenced by additional factors, which are not well represented with the model inputs. The model's generalizability across multiple sites at the ARM SGP site demonstrates its potential for transfer to greater distances, offering a scalable approach for enhancing boundary layer parameterization in atmospheric models.

Chu, Yufei [Stony Brook Univ., NY (United States)]↗

Turbulence Properties at the ENA site

The data set includes the turbulence properties as derived from the data collected by the Doppler lidar (DL) at the ARM Eastern North Atlantic (ENA) site. Among the reported turbulence properties included are profiles of variance of vertical air motion, skewness of vertical air motion, updraft strength, downdraft strength, updraft fraction, and downdraft fraction. Due to the persistent presence of drizzle at the location, the DL echoes were classified as those corresponding to aerosols and drizzle using the backscatter cross-section threshold. This daily backscatter cross-section threshold is also reported. The data set begins on January 1, 2015 and ends January 1, 2021. The gaps in the data set correspond to the periods when the instrument was not operational. The header of the NetCDF file is as below. netcdf dl.turb.stats { dimensions: time = 52609; height = 241; threshold_time = 1828; variables: double time(time); time:Long_Name = "seconds since 1970-1-1"; time:units = "seconds"; double height(height); height:Long_Name = "height of the center of the range gate above ground"; height:units = "meters"; double echo_fraction(time, height); echo_fraction:Long_Name = "Fraction of hour with valid Doppler lidar echoes" ; echo_fraction:units = "N/A"; double wvar(time, height); wvar:Long_Name = "Variance of vertical air motion"; wvar:units = "m^2/s^2"; double wskew(time, height); wskew:Long_Name = "Skewness of vertical air motion"; wskew:units = "N/A"; double updraft_velocity(time, height); updraft_velocity:Long_Name = "Hourly averaged updraft velocity"; updraft_velocity:units = "m/s"; double downdraft_velocity(time, height); downdraft_velocity:Long_Name = "Hourly averaged downdraft velocity"; downdraft_velocity:units = "m/s"; double updraft_fraction(time, height); updraft_fraction:Long_Name = "Hourly updraft fraction"; updraft_fraction:units = "N/A"; double downdraft_fraction(time, height); downdraft_fraction:Long_Name = "Hourly downdraft fraction"; downdraft_fraction:units = "N/A"; double threshold_time(threshold_time); threshold_time:Long_Name = "Seconds since 1970-1-1 for the day"; threshold_time:units = "seconds"; double bscat_threshold(threshold_time); bscat_threshold:Long_Name = "Backscatter threshold used for discerning lidar returns from hydrometeors and aerosols."; bscat_threshold:units = "1/(m*Sr)"; }

54 ENVIRONMENTAL SCIENCES↗

ARM SGP PBLH and MLH datasets from Raman lidar and Doppler lidar

The planetary boundary layer (PBL) plays a critical role in the atmosphere by transferring heat, moisture, and momentum. The warm PBL has a distinct diurnal cycle including the daytime convective mixing layer (ML) and nighttime residual layer developments. Thus, simultaneous determinations of PBL height (PBLH) and ML height (MLH) are necessary for studying PBL characterization and processes. Here, new approaches are developed to provide reliable PBLH and MLH estimates to characterize warm PBL evolution. The approaches use Raman lidar (RL) water vapor mixing ratio (WVMR) and Doppler lidar (DL) vertical velocity measurements at the Southern Great Plains (SGP) atmospheric observatory, which was established by the Atmospheric Radiation Measurement (ARM) User Facility. Compared to widely used lidar aerosol measurements for PBLH, WVMR is a better tracer for PBL vertical mixing. For PBLH, the approach classifies PBL water vapor structures into a few general patterns, then uses a slope method and dynamic threshold method to determine PBLH. For MLH, wavelet analysis is used to reconstruct 2D variance from DL vertical wind velocity measurements according to the turbulence eddy size to minimize the impacts of gravity wave and eddy size on variance calculations; then, a dynamic threshold method is used to determine MLH. Remotely-sensed PBLHs and MLHs are compared with radiosonde measurements based on the Richardson number method. Good agreements between them confirm that the proposed new algorithms are reliable for PBLH and MLH characterization. The algorithms are applied to warm-season RL and ML measurements at the SGP site for five years to study warm-season PBL structure and processes. The weekly composited diurnal evolutions of PBLHs and MLHs in a warm climate were provided to illustrate diurnal and seasonal PBL evolutions. This reliable data set of PBLH and MLH values will be valuable for studying PBL processes, model evolution, and PBL parameterization improvements. The MLH dataset includes the MLH in values of km above ground level. The PBLH dataset includes the PBLH in values of km above ground level, along with a flag ("situation_PBLH") to determine the state of the PBL (1 = Cloudy Condition, 2 = Stable Layer, 3 = Multi-layer WVMR structure, 4 = Well-Mixed PBL, 5 = A de-coupled layer, 6 = Other).

mixing layer height↗

Globally Gridded Groundwater Extraction Volumes and Costs under Six Depletion and Ponded Depth Targets

This repository contains simulated outputs from superwell – a hydro-economic tool for long-term assessment of groundwater cost and supply – providing globally gridded groundwater extractable volumes and associated unit costs ($/km³) for accessible groundwater production, based on a variety of user-defined depletion and ponded depth scenarios. Key model documentation: Niazi, H., Ferencz, S. B., Graham, N. T., Yoon, J., Wild, T. B., Hejazi, M., Watson, D. J., & Vernon, C. R. (2025). Long-term hydro-economic analysis tool for evaluating global groundwater cost and supply: Superwell v1.1. Geoscientific Model Development, 18(5), 1737-1767. https://doi.org/10.5194/gmd-18-1737-2025 Find the source code of the superwell model on GitHub: https://github.com/JGCRI/superwell Repository Overview Main output: superwell_outputs.7z contains 6 files (4.5 GB) named as superwell_py_deep_all_0.*PD_0.*DL.csv. These files present superwell outputs of global groundwater extraction volumes and cost estimates on a 0.5° scale for six scenarios with different Ponded Depth (PD; 0.3 and 0.6 m) and Depletion Limit (DL; 5%, 25%, and 40% of available volume) targets over the entire pumping lifetime of a grid cell superwell_py_deep_all_0.3PD_0.25DL_sample_100.csv contains superwell outputs for 100 data points sampled to match the global inputs' distribution superwell_py_deep_all_0.3PD_0.25DL_Grid_72548.csv contains superwell output for a single grid cell concept_v5.png provides an overview of the superwell workflow Outputs Description year_number: year of pumping depletion_limit: set depletion limit (DL) as a volume fraction of total available groundwater Mappings: continent, country, gcam_basin_id, Basin_long_name, grid_id: geographic identifiers and basin information Inputs: grid_area (km²): area of the grid cell whyclass: hydrogeological classification of the aquifer permeability (m/day), porosity (%), total_thickness (m), depth_to_water (m): aquifer properties. The geo-processed input data has been published separately: https://doi.org/10.57931/2307831 Model outputs: orig_aqfr_sat_thickness (m), aqfr_sat_thickness (m): original and remaining/instantaneous saturated thickness of the aquifer hydraulic_conductivity (m/day), transmissivity (m²/day): hydraulic properties of the aquifer radius_of_influence (m), areal_extent (km²): well radius and area of influence from the center of the well number_of_wells (-): number of wells in a grid cell determined by a ratio of well area and grid area max_drawdown (m), drawdown (m), drawdown_interference (m): well and aquifer drawdown during extraction total_head (m): total lift for the groundwater (depth to water plus drawdown) total_well_length (m): total depth of wells drilled well_yield (m³/day): pumping rate or well yield power (kW), energy (kWh): power and energy required for pumping groundwater Volume Outputs: volume_produced_perwell (m³), cumulative_vol_produced_perwell (m³): production volume metrics per well volume_produced_allwells (m³), cumulative_vol_produced_allwells (m³): aggregate extraction volumes for all wells in a grid cell available_volume (m³): available groundwater in storage for the grid cell as determined by aquifer properties depleted_vol_fraction: fraction of total volume pumped over available volumes in a grid cell (same as depletion limit) Cost Outputs: well_installation_cost ($): well installation cost based on the hydrogeological complexity of the aquifer annual_capital_cost, maintenance_cost, nonenergy_cost ($): nonenergy costs energy_cost_rate ($/kWh): electricity rate energy_cost ($): energy cost of pumping groundwater total_cost_perwell ($), total_cost_allwells ($): total annual energy and non-energy cost for each and all wells in a grid cell a unit_cost ($/m³), unit_cost_per_km3 ($/km³), unit_cost_per_acreft ($/acre-ft): total cost of pumping a unit of groundwater, indicated for different spatial units Key Resources Model documentation: Niazi, H., Ferencz, S., Graham, N., Yoon, J., Wild, T., Hejazi, M., Watson, D., & Vernon, C. (2024; In-prep). Long-term Hydro-economic Assessment Tool for Evaluating Global Groundwater Cost and Supply: Superwell v1. Geoscientific Model Development. Input data: Niazi, H., Watson, D., Hejazi, M., Yonkofski, C., Ferencz, S., Vernon, C., Graham, N., Wild, T., & Yoon, J. (2024). Global Geo-processed Data of Aquifer Properties by 0.5° Grid, Country and Water Basins. MSD-LIVE Data repository. https://doi.org/10.57931/2307831 superwell source code: https://github.com/JGCRI/superwell Cite as Niazi, H., Ferencz, S., Yoon, J., Graham, N., Wild, T., Hejazi, M., Watson, D., & Vernon, C. (2024). Globally Gridded Groundwater Extraction Volumes and Costs under Six Depletion and Ponded Depth Targets. MSD-LIVE Data repository. https://doi.org/10.57931/2307832 Contact Reach out to Hassan Niazi or Stephen Ferencz or open an issue in the superwell repository for questions or suggestions.

Earth Systems↗

The radiolysis of tryptophan and leucine with P-32 beta-radiation

The paper extends earlier experiments on the radiolysis of DL-tryptophan using P-32 beta-radiation to longer reaction times, observing complete destruction of tryptophan by secondary, nonradiolytic processes. In addition DL-leucine is irradiated with P-32 beta-irradiation at -196 C, leading to radiolyses to the extents of about 20-30%, but observing no concomitant asymmetric bias. The complete absence of asymmetric bias in the present and earlier (Bonner et al., 1979) radiolyses of aqueous tryptophan at -25 C and the present radiolyses of water-free leucine at -196 C using P-32 beta-radiation and its accompanying bremsstrahlung leave it an open question whether or not the Vester-Ulbricht beta-decay/bremsstrahlung mechanism for the origin of optical activity is a viable one.

Blair, N. E.↗

The radiolysis and radioracemization of amino acids on silica surfaces

Results are presented of experiments on the radioracemization of amino acids in the presence of silica surfaces such as may have been found on the prebiotic earth. L-leucine and a DL-leucine mixture deposited on samples of 1-quartz and an amorphous silica preparation (Syloid 63) was subjected to Co-60 gamma-ray irradiation, then analyzed by gas chromatography to determine the radiolysis and racemization rates. The quartz surface is found to have a marginal efficacy in enhancing radiolysis when compared with a crystalline L-leucine control, although enhancing radioracemization symmetrically by a factor of two. Both the radiolysis and radioracemization of L-leucine and DL-leucine on a Syloid-63 silica surface are observed to increase with increasing radiation dose, and to be substantially greater than in the crystalline controls. Additional experiments with the nonprotein amino acid isovaline deposited on Syloid 63 confirm the greater radiolysis susceptibility of amino acids deposited on silica with respect to the crystalline state, although racemization is not observed. The observations suggest that the presence of a silica surface would have a deleterious effect on any mechanism for the origin of molecular chirality relying on stereoselective beta-radiolysis.

Bonner, W. A.↗

Search for effect of longitudinally polarized protons on optically active amino acids

The influence of irradiation by longitudinally polarized protons on the differential decomposition of amino acid is investigated in a test of the Vester Ulbricht hypothesis that L-amino acids were produced in preference to the D isomers on the primitive earth due to the effects of parity violation in beta decay. Samples of DL-leucine were irradiated with protons of both positive and negative longitudinal polarization, then analyzed by gas chromatography. Despite advantages of higher polarization, lower velocity and higher ionization density of protons with respect to electrons, proton irradation is found to lead to no detectable asymmetries in DL-leucine degradation, even at 50% gross degradation.

Lemmon, R. M.↗

Determination of hydrogen permeability in commercial and modified superalloys

The results of hydrogen permeability measurements on several iron- and cobalt-base alloys as well as on two long-ranged ordered alloys over the range of 705 to 870 C (1300 to 1600 F) are summarized. The test alloys included wrought alloys N-155, IN 800, A-286, 19-9DL, and 19-9DL modifications with aluminum, niobium, and misch metal. In addition, XF-818, CRM-6D, SA-F11, and HS-31 were evaluated. Two wrought long-range ordered alloys, Ni3Al and (Fe,Ni)3(V,Al) were also evaluated. All tests were conducted at 20.7 MPa pressure in either pure and/or 1% CO2-doped H2 for test periods as long as 133 h. Detailed analyses were conducted to determine the relative permeability rankings of these alloys and the effect of doping, exit surface oxidation, specimen design variations, and test duration on permeability coefficient, and permeation activation energies were determined. The two long-range ordered alloys had the lowest permeability coefficients in pure H2 when compared with the eight commercial alloys and their modifications. With CO2 doping, significant decrease in permeability was observed in commercial alloys--no doped tests were conducted with the long-range ordered alloys.

Bhattacharyya, S.↗

CO emission structure in the Taurus molecular cloud complex

A survey is reported of the J = 1 - 0 2.6-mm transition of CO, covering 11.5 square degrees toward the Taurus dark cloud complex, including Heiles Cloud 2 and DL 2. The main finding is the presence of systematic velocity structure in the Barnard 18 region of DL 2. This structure is tentatively interpreted as an expanding shell of molecular gas of diameter about 3 pc centered on a cluster of young T Tauri stars. The estimated momentum of the shell is consistent with that expected from winds from the young stars seen in the shell.

Murphy, D. C.↗

Simulation of double layers in a model auroral circuit with nonlinear impedance

A reduced circuit description of the U-shaped potential structure of a discrete auroral arc, consisting of the flank transmission line plus parallel-electric-field region, is used to provide the boundary condition for one-dimensional simulations of the double-layer evolution. The model yields asymptotic scalings of the double-layer potential, as a function of an anomalous transport coefficient alpha and of the perpendicular length scale l(a) of the arc. The arc potential phi(DL) scales approximately linearly with alpha, and for alpha fixed phi (DL) about l(a) to the z power. Using parameters appropriate to the auroral zone acceleration region, potentials of phi (DPL) 10 kV scale to projected ionospheric dimensions of about 1 km, with power flows of the order of magnitude of substorm dissipation rates.

Smith, R. A.↗

Application of boost guidance to NASA sounding rocket launch operations at the White Sands Missile Range

This paper addresses the unique problems associated with launching the Black Brant V, VIII, and IX sounding rocket vehicles at White Sands Missile Range (WSMR) and the significance of the introduction of the S19 to the NASA Goddard Space Flight Center Wallops Flight Facility sounding rocket program in terms of launch flexibility, improved impact dispersion, higher flight reliability, and reduced program costs. This paper also discusses salient flight results from NASA 36.011UL (the first S19 guided Black Brant launched at WSMR) and the NASA Comet Halley missions (36.010DL and 36.017DL).

Montag, W. H.↗

Electric fields and double layers in plasmas

Various mechanisms for driving double layers (DLs) in plasmas are described, including applied potential drops, currents, contact potentials, and plasma expansions. Somne dynamic features of the DLs are discussed; and it is demonstrated that DLs and the currents through them undergo slow oscillations, determined by the ion transit time across an effective length of the system in which the DLs form. It is shown that a localized potential dip forms at the low potential end of a DL, which interrupts the electron current through it according to the Langmuir criterion whenever the ion flux into the DL is disrupted. Also considered is the generation of electric fields perpendicular to the ambient magnetic field by contact potentials.

Singh, Nagendra↗

Effect of double layers on magnetosphere-ionosphere coupling

The dynamic aspects of auroral current structures are reviewed with emphasis on consequences for models of microscopic turbulence (MT). A number of models of MT are introduced into a large-scale model of Alfven wave propagation to determine the effect of various models on the overall structure of auroral currents. The effect of a double layer (DL) electric field which scales with the plasma temperature and the Debye length is compared with the effect of anomalous resistivity due to electrostatic ion cyclotron turbulence in which the electric field scales with the magnetic field strength. It is shown that the DL model is less diffusive than the resistive model, indicating the possibility of narrow intense current structures.

Lysak, Robert L.↗

Anomalous transport in discrete arcs and simulation of double layers in a model auroral circuit

The evolution and long-time stability of a double layer in a discrete auroral arc requires that the parallel current in the arc, which may be considered uniform at the source, be diverted within the arc to charge the flanks of the U-shaped double-layer potential structure. A simple model is presented in which this current re-distribution is effected by anomalous transport based on electrostatic lower hybrid waves driven by the flank structure itself. This process provides the limiting constraint on the double-layer potential. The flank charging may be represented as that of a nonlinear transmission. A simplified model circuit, in which the transmission line is represented by a nonlinear impedance in parallel with a variable resistor, is incorporated in a 1-d simulation model to give the current density at the DL boundaries. Results are presented for the scaling of the DL potential as a function of the width of the arc and the saturation efficiency of the lower hybrid instability mechanism.

Smith, Robert A.↗

Martian impact craters: Continuing analysis of lobate ejecta sinuosity

The lobate ejecta morphology surrounding most fresh Martian impact craters can be quantitatively analyzed to determine variations in ejecta sinuosity with diameter, latitude, longitude, and terrain. The results of such studies provide another clue to the question of how these morphologies formed: are they the results of vaporization of subsurface volatiles or caused by ejecta entrainment in atmospheric gases. Kargel provided a simple expression to determine the degree of non-circularity of an ejecta blanket. This measure of sinuosity, called 'lobateness', is given by the ratio of the ejecta perimeter to the perimeter of a circle with the same area as that of the ejecta. The Kargel study of 538 rampart craters in selected areas of Mars led to the suggestion that lobateness increased with increasing diameter, decreased at higher latitude, and showed no dependence on elevation or geologic unit. Major problems with the Kargel analysis are the limited size and distribution of the data set and the lack of discrimination among the different types of lobate ejecta morphologies. Bridges and Barlow undertook a new lobateness study of 1582 single lobe (SL) and 251 double lobe (DL) craters. The results are summarized. These results agree with the finding of Kargel that lobateness increases with increasing diameter, but found no indication of a latitude dependence for SL craters. The Bridges and Barlow study has now been extended to multiple lobe (ML) craters. Three hundred and eighty ML craters located across the entire Martian surface were studied. ML craters provide more complications to lobateness studies than do SL and DL craters - in particular, the ejecta lobes surrounding the crater are often incomplete. Since the lobateness formula compares the perimeter of the ejecta lobe to that of a circle, the analysis was restricted only to complete lobes. The lobes are defined sequentially starting with the outermost lobe and moving inward.

Barlow, Nadine G.↗

Hubble space telescope observations of young star clusters in NGC-4038/4039, 'the antennae' galaxies

New, high-resolution images of the disks of NGC 4038/4039 obtained with the Wide Field Camera of the Hubble Space Telescope (HST) are presented. NGC 4038/4039, nicknamed The Antennae, is a prototypical example of a pair of colliding galaxies believed to be at an early stage of a merger. Down to the limiting magnitude of V approximately 23 mag, the HST images reveal a population of over 700 blue pointlike objects within the disks. The mean absolute magnitude of these objects is M(sub V) = -11 mag, with the brightest objects reaching M(sub V) approximately -15. Their mean apparent color indices ar U - V = -0.7 mag and V - 1 = 0.8 mag on the Johnson UVI passband system, while their mean indices corrected for internal reddening are (u - v)(sub 0) = -1.0 mag and (V - I(sub 0) = 0.5. Their mean effective radius, determined from slightly resolved images, is 18 pc (for H(sub 0) = 50 km/s /Mpc). Based on their luminosities and resolution, most of these objects cannot be individual stars, but are likely young compact star clusters. The brighter ones are similar to the objects found in NGC 1275 and NGC 7252, which appear to be young globular clusters formed during recent galazy mergers. Based on their U - V and V - I colors, the brightest, bluest clusters of NGC 4038/4039 appear to be less than 10 Myr old. Most of these bright clusters are relatively tightly clustered themselves, with typically a dozen individual clusters belonging to a complex identified as a giant H II region from ground-based observations. The cluster luminosity function (LF) is approximately a power law, phi(L)dL proportional to L(exp -1.78+/-0.05)dL, with no hint of a turnover at fainter magnitudes. This power-law shape agrees with the LF of Magellanic Cloud clusters and Galactic open clusters, but differs from the LF of old globular cluster systems that is typically Gaussian with a Full Width at Half Maximum (FWHM) of approximately 3 mag. Besides the blue clusters, we also find about a dozen extremely red objects with V - I greater than 3.0. The highest number density of these red objects is found in the SE quadrant, where star formation appears to be most recent. We propose that these objects may be very young star clusters still embedded in their placental dust cocoons.

Whitmore, Bradley C.↗