Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “genetic algorithms”

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 775 records · Page 43

Machine learning a bond order potential model to study thermal transport in WSe 2 nanostructures

Here, nanostructures of transition metal di-chalcogenides (TMDCs) exhibit exotic thermal, chemical and electronic properties, enabling diverse applications from thermoelectrics and catalysis to nanoelectronics. The thermal properties of these nanoscale TMDCs are of particular interest for thermoelectric applications. Thermal transport studies on nanotubes and nanoribbons remain intractable to first principles calculations whereas existing classical molecular models treat the two chalcogen layers in a monolayer with different atom types; this imposes serious limitations in studying multi-layered TMDCs and dynamical phenomena such as nucleation and growth. Here, we overcome these limitations using machine learning (ML) and introduce a bond order potential (BOP) trained against first principles training data to capture the structure, dynamics, and thermal transport properties of a model TMDC such as WSe 2 . The training is performed using a hierarchical objective genetic algorithm workflow to accurately describe the energetics, as well as thermal and mechanical properties of a free-standing sheet. As a representative case study, we perform molecular dynamics simulations using the ML-BOP model to study the structure and temperature-dependent thermal conductivity of WSe 2 tubes and ribbons of different chiralities. We observe slightly higher thermal conductivities along the armchair direction than zigzag for WSe 2 monolayers but the opposite effect for nanotubes, especially of smaller diameters. We trace the origin of these differences to the anisotropy in thermal transport and the restricted momentum selection rules for phonon–phonon Umpklapp scattering. The developed ML-BOP model is of broad interest and will facilitate studies on nucleation and growth of low dimensional WSe 2 structures as well as their transport properties for thermoelectric and thermal management applications.

36 MATERIALS SCIENCE↗

Combining artificial intelligence and physics-based modeling to directly assess atomic site stabilities: from sub-nanometer clusters to extended surfaces

The performance of functional materials is dictated by chemical and structural properties of individual atomic sites. In catalysts, for instance, the thermodynamic stability of constituting atomic sites is a key descriptor from which more complex properties, such as molecular adsorption energies and reaction rates, can be derived. In this study, we present a widely applicable machine learning (ML) approach to instantaneously compute the stability of individual atomic sites in structurally and electronically complex nano-materials. Conventionally, we determine such site stabilities using computationally intensive first-principles calculations. With our approach, we predict the stability of atomic sites in sub-nanometer metal clusters of 3–55 atoms with mean absolute errors in the range of 0.11–0.14 eV. To extract physical insights from the ML model, we introduce a genetic algorithm (GA) for feature selection. This algorithm distills the key structural and chemical properties governing the stability of atomic sites in size-selected nanoparticles, allowing for physical interpretability of the models and revealing structure–property relationships. The results of the GA are generally model and materials specific. In the limit of large nanoparticles, the GA identifies features consistent with physics-based models for metal–metal interactions. By combining the ML model with the physics-based model, we predict atomic site stabilities in real time for structures ranging from sub-nanometer metal clusters (3–55 atom) to larger nanoparticles (147 to 309 atoms) to extended surfaces using a physically interpretable framework. Finally, we present a proof of principle showcasing how our approach can determine stable and active nanocatalysts across a generic materials space of structure and composition.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Thermodynamic stability versus kinetic accessibility: Pareto fronts for programmable self-assembly

A challenge in designing self-assembling building blocks is to ensure the target state is both thermodynamically stable and kinetically accessible. These two objectives are known to be typically in competition, but it is not known how to simultaneously optimize them. We consider this problem through the lens of multi-objective optimization theory: we develop a genetic algorithm to compute the Pareto fronts characterizing the tradeoff between equilibrium probability and folding rate, for a model system of small polymers of colloids with tunable short-ranged interaction energies. We use a coarse-grained model for the particles' dynamics that allows us to efficiently search over parameters, for systems small enough to be enumerated. For most target states there is a tradeoff when the number of types of particles is small, with medium-weak bonds favouring fast folding, and strong bonds favouring high equilibrium probability. Additionally, the tradeoff disappears when the number of particle types reaches a value m *, that is usually much less than the total number of particles. This general approach of computing Pareto fronts allows one to identify the minimum number of design parameters to avoid a thermodynamic–kinetic tradeoff. However, we argue, by contrasting our coarse-grained model's predictions with those of Brownian dynamics simulations, that particles with short-ranged isotropic interactions should generically have a tradeoff, and avoiding it in larger systems will require orientation-dependent interactions.

