SEARCH · Engineering Papers

Results for “Constants”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Effect of diffusion constant on the morphology of dendrite growth in lithium metal batteries

Lithium dendrites can lead to a short circuit and battery failure, and developing strategies for their suppression is of considerable importance. In this work, we study the growth of dendrites in a simple model system where the solvent is a continuum and the lithium ions are hard spheres that can deposit by sticking to existing spheres or the electrode surface. Using stochastic dynamics simulations, we investigate the effect of applied voltage and diffusion constant on the growth of dendrites. We find that the diffusion constant is the most significant factor, and the inhomogeneity of the electric field does not play a significant role. Additionally, the growth is most pronounced when the applied voltage and diffusion constant are both low. We observe a structural change from broccoli to cauliflower shape as the diffusion constant is increased. The simulations suggest that a control of electrolyte parameters that impact lithium diffusion might be an attractive route to controlling dendrite growth.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

A machine learning framework for elastic constants predictions in multi-principal element alloys

On the one hand, multi-principal element alloys (MPEAs) have created a paradigm shift in alloy design due to large compositional space, whereas on the other, they have presented enormous computational challenges for theory-based materials design, especially density functional theory (DFT), which is inherently computationally expensive even for traditional dilute alloys. In this paper, we present a machine learning framework, namely PREDICT (PRedict properties from Existing Database In Complex alloys Territory), that opens a pathway to predict elastic constants in large compositional space with little computational expense. The framework only relies on the DFT database of binary alloys and predicts Voigt–Reuss–Hill Young’s modulus, shear modulus, bulk modulus, elastic constants, and Poisson’s ratio in MPEAs. We show that the key descriptors of elastic constants are the A–B bond length and cohesive energy. The framework can predict elastic constants in hypothetical compositions as long as the constituent elements are present in the database, thereby enabling property exploration in multi-compositional systems. We illustrate predictions in a FCC Ni-Cu-Au-Pd-Pt system.

Linton, Nathan (ORCID:0000000315485613)↗

DOE OSTI2023

Temperature-dependent elastic constants of thorium dioxide probed using time-domain Brillouin scattering

Herein we report the adiabatic elastic constants of single-crystal thorium dioxide over a temperature range of 77–350 K. Time-domain Brillouin scattering, an all-optical, non-contact picosecond ultrasonic technique, is used to generate and detect coherent acoustic phonons that propagate in the bulk perpendicular to the surface of the crystal. These coherent acoustic lattice vibrations have been monitored in two hydrothermally grown single-crystal thorium dioxide samples along the (100) and (311) crystallographic directions. The three independent elastic constants of the cubic crystal (C 11 , C 12 , and C 44 ) are determined from the measured bulk acoustic velocities. The longitudinal wave along the (100) orientation provided a direct measurement of C 11 . Measurement of C 44 and C 12 was achieved by enhancing the intensity of quasi-shear mode in a (311) oriented crystal by adjusting the polarization angle relative to the crystal axes. We find the magnitude of softening of the three elastic constants to be ~2.5% over the measured temperature range. Good agreement is found between the measured elastic constants with previously reported values at room temperature, and between the measured temperature-dependent bulk modulus with calculated values. We find that semi-empirical models capturing lattice anharmonicity adequately reproduce the observed trend. We also determine the acoustic Grüneisen anharmonicity parameter from the experimentally derived temperature-dependent bulk modulus and previously reported temperature-dependent values of volumetric thermal expansion coefficient and heat capacity. This work presents measurements of the temperature-dependent elasticity in single-crystal thorium dioxide at cryogenic temperature and provides a basis for testing ab initio theoretical models and evaluating the impact of anharmonicity on thermophysical properties.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2023

Solvent Phase Optimizations Improve Correlations with Experimental Stability Constants for Aqueous Lanthanide Complexes

