Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Model fitting”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Candidate massive galaxies at z ~ 4 in the Dark Energy Survey

Using stellar population models, we predicted that the Dark Energy Survey (DES) - due to its special combination of area (5000 deg. sq.) and depth ($i = 24.3$) - would be in the position to detect massive ($\gtrsim 10^{11}$ M$_{\odot}$) galaxies at $z \sim 4$. We confront those theoretical calculations with the first $\sim 150$ deg. sq. of DES data reaching nominal depth. From a catalogue containing $\sim 5$ million sources, $\sim26000$ were found to have observed-frame $g-r$ vs $r-i$ colours within the locus predicted for $z \sim 4$ massive galaxies. We further removed contamination by stars and artefacts, obtaining 606 galaxies lining up by the model selection box. We obtained their photometric redshifts and physical properties by fitting model templates spanning a wide range of star formation histories, reddening and redshift. Key to constrain the models is the addition, to the optical DES bands $g$, $r$, $i$, $z$, and $Y$, of near-IR $J$, $H$, $K_{s}$ data from the Vista Hemisphere Survey. We further applied several quality cuts to the fitting results, including goodness of fit and a unimodal redshift probability distribution. We finally select 233 candidates whose photometric redshift probability distribution function peaks around $z\sim4$, have high stellar masses ($\log($M$^{*}$/M$_{\odot})\sim 11.7$ for a Salpeter IMF) and ages around 0.1 Gyr, i.e. formation redshift around 5. These properties match those of the progenitors of the most massive galaxies in the local universe. This is an ideal sample for spectroscopic follow-up to select the fraction of galaxies which is truly at high redshift. These initial results and those at the survey completion, which we shall push to higher redshifts, will set unprecedented constraints on galaxy formation, evolution, and the re-ionisation epoch.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Sorption testing and generalized composite surface complexation models for determining uranium sorption parameters at a proposed in-situ recovery site

Solid-phase iron concentrations and generalized composite surface complexation models were used to evaluate procedures in determining uranium sorption on oxidized aquifer material at a proposed U in situ recovery (ISR) site. At the proposed Dewey Burdock ISR site in South Dakota, USA, oxidized aquifer material occurs downgradient of the U ore zones. Solid-phase Fe concentrations did not explain our batch sorption test results,though total extracted Fe appeared to be positively correlated with overall measured U sorption. Batch sorption test results were used to develop generalized composite surface complexation models that incorporated the full genericsorption potential of each sample, without detailed mineralogiccharacterization. The resultant models provide U sorption parameters (site densities and equilibrium constants) for reactive transport modeling. The generalized composite surface complexation sorption models were calibrated to batch sorption data from three oxidized core samples using inverse modeling, and gave larger sorption parameters than just U sorption on the measured solidphase Fe. These larger sorption parameters can significantly influence reactive transport modeling, potentially increasing U attenuation. Because of the limited number of calibration points, inverse modeling required the reduction of estimated parameters by fixing two parameters. The best-fit models used fixed values for equilibrium constants, with the sorption site densities being estimated by the inversion process. While these inverse routines did provide best-fit sorption parameters, local minima and correlated parameters might require further evaluation. Despite our limited number of proxy samples, the procedures presented provide a valuable methodology to consider for sites where metal sorption parameters are required. Furthermore, these sorption parameters can be used in reactive transport modeling to assess downgradient metal attenuation, especially when no other calibration data are available, such as at proposed U ISR sites.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Stellar population analysis of MaNGA early-type galaxies: IMF dependence and systematic effects

We study systematics associated with estimating simple stellar population (SSP) parameters – age, metallicity [M/H], α-enhancement [α/Fe], and initial mass function (IMF) shape – and associated M * /L gradients, of elliptical slow rotators (E-SRs), fast rotators (E-FRs), and S0s from stacked spectra of galaxies in the MaNGA survey. These systematics arise from (i) how one normalizes the spectra when stacking; (ii) having to subtract emission before estimating absorption line strengths; (iii) the decision to fit the whole spectrum or just a few absorption lines; (iv) SSP model differences (e.g. isochrones, enrichment, IMF). The MILES+Padova SSP models, fit to the H β , 〈Fe〉, TiO 2SDSS , and [MgFe] Lick indices in the stacks, indicate that out to the half-light radius R e : (a) ages are younger and [α/Fe] values are lower in the central regions but the opposite is true of [M/H]; (b) the IMF is more bottom-heavy in the center, but is close to Kroupa beyond about R e /2; (c) this makes M * /L about 2 × larger in the central regions than beyond R e /2. While the models of Conroy et al. return similar [M/H] and [α/Fe] profiles, the age and (hence) M * /L profiles can differ significantly even for solar abundances and a Kroupa IMF; different responses to non-solar abundances and IMF parametrization further compound these differences. There are clear (model independent) differences between E-SRs, E-FRs, and S0s: younger ages and less enhanced [α/Fe] values suggest that E-FRs and S0s are not SSPs, but relaxing this assumption is unlikely to change their inferred M * /L gradients significantly.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

