Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Functional coupling”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Revisiting the electrochemical reduction of CO 2 on Au 25 (SR)$\overline{_{18}}$ nanocluster

Experiments show that Au25(SR)$\overline{_{18}}$ is a potential electrocatalyst for the electrochemical reduction of CO 2 , but elucidation of the reaction mechanisms remain challenging experimentally. We evaluate key steps in the active site formation, CO 2 reduction, and H 2 evolution on the nanocluster employing density functional theory coupled with an explicit solvent model of the electrochemical interface to calculate activation barriers. The predicted preferred pathway for activating the nanocluster involves dethiolation, resulting in an exposed Au site that is both active and selective for CO 2 reduction. Finally, the alternative S site is found to be kinetically prohibitive and does not facilitate CO 2 reduction.

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

Enabling complete multichannel nonadiabatic dynamics: A global representation of the two-channel coupled, 1,2 1 A and 1 3 A states of NH 3 using neural networks

Global coupled three-state two-channel potential energy and property/interaction (dipole and spin–orbit coupling) surfaces for the dissociation of NH 3 (Ã) into NH + H 2 and NH 2 + H are reported. The permutational invariant polynomial-neural network approach is used to simultaneously fit and diabatize the electronic Hamiltonian by fitting the energies, energy gradients, and derivative couplings of the two coupled lowest-lying singlet states as well as fitting the energy and energy gradients of the lowest-lying triplet state. The key issue in fitting property matrix elements in the diabatic basis is that the diabatic surfaces must be smooth, that is, the diabatization must remove spikes in the original adiabatic property surfaces attributable to the switch of electronic wavefunctions at the conical intersection seam. Here, we employ the fit potential energy matrix to transform properties in the adiabatic representation to a quasi-diabatic representation and remove the discontinuity near the conical intersection seam. The property matrix elements can then be fit with smooth neural network functions. The coupled potential energy surfaces along with the dipole and spin–orbit coupling surfaces will enable more accurate and complete treatment of optical transitions, as well as nonadiabatic internal conversion and intersystem crossing.

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

Feasibility of a Novel Density Functional Method Outside the Kohn-Sham Framework for Modeling Global Potential Energy Surfaces of Molecular Chemical Reactions (Final Technical Report)

The project aimed towards the construction of computational methods capable of modeling the global potential energy surface (PES) of small molecules, molecular ions, and radicals — including the parts of the PES which correspond to chemical reactions, and the reaction paths connecting reactants, intermediates, and reaction products. This goal may seem humble at first glance, but for chemical systems with more than about six atoms in total, at the time the project was proposed, there were no established theoretical methods capable of simulating such systems reliably, not even for small molecules in the gas phase which are electronically benign. This restriction severely hampers our ability to model and control chemical processes under harsh conditions. With the goal of constructing a method capable of modeling such chemical systems, we proposed to pursue a novel approach towards a Multi-Configuration (MC) DFT outside the traditional frameworks of Kohn-Sham theory and other methods of coupling wave function theory with DFT. Rather than being a complete active space (CAS) method, the proposed DOCI-DFT would employ a special restricted form of the active space wave function, called Doubly-Occupied CI (DOCI)—this wave function form is sufficient to describe not only heterolytic, but also homolytic bond dissociation processes at the zeroth order (i.e., as active space wave function). There are no other standard mean field methods which can do so. In the original proposal, we also outlined strong formal and practical arguments speaking for this method. The project proposal was accepted by the Department of Energy and provided two years of seed funding for one graduate student, as well as two years of two weeks PI summer-salary for the PI. However, ultimately the project could not be effectively pursued due to severe interferences outside the context scientific problems and was cancelled. Apart from two weeks of PI salary (and associated fringe & overhead costs issues by the performing institution), no costs were charged against the award, with the entire rest of funding returned to the Department of Energy. This report reiterates the primary background information regarding the project, its original goals, and outlines the preliminary work performed during the two weeks of DOE funding charged.

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

Ba 6 (Cu x Z y )Sn 4 S 16 (Z = Mg, Mn, Zn, Cd, In, Bi, Sn): High Chemical Flexibility Resulting in Good Nonlinear-Optical Properties