Stability constants provide insight into ion complexation in water. While computational studies have been shown to model the energy of the complexation successfully using a thermodynamic cycle approach, it does not extend to calculating the stability constants for 1:1 lanthanide to ligand complexes in solution. Using B3LYP and 6-31+G* Pople basis with small core effective core potential (ECP) on the lanthanum ion, and a solvent model based on the full solute electron density (SMD) solvation model we computed and compared with previously published stability constants of the ligands: acetate, acetohydroximate, acetylacetonate, methanoate, tropolonate, hydroxide, catecholate, malonate, oxalate, phthalate, and sulfate. The best R 2 values for the thermodynamic cycle can only be determined by separating the mono and divalent ions to achieve an R 2 value of 0.86 and 0.74 for mono and divalent ions, respectively. We show that by optimizing the lanthanide-ligand structures in implicit solvent, we achieve an improved correlation between experimental and computed stability constants of R 2 value of 0.89 for the combined mono and divalent ions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

The well-tempered cosmological constant: fugue in B$^\flat$

Zero point fluctuations of quantum fields should generate a large cosmological constant energy density in any spacetime. How then can we have anything other than de Sitter space without fine tuning? Well tempering—dynamical cancellation of the cosmological constant using degeneracy within the field equations— can replace a large cosmological constant with a much lower energy state. Here we give an explicit mechanism to obtain a Minkowski solution, replacing the cosmological constant with zero, and testing its attractor nature and persistence through a vacuum phase transition. We derive the general conditions that Horndeski scalar-tensor gravity must possess, and evolve in a fugue of functions, to deliver nothing and make the universe be flat.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Beyond leading twist: 𝜌 meson decay constants and distribution amplitudes in a self-consistent light-front quark model

In this study, we present a comprehensive analysis of decay constants and chiral-even and chiral-odd distribution amplitudes (DAs) up to twist 4 for the 𝜌 meson in the standard light-front quark model (LFQM) based on the Bakamjian-Thomas (BT) construction. For the 𝜌 meson, which possesses both longitudinal (ℎ=0) and transverse (ℎ =±1) polarizations, two types of decay constants, 𝑓$^{∥}_{𝜌}$ and 𝑓$^{⊥}_{𝜌}$, arises accordingly. We demonstrate that these decay constants can be self-consistently extracted from both local (𝑧𝜇 =0) and nonlocal (𝑧𝜇 ≠0) matrix elements ⟨0⁢|$\overline{𝑞}$(𝑧)⁢Γ⁢𝑞⁡(−𝑧)|⁢𝜌⁡(𝑃,ℎ)⟩, with Γ=(𝛾 𝜇 ,𝜎 𝜇⁢𝜈 ,𝛾 𝜇 ⁢𝛾 5 ,𝟏), in a manner independent of current components, polarizations, and reference frames. In particular, we emphasize the role of nonlocal matrix elements involving axial-vector and scalar currents, where mixing between 𝑓$^{∥}_{𝜌}$ and 𝑓$^{⊥}_{𝜌}$ occurs. We show that this mixing is consistently resolved through the BT construction, ensuring the proper extraction of these decay constants.

Electroweak interaction↗

DOE OSTI2024

Constant-Depth Preparation of Matrix Product States with Adaptive Quantum Circuits

Adaptive quantum circuits, which combine local unitary gates, midcircuit measurements, and feedforward operations, have recently emerged as a promising avenue for efficient state preparation, particularly on near-term quantum devices limited to shallow-depth circuits. Matrix product states (MPS) comprise a significant class of many-body entangled states, efficiently describing the ground states of one-dimensional gapped local Hamiltonians and finding applications in a number of recent quantum algorithms. Recently, it has been shown that the Affleck-Kennedy-Lieb-Tasaki state—a paradigmatic example of an MPS—can be exactly prepared with an adaptive quantum circuit of constant depth, an impossible feat with local unitary gates alone due to its nonzero correlation length [Smith , PRX Quantum 4, 020315 (2023)]. In this work, we broaden the scope of this approach and demonstrate that a diverse class of MPS can be exactly prepared using constant-depth adaptive quantum circuits, outperforming theoretically optimal preparation with unitary circuits. We show that this class includes short- and long-ranged entangled MPS, symmetry-protected topological (SPT) and symmetry-broken states, MPS with finite Abelian, non-Abelian, and continuous symmetries, resource states for MBQC, and families of states with tunable correlation length. Moreover, we illustrate the utility of our framework for designing constant-depth sampling protocols, such as for random MPS or for generating MPS in a particular SPT phase. We present sufficient conditions for particular MPS to be preparable in constant time, with global on-site symmetry playing a pivotal role. Altogether, this work demonstrates the immense promise of adaptive quantum circuits for efficiently preparing many-body entangled states and provides explicit algorithms that outperform known protocols to prepare an essential class of states. Published by the American Physical Society 2024