RUPRECHT 147: THE OLDEST NEARBY OPEN CLUSTER AS A NEW BENCHMARK FOR STELLAR ASTROPHYSICS

Ruprecht 147 is a hitherto unappreciated open cluster that holds great promise as a standard in fundamental stellar astrophysics. We have conducted a radial velocity survey of astrometric candidates with Lick, Palomar, and MMT observatories and have identified over 100 members, including 5 blue stragglers, 11 red giants, and 5 double-lined spectroscopic binaries (SB2s). We estimate the cluster metallicity from spectroscopic analysis, using Spectroscopy Made Easy (SME), and find it to be [M/H] = +0.07 {+-} 0.03. We have obtained deep CFHT/MegaCam g'r'i'z' photometry and fit Padova isochrones to the (g' - i') and Two Micron All Sky Survey (J - K{sub S} ) color-magnitude diagrams, using the {tau}{sup 2} maximum-likelihood procedure of Naylor, and an alternative method using two-dimensional cross-correlations developed in this work. We find best fits for Padova isochrones at age t = 2.5 {+-} 0.25 Gyr, m - M = 7.35 {+-} 0.1, and A{sub V} = 0.25 {+-} 0.05, with additional uncertainty from the unresolved binary population and possibility of differential extinction across this large cluster. The inferred age is heavily dependent on our choice of stellar evolution model: fitting Dartmouth and PARSEC models yield age parameters of 3 Gyr and 3.25 Gyr, respectively. At {approx}300 pc and {approx}3 Gyr, Ruprecht 147 is by far the oldest nearby star cluster.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

A predictive 0-D HCCI combustion model for ethanol, natural gas, gasoline, and primary reference fuel blends

Homogeneous Charge Compression Ignition (HCCI) is a promising advanced combustion concept with high thermal efficiency and low exhaust emissions. This work purposes a computationally-efficient, zero-dimensional (0-D) HCCI combustion model that does not need to solve any differential equations. To develop the burn rate correlations, experimental HCCI data of ethanol, natural gas, E10-gasoline, and Primary Reference Fuel (PRF) blends was collected on a CFR engine. The burn rate model is built based on the individual cycle mass fraction burned (MFB) curves calculated from the experimental data. CA0 can be predicted from an ignition delay correlation. Once CA0 is known, CA50 and CA90 can be predicted based on CA0 and the charge-mass equivalence ratio Φ'. Then, with CA0, 50 and 90 known, a Wiebe function can be constructed to represent the MFB curve and burn rate. For PRF blends, low-temperature and high-temperature heat releases (LTHR and HTHR) are modeled by two Wiebe functions with Φ' and PRF number dependence. The fitted model has high accuracy compared to the experimental data with R 2 values generally greater than 0.97. Finally, various ignition delay correlations from the literature are also tested against the experimentally collected data and some modifications to the correlations are suggested.

42 ENGINEERING↗
DOE OSTI2020

Relationship between ductility and the porosity of additively manufactured AlSi10Mg

Additive manufacturing via selective laser melting can result in variable levels of internal porosity both between build plates and within components from the same build. In this work, sample porosity levels were compared to tensile properties for 176 samples spanning eight different build plates. Sample porosity was measured both by Archimedes density, which provided an estimation of overall porosity, and by observation of voids in the fracture surface, which provided an estimation of the porosity at the failure plane. The porosity observed at the fracture surface consistently demonstrated higher porosity than that suggested by Archimedes density. The porosity values obtained from both methods were compared against the mechanical results. Sample porosity appears to have some correlation to the ultimate tensile strength, yield strength, and modulus, but the strongest relationship is observed between porosity and ductility. Three different models were used to relate the fracture surface porosity to the ductility. The first method was a simple linear regression analysis, while the other two models have been used to relate porosity to ductility in cast alloys. It is shown that all three models fit the data well over the observed porosity ranges, suggesting that the models taken from casting theory can extend to additively manufactured metals. Finally, it is proposed that the non-destructive Archimedes method could be used to estimate an approximate sample ductility through the use of correlations realized here. Such a relationship could prove useful for design and for a deeper understanding of the impact of pores on tensile behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Direct identification of predator-prey dynamics in gyrokinetic simulations