Seven acentric sulfides Ba 6 (Cu x Z y )Sn 4 S 16 (Z = Mg, Mn, Zn, Cd, In, Bi, Sn) were grown by a high-temperature salt flux method. The crystal structures of the Ba 6 (Cu x Z y )Sn 4 S 16 (Z = Mg, Mn, Zn, Cd, In, Bi, Sn) compounds were determined by single-crystal X-ray diffraction with the aid of solid-state NMR spectroscopy. The Ba 6 (Cu x Z y )Sn 4 S 16 (Z = Mg, Mn, Zn, Cd, In, Bi) compounds are isostructural and crystallize in the Ba 6 Ag 4 Sn 4 S 16 structure type. The Sn-containing compound exhibits high structural similarity to Ba 6 (Cu x Z y )Sn 4 S 16 (Z = Mg, Mn, Zn, Cd, In, Bi) with the presence of an interstitial atomic position partially occupied by Sn atoms. The chemical bonding characteristics of Ba 6 (Cu 2.9 Sn 0.4 )Sn 4 S 16 were understood with electron localization function calculations coupled with crystal orbital Hamilton population calculations. The Ba–S and Cu–S interactions are dominantly ionic, but the Sn–S interactions consist of strong covalent bonding characteristics in Ba 6 (Cu 2.9 Sn 0.4 )Sn 4 S 16 . The monovalent Cu atoms, mixed with certain metals with various oxidation states, significantly shift the optical properties of the Ba 6 (Cu x Z y )Sn 4 S 16 (Z = Mg, Mn, Zn, Cd, In, Bi) compounds. This results in a good balance between the second-harmonic-generation (SHG) response and laser damage threshold (LDT). Ba 6 (Cu 1.9 Zn 1.1 )Sn 4 S 16 possesses a high SHG response and a high LDT of 2.8 × AGS and 3 × AGS, respectively. Here, a density functional theory calculation revealed that CuS 4 and SnS 4 tetrahedra significantly contribute to the SHG response in Ba 6 (Cu 2 Mg)Sn 4 S 16 , which also confirmed that CuS 4 tetrahedra are crucial for the stability and optical properties of the Ba 6 (Cu x Z y )Sn 4 S 16 (Z = Mg, Mn, Zn, Cd, In, Bi, Sn) compounds revealed by electronic structure analysis.

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

Experimental Insights into the Coupling of Methane Combustion and Steam Reforming in a Catalytic Plate Reactor in Transient Mode

The microstructured reactor concept is very promising technology to develop a compact reformer for distributed hydrogen generation. In this work, a catalytic plate reactor (CPR) is developed and investigated for the coupling of methane combustion (MC) and methane steam reforming (MSR) over Pt/Al 2 O 3 -coated microchannels in cocurrent and counter-current modes in transient experiments during start-up. A three-dimensional (3D) computational fluid dynamics (CFD) simulation shows uniform velocity and pressure distribution profiles in microchannels. For a channel velocity from 5.1 to 57.3 m/s in the combustor, the oxidation of methane is complete and self-sustainable without explosion, blow-off, or extinction; nevertheless, flashbacks are observed in counter-current mode. In the reformer, the maximum methane conversion is 84.9% in cocurrent mode, slightly higher than that of 80.2% in counter-current mode at a residence time of 33 ms, but at the cost of three times higher energy input in the combustor operating at ~1000 °C. Nitric oxide (NO) is not identified in combustion products, but nitrous oxide (N 2 O) is a function of coupling mode and forms significantly in cocurrent mode. Furthermore, this research would be helpful to establish the start-up strategy and environmental impact of compact reformers on a small scale.

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

Enhancing the accuracy of density functional tight binding models through ChIMES many-body interaction potentials

Semi-empirical quantum models such as Density Functional Tight Binding (DFTB) are attractive methods for obtaining quantum simulation data at longer time and length scales than possible with standard approaches. However, application of these models can require lengthy effort due to the lack of a systematic approach for their development. In this work, we discuss the use of the Chebyshev Interaction Model for Efficient Simulation (ChIMES) to create rapidly parameterized DFTB models, which exhibit strong transferability due to the inclusion of many-body interactions that might otherwise be inaccurate. We apply our modeling approach to silicon polymorphs and review previous work on titanium hydride. We also review the creation of a general purpose DFTB/ChIMES model for organic molecules and compounds that approaches hybrid functional and coupled cluster accuracy with two orders of magnitude fewer parameters than similar neural network approaches. In all cases, DFTB/ChIMES yields similar accuracy to the underlying quantum method with orders of magnitude improvement in computational cost. In conclusion, our developments provide a way to create computationally efficient and highly accurate simulations over varying extreme thermodynamic conditions, where physical and chemical properties can be difficult to interrogate directly, and there is historically a significant reliance on theoretical approaches for interpretation and validation of experimental results.

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

