Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “STRENGTH CALCULATIONS”

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 901 records · Page 50
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Measuring the Debye energy in superconductors via two-electron photoemission spectroscopy

In this work, we demonstrate theoretically that double angle resolved photoemission spectroscopy (2eARPES) can directly probe the existence of Cooper pairs away from the Fermi surface, and can thus provide insight into the characteristic energy scale around the Fermi surface, the Debye energy, in which electrons are bound into Cooper pairs. To this end, we compute the photoelectron counting rate 𝑃 (2) in two different types of unconventional superconductors, a 𝑑 𝑥 2 −𝑦 2 -wave superconductor, and a topological superconductor with a broken time-reversal symmetry. We show that 𝑃 (2) provides insight into the relative strength of intra- and interband pairing in multiband systems, as well as into the spin polarization of the bands.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Unraveling excitation energy transfer assisted by collective behaviors of vibrations

Abstract We investigate how collective behaviors of vibrations such as cooperativity and interference can enhance energy transfer in a nontrivial way, focusing on an example of a donor–bridge–acceptor trimeric chromophore system coupled to two vibrational degrees of freedom. Employing parameters selected to provide an overall uphill energy transfer from donor to acceptor, we use numerical calculations of dynamics in a coupled exciton–vibration basis, together with perturbation-based analytics and calculation of vibronic spectra, to identify clear spectral features of single- and multi-phonon vibrationally-assisted energy transfer (VAET) dynamics, where the latter include up to six-phonon contributions. We identify signatures of vibrational cooperation and interference that provide enhancement of energy transfer relative to that obtained from VAET with a single vibrational mode. We observe a phononic analogue of two-photon absorption, as well as a novel heteroexcitation mechanism in which a single phonon gives rise to simultaneous excitation of both the trimeric system and the vibrational degrees of freedom. The impacts of vibrations and of the one- and two-phonon VAET processes on the energy transfer are seen to be quite different in the weak and strong site–vibration coupling regimes. In the weak coupling regime, two-phonon processes dominate, whereas in the strong coupling regime up to six-phonon VAET processes can be induced. The VAET features are seen to be enhanced with increasing temperature and site–vibration coupling strength, and are reduced in the presence of dissipation. We analyze the dependence of these phenomena on the explicit form of the chromophore–vibration couplings, with comparison of VAET spectra for local and non-local couplings.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Effect of Si impurities on microstructure and tensile properties of a cast Al-Mg-Fe alloy

Al-Mg alloys are attractive for structural castings owing to their superior strength and ductility in the as-cast state. However, given the tight tolerance for impurities, Al-Mg alloys produced from secondary sources still face multiple challenges. Here, we report the effect of increased Si impurity (0.05–1.6 wt%), which is commonly found in secondary Al sources, on microstructure and tensile properties of a cast Al-4.3Mg-1.6Fe (wt%) alloy, commercially referred to as Castaduct-42 alloy. Microstructural characterization revealed that Si addition increased the volume fraction, size, and aspect ratio of primary Al13Fe4 intermetallic particles as well as the volume fraction of other binary and ternary eutectic phases. Tensile testing results demonstrated that increasing Si impurities from 0.05 to 1.6 wt% reduced ductility from 14.3 ± 1.4 % to 2.2 ± 1.0 %. A particle cracking damage accumulation model coupled with failure analysis indicated primary Al13Fe4 intermetallic particles to be the major contributing factor to the deterioration in ductility with increasing Si content. Additionally, the stabilizing effect of Si on primary Al 13 Fe 4 was inconsistent with the CALPHAD calculation results based on existing CALPHAD databases which predict a slightly decreasing primary Al13Fe4 phase fraction with increasing Si concentration. This work provides new insights into the phase stability and mechanical behavior of the lesser studied Al-Mg-Fe-Si alloy system that will contribute to the development of sustainable cast Al-Mg based alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Deprotonated sulfamic acid and its homodimers: Does sulfamic acid adopt zwitterion during cluster growth?