The interaction between spontaneously formed zonal flows and small-scale turbulence in nonlinear gyrokinetic simulations is explored in a shearless closed field line geometry. It is found that when clear limit cycle oscillations prevail, the observed turbulent dynamics can be quantitatively captured by a simple Lotka-Volterra type predator-prey model. Fitting the time traces of full gyrokinetic simulations by such a reduced model allows extraction of the model coefficients. Scanning physical plasma parameters, such as collisionality and density gradient, it was observed that the effective growth rates of turbulence (i.e., the prey) remain roughly constant, in spite of the higher and varying level of primary mode linear growth rates. The effective growth rate that was extracted corresponds roughly to the zonal-flow-modified primary mode growth rate. It was also observed that the effective damping of zonal flows (i.e., the predator) in the parameter range, where clear predator-prey dynamics is observed, (i.e., near marginal stability) agrees with the collisional damping expected in these simulations. This implies that the Kelvin-Helmholtz-like instability may be negligible in this range. The results imply that when the tertiary instability plays a role, the dynamics becomes more complex than a simple Lotka-Volterra predator prey.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Characterizing Baselines and Change in Gas Hydrate Systems using EM Methods

The objective of this project was to advance our understanding of gas hydrate systems in nature by characterizing their electrical properties in the field and in the laboratory. In the laboratory measurements, methane hydrate was synthesized from pure water ice and flash frozen seawater, with varying amounts of sand or silt added. Electrical conductivity was determined by impedance spectroscopy, using equivalent circuit modeling to separate the effects of electrodes and to gain insight into conduction mechanisms. Silt and sand increase the conductivity of pure hydrate, inferred to be contaminant NaCl contributing to conduction in hydrate, to a peak conductivity in agreement with peak resistivities observed in well logs through massive hydrate (3,000--10,000 Ωm). The addition of silt and sand lowers the conductivity of hydrate synthesized from seawater, by an amount consistent with Archie's Law. All samples were characterized using cryogenic scanning electron microscopy and energy dispersive spectroscopy, which shows good connectivity of salt and brine phases. Electrical conductivity measurements of pure hydrate and hydrate mixed with silt during pressure-induced dissociation supports previous conclusions that sediment increases dissociation rate. In order to characterize gas hydrate systems in the field, we collected 360 line kilometers of controlled-source electromagnetic data on Walker Ridge 313, Orca Basin (WR100), Mad Dog (GC781), and Green Canyon 955 in the Gulf of Mexico, all areas with known or seismically inferred gas hydrate deposits and which have be drilled or targeted for future drilling. We deep-towed an EM transmitter that generates an alternating electric field which propagates through the seafloor geology. Data were recorded on 6 receivers towed behind the transmitter at distances between 550 and 1550 m. In the presence of conductive geology, the electric fields will be attenuated, and conversely, in resistive geology the fields will be preserved. Our data were inverted using a 2D inversion method that first optimizes the model-data misfit, then finds the smoothest model fitting the data. This ensures that resistivity structures present in the final model are likely necessary. At each of the proposed drilling sites we found increased resistivity, interpreted as increased hydrate concentrations. However, not only were the primary drilling sites not always more resistive than the alternate sites, at WR313 the strongest resistors were not at the locations targeted for drilling.

03 NATURAL GAS↗
DOE OSTI2019

GD358: Three Decades of Observations for the In-depth Asteroseismology of a DBV Star