First-Principles Study on the Role of Cu and Cl-Based Dopants in NiO

Utilization of wide band gap oxide-based materials in thin-film solar energy technologies has increased in recent years. Among the numerous candidate oxide materials, NiO has shown many desirable optoelectronic properties that are applicable to thin-film PV technologies such as cadmium telluride PV. However, one critical factor requiring further investigation is the p-type doping behavior of NiO, specifically when the cadmium telluride solar cell undergoes conventional processes such as copper doping and chlorine-based activation treatment. The previous literature has shown a large degree of variability in hole concentrations in NiO when copper is used as the primary dopant. This study uses first-principles computational modeling based on density functional theory coupled with defect equilibria calculations to quantitatively explore the role of copper and chlorine-based extrinsic dopants in the p-type doping activity of NiO. The study reveals the importance of extrinsic dopants and their binding interactions with nickel vacancies to effectively p-dope NiO. It is suggested that both the formation of V Ni + Cu Ni and V Ni + Cl O defect pairs under a supersaturated state of NiO are potential mechanisms for increasing hole densities. On the other hand, the production of 2Cu Ni severely limits the effectiveness of p-doping in NiO, even in the presence of the aforementioned defect pairs. Furthermore, the study provides a guideline for experimentalists interested in using copper or chlorine species to understand how to controllably p-dope NiO during thin-film synthesis.

chlorine↗
DOE OSTI2022

Probing eddy size and its effective mixing length in stably stratified roughness sublayer flows

Stably stratified roughness sublayer flows are ubiquitous yet remain difficult to represent in models and to interpret using field experiments. Here, continuous high-frequency potential temperature profiles from the forest floor up to 6.5 times the canopy height observed with distributed temperature sensing (DTS) are used to link eddy topology to roughness sublayer stability correction functions and coupling between air layers within and above the canopy. The experiments are conducted at two forest stands classified as hydrodynamically sparse and dense. Near-continuous profiles of eddy sizes (length scales) and effective mixing lengths for heat are derived from the observed profiles using a novel conditional sampling approach. The approach utilizes potential temperature isoline fluctuations from a statically stable background state. The transport of potential temperature by an observed eddy is assumed to be conserved (adiabatic movement) and we assume that irreversible heat exchange between the eddy and the surrounding background occurs along the (vertical) periphery of the eddy. This assumption is analogous to Prandtl's mixing-length concept, where momentum is transported rapidly vertically and then equilibrated with the local mean velocity gradient. A distinct dependence of the derived length scales on background stratification, height above ground, and canopy characteristics emerges from the observed profiles. Implications of these findings for (1) the failure of Monin–Obukhov similarity in the roughness sublayer and (2) above-canopy flow coupling to the forest floor are examined. The findings have practical applications in terms of analysing similar DTS data sets with the proposed approach, modelling roughness sublayer flows, and interpreting nocturnal eddy covariance measurements above tall forested canopies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

New Insights into the Bulk and Surface Defect Structures of Ceria Nanocrystals from Neutron Scattering Study

Neutron diffraction and pair distribution function studies coupled with Raman spectroscopy have successfully unraveled the detailed oxygen defect structures of ceria nanocubes and nanorods. Two types of defect sites are revealed for the ceria nanocrystals: surface and bulk defects. It is proposed that the surface oxygen defects in both types of CeO 2 nanocrystals are predominantly the partially reduced Ce 3 O 5+x , with the bulk defect structures dominated by interstitial Frenkel-type oxygen vacancies. Ceria nanorods possess much higher concentration of surface oxygen defects relative to the nanocubes, albeit with only slightly higher concentration of bulk Frenkel-type oxygen vacancies. Upon annealing the nanorod sample at 600 °C under vacuum (~10 –4 to 10 –5 mbar), a partially reduced ceria phase with long-range oxygen vacancy ordering (Ce 3 O 5+x ) has been observed experimentally for the first time. This intriguing observation that surface defect phases can take on ordered defect sublattices under certain conditions is of great value in understanding the temperature-dependent catalytic performance of ceria nanocrystals. Furthermore, a drastic decrease of the surface vacancies in the ceria nanocrystals is observed upon exposure to SO 2 , especially for the nanorods, a likely origin for the sulfur poisoning effect on ceria-based materials. This study suggests that tailoring surface morphology is a promising strategy to control defect properties of ceria nanomaterials. It also provides fundamental insights to stabilize surface oxygen defects in CeO 2 nanocrystals to achieve high redox performance under corrosive environments such as under SO 2 /SO x exposure.