Here we present a joint experimental and computational study on the geometric and electronic structures of deprotonated sulfamic acid (SA) clusters [(SA) n –H] – (n = 1, 2) employing negative ion photoelectron spectroscopy and high-level ab initio calculations. The photoelectron spectra provide the vertical/adiabatic detachment energy (VDE/ADE) of the sulfamate anion (SM – ) H 2 N•SO 3 – at 4.85 ± 0.05 and 4.58 ± 0.08 eV, respectively, and the VDE and ADE of the SM – •SA dimer at 6.41 ± 0.05 and 5.87 ± 0.08 eV, respectively. The significantly increased electron binding energies of the dimer confirm the enhanced electronic stability upon the addition of one SA molecule. The CCSD(T)-predicted VDEs/ADEs agree excellently with the experimental data, confirming the identified structures as the most stable ones. Two types of dimer isomers possessing different hydrogen bonding (HB) motifs are identified, corresponding to SM – binding to a zwitterionic SA (SM – •SA z ) and a canonical SA (SM – •SA c ), respectively. Two N–H$\cdots$O HBs and one superior O–H$\cdots$O HB are formed in the lowest-lying SM – •SA c , while SM – •SA z has three moderate N–H$\cdots$O HBs, with the former being 4.71 kcal/mol more stable. Further theoretical analyses reveal that the binding strength advantage of SM – •SA c over SM – •SA z arises from its significant contributions of orbital interactions between fragments, illustrating that sulfamate strongly interacts with its parent SA acid and preferably chooses the canonical SA in the subsequent cluster formations. Given the prominent presence of SA, this study provides the first evidence that the canonical dimer model of sulfamic acid should exist as a superior configuration during cluster growth.

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

Dipole-lattice nanoparticle resonances in finite arrays

We investigate how the periodic lattices define the collective optical characteristics of the silicon and titanium nanoparticle arrays. We examine the effects of dipole lattice on the resonances of optical nanostructures, including those made of lossy materials, such as titanium. Our approach involves employing coupled-electric-magnetic-dipole calculations for finite-size arrays, as well as lattice sums for effectively infinite arrays. Our model shows that the convergence to the infinite-lattice limit is faster when the resonance is broad, requiring fewer array particles. Our approach differs from previous works by altering the lattice resonance through modifications in the array period. We observed that a higher number of nanoparticles is necessary to achieve convergence to the infinite-array limit. Additionally, we observe that the lattice resonances excited next to higher diffraction orders (such as second order) converge more quickly toward the ideal case of an infinite array than the lattice resonances related to the first diffraction order. This work reports on the significant advantages of using a periodic arrangement of lossy nanoparticles and the role of collective excitation in enhancing response from transition metals, such as titanium, nickel, tungsten, and so on. The periodic arrangement of nanoscatterers allows for the excitation of strong dipoles, boosting the performance of nanophotonic devices and sensors by improving the strength of localized resonances.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Distorted dislocation cores and asymmetric glide resistances in titanium

The determination of Critical Resolved Shear Stress (CRSS) in titanium for basal, prismatic, and pyramidal slip-planes without empirical constants is presented by combining Density Functional Theory (DFT) and anisotropic elasticity. A new mechanism leading to tension-compression (T-C) asymmetry in the CRSS levels has been revealed for the first time. The conditions for this asymmetry are established, involving a complex interplay between the dislocation core-structure and core-advance behavior. The three conditions for T-C asymmetry that need to be simultaneously satisfied can be summarized as: (1) a medium stacking fault width, d, (3 < d/b F < 10, where b F is the magnitude of the Burgers vector of the full dislocation), (2) an asymmetric core-structure of the extended dislocation (ξ 1 ≠ ξ 2 , where ξ 1 and ξ 2 are the core-widths of the first and second partials, respectively), and (3) intermittent motion of the partials (Δd/b F ≠ 0, where Δd is the magnitude of fluctuation in stacking fault width during intermittent motion). Pyramidal-slip in titanium satisfies all three conditions, resulting in significant T-C asymmetry. The CRSS theory considers a Wigner-Seitz (W-S) cell based domain area assigned to each lattice site for the calculations of core-energies accurately capturing the slip-plane lattice. The W-S based approach is essential due to the lower symmetry of the HCP crystal circumventing potential errors associated with the one-dimensional atomic-row approximation. Dislocation core structures are obtained for all the slip-systems in titanium showing significant distortion of the disregistry profile governing the core-advance behavior. The CRSS values predicted from the theory show agreement with the experimental CRSS levels reported in the literature.

Dislocation core↗
DOE OSTI2021

Vapor-Phase Intercalation of Cesium into Black Phosphorous