36 MATERIALS SCIENCE↗

Experimentally informed structure optimization of amorphous TiO 2 films grown by atomic layer deposition

Amorphous titanium dioxide TiO 2 (a-TiO 2 ) has been widely studied, particularly as a protective coating layer on semiconductors to prevent corrosion and promote electron–hole conduction in photoelectrochemical reactions. The stability and longevity of a-TiO 2 is strongly affected by the thickness and structural heterogeneity, implying that understanding the structure properties of a-TiO 2 is crucial for improving the performance. This study characterized the structural and electronic properties of a-TiO 2 thin films (~17 nm) grown on Si by atomic layer deposition (ALD). Fluctuation spectra V(k) and angular correlation functions were determined with 4-dimensional scanning transmission electron microscopy (4D-STEM), which revealed the distinctive medium-range ordering in the a-TiO 2 film. A realistic atomic model of a-TiO 2 was established guided by the medium-range ordering and the previously reported short-range ordering of a-TiO 2 film, as well as the interatomic potential. The structure was optimized by the StructOpt code using a genetic algorithm that simultaneously minimizes energy and maximizes the match to experimental short- and medium-range ordering. The StructOpt a-TiO 2 model presents improved agreements with the medium-range ordering and the k-space location of the dominant 2-fold angular correlations compared with a traditional melt-quenched model. The electronic structure of the StructOpt a-TiO 2 model was studied by ab initio calculations and compared to the crystalline phases and experimental results. Finally, this work uncovered the medium-range ordering in a-TiO 2 thin films and provided a realistic a-TiO 2 structure model for further investigation of structure–property relationships and materials design. In addition, the improved multi-objective optimization package StructOpt was provided for structure determination of complex materials guided by experiments and simulations.

36 MATERIALS SCIENCE↗

Finding the Pareto front for high-entropy-alloy catalysts

Finding catalysts that have both high activity and high stability presents a long-standing challenge. Since optimizing activity and stability are conflicting objectives, the best one can do is find the Pareto front that yields optimal tradeoffs between these features. On the Pareto front, there is a trade-off where a portion of catalytic activity must be sacrificed to gain further stability and vice versa . Here, we provide a method to optimize the front by designing a multi-objective genetic algorithm that combines machine learning, graph neural network calculations, and density functional calculations. The application considered is the oxygen evolution reaction catalyzed by high-entropy alloys. We find that the Pareto front generally contains alloys with diverse elements, but that enhancing stability inevitably inflicts a toll on activity. We compare the general conclusions of our work to a survey of 545 experiments.

Zhang, Chengyi [Univ. of Auckland (New Zealand)]↗

Identification of an advanced constitutive model of Magnesium alloy AZ31B

The main aim of this paper is to study the flow behavior of the AZ31B magnesium alloy by means of tensile tests performed in extended ranges of temperature and strain rates. The flow stress-strain curves analyzed by power law type constitutive equation can only fit well with experimental curves at the work-hardening stage. A new mathematical model is studied to describe the softening behavior of material based on tensile experiments. The relative parameters are obtained by fitting the equation with the experimental data. The genetic algorithm has been used to obtain the global optimal fitting parameters. The comparison between the fitted and experimental data proves the effectiveness of the model. The results indicate that this model leads to a better simulation of the flow stress during the softening stage than that of the power law equation. Based on this model, the deep drawing process has been simulated with the commercial finite element code FORGE registered. The punch load and thickness distribution of AZ31 sheet have been studied. The study of the results is helpful to the application of the stamping technology for the magnesium alloy sheet.