36 MATERIALS SCIENCE↗
DOE OSTI2022

PDFDataExtractor: A Tool for Reading Scientific Text and Interpreting Metadata from the Typeset Literature in the Portable Document Format

The layout of portable document format (PDF) files is constant to any screen, and the metadata therein are latent, compared to mark-up languages such as HTML and XML. No semantic tags are usually provided, and a PDF file is not designed to be edited or its data interpreted by software. However, data held in PDF files need to be extracted in order to comply with opensource data requirements that are now government-regulated. In the chemical domain, related chemical and property data also need to be found, and their correlations need to be exploited to enable data science in areas such as data-driven materials discovery. Such relationships may be realized using text-mining software such as the “chemistry-aware” natural-language-processing tool, ChemDataExtractor; however, this tool has limited data-extraction capabilities from PDF files. This study presents the PDFDataExtractor tool, which can act as a plug-in to ChemDataExtractor. It outperforms other PDF-extraction tools for the chemical literature by coupling its functionalities to the chemical-named entityrecognition capabilities of ChemDataExtractor. The intrinsic PDF-reading abilities of ChemDataExtractor are much improved. The system features a template-based architecture. This enables semantic information to be extracted from the PDF files of scientific articles in order to reconstruct the logical structure of articles. While other existing PDF-extracting tools focus on quantity mining, this template-based system is more focused on quality mining on different layouts. PDFDataExtractor outputs information in JSON and plain text, including the metadata of a PDF file, such as paper title, authors, affiliation, email, abstract, keywords, journal, year, document object identifier (DOI), reference, and issue number. With a self-created evaluation article set, PDFDataExtractor achieved promising precision for all key assessed metadata areas of the document text.

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

Quantum Vibronic Effects on the Electronic Properties of Molecular Crystals

We present a study of molecular crystals, focused on the effect of nuclear quantum motion and anharmonicity on their electronic properties. We consider a system composed of relatively rigid molecules, a diamondoid crystal, and one composed of floppier molecules, NAI-DMAC, a thermally activated delayed fluorescence compound. We compute fundamental electronic gaps at the density functional theory (DFT) level of theory, with the Perdew–Burke–Erzenhof (PBE) and strongly constrained and approximately normed (SCAN) functionals, by coupling first-principles molecular dynamics with a nuclear quantum thermostat. We find a sizable zero-point renormalization (ZPR) of the band gaps, which is much larger in the case of diamondoids (0.6 eV) than for NAI-DMAC (0.22 eV). We show that the frozen phonon (FP) approximation, which neglects intermolecular anharmonic effects, leads to a large error (~50%) in the calculation of the band gap ZPR. Instead, when using a stochastic method, we obtain results in good agreement with those of our quantum simulations for the diamondoid crystal. However, the agreement is worse for NAI-DMAC where intramolecular anharmonicities contribute to the ZPR. Our results highlight the importance of accurately including nuclear and anharmonic quantum effects to predict the electronic properties of molecular crystals.

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

A Comparison of Electronic Structure Methods for Predicting the Hydrogenation Energies of Candidate Molecules for Hydrogen Storage

The development of novel energy materials and fuels is required to expand current available energy sources. Aiming to reach this goal, there is growing interest in using molecular hydrogen as an energy carrier due to its abundance and high energy density. Liquid organic hydrogen carriers (LOHCs) are a promising route to the large-scale storage and transport of hydrogen for use in the energy economy. The search for thermodynamically viable LOHC molecules for real world use has led to a set of constraints on the dehydrogenation enthalpy and the minimum gravimetric hydrogen capacity. These constraints allow one to formulate the search for an ideal LOHC candidate molecule as an optimization problem well suited to the strengths of machine learning and artificial intelligence computational approaches. A critical barrier to a large-scale, high-throughput screening of LOHC candidate molecules is the lack of reliable training data. Computational electronic structure methods including density functional theory, coupled cluster approximations, and diffusion Monte Carlo can be used to provide training data where experimental data are either unreliable or do not exist. In this work, we use these methods to calculate the dehydrogenation energies and enthalpies of candidate LOHC molecules.

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