Smith, Kevin C. (ORCID:0000000223971518)↗

DOE OSTI2022

Constant-depth circuits for dynamic simulations of materials on quantum computers

Abstract Dynamic simulation of materials is a promising application for near-term quantum computers. Current algorithms for Hamiltonian simulation, however, produce circuits that grow in depth with increasing simulation time, limiting feasible simulations to short-time dynamics. Here, we present a method for generating circuits that are constant in depth with increasing simulation time for a specific subset of one-dimensional (1D) materials Hamiltonians, thereby enabling simulations out to arbitrarily long times. Furthermore, by removing the effective limit on the number of feasibly simulatable time-steps, the constant-depth circuits enable Trotter error to be made negligibly small by allowing simulations to be broken into arbitrarily many time-steps. For an N -spin system, the constant-depth circuit contains only $\mathcal {O}(N^{2})$ O ( N 2 ) CNOT gates. Such compact circuits enable us to successfully execute long-time dynamic simulation of ubiquitous models, such as the transverse field Ising and XY models, on current quantum hardware for systems of up to 5 qubits without the need for complex error mitigation techniques. Aside from enabling long-time dynamic simulations with minimal Trotter error for a specific subset of 1D Hamiltonians, our constant-depth circuits can advance materials simulations on quantum computers more broadly in a number of indirect ways.

Materials simulation↗

DOE OSTI2022

Identification of Distribution of Elastic Constants of SiC/SiC Braided Tubes

There has been considerable interest in chemical vapor infiltration (CVI) manufactured silicon carbide fiber and silicon carbide matrix (SiCf/SiCm) composite tubes due to their superior mechanical properties. There are many parameters in the manufacturing of SiCf/SiCm composite tubes such as elastic constants of constituent SiC fiber and SiC matrix, braiding angle, porosity, and others whose values vary over a range. A sensitivity analysis will be carried out to determine the dominant parameters among the abovementioned ones. We intend to determine the variability in elastic constants of SiCf/SiCm composite tubes considering the variability in dominant parameters. A chosen number of combinations of dominant variables will be sampled through Latin hypercube sampling (LHS) method. The elastic constants of composite tube namely Young's moduli in circumferential and longitudinal direction, Poisson's ratio and shear modulus will be calculated at sampled points through finite element (FE) analysis. A polynomial response surface (PRS) will be built for each of the elastic constants to be used a surrogate to FE analysis. once the predictive accuracy of PRS is verified at the validation set of data points, the PRS will be invoked 105 times in Monte Carlo simulations (MCS). The variability can be calculated based on MCS-generated data which will be presented at the conference.

42 ENGINEERING↗

DOE OSTI2022

Inefficient intramolecular vibrational energy redistribution for the H + HO2 reaction and negative internal energy dependence for its rate constant

Quasiclassical trajectories (QCT) and newly constructed global potential energy surfaces are used to compute thermal and nonthermal rate constants for the H + HO 2 reaction. The thermal QCTs rate constants are up to 50% smaller than transition state theory (TST) rate constants based on the same level of electronic structure theory. This reduction is demonstrated to result from inefficient intramolecular vibrational energy redistribution (IVR) in the transient H 2 O 2 well, with a significant fraction of trajectories that reach the H 2 O 2 well promptly dissociating back to reactants instead of via the heavily statistically favored 2OH channel. The nonstatistical reduction factor, κ IVR , that quantifies this effect is shown to increase in importance with temperature, with κ IVR = 0.81 at 300 K and 0.47 at 2500 K. Finally, we show that inefficient IVR causes H + HO 2 rate constants mediated by H 2 O 2 to depend inversely on the initial vibrational excitation of HO 2 .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Preliminary parameter studies to form pellets of hard materials for optical constants derived via single-angle reflectance spectroscopy