36 MATERIALS SCIENCE↗

Atomic structure and magnetic properties of Fe 1– x Co x alloys

Using genetic algorithm with first-principle calculations, we searched for low-energy crystal structures of Fe 1−x Co x alloys. We found that Fe 1−x Co x alloys are highly configurationally degenerate with many additional off-stoichiometric stable structures to the well-known B 2 structure. The average magnetic moment of Fe atom increases with concentration of Co in the alloy, while that of Co atom is almost constant, which are consistent with experiments and earlier studies. The magnetic moment of Fe atom is strongly dependent on the number of Co nearest neighbor and it increases with this number.

36 MATERIALS SCIENCE↗

Time and space resolved measurement of the electron temperature, mass density and ionization state in the ablation plasma between two exploding Al wires

We have determined the properties of plasma around and between two exploding wires using high-resolution x-ray absorption spectroscopy. Plasma densities and temperatures ranging from Greater-Than-Or-Equivalent-To 0.1g/cm{sup 3} and a few eV to less than 0.01 g/cm{sup 3} and 30 eV have been measured in experiments at Cornell University with two 40 {mu}m aluminum (Al) wires spaced 1 mm apart driven by {approx}150 kA peak current pulses with 100 ns rise time. The wire plasma was backlit by the 1.4-1.6 keV continuum radiation produced by a Mo wire X-pinch. The spectrometer employed two spherically bent quartz crystals to record the absorption and backlighter spectra simultaneously. The transition between the dense Al wire core and the coronal plasma is seen as a transition from cold K-edge absorption to Mg-, Na-, and finally Ne-like absorption at the boundary. In the plasma that accumulates between the wires, ionization states up to C-Like Al are observed. The spectrometer geometry and {approx}2{mu}m X-pinch source size provide 0.3 eV spectral resolution and 20 {mu}m spatial resolution enabling us to see 1s{yields} 2p satellite transitions as separate lines as well as O-, F-, and Ne-like 1s{yields} 3p transitions that have not been seen before. A step wedge was used to calibrate the transmission, enabling density to be measured within a factor of two and temperature to be measured within {+-}25%. A genetic algorithm was developed to fit synthetic spectra calculated using the collisional-radiative code SCRAM to the experimental spectra. In order to obtain agreement it was necessary to assume multiple plasma regions with variable thicknesses, thereby allowing the inferred plasma conditions to vary along the absorption path.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Optimal design and fabrication of three-dimensional calibration specimens for scanning probe microscopy

Micro-/nano-scale roughness specimens are highly demanded to synthetically calibrate the scanning probe microscopy (SPM) instrument. In this study, three-dimensional (3D) specimens with controllable main surface evaluation parameters were designed. In order to improve the design accuracy, the genetic algorithm was introduced into the conventional digital filter method. A primary 3D calibration specimen with the dimension of 10 {mu}m x 10 {mu}m was fabricated by electron beam lithography. Atomic force microscopy characterizations demonstrated that the statistical and spectral parameters of the fabricated specimen match well with the designed values. Such a kind of 3D specimens has the potential to calibrate the SPM for applications in quantitative surface evaluations.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

O{sub 2}(a{sup 1}{Delta}{sub g}) + Mg, Fe, and Ca: Experimental kinetics and formulation of a weak collision, multiwell master equation with spin-hopping