Cesium vapors were charged into black phosphorous (BP) flakes at varied times and at a temperature gradient of 150 °C. The X-ray diffraction (XRD) measurements of these samples suggest a reduction in the strength of van der Waal interactions between BP layers leading to the loss of coherence of out-of-plane peaks. At the same time, the three main Raman modes of BP (A g 1 , B 2 g , and A g 2 ) steadily redshifted as exposure times were increased, with modes B 2 g and A g 2 shifting faster than A g 1 . After initial rapid downshifts of active BP phonon modes, this intercalation strategy showed its limits following prolonged exposure times. Saturation of BP flakes by Cs vapors ensued and the kinetics was fitted with an exponential decay function. Furthermore, the thermoelectric power (TEP) of cesiated BP exhibited an inversion in sign from positive to negative around 400 K, lending credence to the transformation of as-prepared BP which is a p-type semiconductor to an n-type equivalent due to Cs atom intercalation driven shifting of the Fermi level toward the conduction band of BP and the donation of electrons from Cs. Furthermore, density functional theory (DFT) calculations were used to delve deeper into understanding Cs intercalation on the structural evolution of BP.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Adhesion-force micro-scale study of desiccating granular material

Experiments on five-, four-, three- and two-wet-hydrophilic-grain clusters were performed to investigate evolution of adhesion of granular media during drying on the micro-scale. The experiments show that the adhesion-force of a cluster initially grows at most to three times the original value before decreasing to zero by the end of evaporation. The adhesion-force is composed of capillary pressure force acting over the liquid/solid contact surface area, and surface tension forces acting over the three-phase contact perimeter length. This is in contrast with most macro-scale phenomenological models, in which the only desaturation process variables affecting strength are suction and saturation. Both the contact surface area and contact perimeter length are reduced to zero upon complete liquid evaporation. The morphology of an evaporating water body evolves through slow flow controlled by evaporation rate, interrupted by various modes of fast air entry, which are non-equilibrium jumps of liquid/gas interfaces (Haines jumps). The instabilities involve large adhesion force discontinuities and substantial water mass reconfiguration with water flow in an extremely short time, which makes the process transient. The reconfigurations can reduce the original multi-grain water clusters to four-, three- and two-grain clusters by way of three different instability modes: of thin-sheet instability, or meniscus snap-through instability, depending on the sign of the Gauss curvature of the liquid surface, or finally, for two-grain bridges only, a liquid wire pinch-off. For larger meso-scale assemblies, however, the global adhesion-force evolution is little affected by the jumps. The air entries are potential sites for drying cracks. The (approximately) calculated capillary pressure for two- and three-grain clusters, in no cases is seen to reach high values, predicted from water retention curves.

Engineering↗
DOE OSTI2023

Plasma production and ion heating systems for the Material Plasma Exposure eXperiment

Plasma production and ion heating in the Material Plasma Exposure eXperiment (MPEX), whose design is nearing completion, is accomplished using continuous wave RF power with average power density up to 1.6 MW/m2 at the plasma interface. Plasma is produced using helicon waves coupled through a single helical antenna at 13.56 MHz with power supplied by three 100 kW fixed-frequency RF generators feeding a power combiner network followed by the matching network and launcher. The helicon source can utilize various gasses including hydrogen, deuterium, and helium, with a magnetic field strength in the source region up to 0.2 T, and maximum device |B| of 2.5 T. Power is coupled to ions via ion cyclotron heating at the fundamental resonance using a pair of phased helical antennas operating in the frequency range 4-9 MHz, that launch waves towards the resonance from the antenna region where ω > ωci. ICH power is supplied by a single 500 kW tunable RF transmitter through a 90° power splitter and matching/decoupling network. In the case of both the helicon and ICH systems the antennas are located external to the vacuum, with power transferred through novel water-cooled coaxial vacuum windows consisting of fused quartz outer cylinders and silicon nitride inner cylinders with forced convection water cooling between them. The antenna enclosures are pressurized with dry air to 3 bar absolute for voltage standoff.Several 3-D COMSOL models have been created to simulate the two systems. A model of the helicon region utilizing a cold plasma dielectric tensor with accurate magnetic field and realistic plasma density profiles has been used to calculate the plasma loading/complex antenna input impedance at the launcher feed for various ne and |B| values, for the purpose of estimating power handling. It also incorporates the geometry of all launcher structures relevant to this determination. A still more detailed model of the launcher together with a lossy dielectric plasma surrogate has been used to determine RF electric field values and power losses in the device components.Similar models have been produced to predict the performance and power handling of the ICH launcher. For this device the impedance matrix of the two-element antenna array is calculated using a warm plasma model, necessary to properly determine the wave propagation and absorption.

Goulding, Richard↗
DOE OSTI2025

Metrics and extrapolation of resonant magnetic perturbation thresholds for ELM suppression