Infrared reflectance spectra can be influenced by many factors: the substrate and the thickness of the layer for liquids and the surface micromorphology, the form (powder, crystal) and the particle size for solids. All these parameters can play a role in the appearance of the measured spectrum. To avoid collecting multitudes of spectra to cover all such scenarios, the optical constants n and k, which are intrinsic properties of a material, can instead be used to model the reflectance spectrum. For solids, two techniques are often used to derive optical constants: ellipsometry and single-angle reflectance spectroscopy. For both methods, best results are usually obtained from single crystals. We have recently demonstrated for ammonium sulfate (a relatively soft material) that by optimizing certain conditions, high quality pellets with specularly reflective surfaces can be used in lieu of crystals. This was confirmed by the excellent agreement between the optical constants derived by these two methods. This work focuses on the possible extension of these methods to harder materials, starting with sodium sulfate. The first goal is to see if high quality pellet surfaces can be obtained as for ammonium sulfate. The reflectance values and the associated optical constants can also be obtained.

Diaz, Emmanuela↗

DOE OSTI2023

Challenges for the statistical gravitational-wave method to measure the Hubble constant

Gravitational waves (GW) can be employed as standard sirens that will soon measure the Hubble constant with sufficient precision to weigh in on the ~5σ Hubble tension. Most GW sources will have no identified electromagnetic counterpart, leading to uncertainty in the redshift of the source, and in turn a degeneracy between host galaxy distance, redshift, and H 0 . In the case where no electromagnetic counterparts are identified, it has been proposed that a statistical canvassing of candidate GW hosts, found in a large galaxy survey for example, can be used to accurately constrain the Hubble constant. We study and simulate this "galaxy voting" method to compute H 0 . Here, we find that the Hubble constant posterior is in general biased relative to the true value even when making optimistic assumptions about the statistical properties of the sample. Using the MICECAT light-cone catalog, we find that the bias in the H0 posteriors depends on the realization of the underlying galaxy sample and the precision of the GW source distance measurement.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Modeling Wavelength Dependent Mid-infrared (5.5–25 μm) Optical Constants of Silicate Glasses: A Genetic Algorithm Approach

Wavelength-dependent mid-infrared (400–1,800 cm -1 ; 5.56–25 μm) optical constants (real and imaginary indices of refraction; n and k) are determined using reflectance spectra at a spectral sampling of 4 cm -1 for several silicate glasses of varying SiO 2 wt% which include (a) basaltic volcanic glass from Kilauea, Hawaii, (b) synthetic andesite, (c) two synthetic dacites, (d) obsidian volcanic glass from Mount Lassen, California, (e) synthetic rhyodacite, and (f) rhyolitic volcanic glass from Mexico. Because glasses are optically isotropic, no specific orientation was required for spectral measurements, and polished glass samples were measured at random orientations using micro-FTIR spectrometer. Lorentz-Lorenz dispersion theory and Fresnel reflectance model for high symmetry materials were used to model the optical constants by optimizing oscillator parameters in modeled spectra to match laboratory spectra. A genetic algorithm (GA) approach automatically finds the natural oscillators and their parameters for each spectrum, and then uses these parameters in a non-linear least squares optimization routine. The study compared spectral parameters such as Christiansen feature, reststrahlen bands (RBs), and the peaks centered around 860–1,100 (peak 1) and 400–480 cm -1 (peak 2) of n and k to their respective SiO 2 wt% of the glasses. CF, RBs, peak 1 of n and k, and peak 2 of n shift to higher wavenumbers with increased SiO 2 wt%, whereas peak 2 of k shifts to lower wavenumbers with increased SiO 2 wt%. Derived optical constants of these glasses will improve quantitative abundance mapping of volcanic materials on the surfaces of silicate targets in the Solar System.