Here, we report on the analysis of 34 years of photometric observations of the pulsating helium atmosphere white dwarf GD358. The complete data set includes archival data from 1982 to 2006, and 1195.2 hr of new observations from 2007 to 2016. From this data set, we extract 15 frequencies representing g-mode pulsation modes, adding 4 modes to the 11 modes known previously. We present evidence that these 15 modes are ℓ = 1 modes, 13 of which belong to a consecutive sequence in radial overtone k. We perform a detailed asteroseismic analysis using models that include parameterized, complex, carbon and oxygen core composition profiles to fit the periods. Recent spectroscopic analyses place GD358 near the red edge of the DBV instability strip, at 24,000 ± 500 K and a $\mathrm{log}\,{\text{}}g$ of 7.8 ± 0.08 dex. The surface gravity translates to a mass range of 0.455–0.540 ${M}_{\odot }$. Our best-fit model has a temperature of 23,650 K and a mass of 0.5706 ${M}_{\odot }$. That is slightly more massive than what is suggested by the most recent spectroscopy. We find a pure helium layer mass of 10–5.50, consistent with the result of previous studies and the outward diffusion of helium over time.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Photometric redshifts and intrinsic alignments: Degeneracies and biases in the 3 × 2 pt analysis

We present a systematic study of cosmological parameter bias in weak lensing and large-scale structure analyses for upcoming imaging surveys induced by the interplay of intrinsic alignments (IA) and photometric redshift (photo-z) model misspecification error. We first examine the degeneracies between the parameters of the tidal alignment–tidal torquing (TATT) model for IA and of a photo-z model including a mean shift ( Δ z ¯ ) and variance ( σ z ) for each tomographic bin of lenses and sources, under a variety of underlying true IA behaviors. We identify strong degeneracies between: (1) the redshift scaling of the tidal alignment amplitude and the mean shift and variances of source bins, (2) the redshift scaling of the tidal torquing amplitude and the variance of the lowest- z source bin, and (3) the IA source density weighting and the mean shift and variance of several source bins. We then use this information to guide our exploration of the level of cosmological parameter bias which can be induced given incorrect modeling of IA, photo-z, or both. We find that marginalizing over all the parameters of TATT is generally sufficient to preclude cosmological parameter bias in the scenarios we consider. However, this does not necessarily mean that IA and photo-z parameters are themselves unbiased, nor does it mean that the best-fit model is a good fit to the data. We also find scenarios where the inferred parameters produce χ d . o . f . 2 values indicative of a good fit but cosmological parameter bias is significant, particularly when the IA source density weighting parameter is not marginalized over. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

A Perspective on Traditional and Data Driven Electrochemical Modeling and Analysis

To understand the behavior of electrochemical systems, we need to reduce the dimensionality of the measured current-voltage-time (I-V-t) data by fitting models, thus enabling us to analyze and compare the governing physics. Traditionally, the process for this is an 'expert first' approach: defining the model and its explicit assumptions based on inductive reasoning or empirical observation, fitting small portions of the I-V-t data where assumptions are most valid or carefully designing experiments to enforce key assumptions, and then interpreting the model parameters. However, modern data-driven methods enable a new paradigm: a 'data first' approach, where the latent behaviors governing the system's measured response are identified directly using machine-learning models that optimize both model structure and parameters from the I-V-t data, guaranteeing that the learned model explains as much of the observed system response as possible. After model identification, the model can then be interrogated by an expert to connect observed behaviors with underlying physics. This talk will review several different types of electrochemical analysis (electrochemical impedance, differential voltage-capacity, electrochemical kinetics) and compare the traditional and data-driven methods for analyzing the data.

42 ENGINEERING↗
DOE OSTI2015

The shape of protein–polymer conjugates in dilute solution

ABSTRACT Protein–polymer conjugation can significantly affect many different aspects of protein behavior, ranging from their solution properties to their ability to form solution and bulk nanostructured materials. An underlying fundamental question is how the molecular design affects the shape of the conjugate and, consequently, its properties. This work measures the molecular configuration of model protein–polymer conjugates in dilute solution using small‐angle neutron scattering (SANS) and uses quantitative model fitting to understand the shape of the molecules. Form factor measurements of four model bioconjugates of the red fluorescent protein mCherry and the polymers poly( N ‐isopropylacrylamide), poly(hydroxypropyl acrylate), poly(oligoethylene glycol acrylate), and poly(ethylene glycol) show that these protein–polymer conjugates are well described by a recently developed scattering function for colloid–polymer conjugates that explicitly incorporates excluded volume interactions in the polymer configuration. In the regime where the protein does not exhibit strong interactions with the polymer, modeling the protein–polymer interactions using a purely repulsive Weeks–Chandler–Andersen potential also leads to a coarse‐grained depiction of the conjugate that agrees well with its scattering behavior. The coarse‐grained model can additionally be used for systems with varying protein–polymer interactions, ranging from purely repulsive to strongly attractive, which may be useful for conjugates with strong electrostatic or hydrophobic attractive interactions. © 2015 Wiley Periodicals, Inc. J. Polym. Sci., Part A: Polym. Chem. 2016, 54, 292–302