The first excited electronic state of molecular oxygen, O{sub 2}(a{sup 1}{Delta}{sub g}), is formed in the upper atmosphere by the photolysis of O{sub 3}. Its lifetime is over 70 min above 75 km, so that during the day its concentration is about 30 times greater than that of O{sub 3}. In order to explore its potential reactivity with atmospheric constituents produced by meteoric ablation, the reactions of Mg, Fe, and Ca with O{sub 2}(a) were studied in a fast flow tube, where the metal atoms were produced either by thermal evaporation (Ca and Mg) or by pulsed laser ablation of a metal target (Fe), and detected by laser induced fluorescence spectroscopy. O{sub 2}(a) was produced by bubbling a flow of Cl{sub 2} through chilled alkaline H{sub 2}O{sub 2}, and its absolute concentration determined from its optical emission at 1270 nm (O{sub 2}(a{sup 1}{Delta}{sub g}- X{sup 3}{Sigma}{sub g}{sup -}). The following results were obtained at 296 K: k(Mg + O{sub 2}(a) + N{sub 2}{yields} MgO{sub 2}+ N{sub 2}) = (1.8 {+-} 0.2) Multiplication-Sign 10{sup -30} cm{sup 6} molecule{sup -2} s{sup -1}; k(Fe + O{sub 2}(a) {yields} FeO + O) = (1.1 {+-} 0.1) Multiplication-Sign 10{sup -13} cm{sup 3} molecule{sup -1} s{sup -1}; k(Ca + O{sub 2}(a) + N{sub 2}{yields} CaO{sub 2}+ N{sub 2}) = (2.9 {+-} 0.2) Multiplication-Sign 10{sup -28} cm{sup 6} molecule{sup -2} s{sup -1}; and k(Ca + O{sub 2}(a) {yields} CaO + O) = (2.7 {+-} 1.0) Multiplication-Sign 10{sup -12} cm{sup 3} molecule{sup -1} s{sup -1}. The total uncertainty in these rate coefficients, which mostly arises from the systematic uncertainty in the O{sub 2}(a) concentration, is estimated to be {+-}40%. Mg + O{sub 2}(a) occurs exclusively by association on the singlet surface, producing MgO{sub 2}({sup 1}A{sub 1}), with a pressure dependent rate coefficient. Fe + O{sub 2}(a), on the other hand, shows pressure independent kinetics. FeO + O is produced with a probability of only {approx}0.1%. There is no evidence for an association complex, suggesting that this reaction proceeds mostly by near-resonant electronic energy transfer to Fe(a{sup 5}F) + O{sub 2}(X). The reaction of Ca + O{sub 2}(a) occurs in an intermediate regime with two competing pressure dependent channels: (1) a recombination to produce CaO{sub 2}({sup 1}A{sub 1}), and (2) a singlet/triplet non-adiabatic hopping channel leading to CaO + O({sup 3}P). In order to interpret the Ca + O{sub 2}(a) results, we utilized density functional theory along with multireference and explicitly correlated CCSD(T)-F12 electronic structure calculations to examine the lowest lying singlet and triplet surfaces. In addition to mapping stationary points, we used a genetic algorithm to locate minimum energy crossing points between the two surfaces. Simulations of the Ca + O{sub 2}(a) kinetics were then carried out using a combination of both standard and non-adiabatic Rice-Ramsperger-Kassel-Marcus (RRKM) theory implemented within a weak collision, multiwell master equation model. In terms of atmospheric significance, only in the case of Ca does reaction with O{sub 2}(a) compete with O{sub 3} during the daytime between 85 and 110 km.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Calibration and analysis of spatially resolved x-ray absorption spectra from a nonuniform plasma

We report here the calibration and analysis techniques used to obtain spatially resolved density and temperature measurements of a pair of imploding aluminum wires from x-ray absorption spectra. A step wedge is used to measure backlighter fluence at the film, allowing transmission through the sample to be measured with an accuracy of {+-}14% or better. A genetic algorithm is used to search the allowed plasma parameter space and fit synthetic spectra with 20 {mu}m spatial resolution to the measured spectra, taking into account that the object plasma nonuniformity must be physically reasonable. The inferred plasma conditions must be allowed to vary along the absorption path in order to obtain a fit to the spectral data. The temperature is estimated to be accurate to within {+-}25% and the density to within a factor of two. This information is used to construct two-dimensional maps of the density and temperature of the object plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Plasma density measurements in tungsten wire-array Z-pinches

Measurements of the plasma density profile near the exploding wires in 1 MA tungsten (W) wire-array Z-pinches have been made using calibrated x-ray absorption. As many as 5 x-ray images per pulse were obtained between 65 and 160 ns after the start of the 100 ns rise time current pulse. Measured W ion densities range from above 10{sup 19}/cm{sup 3} close to the wire to {approx}10{sup 17}/cm{sup 3} about 1 mm away from the wire in the plasma stream. After accurate geometrical registration of the individual wires in each successive image in a pulse using the Genetic Algorithm, the temporal evolution of the axial modulation wavelength distribution of the ablation rate from the wires in each array and the global mass-ablation rate as a function of time are presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Investigation of a polychromatic tomography method for the extraction of the three-dimensional spatial structure of implosion core plasmas

A spectroscopic method has been developed to extract the three-dimentional spatial structure (i.e., spatial distribution of electron temperature and density) of inertial confinement fusion implosion cores based on the analysis of space-resolved spectra from a tracer element recorded along three quasi-orthogonal lines of sight. We discuss a spectral model that computes space-resolved spectra for a given spatial structure. This model is then combined with a multi-objective search and optimization technique driven by a Pareto genetic algorithm to perform the inversion and to extract the spatial structure of the implosion core from a simultaneous and self-consistent analysis of a set of space-resolved spectra. This method is investigated with a series of synthetic data test cases to explore its reliability, requirements, and limitations. We have found a constraint parameter P{sub const} such that the method is robust and the extracted spatial structure is reliable when P{sub const}> 1. The idea of polychromatic tomography is general and has potential to extract the spatial structure of other laboratory high energy-density plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Approximate photochemical dynamics of azobenzene with reactive force fields

We have fitted reactive force fields of the ReaxFF type to the ground and first excited electronic states of azobenzene, using global parameter optimization by genetic algorithms. Upon coupling with a simple energy-gap transition probability model, this setup allows for completely force-field-based simulations of photochemical cis→trans- and trans→cis-isomerizations of azobenzene, with qualitatively acceptable quantum yields. This paves the way towards large-scale dynamics simulations of molecular machines, including bond breaking and formation (via the reactive force field) as well as photochemical engines (presented in this work)

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Optical and electronic properties of mixed Ag-Au tetramer cations

We present experimental and theoretical studies of the optical response of mixed Ag {sub n} Au {sub 4−n}{sup +} (n=1–3) clusters in the photon energy range ℏω = 1.9–3.5 eV. Absorption spectra are recorded by a newly built longitudinal molecular beam depletion spectroscopy apparatus providing lower limits to absolute photodissociation cross sections. The experimental data are compared to optical response calculations in the framework of long-range corrected time-dependent density functional theory with initial cluster geometries obtained by the unbiased Birmingham Cluster Genetic Algorithm coupled with density functional theory. Experiments and excited state calculations shed light on the structural and electronic properties of the mixed Ag-Au tetramer cations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Experimental determination of the complex stiffness tensor and Euler angles in anisotropic media using ultrasonic waves

The aim of this work is to measure the complex elastic tensor and Euler angles in very complex anisotropic media like austenitic steel welds, by inverse problem resolution from experimental data. The obtained experimental characteristics of the anisotropic material will be injected in a FE code developed by EDF enabling the simulation of an actual ultrasonic NDE of welds. The present work aims to provide reliable input data to the 3D future development of the code. In particular, this complex elastic tensor will allow to predict by modeling beam skewing ant attenuation in an austenitic weld. The investigation of such anisotropic media is very complex because of the directional dependency of the elastic stiffness tensor. Then we will discuss the use of a hybrid genetic algorithm to overcome this difficulty. The identification method is based on waveforms spectra reconstruction associated to a physical model describing wave propagation in plates, during underwater measurements. The entire procedure is qualified and validated using simulated data. Moreover, a comparison of the estimated elastic coefficients with literature values and ultrasonic measurements obtained in transmission is also given, at the end of the paper.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