58 GEOSCIENCES↗

DOE OSTI2023

Neutral gas pressure dependence of ion–ion mutual neutralization rate constants using Landau–Zener theory coupled with trajectory simulations

In this computational study, we describe a self-consistent trajectory simulation approach to capture the effect of neutral gas pressure on ion–ion mutual neutralization (MN) reactions. The electron transfer probability estimated using Landau–Zener (LZ) transition state theory is incorporated into classical trajectory simulations to elicit predictions of MN cross sections in vacuum and rate constants at finite neutral gas pressures. Electronic structure calculations with multireference configuration interaction and large correlation consistent basis sets are used to derive inputs to the LZ theory. The key advance of our trajectory simulation approach is the inclusion of the effect of ion-neutral interactions on MN using a Langevin representation of the effect of background gas on ion transport. For H + – H – and Li + – H(D) – , our approach quantitatively agrees with measured speed-dependent cross sections for up to ~10 5 m/s. For the ion pair Ne + – Cl – , our predictions of the MN rate constant at ~1 Torr are a factor of ~2 to 3 higher than the experimentally measured value. Similarly, for Xe + – F – in the pressure range of ~20 000–80 000 Pa, our predictions of the MN rate constant are ~20% lower but are in excellent qualitative agreement with experimental data. Here, the paradigm of using trajectory simulations to self-consistently capture the effect of gas pressure on MN reactions advanced here provides avenues for the inclusion of additional nonclassical effects in future work.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

How to calculate the rate constants for nonradiative transitions between the M S components of spin multiplets?

Predicting the rates of spin-dependent processes characterized by nonradiative transitions between electronic states with different spin multiplicities is important for understanding the mechanisms of many photochemical and catalytic reactions. To calculate these rates, it is necessary to define the spin state representation and the couplings between these states that drives the interstate transitions. In this work, we describe three different approaches to calculating the spin-orbit coupling (SOC), transition probabilities, and rate constants between the MS components of the electronic states with different spin multiplicities. We implemented these approaches in our nonadiabatic statistical theory (NAST) software package, which predicts the transition probabilities and rate constants of spin-dependent processes using information obtained from electronic structure calculations. Here, we discuss the advantage and drawbacks of each approach and, as an example, calculate the rate constants for transitions between the spin states in the active site model of the protein rubredoxin.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Relativistic configuration-interaction and many-body-perturbation-theory calculations of U i hyperfine constants

Neutral uranium (U i) is a very difficult atom for theoretical calculations due to a large number (six) of valence electrons, strong valence-valence and valence-core correlations, high density of states, and relativistic effects. Configuration-interaction many-body perturbation theory (CI-MBPT) can efficiently treat valence-core correlations and relativistic effects, but because the formalism was developed for a Dirac-Hartree-Fock (DHF) starting potential that does not contain valence electrons, quite large CI space is needed to compensate for the + 6 charge of such a potential. Much more efficient is the relativistic configuration-interaction (RCI) approach, which uses a relatively accurate starting DHF potential that includes some valence electrons to make the valence-electron Hamiltonian diagonally dominated for some states. Here we report calculations of U i hyperfine constants of several low-energy states using the RCI method with the starting potential that includes four f valence electrons. With this starting potential, it is possible to use the single-configuration approximation or small basis sets to obtain quite accurate results for hyperfine-structure constants. In fact, by scaling the nuclear magnetic moment, the agreement for five levels was within 5% and a new magnetic moment can be recommended, 0.43(2). When two states are mixed, it is difficult to predict the exact mixing fractions, so in this case we proposed a method in which the mixing fraction is found by optimization of one parameter in the MBPT correction to obtain correct g factors. The resulting hyperfine constants are significantly improved. Finally, the CI-MBPT approach was also tested, with some success for low-energy states, with the limited basis set. The methods investigated here can be further developed to include more extensive data sets to improve accuracy and can be applied to other atoms and for calculations of other properties, for example, relevant to fundamental symmetry tests.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2020