Adjusting the Energy Profile for CH–O Interactions Leads to Improved Stability of RNA Stem-Loop Structures in MD Simulations

The role of ribonucleic acid (RNA) in biology continues to grow, but insight into important aspects of RNA behavior is lacking, such as dynamic structural ensembles in different environments, how flexibility is coupled to function, and how function might be modulated by small molecule binding. In the case of proteins, much progress in these areas has been made by complementing experiments with atomistic simulations, but RNA simulation methods and force fields are less mature. It remains challenging to generate stable RNA simulations, even for small systems where well-defined, thermostable structures have been established by experiments. Further many different aspects of RNA energetics have been adjusted in force fields, seeking improvements that are transferable across a variety of RNA structural motifs. In this work, the role of weak CH···O interactions is explored, which are ubiquitous in RNA structure but have received less attention in RNA force field development. By comparing data extracted from high-resolution RNA crystal structures to energy profiles from quantum mechanics and force field calculations, it is shown that CH···O interactions are overly repulsive in the widely used Amber RNA force fields. A simple, targeted adjustment of CH···O repulsion that leaves the remainder of the force field unchanged was developed. Then, the standard and modified force fields were tested using molecular dynamics (MD) simulations with explicit water and salt, amassing over 300 μs of data for multiple RNA systems containing important features such as the presence of loops, base stacking interactions as well as canonical and noncanonical base pairing. In this work and others, standard force fields lead to reproducible unfolding of the NMR-based structures. Including a targeted CH···O adjustment in an otherwise identical protocol dramatically improves the outcome, leading to stable simulations for all RNA systems tested.

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

Stabilizing Aqueous Three-Dimensional Printed Constructs Using Chitosan-Cellulose Nanocrystal Assemblies

The assembly and binding of nanoparticles at the interfaces of aqueous two-phase systems enable the three-dimensional (3D) printing of all-aqueous naturally occurring materials. When a dispersion of cellulose nanocrystals (CNCs) in an aqueous solution of polyethylene glycol (PEG) is brought into contact with chitosan dissolved in an aqueous solution of dextran, the CNCs and chitosan diffuse to the interface between the two immiscible aqueous solutions, electrostatically interact, and form a solid, membranous layer sufficiently rapidly to 3D print tubules of one liquid in the other. Additionally, the diameter, length, spatial arrangement, and stability of the printed tubules can be broadly controlled. Adsorption and directional diffusion of ionic species across the membranous layer make heavy metal ion removal possible. The results present a platform for fabricating and developing all-aqueous compartmentalized systems where function can be independently coupled to the inherent functionality of the nanoparticles or ligands.

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

Unraveling the CO Oxidation Mechanism over Highly Dispersed Pt Single Atom on Anatase TiO 2 (101)

Catalysts with noble metals deposited as single atoms on metal oxide supports have recently been studied extensively due to their maximized metal utilization and potential for performing difficult chemical conversions owing to their unique electronic properties. Understanding of the reaction mechanisms on supported single-metal atoms is still limited but is highly important for designing more efficient catalysts. In this study, we report the complexity of the CO oxidation reaction mechanism on Pt single atoms supported on anatase TiO 2 (Pt SA /a-TiO 2 ) by coupling density functional theory (DFT) calculations and microkinetic analysis with kinetic measurements, in situ/operando infrared, and X-ray absorption spectroscopies. Starting from the adsorbed Pt SA occupying an O vacancy induced by reductive pretreatment, we show that CO oxidation follows a complex mechanism consisting of initiation steps to reorganize the active site and multibranch reactive cycles, with the Pt SA /a-TiO 2 catalyst not returning to its initial configuration. The initiation step consists of CO and O 2 adsorption healing the O vacancy, followed by CO oxidation using gas-phase CO to form Pt(CO). The reactive cycle alternates O 2 adsorption and dissociation to oxidize the catalyst to Pt(O)(O)(CO) and branching pathways of competing Langmuir–Hinshelwood (LH)- or Eley–Rideal (ER)-type CO oxidation steps to reduce it again to Pt(CO). In situ/operando infrared experiments, including cryogenic CO adsorption and isotopic CO exchange, confirm the combined involvement of strongly adsorbed CO and gas-phase CO in an Eley–Rideal step along the reaction cycle. Microkinetic modeling shows that Pt single atoms are present in a mixture of Pt(CO), Pt(CO)(O 2 ), Pt(O)(CO)(O 2 ), and Pt(CO)(CO 3 ) structures as the main intermediates during steady-state CO oxidation, all having the C–O vibrational stretch close to the experimentally observed value of 2115 cm –1 . Microkinetic modeling also shows that the fractional orders of CO and O 2 measured experimentally originate from multiple steps with a high degree of rate control and not from a simple competitive adsorption. The results demonstrate the complex reaction pathways that even CO oxidation on a simple single-atom system can follow, providing mechanistic insights for designing efficient Pt-based single-atom catalysts. Here, we further show that microkinetic modeling results are sensitive to changes in energies of intermediate and transition states within errors of density functional theory, which can ultimately lead to incorrect conclusions regarding the reaction pathways and most abundant reaction intermediates if not accounted for by experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Estimation of wintertime cloud radiative effects in the Western Arctic, a function of cloud-moisture-coupling and sea ice conditions