Lam, Christopher N.↗
DOE OSTI2016

Nonlinear threshold effect in the Z-scan method of characterizing limiters for high-intensity laser light

A threshold model is described which permits one to determine the properties of limiters for high-powered laser light. It takes into account the threshold characteristics of the nonlinear optical interaction between the laser beam and the limiter working material. The traditional non-threshold model is a particular case of the threshold model when the limiting threshold is zero. The nonlinear characteristics of carbon nanotubes in liquid and solid media are obtained from experimental Z-scan data. Specifically, the nonlinear threshold effect was observed for aqueous dispersions of nanotubes, but not for nanotubes in solid polymethylmethacrylate. The threshold model fits the experimental Z-scan data better than the non-threshold model. Output characteristics were obtained that integrally describe the nonlinear properties of the optical limiters.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

EVOLUTIONARY TRACKS FOR BETELGEUSE

We have constructed a series of nonrotating quasi-hydrostatic evolutionary models for the M2 Iab supergiant Betelgeuse (α Orionis). Our models are constrained by multiple observed values for the temperature, luminosity, surface composition, and mass loss for this star, along with the parallax distance and high-resolution imagery that determines its radius. We have then applied our best-fit models to analyze the observed variations in surface luminosity and the size of detected surface bright spots as the result of up-flowing convective material from regions of high temperature in the surface convective zone. We also attempt to explain the intermittently observed periodic variability in a simple radial linear adiabatic pulsation model. Based on the best fit to all observed data, we suggest a best progenitor mass estimate of 20{sub −3}{sup +5} M{sub ⊙} and a current age from the start of the zero-age main sequence of 8.0–8.5 Myr based on the observed ejected mass while on the giant branch.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Investigating the differential emission measure and energetics of microflares with combined SDO/AIA and RHESSI observations

An important question in solar physics is whether solar microflares, the smallest currently observable flare events in X-rays, possess the same energetic properties as large flares. Recent surveys have suggested that microflares may be less efficient particle accelerators than large flares, and hence contribute less non-thermal energy, which may have implications for coronal heating mechanisms. We therefore explore the energetic properties of microflares by combining EUV and X-ray measurements. We present forward-fitting differential emission measure (DEM) analysis of 10 microflares. The fitting is constrained by combining, for the first time, high-temperature Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) observations and flux data from the Solar Dynamics Observatory (SDO) Atmospheric Imaging Assembly (AIA). Two fitting models are tested for the DEM; a Gaussian distribution and a uniform DEM profile. A Gaussian fit proved unable to explain the observations for any of the studied microflares. However, 8 of 10 events studied were reasonably fit by a uniform DEM profile. Hence microflare plasma can be considered to be significantly multi-thermal, and may not be significantly peaked or contain resolvable fine structure, within the uncertainties of the observational instruments. The thermal and non-thermal energy is estimated for each microflare, comparing the energy budget with an isothermal plasma assumption. From the multi-thermal fits the minimum non-thermal energy content was found to average approximately 30% of the estimated thermal energy. By comparison, under an isothermal model the non-thermal and thermal energy estimates were generally comparable. Hence, multi-thermal plasma is an important consideration for solar microflares that substantially alters their thermal and non-thermal energy content.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Estimation of Fission Product Transport Parameters for Cesium in the AGR-3/4 TRISO Fuel Experiment