This large database study of resonant magnetic perturbation (RMP) edge localized mode (ELM) suppression thresholds in the AUG, DIII-D, EAST, and KSTAR tokamaks details the key strengths and weaknesses of RMP metrics. The RMP ELM suppression database used for this work contains plasma information at the time of transition from ELMing to ELM suppressed states where a clear experimental threshold is identified. The experimental threshold distributions are compared for five metrics: (1) the island overlap width, (2) pedestal top Chirikov overlap, (3) peeling edge displacement, (4) pedestal top resonant drive, and (5) edge dominant mode overlap. The distributions, the regularity of the dependence on RMP coil currents, and the sensitivities of a given metric to equilibrium reconstruction details are compared. The overlap metric proves to be a good compromise between including the appropriate plasma response physics and maintaining a numerical robustness. This quantity does not exhibit clear power-law scalings for projection, but machine learning can assist in predicting thresholds within the existing parameter ranges and providing uncertainty quantification of those predictions. Two new first-principles models, one utilizing a threshold from the non-linear Modified Rutherford equation evaluated at the pedestal top and one utilizing the SLAYER code to calculate the linear tearing threshold from torque balance, offer possible paths to extrapolation beyond the existing database parameter space.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

The effects of heat treatment on the response of WE43 Mg alloy: crystal plasticity finite element simulation and SEM-DIC experiment

The present report addressed the effects of heat treatment on the mechanical response of a WE43 Mg alloy using an integrated framework of SEM-DIC experiment and CPFE simulation. Both macroscopic responses including yield strength, ultimate strength, ductility, and microscopic responses, including local displacement and strain maps, were experimentally investigated. The focus of this work is to use the CPFE simulation as an integrated computational tool to study the effects of heat treatment. The CPFE framework was evaluated using the local fields of displacement and strain obtained from the SEM-DIC experiment rather than the conventional scheme of using macroscopic responses at different loading directions. Subsequently, the information which is available using CPFE, such as the critical resolved shear stress (CRSS) and relative slip activity, was used to study the effects of heat treatment on the response of WE43 Mg alloy. The contributions of different strengthening mechanisms on the CRSS were addressed. The results show that effects of heat treatment can be captured using the predominant mechanisms of the grain size effect and the influence of precipitates. Finally, it has been shown that classical Hall-Petch in which one constant can capture the size effects, should be modified. To do so, each deformation mode should have a unique Hall-Petch constant, which are calculated here for the WE43 Mg alloy.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Impact of control blade insertion on the deformation behavior of SiC-SiC channel boxes in BWRs

This article describes the analysis of distortion of a silicon carbide fiber-reinforced, silicon carbide matrix (SiC-SiC) composite channel box under in-reactor conditions of a boiling water reactor (BWR). The BWR core has significant gradients in the fast neutron flux across the channel box due to the presence of water rods within the fuel assemblies, and these gradients increase further with the insertion of control blades. As a result of the temperature and neutron flux dependent irradiation-induced swelling of SiC, the SiC-SiC composite channel box can undergo distortion. In this work, we evaluate the SiC-SiC channel box distortion for three different control blade positions. This analysis is based on the neutron flux and temperature distributions in the BWR core calculated using the neutronics code MPACT and thermal-hydraulics code CTF. This calculation is coupled through temperature feedback. Subsequently, we have performed structural analysis based on the calculated neutron flux and temperature distributions to determine the deformation and stress development in the channel box. The structural analysis was performed using the fuel performance modeling code BISON and the commercial finite element analysis software Abaqus. The results indicate that large gradients in fast neutron flux (up to 35–40% across a single axial level) will develop across the channel box. Due to these gradients, the channel box will undergo time-dependent bending for all the control blade positions in the assembly. The time-dependent bowing behavior is dominated by the transient swelling of SiC-SiC material under non-uniform neutron flux, and changes with variation in the control blade position. The bending will cause temporary interference between the channel box and control blade, and the interference is expected to be most severe for the fully inserted control blade position. The developed stresses due to differential swelling in the channel box exceed the proportional limit stress of the material, which may cause matrix microcracking in the channel box. However, the stresses remain below the tensile strength of the material, and therefore, development of a full, through-thickness crack in the channel box is not expected. Finally, further work is recommended to explore and evaluate the mitigation strategies.

42 ENGINEERING↗
DOE OSTI2024

Two-dimensional coherent spectrum of high-spin models via a quantum computing approach