Cloud radiative effects are estimated using up- and downward ground-based observations of broadband thermal infrared sensors in cloudy conditions which are provided by the ARM data repository. The ARM Surface Radiation measurement is analyzed. To study the cloud radiative effects, the observed thermal infrared measurements under cloudy conditions are analyzed together with microphysical cloud properties like liquid and ice water path obtained from the CloudNet algorithms applied to the remote-sensing instrumentation suite at the ARM site in Utqiagvik. Cloud observations have been sorted depending on the coupling of the cloud with upwind water vapor transport events that can interact with sea ice upwind. Statistics of surface net thermal fluxes and cloud water/ice properties show a distinguished distribution for coupled and decoupled clouds. It was found that coupled clouds contain more liquid water and are more efficient reduce the cooling of surface LW radiation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Melting point of iron at high pressure: An assessment of uncertainties and effect of electronic temperature

An accurate calculation of the melting point of iron at various pressures in the Earth's core is important for understanding the core structure, geodynamo, and the Earth's history. Previous studies have assessed the melt line of iron at these extreme conditions using various experimental measurement techniques as well as both ab initio and classic molecular dynamics simulations. However, experimental measurements have uncertainties up to several hundred Kelvin, and inconsistencies remain among simulation results. Here in this work, we propose an iterative framework that couples density functional theory (DFT) calculations and molecular dynamics simulations performed using an ensemble of interatomic potentials to assess the effect of electronic temperature on the melting point. We systematically validate the potentials by comparing lattice constants and phonon dispersion curves at 0 K and enthalpy differences between liquid and HCP, FCC, BCC phases of iron close to the melt line at 300 GPa with DFT. Our results show that HCP iron melts at 6144 K (at 300 GPa), BCC phase is thermodynamically unstable, and FCC is metastable at this temperature. The melting points of FCC and BCC phases at 300 GPa are 5858 and 5647 K, respectively.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

R-matrix calculations for opacities: I. Methodology and computations

Abstract An extended version of the R -matrix methodology is presented for calculation of radiative parameters for improved plasma opacities. Contrast and comparisons with existing methods primarily relying on the distorted wave approximation are discussed to verify accuracy and resolve outstanding issues, particularly with reference to the opacity project (OP). Among the improvements incorporated are: (i) large-scale Breit–Pauli R -matrix calculations for complex atomic systems including fine structure, (ii) convergent close coupling wave function expansions for the ( e + ion) system to compute oscillator strengths and photoionization cross sections, (iii) open and closed shell iron ions of interest in astrophysics and experiments, (iv) a treatment for plasma broadening of autoionizing resonances as function of energy-temperature-density dependent cross sections, (v) a ‘top-up’ procedure to compare convergence with R -matrix calculations for highly excited levels, and (vi) spectroscopic identification of resonances and bound ( e + ion) levels. The present R -matrix monochromatic opacity spectra are fundamentally different from OP and lead to enhanced Rosseland and Planck mean opacities. An outline of the work reported in other papers in this series and those in progress is presented. Based on the present re-examination of the OP work, opacities of heavy elements might require revisions in high temperature-density plasma sources.

Pradhan, A. K. (ORCID:0000000187753643)↗