A one-dimensional (1D) finite-element model of fission product transport in the AGR-3/4 experiment has been developed using the Multiphysics Object Oriented Simulation Environment (MOOSE) framework and implemented in the fuel performance code, BISON. The model resolves capsule-specific geometries, materials, and temperature histories and simulates radial migration of fission products from the fuel compact through the inner ring, outer ring, and into the sink ring. Model parameters governing diffusion and sorption were estimated for key fission products – cesium (Cs), and europium (Eu) – by simultaneously fitting modeled isotopic concentration profiles and total ring inventories to a post-irradiation experimental measurement. These data include gamma scanning, liquid scintillation for Sr-90, radial deconsolidation leach-burn-leach analysis, tomographic reconstructions, and destructive physical sampling. A mortar-based interfacial sorption framework was implemented to enforce physically consistent mass transfer and flux conservation across gas gaps. Two classes of parameter sets were derived: a least-squares best-fit, and a safety-oriented conservative-fit, what applies strong penalties for underprediction of sink inventories. Across all twelve capsules, the model successfully reproduces the dominant radial transport trends for Cs, Sr, with decreasing concentrations from the compact outward through successive rings. Cs behavior is captured most consistently, while strontium predictions reveal systematic trade-offs between compact accuracy and conservative sink-ring bounding. The results demonstrate that sink ring weighted calibration provides conservative, safety relevant bounds on low temperature fission product transport, but at the cost of underpredicting compact inventories for Sr isotopes. These discrepancies highlight the need for additional physics, including fast-slow diffusion model, incorporating trapping mechanism in the transport behavior. Overall, this work establishes a robust, capsule-specific modeling framework for AGR-3/4 fission product transport and provides a defensible basis for parameter selection in source-term and fuel performance analyses for high temperature gas-cooled reactors.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Machine learning for impurity charge-state transition levels in semiconductors from elemental properties using multi-fidelity datasets

Quantifying charge-state transition energy levels of impurities in semiconductors is critical to understanding and engineering their optoelectronic properties for applications ranging from solar photovoltaics to infrared lasers. While these transition levels can be measured and calculated accurately, such efforts are time-consuming and more rapid prediction methods would be beneficial. Here, we significantly reduce the time typically required to predict impurity transition levels using multi-fidelity datasets and a machine learning approach employing features based on elemental properties and impurity positions. We use transition levels obtained from low-fidelity (i.e., local-density approximation or generalized gradient approximation) density functional theory (DFT) calculations, corrected using a recently proposed modified band alignment scheme, which well-approximates transition levels from high-fidelity DFT (i.e., hybrid HSE06). Further, the model fit to the large multi-fidelity database shows improved accuracy compared to the models trained on the more limited high-fidelity values. Crucially, in our approach, when using the multi-fidelity data, high-fidelity values are not required for model training, significantly reducing the computational cost required for training the model. Our machine learning model of transition levels has a root mean squared (mean absolute) error of 0.36 (0.27) eV vs high-fidelity hybrid functional values when averaged over 14 semiconductor systems from the II–VI and III–V families. As a guide for use on other systems, we assessed the model on simulated data to show the expected accuracy level as a function of bandgap for new materials of interest. Finally, we use the model to predict a complete space of impurity charge-state transition levels in all zinc blende III–V and II–VI systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2015

WHEN DID ROUND DISK GALAXIES FORM?

When and how galaxy morphology, such as the disk and bulge seen in the present-day universe, emerged is still not clear. In the universe at z ≳ 2, galaxies with various morphologies are seen, and star-forming galaxies at z ∼ 2 show the intrinsic shape of bar-like structures. Then, when did the round disk structure form? Here we take a simple and straightforward approach to see the epoch when a round disk galaxy population emerged by constraining the intrinsic shape statistically based on the apparent axial ratio distribution of galaxies. We derived the distributions of the apparent axial ratios in the rest-frame optical light (∼5000 Å) of star-forming main-sequence galaxies at 2.5 > z > 1.4, 1.4 > z > 0.85, and 0.85 > z > 0.5, and found that their apparent axial ratios show peaky distributions at z ≳ 0.85, while a rather flat distribution at the lower redshift. By using a tri-axial model (A > B > C) for the intrinsic shape, we found that the best-fit models give the peaks of the B/A distribution of 0.81 ± 0.04, 0.84 ± 0.04, and 0.92 ± 0.05 at 2.5 > z > 1.4, 1.4 > z > 0.85, and 0.85 > z > 0.5, respectively. The last value is close to the local value of 0.95. Thickness (C/A) is ∼0.25 at all the redshifts and is close to the local value (0.21). The results indicate that the shape of the star-forming galaxies in the main sequence changes gradually, and that the round disk is established at around z ∼ 0.9. The establishment of the round disk may be due to the cessation of a violent interaction between galaxies or the growth of a bulge and/or a supermassive black hole residing at the center of a galaxy that dissolves the bar structure.

79 ASTRONOMY AND ASTROPHYSICS↗