Here in this work we present and benchmark a quantum computing approach to calculate the two-dimensional coherent spectrum (2DCS) of high-spin models. Our approach is based on simulating their real-time dynamics in the presence of several magnetic field pulses, which are spaced in time. We utilize the adaptive variational quantum dynamics simulation algorithm for the study due to its compact circuits, which enables simulations over sufficiently long times to achieve the required resolution in frequency space. Specifically, we consider an antiferromagnetic quantum spin model that incorporates Dzyaloshinskii-Moriya interactions and single-ion anisotropy. The obtained 2DCS spectra exhibit distinct peaks at multiples of the magnon frequency, arising from transitions between different eigenstates of the unperturbed Hamiltonian. By comparing the one-dimensional coherent spectrum with 2DCS, we demonstrate that 2DCS provides a higher resolution of the energy spectrum. We further investigate how the quantum resources scale with the magnitude of the spin using two different binary encodings of the high-spin operators: the standard binary encoding and the Gray code. At low magnetic fields both encodings require comparable quantum resources, but at larger field strengths the Gray code is advantageous. Numerical simulations for spin models with increasing number of sites indicate a polynomial system-size scaling for quantum resources. Lastly, we compare the numerical 2DCS with experimental results on a rare-earth orthoferrite system. The observed strength of the magnonic high-harmonic generation signals in the 2DCS of the quantum high-spin model aligns well with the experimental data, showing significant improvement over the corresponding mean-field results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Effect of Cooling Rate on Nano-Eutectic Formation in Laser Surface Remelted and Rare Earth Modified Hypereutectic Al-20Si Alloys

Laser Surface Remelting (LSR) was applied to arc-melted Al-20Si-0.2Sr, Al-20Si-0.2Ce, and Al-20Si hypereutectic alloys to refine microstructures. Experiments revealed that microstructures in the melt pool varied from fully eutectic to a mixture of Al dendrites and inter-dendritic eutectic. We calculated cooling rates using the Eagar-Tsai model and correlated cooling rates with characteristic microstructures, revealing that a cooling rate on the order of 104 K/s could lead to maximized fully eutectic microstructure morphology. Due to rapid solidification, the Si composition in the LSR eutectic was measured at 18.2 wt.%, higher than the equilibrium eutectic composition of 12.6 wt.%Si. Compared to Al-20Si, Ce addition had no significant effect on the volume fraction of the fully eutectic structure but refined Si fibers to approximately 30 nm in diameter. Sr addition did not further refine the diameter of eutectic Si fibers compared to Al-20Si but increased the volume fraction of the fully eutectic microstructure morphology. The refinement ratio (φ) of the Si fiber diameter from the bottom of the melt pool to the surface for the three alloys was similar, at around 28%. The established correlation between the cooling rate and the size and morphology of the microstructure within the melt pool will enable tailoring of the microstructure in laser-processed as well as deposited alloys for high strength and plasticity.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Alpha-Quartz Plastic Strength Investigation Via Diffraction Experiments on Novaculite Using A D-Dia and Elastic Plastic Self- Consistent Interpretation [Thesis}

X-Ray Tomography. Porosity is commonly measured using mercury injection (MI) or water immersion porosimeter (WIP). Both MI and WIP utilize pressure to fill open voids with mercury or water respectively. A measure of the volume change of the sample is then used to estimate the percentage of open voids. However, for the purpose of rheological study it is important to get an accurate measure of open and closed voids. X-ray tomography has been selected as it offers a three-dimensional view into the sample that can quantify all pores limited only by the voxel size which is in the micron range for this study. Radiographs for tomography were collected at the Material Science and Technology Division at Los Alamos National Lab using the Carl Zeiss Xradio 520 instrument and Scout-and-Scan version 16.1 operating software. Two samples were imaged, the starting material and the deformed sample SiO2_65. 3001 radiographs were taken of the starting material with a 6 second exposure time using a 4x objective lens. The x-ray beam was set to 60 kilovoltage peak (kVp) and 5 watts. 1901 radiographs were taken of SiO2_65 with a 25 second exposure time using a 10x objective lens. The x-ray beam was set to 80 kVp and 7 watts. Radiograph files were analyzed by Brian Patterson using Avizo. Void and inclusion volumes were output by voxel sized (1.03 μm) slices, used to calculate a total percentage volume for the starting material.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Microstructure, electrical resistivity, and tensile properties of neutron-irradiated Cu–Cr–Nb–Zr