Improved Infrared Optical Constants from Pressed Pellets: II. Ellipsometric n and k Values for Ammonium Sulfate with Variability Analysis

Infrared spectral imaging is facilitated via the comparison of spectra recorded in situ to a databank of actual or synthetic infrared reflectance spectra. It has recently been shown that reference spectra corresponding to the many different morphological forms of the same chemical can be generated synthetically using the imaginary, k, and real, n, components of the complex refractive index, n~ = n + ik. One method to obtain the n and k vectors is infrared ellipsometry, which measures the changes in amplitude, tan ?, and phase, ?, of polarized light reflected from the sample both as a function of wavenumber and angle of incidence. The method requires specularly reflected light so best results are usually obtained with polished planar samples of large surface area. Due to the difficulties of obtaining such samples, however, we investigate the possibility of pressing powders of neat materials and obtaining the corresponding optical constants from the pellets. In this paper, variability in the sample pellet and preparation method is investigated, as is variability in the fitting procedure for the derived optical constants. The n/k vectors are derived from the measured ellipsometric parameters, tan ? and ?, as they are fit by an oscillator model which yield n(?) and k(?) vectors as a function of wavenumber, ?. Construction of the oscillator model is not automatic and depends on significant input from the analyst as well as the sample’s physical characteristics. For the pellet pressing, the experimental variability was found to be minimized for size-selected powdered samples as gauged by the minimal variance in ? and ? for three different pellets; similarly, the analytical precision for multiple measurements of the same pellet was also quite good, suggesting that a pressed pellet is a viable sample preparation method. The greatest variability came in the analytic fitting procedure with differences up to 20% arising from four different fits to the same data set. The final ellipsometric n/k data are compared to literature values obtained from crystalline (NH4)2SO4 samples as well as single-angle reflectance measurements that also used pressed pellets. Comparison with the previous literature values shows generally good agreement although higher k-values are observed for the independent sets of data derived from pressed pellets. These data are suggested as an improved set of optical constants for (NH4)2SO4.

optical constants, ellipsometry, ammonium sulfate,↗

DOE OSTI2020

Infrared Optical Constants from Pressed Pellets of Powders: I. Improved n and k Values of (NH4)(2)SO(4)from Single-Angle Reflectance

In combination with other parameters, the real, n(?), and imaginary, k(?), components of the complex refractive index, = n + ik, can be used to simulate the optical properties of a material in different forms, e.g. its infrared spectra. But obtaining reliable values of the optical constants n/k for solid materials is challenging due to the lack of optical quality specimens, usually crystals, large enough to measure. An alternative to crystals is to press the powder into a uniform disk. We have produced pellets from ammonium sulfate, (NH4)2SO4, powder and derived the pellets’ n and k values via single-angle reflectance using a specular reflectance device in combination with a Fourier transform infrared spectrometer. The single-angle technique measures amplitude of light reflected from the material as a function of wavelength over a wide spectral domain; the optical constants are determined from the reflectance data using the Kramers-Kronig relationship. We investigate several parameters associated with the pellets and pellet formation and their effects upon delivering the most reliable n/k values. Parameters studied include pellet diameter, mass and density (void space), drying, grinding, sieving and particle size in the pellet formation, as well as pressing pressure and duration. Of these parameters, using size-selected mixtures of dried, small (< 100 µm) particles and pressing at = 10 tons for at least 30 minutes were found key to forming highly reflective samples. Comparison to two sets of previous literature n(?) and k(?) values obtained from crystalline (NH4)2SO4 both as crystal reflectance as well as extinction spectra of aerosols measured in a flow tube show reasonable agreement, but suggests the present values, as confirmed from two independent techniques, represent a substantial improvement for n/k values for (NH4)2SO4, also demonstrating promise to measure the optical constants of other materials.

Optical constants, index of refraction, single-ang↗