High strength, high conductivity copper alloys that can resist creep at high temperatures are one of the primary candidates for efficient heat exchangers in fusion reactors. Cu–Cr–Nb–Zr (CCNZ) alloys, which were designed to improve the strength and creep life of ITER Cu–Cr–Zr (CCZ) reference alloys, have been found to have comparable electrical conductivity and tensile properties to CCZ alloys. The measured creep rupture times for these improved alloys is about ten times higher than the ITER reference alloys at 90–125 MPa at 500 °C. However, the effects of neutron irradiation on these alloys, and the ensuing material properties, have not been studied; thus, their utility in a fusion reactor environment is not well understood. This study characterizes the room temperature mechanical and electrical properties of a neutron-irradiated CCNZ alloy and compares them to a neutron-irradiated ITER reference heat sink CCZ alloy. Tensile specimens were neutron irradiated in the High Flux Isotope Reactor (HFIR) to 5 dpa between 250 °C and 325 °C. Post-irradiation characterization included electrical resistivity measurements, hardness, and tensile tests. Microstructural evaluation used scanning electron microscopy, energy dispersive x-ray spectroscopy, and atom probe tomography to characterize the irradiation-produced changes in the microstructure and investigate the mechanistic processes leading to post-irradiation properties. Transmutation calculations were validated with composition measurements from atom probe data and used to calculate contributions to the increased electrical resistivity measured after irradiation. Comparisons with CCZ alloys in the same irradiation heat found that the post-irradiated CCNZ and CCZ alloys had comparable electrical resistivity. Although CCNZ alloys suffered more irradiation hardening than CCZ, the overall tensile behavior deviated very little from non-irradiated values in the temperature range studied.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Magnetic phase diagram of a two-orbital model for bilayer nickelates with varying doping

Motivated by the recently discovered high-T c bilayer nickelate superconductor La 3⁢ Ni 2 ⁢O 7 , we comprehensively research a bilayer 2×2×2 cluster for different electronic densities n by using the Lanczos method. We also employ the random-phase approximation to quantify the first magnetic instability with increasing Hubbard coupling strength, also varying n. Based on the spin structure factor S(q), we have obtained a rich magnetic phase diagram in the plane defined by n and U/W, at fixed Hund coupling, where U is the Hubbard strength and W the bandwidth. We have observed numerous states, such as A-AFM, Stripes, G-AFM, and C-AFM. At half-filling, n=2 (two electrons per Ni site, corresponding to N=16 electrons), the canonical superexchange interaction leads to a robust G-AFM state (π,π,π) with antiferromagnetic couplings both in-plane and between layers. By increasing or decreasing electronic densities, ferromagnetic tendencies emerge from the “half-empty” and “half-full” mechanisms, leading to many other interesting magnetic tendencies. In addition, the spin-spin correlations become weaker both in the hole or electron doping regions compared with half-filling. At n=1.5 (or N=12), density corresponding to La 3 ⁢Ni 2 ⁢O 7 , we obtained the “Stripe 2” ground state (antiferromagnetic coupling in one in-plane direction, ferromagnetic coupling in the other, and antiferromagnetic coupling along the z axis) in the 2×2×2 cluster. In addition, we obtained a much stronger AFM coupling along the z axis than the magnetic coupling in the xy plane. The random-phase approximation calculations with varying n give very similar results as Lanczos, even though both techniques are based on quite different procedures. Additionally, a state with q/π=(0.6,0.6,1) close to the E-phase wavevector is found in our RPA calculations by slightly reducing the filling to n=1.25, possibly responsible for the E-phase SDW recently observed in experiments. In conclusion, our predictions can be tested by chemically doping La 3 ⁢Ni 2 ⁢O 7 .

36 MATERIALS SCIENCE↗
DOE OSTI2020

Semiglobal diabatic potential energy matrix for the N–H photodissociation of methylamine

We constructed an analytic diabatic potential energy matrix (DPEM) that describes the N–H photodissociation of methylamine; the electronic state space includes the ground and first excited singlet states. The input for the fit was calculated by extended multi-state complete active space second-order perturbation theory. The data were diabatized using the dipole–quadrupole diabatization method in which we incorporated a coordinate-dependent weighting scheme for the contribution of the quadrupole moments. Therefore, to make the resulting potential energy surfaces semiglobal, we extended the anchor points reactive potential method, a multiscale approach that assigns the internal coordinates to categories with different levels of computational treatment. Key aspects of the adiabatic potential energy surfaces obtained by diagonalizing the DPEM agree with the available experimental and theoretical data at energies relevant for photochemical studies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